警告消息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)