rm(list=ls(all=T))
options(digits=4, scipen=12)
library(magrittr)
Introduction
議題:用申請假釋者的屬性,預測他會不會違反假釋規定
學習重點:
- 設定隨機種子set.seed(),依比例分割資料
- 從邏輯式回歸的係數推算自變數的邊際效果
- 勝率和機率、勝率比和機率差 之間的換算
- 臨界機率對混淆矩陣(期望報酬)的影響payoff = matrix(c(0,-100,-10,-50),2,2); payoff
- AUC的實質意義payoff = matrix(c(100, -80, -20, 100),2,2); payoff
- 如何(從報酬矩陣)決定臨界機率
- 什麼是抽樣偏差,如何避免、如何修正
1 資料處理 Loading the Dataset
【1.1 讀進資料】How many parolees are contained in the dataset?
【1.2 底線機率】How many of the parolees in the dataset violated the terms of their parole?
2 整理資料框 Creating Our Prediction Model
【2.1 類別變數】Which variables in this dataset are unordered factors with at least three levels?
【2.2 資料框摘要】How does the output of summary() change for a factor variable as compared to a numerical variable?
3 資料分割 Splitting into a Training and Testing Set
【3.1 指定隨機種子、依比例分割資料】Roughly what proportion of parolees have been allocated to the training and testing sets?
【3.2 隨機種子的功用】Now, suppose you re-ran lines [1]-[5] of Problem 3.1. What would you expect? If you instead ONLY re-ran lines [3]-[5], what would you expect? If you instead called set.seed() with a different number and then re-ran lines [3]-[5] of Problem 3.1, what would you expect?
4 建立模型 Building a Logistic Regression Model
【4.1 顯著性】What variables are significant in this model?
【4.2 從回歸係數估計邊際效用】What can we say based on the coefficient of the multiple.offenses variable?
【4.3 從預測值估計勝率和機率】According to the model, what are the odds this individual is a violator? What is the probability this individual is a violator?
5 驗證模型 Evaluating the Model on the Testing Set
【5.1 從測試資料預測機率】What is the maximum predicted probability of a violation?
【5.2 從混淆矩陣計算敏感性、明確性、正確率】What is the model’s sensitivity, specificity, accuracy?
【5.3 底線機率】What is the accuracy of a simple model that predicts that every parolee is a non-violator?
【5.4 根據報償矩陣調整臨界機率】Which of the following most likely describes their preferences and best course of action?
【5.5 正確率 vs 辨識率】Which of the following is the most accurate assessment of the value of the logistic regression model with a cutoff 0.5 to a parole board, based on the model’s accuracy as compared to the simple baseline model?
【5.6 計算辨識率】Using the ROCR package, what is the AUC value for the model?
【5.7 辨識率的定義】Describe the meaning of AUC in this context.
6 抽樣偏差 Identifying Bias in Observational Data
【6.1 如何避免、診斷、修正抽樣偏差】How could we improve our dataset to best address selection bias?
LS0tDQp0aXRsZTogIkFTMy0yIFByZWRpY3RpbmcgUGFyb2xlIFZpb2xhdG9ycyINCmF1dGhvcjogIuWNk+mbjeeEtiBEOTk0MDEwMDAxIg0Kb3V0cHV0OiBodG1sX25vdGVib29rDQotLS0NCg0KYGBge3IgZWNobz1ULCBtZXNzYWdlPUYsIGNhY2hlPUYsIHdhcm5pbmc9Rn0NCnJtKGxpc3Q9bHMoYWxsPVQpKQ0Kb3B0aW9ucyhkaWdpdHM9NCwgc2NpcGVuPTEyKQ0KbGlicmFyeShtYWdyaXR0cikNCmBgYA0KDQotIC0gLQ0KDQojIyMgSW50cm9kdWN0aW9uDQoNCioq6K2w6aGM77ya55So55Sz6KuL5YGH6YeL6ICF55qE5bGs5oCn77yM6aCQ5ris5LuW5pyD5LiN5pyD6YGV5Y+N5YGH6YeL6KaP5a6aKioNCg0KKirlrbjnv5Lph43pu57vvJoqKg0KDQorIOioreWumumaqOapn+eoruWtkHNldC5zZWVkKCnvvIzkvp3mr5TkvovliIblibLos4fmlpkNCisg5b6e6YKP6Lyv5byP5Zue5q2455qE5L+C5pW45o6o566X6Ieq6K6K5pW455qE6YKK6Zqb5pWI5p6cDQorIOWLneeOh+WSjOapn+eOh+OAgeWLneeOh+avlOWSjOapn+eOh+W3riDkuYvplpPnmoTmj5vnrpcgDQorIOiHqOeVjOapn+eOh+Wwjea3t+a3huefqemZoyjmnJ/mnJvloLHphawp55qE5b2x6Z+/cGF5b2ZmID0gbWF0cml4KGMoMCwtMTAwLC0xMCwtNTApLDIsMik7IHBheW9mZg0KKyBBVUPnmoTlr6bos6rmhI/nvqlwYXlvZmYgPSBtYXRyaXgoYygxMDAsIC04MCwgLTIwLCAxMDApLDIsMik7IHBheW9mZg0KKyDlpoLkvZUo5b6e5aCx6YWs55+p6ZmjKeaxuuWumuiHqOeVjOapn+eOhyANCisg5LuA6bq85piv5oq95qij5YGP5beuLOWmguS9lemBv+WFjeOAgeWmguS9leS/ruatow0KDQo8YnI+DQoNCi0gLSAtDQoNCiMjIyMgMSDos4fmlpnomZXnkIYgTG9hZGluZyB0aGUgRGF0YXNldA0KDQrjgJAqKjEuMSDoroDpgLLos4fmlpkqKuOAkUhvdyBtYW55IHBhcm9sZWVzIGFyZSBjb250YWluZWQgaW4gdGhlIGRhdGFzZXQ/DQpgYGB7cn0NCg0KYGBgDQoNCuOAkCoqMS4yIOW6lee3muapn+eOhyoq44CRSG93IG1hbnkgb2YgdGhlIHBhcm9sZWVzIGluIHRoZSBkYXRhc2V0IHZpb2xhdGVkIHRoZSB0ZXJtcyBvZiB0aGVpciBwYXJvbGU/DQpgYGB7cn0NCg0KYGBgDQo8YnI+DQoNCi0gLSAtDQoNCiMjIyMgMiDmlbTnkIbos4fmlpnmoYYgQ3JlYXRpbmcgT3VyIFByZWRpY3Rpb24gTW9kZWwNCg0K44CQKioyLjEg6aGe5Yil6K6K5pW4KirjgJFXaGljaCB2YXJpYWJsZXMgaW4gdGhpcyBkYXRhc2V0IGFyZSB1bm9yZGVyZWQgZmFjdG9ycyB3aXRoIGF0IGxlYXN0IHRocmVlIGxldmVscz8gDQpgYGB7cn0NCg0KYGBgDQoNCuOAkCoqMi4yIOizh+aWmeahhuaRmOimgSoq44CRSG93IGRvZXMgdGhlIG91dHB1dCBvZiBgc3VtbWFyeSgpYCBjaGFuZ2UgZm9yIGEgZmFjdG9yIHZhcmlhYmxlIGFzIGNvbXBhcmVkIHRvIGEgbnVtZXJpY2FsIHZhcmlhYmxlPyANCmBgYHtyfQ0KDQpgYGANCjxicj4NCg0KLSAtIC0NCg0KIyMjIyAzIOizh+aWmeWIhuWJsiBTcGxpdHRpbmcgaW50byBhIFRyYWluaW5nIGFuZCBUZXN0aW5nIFNldA0KDQrjgJAqKjMuMSDmjIflrprpmqjmqZ/nqK7lrZDjgIHkvp3mr5TkvovliIblibLos4fmlpkqKuOAkVJvdWdobHkgd2hhdCBwcm9wb3J0aW9uIG9mIHBhcm9sZWVzIGhhdmUgYmVlbiBhbGxvY2F0ZWQgdG8gdGhlIHRyYWluaW5nIGFuZCB0ZXN0aW5nIHNldHM/DQpgYGB7cn0NCg0KYGBgDQoNCuOAkCoqMy4yIOmaqOapn+eoruWtkOeahOWKn+eUqCoq44CRTm93LCBzdXBwb3NlIHlvdSByZS1yYW4gbGluZXMgWzFdLVs1XSBvZiBQcm9ibGVtIDMuMS4gV2hhdCB3b3VsZCB5b3UgZXhwZWN0PyBJZiB5b3UgaW5zdGVhZCBPTkxZIHJlLXJhbiBsaW5lcyBbM10tWzVdLCB3aGF0IHdvdWxkIHlvdSBleHBlY3Q/IElmIHlvdSBpbnN0ZWFkIGNhbGxlZCBzZXQuc2VlZCgpIHdpdGggYSBkaWZmZXJlbnQgbnVtYmVyIGFuZCB0aGVuIHJlLXJhbiBsaW5lcyBbM10tWzVdIG9mIFByb2JsZW0gMy4xLCB3aGF0IHdvdWxkIHlvdSBleHBlY3Q/DQpgYGB7cn0NCg0KYGBgDQo8YnI+DQoNCi0gLSAtDQoNCiMjIyMgNCDlu7rnq4vmqKHlnosgQnVpbGRpbmcgYSBMb2dpc3RpYyBSZWdyZXNzaW9uIE1vZGVsDQoNCuOAkCoqNC4xIOmhr+iRl+aApyoq44CRV2hhdCB2YXJpYWJsZXMgYXJlIHNpZ25pZmljYW50IGluIHRoaXMgbW9kZWw/DQpgYGB7cn0NCg0KYGBgDQoNCuOAkCoqNC4yIOW+nuWbnuatuOS/guaVuOS8sOioiOmCiumam+aViOeUqCoq44CRV2hhdCBjYW4gd2Ugc2F5IGJhc2VkIG9uIHRoZSBjb2VmZmljaWVudCBvZiB0aGUgYG11bHRpcGxlLm9mZmVuc2VzYCB2YXJpYWJsZT8NCmBgYHtyfQ0KDQpgYGANCg0K44CQKio0LjMg5b6e6aCQ5ris5YC85Lyw6KiI5Yud546H5ZKM5qmf546HKirjgJFBY2NvcmRpbmcgdG8gdGhlIG1vZGVsLCB3aGF0IGFyZSB0aGUgb2RkcyB0aGlzIGluZGl2aWR1YWwgaXMgYSB2aW9sYXRvcj8gIFdoYXQgaXMgdGhlIHByb2JhYmlsaXR5IHRoaXMgaW5kaXZpZHVhbCBpcyBhIHZpb2xhdG9yPw0KYGBge3J9DQoNCmBgYA0KPGJyPg0KDQotIC0gLQ0KDQojIyMjIDUg6amX6K2J5qih5Z6LIEV2YWx1YXRpbmcgdGhlIE1vZGVsIG9uIHRoZSBUZXN0aW5nIFNldA0KDQrjgJAqKjUuMSDlvp7muKzoqabos4fmlpnpoJDmuKzmqZ/njocqKuOAkVdoYXQgaXMgdGhlIG1heGltdW0gcHJlZGljdGVkIHByb2JhYmlsaXR5IG9mIGEgdmlvbGF0aW9uPw0KYGBge3J9DQoNCmBgYA0KDQrjgJAqKjUuMiDlvp7mt7fmt4bnn6npmaPoqIjnrpfmlY/mhJ/mgKfjgIHmmI7norrmgKfjgIHmraPnorrnjocqKuOAkVdoYXQgaXMgdGhlIG1vZGVsJ3MgYHNlbnNpdGl2aXR5YCwgYHNwZWNpZmljaXR5YCwgYGFjY3VyYWN5YD8NCmBgYHtyfQ0KDQpgYGANCg0K44CQKio1LjMg5bqV57ea5qmf546HKirjgJFXaGF0IGlzIHRoZSBhY2N1cmFjeSBvZiBhIHNpbXBsZSBtb2RlbCB0aGF0IHByZWRpY3RzIHRoYXQgZXZlcnkgcGFyb2xlZSBpcyBhIG5vbi12aW9sYXRvcj8NCmBgYHtyfQ0KDQpgYGANCg0K44CQKio1LjQg5qC55pOa5aCx5YSf55+p6Zmj6Kq/5pW06Ieo55WM5qmf546HKirjgJFXaGljaCBvZiB0aGUgZm9sbG93aW5nIG1vc3QgbGlrZWx5IGRlc2NyaWJlcyB0aGVpciBwcmVmZXJlbmNlcyBhbmQgYmVzdCBjb3Vyc2Ugb2YgYWN0aW9uPw0KYGBge3J9DQoNCmBgYA0KDQrjgJAqKjUuNSDmraPnorrnjocgdnMg6L6o6K2Y546HKirjgJFXaGljaCBvZiB0aGUgZm9sbG93aW5nIGlzIHRoZSBtb3N0IGFjY3VyYXRlIGFzc2Vzc21lbnQgb2YgdGhlIHZhbHVlIG9mIHRoZSBsb2dpc3RpYyByZWdyZXNzaW9uIG1vZGVsIHdpdGggYSBjdXRvZmYgMC41IHRvIGEgcGFyb2xlIGJvYXJkLCBiYXNlZCBvbiB0aGUgbW9kZWwncyBhY2N1cmFjeSBhcyBjb21wYXJlZCB0byB0aGUgc2ltcGxlIGJhc2VsaW5lIG1vZGVsPw0KYGBge3J9DQoNCmBgYA0KDQrjgJAqKjUuNiDoqIjnrpfovqjorZjnjocqKuOAkVVzaW5nIHRoZSBgUk9DUmAgcGFja2FnZSwgd2hhdCBpcyB0aGUgQVVDIHZhbHVlIGZvciB0aGUgbW9kZWw/DQpgYGB7cn0NCg0KYGBgDQoNCuOAkCoqNS43IOi+qOitmOeOh+eahOWumue+qSoq44CRRGVzY3JpYmUgdGhlIG1lYW5pbmcgb2YgQVVDIGluIHRoaXMgY29udGV4dC4NCmBgYHtyfQ0KDQpgYGANCjxicj4NCg0KLSAtIC0NCg0KIyMjIyA2IOaKveaoo+WBj+W3riBJZGVudGlmeWluZyBCaWFzIGluIE9ic2VydmF0aW9uYWwgRGF0YQ0KDQrjgJAqKjYuMSDlpoLkvZXpgb/lhY3jgIHoqLrmlrfjgIHkv67mraPmir3mqKPlgY/lt64qKuOAkUhvdyBjb3VsZCB3ZSBpbXByb3ZlIG91ciBkYXRhc2V0IHRvIGJlc3QgYWRkcmVzcyBzZWxlY3Rpb24gYmlhcz8NCmBgYHtyfQ0KDQpgYGANCjxicj4NCg0KLSAtIC0NCg0KPGJyPjxicj48YnI+DQo=