实验设计与分析(第6版,Montgomery)第3章单因子实验:方差分析3.11思考题3.6 R语言解题

本文是实验设计与分析(第6版,Montgomery著,傅珏生译) 第3章单因子实验:方差分析3.11思考题3.6 R语言解题。主要涉及单因子方差分析,正态性假设检验,残差与拟合值的关系图,TukeyHSD法。

(a) Is there evidence to indicate that dosage level affects bioactivity? Use α = 0.05.

X<-c(24,28,37,30,37,44,31,35,42,47,52,38)

A<-factor(rep(1:3, each=4))

miscellany<-data.frame(X,A)

aov.mis<-aov(X~A, data=miscellany)

> summary(aov.mis) Df Sum Sq Mean Sq F value Pr(>F)

A 2 450.7 225.33 7.036 0.0145 *

Residuals 9 288.2 32.03


Signif. codes:

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

There appears to be a different in the dosages.

(b) If it is appropriate to do so, make comparisons between the pairs of means. What conclusions can you draw?

> TukeyHSD(aov.mis)

Tukey multiple comparisons of means

95% family-wise confidence level

Fit: aov(formula = X ~ A, data = miscellany)

$A

diff lwr upr p adj

2-1 7 -4.172869 18.17287 0.2402975

3-1 15 3.827131 26.17287 0.0114434

3-2 8 -3.172869 19.17287 0.1680265

Because there appears to be a difference in the dosages, the comparison of means is appropriate.

opar <- par(mfrow=c(2,2),cex=.8)

plot(aov.mis)

par(opar)

相关推荐
烟锁池塘柳03 小时前
【R语言】R语言中 rbind() 与 merge() 的区别详解
sql·r语言
夏婵语冰1 天前
实用R语言机器学习指南:从数据预处理到模型实战(附配套学习资源)
开发语言·学习·r语言
Yolo566Q2 天前
“R语言+遥感”的水环境综合评价方法实践技术应用
开发语言·r语言
天意生信云2 天前
生信分析自学攻略 | R语言数据筛选和修改
数据结构·经验分享·r语言
烟锁池塘柳04 天前
【R语言】R 语言中 gsub 与正则表达式详解(含 POSIX 与 Perl 风格实例)
正则表达式·r语言·perl
梦想的初衷~5 天前
“R+遥感”的水环境综合评价方法实践技术应用
r语言·水环境·水文
天意生信云6 天前
生信分析自学攻略 | R软件和Rstudio的安装
经验分享·r语言
TS的美梦6 天前
【1:1复刻R版】python版火山图函数一键出图
开发语言·python·r语言·scanpy·火山图
nonono6 天前
深度学习——R-CNN及其变体
深度学习·r语言·cnn
烟锁池塘柳07 天前
【R语言】R 语言中打印含有双引号的字符串时会出现 “\” 的原因解析
r语言