Advertisement
Guest User

Untitled

a guest
Dec 12th, 2017
84
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.29 KB | None | 0 0
  1. library(xlsx)
  2. mydata <- read.xlsx("C:/Users/Nicolai/Downloads/Data_Car_accidents.xlsx", 1)
  3.  
  4. head(mydata)
  5. table(mydata$Accident)
  6.  
  7. accidentTrue <- mydata[which(mydata$Accident == 'YES'),]
  8. accidentFalse <- mydata[which(mydata$Accident == 'NO'),]
  9. accidentTrueTraining <- sample(1:nrow(accidentTrue), 0.7*nrow(accidentTrue))
  10. accidentFalseTraining <- sample(1:nrow(accidentFalse), 0.7*nrow(accidentTrue))
  11.  
  12. training_true <- accidentTrue[accidentTrueTraining,]
  13. training_false <- accidentFalse[accidentFalseTraining,]
  14. trainingData <- rbind(training_true, training_false)
  15.  
  16. test_true <- accidentTrue[-accidentTrueTraining,]
  17. test_false <- accidentFalse[-accidentFalseTraining,]
  18. testData <- rbind(test_true, test_false)
  19.  
  20. setRepositories()
  21.  
  22. library(smbinning)
  23. factor_vars <- c("Gender", "Socioeconomic_status")
  24. continuous_vars <- c("Age", "BAC_â..")
  25.  
  26. adjustedMod <- glm(mydata$Accident ~ mydata$Gender + mydata$Age + mydata$Socioeconomic_status + mydata$BAC_â.., data=trainingData, family=binomial(link="logit"))
  27. adjustedMod
  28. nonadjustedMod <- glm(mydata$Accident ~ mydata$BAC_â.., data=trainingData, family=binomial(link="logit"))
  29.  
  30.  
  31. summary(adjustedMod)
  32. summary(nonadjustedMod)
  33.  
  34. predicted <- plogis(predict(logitMod, testData, type="response"))
  35. predicted
  36.  
  37. Concordance(mydata$Accident, predicted)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement