实验设计与分析(第6版,Montgomery)第4章随机化区组,拉丁方, 及有关设计4.5节思考题4.26~4.27 R语言解题

本文是实验设计与分析(第6版,Montgomery著,傅珏生译) 第章随机化区组,拉丁方, 及有关设计4.5节思考题4.26~4.27 R语言解题。主要涉及方差分析,正交拉丁方。

batch <- c(rep("batch1",5), rep("batch2",5), rep("batch3",5), rep("batch4",5), rep("batch5",5))

acid <- rep(c("oper1", "oper2", "oper3", "oper4", "oper5"),5)

time <- c("A", "B", "C", "D", "E", "B", "C", "D", "E", "A", "C", "D", "E", "A", "B", "D", "E",

"A", "B", "C", "E", "A", "B", "C", "D")

catalyst <- c("a", "b", "y", "d", "e", "y", "d", "e", "a", "b", "e", "a", "b", "y", "d", "b", "y",

"d", "e", "a", "d", "e", "a", "b", "y")

y1 <- c(26,16,19,16,13)

y2 <- c(18,21,18,11,21)

y3 <- c(20,12,16,25,13)

y4 <- c(15,15,22,14,17)

y5 <- c(10,24,17,17,14)

y <- c(y1,y2,y3,y4,y5)

rocket.data <- data.frame(batch, acid,time,catalyst, y)

fit <- lm(y~ acid +batch+ time + catalyst, data=rocket.data)

anova(fit)

> anova(fit)

Analysis of Variance Table

Response: y

Df Sum Sq Mean Sq F value Pr(>F)

acid 4 24.4 6.10 1.0427 0.442543

batch 4 10.0 2.50 0.4274 0.785447

time 4 342.8 85.70 14.6496 0.000941 ***

catalyst 4 12.0 3.00 0.5128 0.728900

Residuals 8 46.8 5.85


Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

order <- c(rep("order 1",4), rep("order 2",4), rep("order 3",4), rep("order 4",4))

oper <- rep(c("oper1", "oper2", "oper3", "oper4"),4)

method <- c("C", "B", "D", "A", "B", "C", "A", "D", "A", "D", "B", "C", "D", "A", "C", "B")

place <- c("b", "y", "d", "a", "a", "d", "y", "b", "d", "a", "b", "y", "y", "b", "a", "d")

y1 <- c(11,10,14,8)

y2 <- c(8,12,10,12)

y3 <- c(9,11,7,15)

y4 <- c(9,8,18,6)

y <- c(y1,y2,y3,y4)

rocket.data <- data.frame(order, oper, method, place, y)

fit <- lm(y~ method + order +oper+ place, data=rocket.data)

anova(fit)

> anova(fit)

Analysis of Variance Table

Response: y

Df Sum Sq Mean Sq F value Pr(>F)

method 3 95.5 31.833 3.4727 0.1669

order 3 0.5 0.167 0.0182 0.9960

oper 3 19.0 6.333 0.6909 0.6157

place 3 7.5 2.500 0.2727 0.8429

Residuals 3 27.5 9.167

相关推荐
猫头虎21 小时前
如何解决 pip install -r requirements.txt 本地轮子路径 ‘./packages/xxx.whl’ 不存在 问题
开发语言·网络·python·r语言·pip·gpu算力·国产
zhangfeng11332 天前
R语言 表达矩阵 count_table 筛选出 行名是 某个 基因的 数据或者某个列中的数据是某个基因的数据
矩阵·r语言·生物信息
zhangfeng11333 天前
geo Counts 数据 ,机器学习 模型的外部验证 ROC外部验证数据处理流程
人工智能·机器学习·r语言·生物信息
zhangfeng11333 天前
生物信息 R语言和 cytoscape 相互沟通的组件RCy3,构建cytoscape网络表 节点类型表 链接边的表,并推送到cytoscape
数据库·r语言·生物信息
zhangfeng11333 天前
R语言 安装老一点的班装包 核心是从CRAN归档(Archive)下载对应版本的安装包
开发语言·r语言
数模加油站4 天前
最新R(4.4.1)及R-studio保姆级安装配置详细教程及常见问题解答
开发语言·windows·数学建模·r语言
czliutz4 天前
R语言绘制股票K线图及布林线
开发语言·r语言
天意生信云7 天前
单细胞转录组:差异基因分析和富集分析
经验分享·数据分析·r语言
zhangfeng11339 天前
wgcna 相关性热图中4个颜色 4个共表达模块 的模块基因是否都要做GO/KEGG分析”,核心取决于你的**研究目标和模块的生物学意义*
开发语言·r语言·生物信息
一口面条一口蒜9 天前
R语言中的获取函数与替换函数
开发语言·r语言