Когда я запускаю GEE, я использую glht () для теста Post ho c, после этого я использую confint (), чтобы получить значение 95% CI. Но что странно, что результат менялся каждый раз, когда я запускаю код снова и снова. Почему?
library(multcomp)
library(geepack)
data(dietox)
#dietox$Weight = as.integer(dietox$Weight)
dietox$Cu <- as.factor(dietox$Cu)
mf <- formula(Weight ~ Cu + Time + I(Time^2) + I(Time^3))
gee1 <- geeglm(mf, data=dietox, id=Pig, corstr="ar1")
CI <- glht(gee1,mcp("Cu" = "Tukey"))
confint(CI)
result:
Estimate lwr upr
Cu035 - Cu000 == 0 -0.4467 -3.7996 2.9061
Cu175 - Cu000 == 0 1.1821 -2.9923 5.3565
Cu175 - Cu035 == 0 1.6288 -2.6616 5.9193
run the code again:
Estimate lwr upr
Cu035 - Cu000 == 0 -0.4467 -3.8002 2.9067
Cu175 - Cu000 == 0 1.1821 -2.9931 5.3573
Cu175 - Cu035 == 0 1.6288 -2.6624 5.9201
again:
Estimate lwr upr
Cu035 - Cu000 == 0 -0.4467 -3.7985 2.9050
Cu175 - Cu000 == 0 1.1821 -2.9910 5.3552
Cu175 - Cu035 == 0 1.6288 -2.6603 5.9179
.......