Displaying 2 results from an estimated 2 matches for "treatmenttx".
2013 Jul 05
1
multcomp on significant interaction in coxme model
...0 1
naked:cocooned:uv - tx 0 0 0 0 -1 1
then I constructed a secon matrix
tmp<-expand.grid(treatment=unique(removal.lapf$treatment),cover=unique(removal.lapf$cover))
X<-model.matrix(~cover*treatment,data=tmp)
X
which gives the following
(Intercept) covernaked treatmenttx treatmentuv covernaked:treatmenttx
covernaked:treatmentuv
1 1 1 1 0
1 0
2 1 1 0 1
0 1
3 1 1 0 0
0...
2012 Apr 16
0
warning message: coxme with package multcomp
...557399
Chisq df p AIC BIC
Integrated loglik 3.04 4 0.55143 -4.96 -1.36
Penalized loglik 3.04 2 0.21879 -0.96 0.84
Model: Surv(day, status) ~ treatment + (1 | r.brood) + (1 | r.worker)
Fixed coefficients
coef exp(coef) se(coef) z p
treatmenttx -22.1625984 2.370859e-10 0.00000 -Inf 0.00
treatmentuv -0.6931777 4.999847e-01 1.22475 -0.57 0.57
Random effects
Group Variable Std Dev Variance
r.brood Intercept 0.0199992574 0.0003999703
r.worker Intercept 0.0199992574 0.0003999703
to see if removal is overall dependend upo...