1

我以前可以在 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
4

0 回答 0