Question

我正在编写函数来输出ANOVA作为输出

我不明白如何从以下信息输出anova对象：

# degrees of freedom 
    repdf = 1
    trtdf = 22
    totaldf = 23
 # sum of square    
    ssrep = 10.3
    sstrt = 14567.2
    sstotal = 14577.2

anova对象数据框或列表还是其他特殊编程类别？

编辑：根据Ben的建议

Source <- c("replication", "Treatments", "Total") 
Df <- c(repdf, trtdf, totaldf)
"Sum Sq" <- c(ssrep, sstrt, sstotal)
anovadf <- data.frame(Source, Df, "Sum Sq")
class(anovadf) <- c("anova","data.frame")

不能告诉我anova对象的str应该是什么样的？任何进一步的帮助

> str(anovadf)
Classes ‘anova’ and 'data.frame':       3 obs. of  3 variables:
 $ Source   : Factor w/ 3 levels "Error","replication",..: 2 3 1
 $ Df       : num  1 22 23
 $ X.Sum.Sq.: Factor w/ 1 level "Sum Sq": 1 1 1

Answer 1

创建一个anova对象，保存，然后在结果上使用str()。来自lm.D9创建的example("lm")对象：

> str(anova(lm.D9))
Classes ‘anova’ and 'data.frame':   2 obs. of  5 variables:
 $ Df     : int  1 18
 $ Sum Sq : num  0.688 8.729
 $ Mean Sq: num  0.688 0.485
 $ F value: num  1.42 NA
 $ Pr(>F) : num  0.249 NA
 - attr(*, "heading")= chr  "Analysis of Variance Table\n" "Response: weight"

所以这是数据框的一个特例。构建数据框a以匹配示例，然后尝试分配类：class(a) <- c("anova","data.frame")。

特别是：

Df <- c(repdf, trtdf, totaldf)
ssq <- c(ssrep, sstrt, sstotal)

anovadf <- data.frame(Df, `Sum Sq`=ssq, `Mean Sq`=ssq/Df, check.names=FALSE)
rownames(anovadf) <- c("replication","treatments","total")
class(anovadf) <- c("anova","data.frame")

anovadf
            Df  Sum Sq Mean Sq
replication  1    10.3   10.30
treatments  22 14567.2  662.15
total       23 14577.2  633.79

你必须对列名有点小心 - 它们必须受反引号保护，你必须使用check.names=FALSE，因为它们不是合法的变量名（它们包含空格）。您可以将F统计值和P值添加到此 - 我没有，因为我不确定适当的错误术语是什么。

构造像对象一样的ANOVA

1 个答案: