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
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
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
3 matches
Mail list logo
