Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- library(xlsx)
- mydata <- read.xlsx("C:/Users/Nicolai/Downloads/Data_Car_accidents.xlsx", 1)
- head(mydata)
- table(mydata$Accident)
- accidentTrue <- mydata[which(mydata$Accident == 'YES'),]
- accidentFalse <- mydata[which(mydata$Accident == 'NO'),]
- accidentTrueTraining <- sample(1:nrow(accidentTrue), 0.7*nrow(accidentTrue))
- accidentFalseTraining <- sample(1:nrow(accidentFalse), 0.7*nrow(accidentTrue))
- training_true <- accidentTrue[accidentTrueTraining,]
- training_false <- accidentFalse[accidentFalseTraining,]
- trainingData <- rbind(training_true, training_false)
- test_true <- accidentTrue[-accidentTrueTraining,]
- test_false <- accidentFalse[-accidentFalseTraining,]
- testData <- rbind(test_true, test_false)
- setRepositories()
- library(smbinning)
- factor_vars <- c("Gender", "Socioeconomic_status")
- continuous_vars <- c("Age", "BAC_â..")
- adjustedMod <- glm(mydata$Accident ~ mydata$Gender + mydata$Age + mydata$Socioeconomic_status + mydata$BAC_â.., data=trainingData, family=binomial(link="logit"))
- adjustedMod
- nonadjustedMod <- glm(mydata$Accident ~ mydata$BAC_â.., data=trainingData, family=binomial(link="logit"))
- summary(adjustedMod)
- summary(nonadjustedMod)
- predicted <- plogis(predict(logitMod, testData, type="response"))
- predicted
- Concordance(mydata$Accident, predicted)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement