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