具有R的二次模型。使用poly（..）和I（..）函數（R語言）

Remi.b

2013-08-01 22:20:29 UTC

view on stackexchange narkive permalink

什麼原因導致以下不同結果？

  var1 = c（0.04875,0.13725,0.28350,0.50975,0.77425,0.94700,0.05325,0.14050,0.29725,0.51525,0.79000,0.95400,0.04625,0.15250,0.29000,0.53300,0.79825,0.95225， 0.05025,0.14625,0.28800,0.52625,0.78200,0.95925,0.04700,0.14225,0.30325,0.53500,0.79325,0.95875,0.04775,0.13850,0.28675,0.54250,0.78300,0.95175,0.05150,0.12725,0.30175,0.54725,0.79475,0.96275,0.05375， 0.14100,0.30050,0.53275,0.78100,0.96175,0.05450,0.15300,0.29650,0.52850,0.80100,0.95675,0.05425,0.13975,0.30875,0.56025,0.80575,0.96100,0.05100,0.15350,0.31175,0.53300,0.78900,0.96000,0.04650,0.13525， 0.29600,0.53625,0.78475,0.96375,0.05375,0.13900,0.29600,0.53725,0.78700,0.95800,0.05075,0.14350,0.29225,0.54525,0.80275,0.95800,0.05050,0.13200,0.29850,0.52700,0.80525,0.96150,0.05150,0.14050,0.29450， 0.52375,0.79450,0.96375,0.05375,0.13525,0.30475,0.55250,0.79425,0.96025,0.04950,0.14500,0.29425,0.52250,0.78475,0.95650,0.05225,0.14425,0.29225,0.53150,0.80425,0.95375）var2 = c（1,2， 3,4,5,6,1,2,3,4,5,6,1,2,3,4,5,6,1,2,3,4,5,6,1,2,3， 4 ，5,6,1,2,3,4,5,6,1,2,3,4,5,6,1,2,3,4,5,6,1,2,3,4,5 ，6,1,2,3,4,5,6,1,2,3,4,5,6,1,2,3,4,5,6,1,2,3,4,5,6 ，1,2,3,4,5,6,1,2,3,4,5,6,1,2,3,4,5,6,1,2,3,4,5,6,1 ，2,3,4,5,6,1,2,3,4,5,6）var3 = c（4,4,4,4,4,4,4,6,6,6,6,6,6 ，8,8,8,8,8,8,10,10,10,10,10,10,12,12,12,12,12,12,12,14,14,14,14,14,14,16 ，16,16,16,16,16,16,18,18,18,18,18,18,20,20,20,20,20,20,20,22,22,22,22,22,22,22,24,24 ，24,24,24,24,26,26,26,26,26,26,28,28,28,28,28,28,28,30,30,30,30,30,30,32,32,32 ，32,32,32,34,34,34,34,34,34,36,36,36,36,36,36,38,38,38,38,38,38,40,40,40,40 ，40,40）summary（aov（var1〜as.factor（var2）* var3））summary（aov（var1〜as.factor（poly（var2,1））* poly（var3,1）））

我可能不太了解poly（..）和I（..）的工作方式。

（我的目的是創建不同程度（二次方，三次方等）的模型並比較它們的BIC（或AIC））。

> poly（1:10，1）1 [1，] -0.49543369 [2，] -0.38533732 [ 3，] -0.27524094 [4，] -0.16514456 [5，] -0.05504819 [6，] 0.05504819 [7，] 0.16514456 [8，] 0.27524094 [9，] 0.38533732 [10，] 0.49543369attr（，“度”）[ 1] 1attr（，“ coefs”）attr（，“ coefs”）$ alpha [1] 5.5attr（，“ coefs”）$ norm2 [1] 1.0 10.0 82.5attr（，“ class”）[1]“ poly” “ matrix”

> poly（1:10，2）1 2 [1，] -0.49543369 0.52223297 [2，] -0.38533732 0.17407766 [3，] -0.27524094 -0.08703883 [4，] -0.16514456 -0.26111648 [5，] -0.05504819 -0.34815531 [6，] 0.05504819 -0.34815531 [7，] 0.16514456 -0.26111648 [8，] 0.27524094 -0.08703883 [9，] 0.38533732 0.17407766 [10，] 0.49543369 0.522232 97attr（，“度”）[1] 1 2attr（，“ coefs”）attr（，“ coefs”）$ alpha [1] 5.5 5.5 attr（，“ coefs”）$ norm2 [1] 1.0 10.0 82.5 528.0attr（，“ class”）[1]“ poly”“ matrix”

> poly（1:10，2，raw = TRUE）的二次方的期望。1 2 [1，] 1 1 [2，] 2 4 [3，] 3 9 [4，] 4 16 [5，] 5 25 [6，] 6 36 [7，] 7 49 [8，] 8 64 [9 ，] 9 81 [10，] 10 100 attr（，“ degree”）[1] 1 2attr（，“ class”）[1]“ poly”“ matrix”