[R] Network meta-analysis, varConstPower in nlme - Thank you!

[Christian Gold]

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?)

Best,
Christian


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!?
> 
> Christian
>