前言

正交实验设计分析很简单，不用那么复杂（分析交互作用）。正交实验设计的目的有两个，1是找到最佳组合方案，可以直观找最大的观测值，也可以通过极差法找到每个因素的最佳水平进行组合。2是对各因素的影响程度进行排序，可以用极差法或者方差分析法，也很简单。

一篇教程就看明白：https://wenku.baidu.com/view/d46a06aab8f3f90f76c66137ee06eff9aef849b8.html

导入数据

library("readxl")       # necessary to import the data from Excel file
library("dplyr")        # filter and reformat data frames
data<- read_excel("lwx.xlsx")%>%mutate_at(vars(T,pH,H,C),as.factor)

方差分析

fit<-lm(Y~T*H*pH*C,data)
summary.aov(fit)
##             Df Sum Sq Mean Sq F value   Pr(>F)    
## T            2   8558    4279    89.1 4.60e-10 ***
## H            2  32387   16194   337.2 5.42e-15 ***
## pH           2  42554   21277   443.1 4.91e-16 ***
## C            2  16711    8356   174.0 1.68e-12 ***
## Residuals   18    864      48                     
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

根据F value比较各因素的重要性，pH(443.1)>H(337.2)>C(174.0)>T(89.1)。