建立多分類logistic模型的R語言代碼

setwd("C:/Users/lenovo/Desktop/迴歸12月1日成果") #設定當前的工作目錄

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

cere

#使用multinom做多類別logistic迴歸

library(nnet)

mult.cere<-multinom(y~x1+x2+x3,data=cere)

summary(mult.cere)

mult.cere1<-update(mult.cere,~.-1)#做係數的顯著性檢驗

mult.cere2<-update(mult.cere,~.-x1)

mult.cere3<-update(mult.cere,~.-x2)

mult.cere4<-update(mult.cere,~.-x3)

anova(mult.cere,mult.cere1)

anova(mult.cere,mult.cere2)

anova(mult.cere,mult.cere3)

anova(mult.cere,mult.cere4)

step.cere<-step(mult.cere)  #逐步迴歸選元

summary(step.cere)

exp(coef(step.cere))#用以解釋模型

cere.pred<-predict(step.cere) #做預測

cere.pred

table(cere$y,cere.pred)

n<-table(cere$y,cere.pred);n  #預測正確的百分比

Category<-levels(cere$y)

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

rbind(Category,Percantage)

cere.pred2<-predict(step.cere,type="p")

cere.pred2