练习R:正交试验极差分析、方差分析与统计图形

原创 数据小兵  2021-06-23 08:38  阅读 787 次
数据小兵成长记

小兵尝试用R来对一个简单正交试验进行数据处理分析,绘制相应的可视化图形。

研究某工艺最优条件,重点考虑4个试验因素ABCD,同时考察AB交互作用,各因素均取2水平。

本例选择L8.2.7正交表安排试验因素与交互作用。

试验方案及结果(部分)表格如下,数据文件命名为design2。

先基于aov()计算方差分析基础数据分析结果,

aov.mod <- aov(results~tem+oxy+wat+ph+tem:oxy,data = design2)

利用FrF2包绘制主效应图以及输出简要描述统计,

summary(MEPlot(aov.mod))

ABCD四因素极差依次为2、3、5.5、1.5,得四因素主次顺序依次为C>B>A>D。

四因素主效应图如下:

可见各因素的优水平A1B1C2D2。(这个结论有瑕疵,注意应当结合交互作用来讨论)。

接下来看看方差分析表结果:

summary(aov.mod)

可见C因素和AB交互显著(P<0.05),另外三个因素不显著。C为主要因素,其他为次要因素。

利用FrF2包绘制AB交互作用图,

IAPlot(aov.mod,select = 1:2)

两条直线明显相交,可辅助判断存在交互作用(仍以显著性检验为准)。

本文完
文/图=数据小兵
参考资料:

郑杰. 试验设计与数据分析: 基于R语言应用[M]. 华南理工大学出版社, 2016.

 

更多R统计文章
 
R语言单一样本t检验案例实现
R语言的正态密度曲线很美
R自带pairs函数矩阵散点图
car包spm函数矩阵散点图
用R语言pie函数做饼图
ggplot2统计图形:常见的4种箱线图
ggplot2统计图形:常见的4种直方图
按列索引按列名称删除指定的列数据
R语言Levene方差齐次检验
table函数:分类数据的频数与频率统计
jiebaR包中文分词及词云制作
R语言带文字标签的散点图
如何获取R自带数据集与R包数据集说明文档?
如何修改R数据框的列名称?
二元正态分布及双变量相关分析简单案例演示
R语言相关系数、显著性检验及可视化的尝试
ggplot2统计图形:常见的4种散点图
R语言scale()函数实现数据标准化
用R语言做单因素方差分析及多重比较
用R语言自动智能化创建时间序列ARIMA模型
用R对连续数据做描述统计
R语言主成分分析
shapiro.test()与ggqqplot()做正态分布检验
Fleiss' Kappa系数处理多个观察者一致性检验
练习R:用ggpubr包ggline绘制线图
双因素方差分析如何做非参数检验替代方法?
练习R:dplyr包实现分组描述统计
练习R:单个连续数据的探索性统计可视化函数EDA
练习R:用lm.ridge()做岭回归分析,可惜无法输出R平方
练习R:lm+plot+abline+text四函数绘制线性拟合散点图
练习R:用lm.ridge()做岭回归分析,可惜无法输出R平方
练习R:分类自变量的线性回归
练习R:factor函数为有序分类变量创建值标签
学习R:识别缺失值与将指定数据编码为缺失值
练习R:car包recode函数多分类变量的重新编码
如何在Rstudio中退出/取消par分区绘图效果?
如何删除environment中的数据集?
练习R:stepAIC多元逐步回归
练习R:用3d空间图展示多元线性回归模型
为什么要学R语言?如何学习R语言?
练习R:用3d空间图展示多元线性回归模型
练习R:dplyr包arrange函数排序
练习R:DoE.base包生成标准L9-3-4正交表
本文地址:http://www.datasoldier.net/archives/3077
版权声明:本文为原创文章,版权归 数据小兵 所有,欢迎分享本文,转载请保留出处!
视频课程《SPSS统计分析:从入门到实践提高》
欢迎订阅SPSS训练营微信公众号

评论已关闭!