为r中的多个级别创建置信区间图
我有这个数据框:为r中的多个级别创建置信区间图,r,regression,gee,R,Regression,Gee,我有这个数据框: mac<-c( 0.125150, 0.045800, -0.955299, -0.232007, 0.120880, -0.041525, 0.290473, -0.648752, 0.113264, -0.403685) narc<-c(-0.634753, 0.000492, -0.178591, -0.202462, -0.592054, -0.583173, -0.632375, -0.176673, -0.680557, -0.062127)
mac<-c( 0.125150, 0.045800, -0.955299, -0.232007, 0.120880, -0.041525, 0.290473, -0.648752, 0.113264, -0.403685)
narc<-c(-0.634753, 0.000492, -0.178591, -0.202462, -0.592054, -0.583173, -0.632375, -0.176673, -0.680557, -0.062127)
darc<-c(-0.434753, 0.000592, -0.278591, -0.402462, -0.692054, -0.783173, -0.732375, -0.576673, -0.880557, -0.162127)
ideo<-c(1,2,3,2,3,2,3,1,2,1)
ex<-data.frame(mac,darc,narc,ideo)
macgeeglm(ideo~…)
给出了一个错误。我编辑了。现在检查plzOk现在没有错误,但是没有ideo
的级别。级别应该是1,2,3,它们不会出现在spss
数据帧中。而ideo
是响应,为什么要对每个级别的响应进行估计?听起来很奇怪。。。
library(geepack)
spss<-coef(summary(geeglm(as.integer(ideo) ~ mac+narc+darc,data = ex, id = ideo,
corstr = "independence"))) %>%
mutate(lowerWald = Estimate-1.96*Std.err, # Lower Wald CI
upperWald=Estimate+1.96*Std.err, # Upper Wald CI
df=1,
ExpBeta = exp(Estimate)) %>% # Transformed estimate
mutate(lWald=exp(lowerWald), # Upper transformed
uWald=exp(upperWald)) # Lower transformed
spss
#subsets for ever ideo
ex1<-ex[ex$ideo%in% as.factor(1),] #right men
ex2<-ex[ex$ideo%in% as.factor(2),] #center men
ex3<-ex[ex$ideo%in% as.factor(3),] #far men