以前は、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)
ほうき::tidy エラー
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
ほうき::glance エラー
Error in summary.rms(x) :
adjustment values not defined here or with datadist for mpg