如何解决如何从 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 举报,一经查实,本站将立刻删除。