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

编程问答 2022-04-23

问题描述

示例：

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 (将#修改为@）

lmrr