[R] Network meta-analysis, varConstPower in nlme - Thank you!
c.gold at magnet.at
Thu Oct 16 14:52:39 CEST 2008
Excellent - now it works. Thank you, Christian!
(I still seem to be unable to replicate the results as reported in the
paper, but that is unrelated to R-help. Thomas, do you have any idea
what I may have entered wrongly?)
Christian Ritz wrote:
> Hi Christian,
> I believe that the argument "var" has changed name to "weights".
> The following lines work for me:
> sigma <- rep(sqrt(.5), nrow(lumley1)) # not nrow=
> lme1 <- lme(Y1 ~ trt.B + trt.C + trt.D + trt.E, random = ~ 1 | trtpair,
> data=lumley1, weights = varConstPower(form=~sigma, fixed=list(power=1)))
> lme2 <- lme(Y2 ~ trt.B + trt.C + trt.D + trt.E, random = ~ 1 | trtpair,
> data=lumley1), weights = varConstPower(form=~sigma, fixed=list(power=1)))
> I hope you can now make it work!?
More information about the R-help