问题描述
我有一个函数(患病率软件包),在给定的数量为正数(x)和采样数为(n)的情况下,可以估算病原体的真实患病率。该功能需要将锯齿加载到计算机上以进行贝叶斯估计。这是单个条目以及输出的功能:
truePrev(x=432,n=1122,SE= ~ dbeta(13.67,3.95),SP= ~ dbeta(15.79,1.30))
mean median mode sd 2.5% 97.5%
TP 0.444 0.443 0.434 0.102 0.233 0.651
SE 0.766 0.776 0.796 0.098 0.559 0.926
SP 0.915 0.933 0.972 0.068 0.740 0.995
Multivariate BGR statistic = 1.0014
BGR values substantially above 1 indicate lack of convergence
现在,我正在尝试使用purrr map2在500行以上的数据框中迭代该函数:
df<- data.frame("ID" = c(1,2,3),"x" = c(5,10,15),"n" = c(20,25,30))
df_out<-map2(.x = x,.y= n,truePrev,SE=~ dbeta(13.67,SP=~
dbeta(15.79,1.30))
我可以获取输出,但是想要将其绑定到原始数据帧,以便它对应于正确的ID行。我一直在研究Stack Overflow,但是此输出的格式保留了我从工作中读取的其他一些建议。当我尝试将输出从列表更改为数据框时,出现此错误:
as.data.frame(df_out)
Error in as.data.frame.default(x[[i]],optional = TRUE) :
cannot coerce class ‘structure("prev",package = "prevalence")’ to a
data.frame
解决方法
暂无找到可以解决该程序问题的有效方法,小编努力寻找整理中!
如果你已经找到好的解决方法,欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@)