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

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

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

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

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

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

y1 <- c(8,7,1,7,3)

y2 <- c(11,2,7,3,8)

y3 <- c(4,9,10,1,5)

y4 <- c(6,8,6,6,10)

y5 <- c(4,2,3,8,8)

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

rocket.data <- data.frame(batch, oper, seed, y)

matrix(rocket.data$seed, nrow=5, byrow=T)

rocket.data

fit <- lm(y~seed+batch+oper, data=rocket.data)

anova(fit)

> anova(fit)Analysis of Variance Table

Response: y

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

seed 4 141.44 35.360 11.3092 0.0004877 ***

batch 4 15.44 3.860 1.2345 0.3476182

oper 4 12.24 3.060 0.9787 0.4550143

Residuals 12 37.52 3.127


Signif. codes:

0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

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

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

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

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

y2 <- c(7,18,11,8)

y3 <- c(5,10,11,9)

y4 <- c(10,10,12,14)

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

rocket.data <- data.frame(batch, oper, seed, y)

matrix(rocket.data$seed, nrow=4, byrow=T)

rocket.data

fit <- lm(y~seed+batch+oper, data=rocket.data)

anova(fit)

> anova(fit)Analysis of Variance Table

Response: y

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

seed 3 72.5 24.1667 13.8095 0.004213 **

batch 3 18.5 6.1667 3.5238 0.088519 .

oper 3 51.5 17.1667 9.8095 0.009926 **

Residuals 6 10.5 1.7500


Signif. codes:

0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

相关推荐
love530love15 小时前
怎么更新 cargo.exe ?(Rust 工具链)
人工智能·windows·python·rust·r语言
烟锁池塘柳04 天前
【R语言】R语言中 rbind() 与 merge() 的区别详解
sql·r语言
夏婵语冰5 天前
实用R语言机器学习指南:从数据预处理到模型实战(附配套学习资源)
开发语言·学习·r语言
Yolo566Q6 天前
“R语言+遥感”的水环境综合评价方法实践技术应用
开发语言·r语言
天意生信云6 天前
生信分析自学攻略 | R语言数据筛选和修改
数据结构·经验分享·r语言
烟锁池塘柳08 天前
【R语言】R 语言中 gsub 与正则表达式详解(含 POSIX 与 Perl 风格实例)
正则表达式·r语言·perl
梦想的初衷~9 天前
“R+遥感”的水环境综合评价方法实践技术应用
r语言·水环境·水文
天意生信云10 天前
生信分析自学攻略 | R软件和Rstudio的安装
经验分享·r语言
TS的美梦10 天前
【1:1复刻R版】python版火山图函数一键出图
开发语言·python·r语言·scanpy·火山图
nonono10 天前
深度学习——R-CNN及其变体
深度学习·r语言·cnn