微信公众号搜"智元新知"关注
微信扫一扫可直接关注哦!

如何从 lm() + anova() 但使用 RcppEigen::fastLm() 获得输出?

如何解决如何从 lm() + anova() 但使用 RcppEigen::fastLm() 获得输出?

示例:

df <- cbind.data.frame(y = rnorm(100),x1 = factor(sample(1:2,100,replace = T)))

mod1 <- lm(y ~ x1,df)
anova(mod1)

输出

Analysis of Variance Table

Response: y
          Df Sum Sq Mean Sq F value Pr(>F)
x1         1  0.388 0.38766  0.3879 0.5349
Residuals 98 97.940 0.99939

使用 RcppEigen:

mod2 <- RcppEigen::fastLm(y = df$y,X = model.matrix(~ x1,df))
anova(mod2)

错误

Error in UseMethod("anova") : 
  no applicable method for 'anova' applied to an object of class "fastLm"

有没有办法在使用 RcppEigen::fastLm() 时复制通过基本 R 函数获得的输出

版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。