问题描述
以下是我所做的以及每次收到的错误消息
gee_fit <- geeglm(rslt~ ns(day,3),data = epidem,id =patid,family = binomial,corstr = 'exchangeable')
获取所选日期值的预测概率
predict(gee_fit,newdata = data.frame(day = c(20,30)),type="response")
day_grid<-seq(from=min(day),to=max(day))
preds<-predict(gee_fit,newdata=list(day_grid),se=T)
XRinv ^ 2％*％rep（res.var，p）中的错误： 需要数字/复杂矩阵/矢量参数
我要去哪里错了？
数据
structure(list(day = structure(c(2,8,13,15,12,29,33,22,17,21,44,18,26,19,5,10,16,41,2,20,28,34,30,37,14,9,24,60,42,48,52,58,51,57,64,25,11,27,32,38,7,35,50,3,59,47,53,61,56,84,4,65,40,74,6,23,49,43,54,68,36,46,55,63,45,72,73,78,75,7),label = "day",format.stata = "%9.0g"),rslt = structure(c(0,1,1),label = "test rslts",patid = structure(c(147,147,148,149,150,151,152,153,154,155,156,157,158,159,160,161,162,163,164,165,166,167,168,169,170,171,172,173,174,175,176,177,178,179,181,182,183,184,185,186,187,188,189,190,191,192,193,194,195,196,197,198,199,200,77,83,108,81,71,66,111,140,107,139,116,117,82,104,96,138,110,86,114,130,67,106,85,133,91,144,105,100,124,103,128,143,95,123,137,134,97,118,119,120,31,126,76,94,109,93,115,141,99,62,39,113,146,121,92,80,87,102,136,131,89,127,145,125,122,129,90,135,112,142,69,98,101,132,88,70),label = "group(record)",format.stata = "%9.0g")),row.names = c(NA,-316L),class = c("tbl_df","tbl","data.frame"))
解决方法
暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！
如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@）