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

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

chemical<-data.frame(

X=c(73,68,74,71,67,73,67,75,72,70,75,68,78,73,68,73,71,75,75,69),

A = gl(4,5,20), #注释1 主因子chemical types

B = gl(5,1,20)) #区组

chemical.aov<-aov(X~A+B, data= chemical)

summary(chemical.aov)

> summary(chemical.aov)

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

A 3 12.95 4.32 2.376 0.121

B 4 157.00 39.25 21.606 2.06e-05 ***

Residuals 12 21.80 1.82


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

There is no difference among the chemical types at α = 0.05 level.

baterial<-data.frame(X=c(13,22,18,39,16,24,17,44,5,4,1,22),

A = gl(3,4,12), #注释1 主因子chemical types

B = gl(4,1,12)) #区组

baterial.aov<-aov(X~A+B, data= baterial)

summary(baterial.aov)

> summary(baterial.aov)

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

A 2 703.5 351.8 40.72 0.000323 ***

B 3 1106.9 369.0 42.71 0.000192 ***

Residuals 6 51.8 8.6


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

There is a difference between the means of the three solutions. The Fisher LSD procedure indicates that solution 3 is significantly different than the other two.

相关推荐
qianmo20212 天前
gpt-4o+deepseek+R生成热力图表
java·数据库·r语言
陈天白2 天前
RNA-seq分析之最佳cutoff(TCGA版)
r语言·生物信息学·rna-seq
Biomamba生信基地3 天前
挑战用R语言硬干一百万单细胞数据分析
开发语言·数据分析·r语言·生信·医药
czhc114007566312 天前
LINUX913 shell:set ip [lindex $argv 0],\r,send_user,spawn ssh root@ip “cat “
tcp/ip·r语言·ssh
zhangfeng113312 天前
win7 R 4.4.0和RStudio1.25的版本兼容性以及系统区域设置有关 导致Plots绘图面板被禁用,但是单独页面显示
开发语言·人工智能·r语言·生物信息
zhangfeng113313 天前
在 R 语言里,`$` 只有一个作用 按名字提取“列表型”对象里的单个元素 对象 $ 名字
开发语言·windows·r语言
高-老师13 天前
R语言生物群落(生态)数据统计分析与绘图实践技术应用
开发语言·r语言·生物群落
大千AI助手13 天前
残差:从统计学到深度学习的核心概念
人工智能·深度学习·resnet·统计学·方差分析·残差·残差分析
WangYan202213 天前
R语言:数据读取与重构、试验设计(RCB/BIB/正交/析因)、ggplot2高级绘图与统计检验(t检验/方差分析/PCA/聚类)
r语言·ggplot2·dplyr
zhangfeng113314 天前
错误于make.names(vnames, unique = TRUE): invalid multibyte string 9 使用 R 语言进行数据处理时
开发语言·r语言·生物信息