Displaying 1 result from an estimated 1 matches for "rats1".
Did you mean:
rats
2012 Sep 14
1
Correlation between random effects in the package coxme
Hello,
Why the correlation between the random effects is negative?
library(coxme)
rats1 <- coxme(Surv(time, status) ~ (1|litter), rats)
random.effects(rats1)[[1]] #one value for each of the 50 litters
print(rats1)
rats2 <- lmekin(time ~ (1|litter), rats)
fixed.effects(rats2)
random.effects(rats2)[[1]] #one value for each of the 50 litters
print(rats2)
cor(random.effects(rats1)...