建立有序多分类logistic模型的R语言代码

setwd("C:/Users/lenovo/Desktop/回归12月1日成果") #设定当前的工作目录

cred=read.table("shuju.txt",header=T)

cred  #读取数据

library(MASS)   #做序列多分类回归

polr.cred<-polr(factor(y)~factor(x1)+factor(x2)+x3,data=cred)  

summary(polr.cred)

drop1(polr.cred,test="Chi")  #系数显著性检验

pred.cred<-predict(polr.cred)#得到属于哪个类

#预测值与观测值的交叉列表看分类效果

n<-table(cred$y,pred.cred);n

Percantage<-c(n[1,1]/sum(n[1,]),n[2,2]/sum(n[2,]),n[3,3]/sum(n[3,]))

rbind(Category,Percantage)