library(openxlsx)
d5.7 = read.xlsx('C:/Users/12290/Desktop/多元统计表格/mvexer5.xlsx',sheet = 'E5.7')
d5.7## 实验号 充磁量 定位角度 线圈匝数 输出力矩
## 1 1 900 10 70 160
## 2 2 900 11 80 215
## 3 3 900 12 90 180
## 4 4 1100 10 70 168
## 5 5 1100 11 80 236
## 6 6 1100 12 90 190
## 7 7 1300 10 70 157
## 8 8 1300 11 80 205
## 9 9 1300 12 90 140anova(lm(输出力矩~充磁量+定位角度+线圈匝数+线圈匝数*充磁量*定位角度+充磁量*定位角度+充磁量*线圈匝数+定位角度*线圈匝数,data = d5.7))## Analysis of Variance Table
##
## Response: 输出力矩
## Df Sum Sq Mean Sq F value Pr(>F)
## 充磁量 1 468.2 468.2 1.2644 0.34268
## 定位角度 1 104.2 104.2 0.2813 0.63259
## 充磁量:线圈匝数 1 342.2 342.2 0.9243 0.40727
## 定位角度:线圈匝数 1 5582.7 5582.7 15.0771 0.03026 *
## 充磁量:定位角度:线圈匝数 1 44.1 44.1 0.1191 0.75284
## Residuals 3 1110.8 370.3
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1P（充磁量）>0.05，说明充磁量对输出力矩没有显著影响
P（定位角度）>0.05，说明定位角度对输出力矩没有显著影响
P（充磁量：线圈匝数）>0.05，说明充磁量和线圈匝数对输出力矩没有显著影响
P（定位角度：线圈匝数）<0.05，说明定位角度和线圈匝数对输出力矩有显著影响
P（充磁量：线圈匝数：定位角度）>0.05，说明充磁量、定位角度和线圈匝数对输出力矩没有显著影响