我以前能够使用broom :: tidy和broom ::浏览rms :: lrm对象。但是,尝试这样做会导致错误。
以前,bms :: tidy对rms:lrm对象提供以下内容(值已组成):
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
但现在我得到了错误。我尝试安装以前版本的扫帚和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