2010-06-10 71 views
3

以我previous post,我尋找相關比(η或η )例程在R.我被無人使用η線性度在GLM程序檢查的事實感到驚訝。非圖形線性估計

讓我們從一個簡單的例子開始:如何檢查雙變量相關的線性?單憑散點圖?

有幾種方法可以做到這一點,一種方法是比較線性和非線性模型,然後應用F檢驗來尋求它們之間的顯着差異。

最後,問題是:如何檢查線性,「非基本」的方式?

+2

一個更重要的問題是_why_你在做這個,爲什麼你不想使用圖形? – hadley 2010-06-10 13:52:26

+0

我想得到一些關於非圖形化程序的信息,只是出於好奇,但這並不意味着我不想使用圖形...相反,我總是看數據以獲得一個關於底層結構的印象(這聽起來有點不妥,但爲了這些目的,我使用了'GGobi','rggobi'和'ggplot2')。所以有... – aL3xa 2010-06-10 14:52:47

回答

5

答案就是你所說的(比較線性和非線性模型)。 例如

model1<-lm(yv~xv) 
model2<-lm(yv~xv+I(xv^2)) #Even if we restrict ourselves to the inclusion of a quadratic term, there are many curves we can describe, depending upon the signs of the linear and quadratic terms 

anova(model1,model2) 

Analysis of Variance Table 

Model 1: yv ~ xv 
Model 2: yv ~ xv + I(xv^2) 
    Res.Df RSS Df Sum of Sq  F Pr(>F) 
1  16 91.057        
2  15 68.143 1 22.915 5.0441 0.0402 * 

越複雜彎曲模式是在所述線性模型(P = 0.04),因此,在這種情況下,我們承認存在在數據曲率證據顯著改善。

+0

只是爲了澄清:這個程序是否與Mandel測試相同? – Beasterfield 2011-12-02 14:43:22

4

RESET(迴歸方程式規格誤差測試)是爲丟失迴歸器設計的,但它經常用於測試非線性。可以在LMTEST包中找到 - 在許多其他有用的測試中。這與你已經做的非常相似。或者,您可以設計一個遞歸殘差測試,以利用輸入非線性變量進行排序時它們可能全部爲正/負的事實。

+1

+1將其應用於我的測試數據:'resettest(yv〜xv,power = 2,type =「regressor」)''給出了'RESET = 5.0441,df1 = 1,df2 = 15,p-值= 0.0402' – 2010-06-10 13:29:28

+1

'LDdiag'軟件包中還有'ramsey'功能 – aL3xa 2010-06-10 15:23:02