2014-04-02 33 views
1

我的當前線性模型是:fit<-lm(ES~Area+Anear+Dist+DistSC+Elevation)- [R線性模型允許二次和一階,但不高於

我已被要求通過進一步此:

使用五個說明適合ES的線性模型變量和 包括最多二次項和一階交互(即允許 區域^ 2和區域*高程,但不允許區域^ 3或 區域*高程*距離)。

從我的研究中,我可以做+I(Area^2)+(Area*Elevation)但這會形成一個巨大的名單。

假設我正確地理解了這個問題,我會加5個平方項和10個項,共計20個項。或者我不需要所有這些?

這真的是最有效的方法嗎?

編輯:

請注意,我在執行逐步迴歸爲空模型和之後的完整模型規劃。使用poly時,我似乎遇到了麻煩。

+1

看看'?poly' – Dason

+0

我做到了,感謝加入平方項。我仍然不確定如何使用它。我如何定義x以允許它給出五個變量平方和時間彼此? – JRSR

回答

2

?formula進一步您的教育:

fit<-lm(ES~ (Area+Anear+Dist+DistSC+Elevation)^2) 

那些不會被平方項,而是部分的你被要求提供什麼...所有的雙向交互(和主要影響)。公式「數學」與常規使用權力不同。要的方式,允許適當的統計解釋使用poly

fit<-lm(ES~ (Area+Anear+Dist+DistSC+Elevation)^2 + 
      poly(Area,2) +poly(Anear,2)+ poly(Dist,2)+ poly(DistSC,2)+ poly(Elevation,2)) 
+2

'poly(Area,Anear,Dist,DistSC,Elevation,degree = 2)'處理所有感興趣的術語。 – Dason

+0

感謝您的評論。我假設我的輸出「poly(Area,Anear,Dist,DistSC,Elevation,degree = 2)1.0.0.0.0」表示面積,而「poly(Area,Anear,Dist,DistSC,Elevation,degree = 2)0.0.0.2 .0「表示DistSC^2。 – JRSR

+0

我的逐步迴歸似乎不喜歡使用聚,而如果我只是擁有所有條款似乎運行良好,我可以解決這個問題嗎? – JRSR