[R] Monotone Transformation

2009-04-17 Thread Feng Jingyu
Hi, I am trying to use R to mimic what I did in SAS. proc transreg data=x ; model identity(GSI)=monotone(group1); output out=d2 pprefix=M; run; Accroding to SAS documentation, the MONOTONE transfomation algorithm comes from (Kruskal 1964, secondary approach to ties). I have

Re: [R] Fit unequal variance model in R

2009-04-04 Thread Feng Jingyu
Yes, I am considering I am a lucky man because you answered my question. Have a great weekend! Dieter Menne wrote: > > > Feng Jingyu wrote: >> >> Hi For my purpose, I need to match variance estimates for each group from >> R and SAS. They do match now. >>

Re: [R] Fit unequal variance model in R

2009-04-04 Thread Feng Jingyu
Hi For my purpose, I need to match variance estimates for each group from R and SAS. They do match now. I notice there are still some difference between those two. For example the standard errors of coeffieicints are different, but I don't care those. Dieter Menne wrote: > > > > To conserve y

Re: [R] Fit unequal variance model in R

2009-04-04 Thread Feng Jingyu
x -1.9580366 -0.7219450 -0.1989222 0.9096203 1.8319336 Residual standard error: 0.2354039 Degrees of freedom: 15 total; 11 residual Dieter Menne wrote: > > > Feng, Jingyu wrote: >> >> I'am trying to develop some code if R, which would correspond to what

Re: [R] Fit unequal variance model in R

2009-04-04 Thread Feng Jingyu
Thanks a lot. The problem is solved. It took me a while to understand the output from the R. With little calculation, I am able to match results from R to SAS. Dieter Menne wrote: > > > Feng Jingyu wrote: >> >> I used gls and it still does not provide me different est

[R] Fit unequal variance model in R

2009-04-01 Thread Feng, Jingyu
I'am trying to develop some code if R, which would correspond to what I did in SAS. The data look like: TreatmentReplicategroup1 GSI Control A 1 0.81301 Control B 1 1.06061 Control C 1 1.26350 Control