logistics regression

# solve the problem of 13 chapter
rm(list = ls())
data <- read.csv("C:/users/mike1/desktop/data/RData/chapter13_oneExercise.csv", header = T, sep = ",")
dim(data)
names(data) # show the tabel 
data[1:10,]# show 10 lines

result <- glm(y~目标速度, data = data, family = binomial(link = probit))
summary(result)

result_fit <- result$fit 
result_fit

result_prd <- result$linear.predictors
result_prd

result_res <- residuals(result, c = "deviance")
result_res

qqnorm(result_res)
qqline(result_res)

plot(result_prd, result_res)
plot(result_fit, result_res)

 

 

 预测值与残差图：