前言
这是一个非常典型的区组设计,可考虑随机区组设计的方差分析。随机区组设计的方差分析在统计上常常归在单因素方差分析里面,“区组”并不是严格意义的“因素”(因素字面是相对处理或者叫干预因素而言的),但在进行统计分析的时候,模型是将其作为因素或协变量来处理的。其分析方法跟单因素的方差分析采用的几乎是相同的函数,只是多了一个自变量而已。
本文涉及到命令有read_excel{readxl}、lm {stats}、Anova {car}、aov {stats}等。
【1】数据载入
library(readxl)
bFd <- read_excel("D:/Temp/rbFdata.xlsx")
【2】方差分析
[2.1]回归法:lm {stats}
Block<-factor(bFd$block)
blrfit=lm(bFd$HBia~Block+bFd$group)