我以前可以在 rms::lrm 对象上使用 broom::tidy 和 broom::glance 。但是,现在尝试这样做会导致错误。
以前,rms:lrm 对象上的 broom::tidy 将提供以下内容(值是由组成的):
term estimate std.error statistic p.value
y>=six -2.1 0.76 -2.78 0.01
y>=eight -4.98 1.16 -4.31 0
mpg 0.08 0.02 3.97 0
和 broom::glance 将输出以下值:
obs max.deriv model.lr df p.value concordance somers.d gamma tau.a r.squared brier g gr gp
但现在我得到了错误。我尝试安装 broom 和 rms 的早期版本,但这并没有解决错误。我目前正在使用R 3.4.2
library(dplyr)
library(rms)
library(broom)
data("mtcars")
df <- mtcars %>% select(mpg, cyl)
df$CYL <- factor(df$cyl, levels = c(6,4,8), labels=c("six", "four", "eight")) #ordered factor for ordinal logistic
fit <- lrm(CYL~mpg, data=df)
tidy(fit)
扫帚::整洁的错误
Error in as.data.frame.default(x) :
cannot coerce class "c("lrm", "rms", "glm")" to a data.frame
In addition: Warning message:
In tidy.default(fit) :
No method for tidying an S3 object of class lrm , using as.data.frameNo method for tidying an S3 object of class rms , using as.data.frameNo method for tidying an S3 object of class glm , using as.data.frame
扫帚::一瞥错误
Error in summary.rms(x) :
adjustment values not defined here or with datadist for mpg