实验设计与分析(第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)

相关推荐
魔力之心10 小时前
R study notes[1]
r语言
Chef_Chen20 小时前
从0开始学习R语言--Day57--SCAD模型
开发语言·学习·r语言
医工交叉实验工坊20 小时前
R 语言绘制六种精美热图:转录组数据可视化实践(基于 pheatmap 包)
开发语言·信息可视化·r语言
AAIshangyanxiu1 天前
最新基于R语言结构方程模型分析与实践技术应用
开发语言·r语言·结构方程模型·生态统计学
biomooc2 天前
R拟合 | 一个分布能看到三个峰,怎么拟合出这三个正态分布的参数? | 高斯混合模型 与 EM算法
r语言
请你喝好果汁6412 天前
R中匹配函数
开发语言·r语言
Tiger Z2 天前
R 语言科研配色 --- 第 81 期 (附免费下载的配色绘图PPT)
r语言·科研·配色
恩喜玛生物5 天前
100个GEO基因表达芯片或转录组数据处理27 GSE83456
深度学习·r语言
Yolo566Q7 天前
基于现代R语言【Tidyverse、Tidymodel】的机器学习方法
开发语言·经验分享·机器学习·r语言
Tiger Z7 天前
《R for Data Science (2e)》免费中文翻译 (第0章) --- Introduction
r语言·数据科学·中文翻译