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

2008-10-16 Thread 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

[R] Network meta-analysis, varConstPower in nlme

2008-10-15 Thread Christian Gold
Dear Thomas Lumley, and R-help list members, I have read your article Network meta-analysis for indirect treatment comparisons (Statist Med, 2002) with great interest. I found it very helpful that you included the R code to replicate your analysis; however, I have had a problem replicating your

Re: [R] Network meta-analysis, varConstPower in nlme

2008-10-15 Thread Christian Ritz
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,