如何计算“鹅”对象的置信区间

Question

我想知道如何计算“鹅”对象的置信区间，即序数广义估计方程的回归输出。 tidy() 或 conf.int() 不接受这种类型的对象。

请参阅下面的可重现示例。

n2 <- 3

n1 <- 7

data <- data.frame(id= c("101", "102", "103", "104", "105", "106", "107", "108", "109", "110",

"101", "102", "103", "104", "105", "106", "107", "108", "109", "110",

"101", "102", "103", "104", "105", "106", "107", "108", "109", "110",

"101", "102", "103", "104", "105", "106", "107", "108", "109", "110",

"101", "102", "103", "104", "105", "106", "107", "108", "109", "110"),

exposure= c("group1", "group3", "group6", "group2", "group4", "group5", "group1", "group3", "group6", "group2",

"group4", "group5", "group7", "group7", "group1", "group3", "group4", "group5", "group7", "group7",

"group6", "group2", "group6", "group2", "group6", "group2", "group6", "group2", "group6", "group2",

"group6", "group2","group6", "group2" , "group6", "group2", "group6", "group2","group6", "group2",

"group6", "group2", "group6", "group2", "group6", "group2", "group6", "group2", "group6", "group2"),

outcome=c("1", "2", "3", "2", "1", "3", "1", "1", "2", "3",

"1", "3", "1", "1", "3", "1", "1", "3", "1", "1",

"1", "3", "1", "1", "3", "1", "1", "3", "1", "1",

"1", "3", "1", "1", "3", "1", "1", "3", "1", "1",

"1", "3", "1", "1", "3", "1", "1", "3", "1", "1"))

data$outcome <- factor(data$outcome, ordered = TRUE)

data <- data %>%

arrange(id)

model1 <- ordgee(outcome ~  exposure,

id = id,

data = data,

corstr= "independence")

summary(model1)

我尝试使用conf.int()、tidy()，但它们不接受“鹅”对象。

我也尝试了here

建议的解决方案

但是当我运行时出现错误。

Answer 1

也在想办法解决这个问题！如果您得到答案请告诉我。