练习R：DoE.base包生成标准L9-3-4正交表

R语言可通过DoE.base包的oa.design()函数实现正交表的输出，但是我发现它默认输出的正交表和标准正交表有些差异，还需要简单调整行列次序，最终转换为标准表。

以L9-3-4标准表举例。

1.直接生成正交表

design1 <- oa.design(L9.3.4)

结果如下：

这9个组合是标准表的组合，但是行没有排序，另外因素列的次序也不对。为了和同行交流嘛，我们的正交表尽量是标准的，不然人家以为你搞错了设计。因此我们有必要调整行排序和列次序。

2.调整列次序

仔细对比发现，标准表的四列依次为DCAB，咱们需要把现在的ABCD次序修改为DCAB。这里就使用R的数据框操作即可。

mydesign1 <- design1[,c(4,3,1,2)]

3.调整行排序

依次按第一列和第二列升序排序，

arrange(mydesign1,D,C)

上表就是一个标准的L9-3-4正交表格式。

同类文章推荐阅读：

做论文要用正交表？我打包送给你

本文完

文/图=数据小兵

更多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函数排序