Displaying 1 result from an estimated 1 matches for "lumley1".
Did you mean:
lumley
2008 Oct 15
2
Network meta-analysis, varConstPower in nlme
...ound it very
helpful that you included the R code to replicate your analysis;
however, I have had a problem replicating your example and wondered if
you are able to give me a hint. When I use the code from the article:
lme1 <- lme(Y1 ~ trt.B + trt.C + trt.D + trt.E, random = ~ 1 | trtpair,
data=lumley1, var = varConstPower(form=~sigma, fixed=list(power=1)))
I get an error message:
Error in lme(Y1 ~ trt.B + trt.C + trt.D + trt.E, random = ~1 | trtpair, :
unused argument(s) (var = list(const = numeric(0), power = numeric(0)))
The problem seems to be in the varConstPower component, but I don&...