警告消息R中的有序逻辑回归

警告消息R中的有序逻辑回归,r,warnings,ordinal,R,Warnings,Ordinal,这是我的数据头(两个): 我想做逻辑回归,其中分数是响应变量。我想做四个独立的 在这里: 我所做的: library(MASS) y <- factor(both$betyg) mod.fit <- polr(y ~ prog + gender + age + season, data=both, Hess=TRUE) summary(mod.fit) 库(MASS) y因为你的结果是有序的,你可能会用ordinal做得更好,但可能需要检查比例优势假设。您所描述的模型与p

这是我的数据头(两个):

我想做逻辑回归,其中分数是响应变量。我想做四个独立的

在这里:

我所做的:

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

y因为你的结果是有序的,你可能会用ordinal做得更好,但可能需要检查比例优势假设。您所描述的模型与
polr
所做的差不多,尽管它们并不像您所说的那样独立。加州大学洛杉矶分校在这方面有很好的教程

至于确定哪种模型是最好的,当处理这些根本不同类型的模型时,我建议交叉验证。预测的准确性并不存在,任何伪R^2指标在不同模型的解释上都会有所不同


此外,由于这个问题更多地涉及统计信息而不是R编码/实现,我建议使用CrossValidated(stats StackExchange站点)。

这确实非常有用。当我使用R代码时收到警告消息。不幸的是,我不知道如何解释。应该是有用的。你说你想做逻辑回归,其中
grade
是顺序响应,但在你的模型代码中,你使用了一些称为
betyg
的未引用变量。如果我没弄错的话,那是瑞典语中“等级”的意思。为了保持一致性,我建议您选择一个名称。
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)