实验设计与分析(第6版,Montgomery)第5章析因设计引导5.7节思考题5.13 R语言解题

本文是实验设计与分析(第6版,Montgomery著,傅珏生译) 第5章析因设计引导5.7节思考题5.13 R语言解题。主要涉及方差分析,正态假设检验,残差分析,交互作用图。

dataframe<-data.frame(

y=c(36,18,30,39,20,37,36,22,33,32,20,34),

column=gl(4,3,12),

row=gl(3,1,12))

summary (dataframe)

dataframe.aov2 <- aov(y~column+row,data=dataframe)

summary (dataframe.aov2)

> summary (dataframe.aov2)

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

column 3 28.9 9.64 2.006 0.214717

row 2 580.5 290.25 60.399 0.000106 ***

Residuals 6 28.8 4.81


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

with(dataframe,interaction.plot(column,row,y,type="b",pch=19,fixed=T,xlab="column",ylab="y"))

plot.design(y~column*row,data=dataframe)

fit <-lm(y~column+row,data=dataframe)

anova(fit)

summary(fit)

par(mfrow=c(2,2))

plot(fit)

par(mfrow=c(2,2))

plot(as.numeric(dataframecolumn), fitresiduals, xlab="column", ylab="Residuals", type="p", pch=16)

plot(as.numeric(dataframerow), fitresiduals, xlab="row", ylab="Residuals", pch=16)

相关推荐
shootero@126.com2 小时前
R语言开发记录,二(创建R包)
r语言
shootero@126.com9 小时前
R语言开发记录,一
开发语言·r语言
十三画者11 天前
【科研绘图系列】R语言绘制circos图形(circos plot)
数据挖掘·数据分析·r语言·数据可视化
Morpheon12 天前
使用 R 处理图像
开发语言·计算机视觉·r语言
WW、forever13 天前
服务器手动安装并编译R环境库包:PROJ→RGDAL
运维·服务器·r语言
WW、forever13 天前
【服务器R环境架构】基于 micromamba下载 R 库包
运维·服务器·r语言
Jet450516 天前
第100+42步 ChatGPT学习:R语言实现阈值调整
开发语言·学习·chatgpt·r语言
Smiling63916 天前
【R】基于R实现贝叶斯分析(一)
r语言·贝叶斯分析·bayes
Chef_Chen16 天前
从0开始学习R语言--Day23--稳健回归
学习·回归·r语言
Morpheon19 天前
R语言非结构化文本挖掘入门指南
开发语言·r语言