我需要通过dyads对标准错误进行聚类,我正在尝试使用rms包。为此,我需要使用lrm instrad of glm来拟合我的模型。当我输入library(rms)
时,我收到以下错误:
Error in loadNamespace(i, c(lib.loc, .libPaths()), versionCheck = vI[[i]]) :
there is no package called ‘mvtnorm’
Error: package or namespace load failed for ‘rms’,
当我尝试使用以下命令通过dyads聚类标准错误时,我得到函数lrm不存在的错误。
fit=lrm(disease ~ age + study + rcs(bmi,3), x=T, y=T, data=dataf)
fit
robcov(fit, cluster=dataf$id)
bootcov(fit,cluster=dataf$id)
也有人知道x=T
和y=T
是什么?我刚开始在R工作并且总体上有统计数据,所以这对我来说都不是很明显。