R中的警告消息順序邏輯回歸

Question

這是我的數據頭（兩者）：

 season  gender age   prog     grade 
    fall    woman  old  FRIST       B
    fall    woman  old  FRIST       A 
    spring  woman  old  FRIST       E 
    spring    man  old  NMATK       C 
    spring  woman  old  NFYSK       A 
    fall    woman  old  FRIST       E 

我想進行邏輯回歸，其中成績是響應變量。 我要使其中四個是獨立的。

這里：

E/A+B+C+D=alpha_1+beta^x_1+beta^y_1+...

D+E/A+B+C=alpha_2+beta^x_2+beta^y_2+...

C+D+E/A+B=alpha_3+beta^x_3+beta^y_3+...

B+C+D+E/A=alpha_4+beta^x_4+beta^y_4+...

我做了什么：

    library(MASS)
y <- factor(both$betyg)
mod.fit <- polr(y ~ prog + gender + age + season, data=both, Hess=TRUE)
summary(mod.fit) 

然后我收到此消息：

警告消息：在polr中（y〜prog +性別+年齡+季節，數據=兩者，Hess = TRUE）：設計似乎是等級不足的，因此要降低一些系數

我知道這不是錯誤，只是警告。 我不知道如何解釋它，或者采取其他措施來避免此消息？

Answer 1

由於您的結果是有序的，因此使用序數可能會更好，但是可能要檢查比例幾率假設。 您所描述的模型幾乎是polr工作，盡管它們並不是您所說的獨立。 UCLA在這方面有一個很好的教程。

至於確定哪種模型最好，在處理根本不同類型的模型時，我建議進行交叉驗證。 預測准確性並不重要，並且任何偽R ^ 2指標在模型之間的解釋都將有所不同。

此外，由於此問題比R編碼/實現更關注統計信息，因此我建議使用CrossValidated（Status StackExchange網站）。