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

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

dataframe <-data.frame(

strength=c(109,110,110,112,116,114,110,115,110,111,112,115,108,109,111,109,114,119,110,108,114,112,120,117),

machine=gl(4,6,24),

operator=gl(3,2,24))

summary (dataframe)

dataframe.aov2 <- aov(strength~operator*machine,data=dataframe)

summary (dataframe.aov2)

> summary (dataframe.aov2)

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

operator 2 160.33 80.17 21.143 0.000117 ***

machine 3 12.46 4.15 1.095 0.388753

operator:machine 6 44.67 7.44 1.963 0.150681

Residuals 12 45.50 3.79


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

with(dataframe,interaction.plot(machine,operator,strength,type="b",pch=19,fixed=T,xlab="operator",ylab="strength"))

plot.design(strength~machine*operator,data=dataframe)

fit <-lm(strength~operator*machine,data=dataframe)

anova(fit)

> anova(fit)

Analysis of Variance Table

Response: strength

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

operator 2 160.333 80.167 21.1429 0.0001167 ***

machine 3 12.458 4.153 1.0952 0.3887526

operator:machine 6 44.667 7.444 1.9634 0.1506807

Residuals 12 45.500 3.792


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

summary(fit)

> summary(fit)

Call:

lm(formula = strength ~ operator * machine, data = dataframe)

Residuals:

Min 1Q Median 3Q Max

-2.5 -1.0 0.0 1.0 2.5

Coefficients:

Estimate Std. Error t value Pr(>|t|)

(Intercept) 1.095e+02 1.377e+00 79.527 <2e-16 ***

operator2 1.500e+00 1.947e+00 0.770 0.4560

operator3 5.500e+00 1.947e+00 2.825 0.0153 *

machine2 3.000e+00 1.947e+00 1.541 0.1493

machine3 -1.000e+00 1.947e+00 -0.514 0.6169

machine4 -5.000e-01 1.947e+00 -0.257 0.8017

operator2:machine2 -3.500e+00 2.754e+00 -1.271 0.2278

operator3:machine2 -4.500e+00 2.754e+00 -1.634 0.1282

operator2:machine3 8.808e-14 2.754e+00 0.000 1.0000

operator3:machine3 2.500e+00 2.754e+00 0.908 0.3818

operator2:machine4 2.500e+00 2.754e+00 0.908 0.3818

operator3:machine4 4.000e+00 2.754e+00 1.453 0.1720


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

Residual standard error: 1.947 on 12 degrees of freedom

Multiple R-squared: 0.827, Adjusted R-squared: 0.6684

F-statistic: 5.214 on 11 and 12 DF, p-value: 0.004136

par(mfrow=c(2,2))

plot(fit)

par(mfrow=c(2,2))

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

plot(as.numeric(dataframeoperator), fitresiduals, xlab="operator", ylab="Residuals", pch=16)

相关推荐
zhangfeng11337 小时前
亲测可用,R语言 ggplot2 箱线图线条控制参数详解,箱线图离散数值控制
开发语言·python·r语言·生物信息
zhangfeng113318 小时前
R 导出 PDF 时中文不显示 不依赖 showtext** 的最简方案(用 extrafont 把系统 TTF 真正灌进 PDF 内核)
开发语言·r语言·pdf·生物信息
生物小卡拉19 小时前
R脚本--PCA分析系列1_v1.0
开发语言·r语言
生物小卡拉19 小时前
R脚本——LefSE分析与可视化-v1
开发语言·r语言
猫头虎2 天前
如何解决 pip install -r requirements.txt 本地轮子路径 ‘./packages/xxx.whl’ 不存在 问题
开发语言·网络·python·r语言·pip·gpu算力·国产
zhangfeng11333 天前
R语言 表达矩阵 count_table 筛选出 行名是 某个 基因的 数据或者某个列中的数据是某个基因的数据
矩阵·r语言·生物信息
zhangfeng11334 天前
geo Counts 数据 ,机器学习 模型的外部验证 ROC外部验证数据处理流程
人工智能·机器学习·r语言·生物信息
zhangfeng11334 天前
生物信息 R语言和 cytoscape 相互沟通的组件RCy3,构建cytoscape网络表 节点类型表 链接边的表,并推送到cytoscape
数据库·r语言·生物信息
zhangfeng11334 天前
R语言 安装老一点的班装包 核心是从CRAN归档(Archive)下载对应版本的安装包
开发语言·r语言
数模加油站5 天前
最新R(4.4.1)及R-studio保姆级安装配置详细教程及常见问题解答
开发语言·windows·数学建模·r语言