问题描述
我在R中具有以下数据帧(my_df):
race trauma y1 y2 y3
0 0 1 0 0
0 0 1 0 0
0 0 1 0 0
0 0 1 0 0
0 0 1 0 0
0 0 1 0 0
0 0 1 0 0
0 0 0 1 0
0 0 0 1 0
0 0 0 1 0
0 0 0 1 0
0 0 0 1 0
0 0 0 1 0
0 0 0 1 0
0 0 0 1 0
0 0 0 1 0
0 0 0 1 0
0 0 0 1 0
0 0 0 1 0
0 0 0 1 0
0 0 0 1 0
0 0 0 1 0
0 0 0 0 1
0 1 1 0 0
0 1 1 0 0
0 1 1 0 0
0 1 1 0 0
0 1 1 0 0
0 1 1 0 0
0 1 1 0 0
0 1 1 0 0
0 1 0 1 0
0 1 0 1 0
0 1 0 1 0
0 1 0 1 0
0 1 0 1 0
0 1 0 1 0
0 1 0 1 0
0 1 0 1 0
0 1 0 1 0
0 1 0 1 0
0 1 0 1 0
0 1 0 1 0
0 1 0 0 1
0 2 1 0 0
0 2 1 0 0
0 2 1 0 0
0 2 1 0 0
0 2 1 0 0
0 2 0 1 0
0 2 0 1 0
0 2 0 1 0
0 2 0 1 0
0 2 0 1 0
0 2 0 1 0
0 2 0 1 0
0 2 0 1 0
0 2 0 1 0
0 2 0 1 0
0 2 0 1 0
0 2 0 1 0
0 2 0 1 0
0 2 0 1 0
0 2 0 1 0
0 2 0 1 0
0 2 0 0 1
0 3 1 0 0
0 3 0 1 0
0 3 0 1 0
0 3 0 1 0
0 3 0 1 0
0 3 0 1 0
0 3 0 1 0
0 3 0 1 0
0 3 0 1 0
0 3 0 1 0
0 3 0 0 1
0 4 1 0 0
0 4 0 1 0
0 4 0 1 0
0 4 0 1 0
0 4 0 1 0
0 5 0 0 1
0 5 0 0 1
1 0 0 1 0
1 0 0 0 1
1 1 0 1 0
1 1 0 1 0
1 1 0 1 0
1 1 0 0 1
1 2 0 1 0
1 2 0 1 0
1 2 0 1 0
1 2 0 0 1
1 3 0 1 0
1 3 0 1 0
1 3 0 0 1
输出dput(my_df)给出
structure(list(race = c(0L,0L,1L,1L),trauma = c(0L,2L,3L,4L,5L,3L),happy = c(1L,3L)),class = "data.frame",row.names = c(NA,-97L))
我使用vgam软件包来拟合累积物流模型,并通过pROC软件包计算一致性指数:
> fit = vglm(cbind(y1,y2,y3) ~ race + trauma,family=cumulative(parallel=TRUE),data=my_df)
> multiclass.roc(as.numeric(cbind(my_df$y1,my_df$y2,my_df$y3)),as.numeric(predict(fit,type="response")),plot=FALSE,legacy.axes=TRUE)
ROC曲线下的面积(即一致性指数)为0.8271。根据我要复制的来源(分类数据分析,Agresti,第3版,第314页),结果应为0.688。拟合的结果与书中介绍的结果一致。我猜我的一致性索引肯定有问题。我不知道是什么。 预先感谢您的帮助。
解决方法
暂无找到可以解决该程序问题的有效方法,小编努力寻找整理中!
如果你已经找到好的解决方法,欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@)