【发布时间】:2021-09-03 01:22:17
【问题描述】:
我在使用 geefirthr 包中的 geefirth 为相关二进制数据运行 Firth 类型 GEE 时遇到问题。主题嵌套在不同的中心。
df <- data.frame(subject_id=c(1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20),
center_id=c(2,2,2,2,2,3,3,3,3,3,3,3,3,4,4,4,4,5,5,5),
x1=c(0,0,0,0,0,0,0,0,0,1,0,1,0,1,1,1,1,0,0,0),
x2=c(0.95,0.95,0.95,0.95,0.95,0.86,0.86,0.86,0.86,0.86,0.86,0.86,0.86,0.90,0.90,0.90,0.90,0.86,0.86,0.86),
x3=c(2.80,2.80,2.8,2.8,2.8,1.5,1.5,1.5,1.5,1.5,1.5,1.5,1.5,2.5,2.5,2.5,2.5,2.6,2.6,2.6),
x4=c(4.72,4.72,4.72,4.72,4.72,3.35,3.35,3.35,3.35,3.35,3.35,3.35,3.35,3.81,3.81,3.81,3.81,4.40,4.40,4.40),
y=c(1,0,0,1,1,0,0,0,0,0,0,0,0,1,1,1,0,0,1,1))
data(df)
geefirthr::geefirth(y ~ x1 + x2 + x3 + x4, id=center_id, data=df, corstr="independence")
我收到以下错误: 矩阵中的错误(x$id,ncol=1):“数据”必须是向量类型,为“NULL”
我尝试了以下方法,is.vector(df) 返回 TRUE。关于如何解决的任何想法?
df <- unlist(df)
df <- as.vector(df)
df <- as.vector(df$center_id)
【问题讨论】:
-
这个包在 CRAN 中是否可用
-
我运行了您的代码,但没有您给出的错误。另外
data(df)不是正确的代码行 -
@akrun 该软件包在 CRAN 和 github 上可用:rdrr.io/github/momenulhaque/geefirthr
-
当我尝试使用 R 4.1.1 使用 install.packages 时,它显示消息包在 CRAN 中不可用
-
@akrun 我使用
remotes::install_github("mhmondol/geefirthr")安装。但是,available::available_on_cran("geefirthr")确认在 R 4.0.2 的 CRAN 上可用。
标签: r