r - 格式化自定义摘要输出以匹配 R 中的 ANOVA 输出

Question

我是 R 的新手。我正在处理一项任务。作业是创建 R 包来模拟方差分析表。我已经创建了任务中规定的所有必要功能。该函数计算出正确的值，但我无法像 R 内置的 anova() 函数那样显示方差分析表。这是我的 summary.oneway 函数

summary.oneway <- function(object, ...){

  #model <- oneway(object)
  fval <- object$FValue
  TAB <- list(t(object$AOV), "Mean Sq."= rbind(object$MSB, object$MSW),
              "F Value" = fval, p.value = object$p.value)

  res <- list(call=object$call, onewayAnova = TAB)
  class(res) <- "summary.oneway"
  res  
}

这是输出:

Analysis of Variance:oneway.formula(formula = coag ~ diet, data = coagdata)[[1]]         Sum of Squares Deg. of Freedomdiet                228               3Residual            112              20$`Mean Sq.`        1[1,] 76.0[2,]  5.6$`F Value`       1 13.57143 $p.value           1 4.658471e-05

实际方差分析输出:

Analysis of Variance TableResponse: coag          Df Sum Sq Mean Sq F value    Pr(>F)    diet       3    228    76.0  13.571 4.658e-05 ***Residuals 20    112     5.6                      ---Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 

如何实现这种格式？我错过了什么？

非常感谢您的帮助。

六合

Answer 1

Anova 输出使用打印方法 print.anova 你可能想看看 methods(print) 尤其是 stats:::print.anova

您很可能想要创建自己的打印功能

 print.oneway <- function(object, ...) {
     foo
     bar
 }