R语言_驾驶员事故预测

# Title     : TODO
# Objective : TODO
# Created by: johnchen
# Created on: 2020/4/24/024
library(openxlsx)
#读取例子5.1数据
d5.1 = read.xlsx('mvstats5.xlsx','d5.1')
logit = glm(y~x1+x2+x3, family = binomial, data = d5.1)#logistic模型，family是二项分布
logit.step = step(logit)
#预测视力正常司机logistic回归结果
pre1 = predict(logit.step, data.frame(x1 = 1))
#预测视力正常司机发生事故概率
p1 = exp(pre1)/(1 + exp(pre1))
##预测视力有问题司机logistic回归结果
pre2 = predict(logit.step,data.frame(x1 = 0))
#预测视力有问题司机发生事故概率
p2 = exp(pre2)/(1 + exp(pre2))
c(p1,p2)#结果显示