For today’s lab, we’ll take a closer look at logistic regression. There’s a huge range of backgrounds in this class, so some of you are probably pretty familiar with logistic regression, and it may be totally new to others! So today, you have two options. The first is to work through another OpenIntro tutorial, and the second is to do a case study on your own. You only need to do one for the lab to be considered complete today.

Option 1: Tutorial

In this tutorial option, you’ll explore:

  • How GPA is related to medical school acceptance
  • Five-year survival of heart transplant patients

If you choose this option, you should include screenshots of the two graphs that best display the results of the model. This is a subjective decision!

Option 2: Case Study on Your Own

If you feel pretty confident with logistic regression (or are just sick of the tutorials) you can choose to extend the demo from class on your own! You should include at least two variables in your model.

Included here is the rough outline of steps you should follow. In your write-up, you should include (1) an EDA graph, (2) A Model Fit Summary Graph, and (3) a short (2-3 sentence) conclusion.

Download Data

# remotes::install_github("OpenIntroStat/openintro") 
library(openintro)
email
## # A tibble: 3,921 x 21
##    spam  to_multiple from     cc sent_email time                image attach
##    <fct> <fct>       <fct> <int> <fct>      <dttm>              <dbl>  <dbl>
##  1 0     0           1         0 0          2012-01-01 01:16:41     0      0
##  2 0     0           1         0 0          2012-01-01 02:03:59     0      0
##  3 0     0           1         0 0          2012-01-01 11:00:32     0      0
##  4 0     0           1         0 0          2012-01-01 04:09:49     0      0
##  5 0     0           1         0 0          2012-01-01 05:00:01     0      0
##  6 0     0           1         0 0          2012-01-01 05:04:46     0      0
##  7 0     1           1         0 1          2012-01-01 12:55:06     0      0
##  8 0     1           1         1 1          2012-01-01 13:45:21     1      1
##  9 0     0           1         0 0          2012-01-01 16:08:59     0      0
## 10 0     0           1         0 0          2012-01-01 13:12:00     0      0
## # … with 3,911 more rows, and 13 more variables: dollar <dbl>, winner <fct>,
## #   inherit <dbl>, viagra <dbl>, password <dbl>, num_char <dbl>,
## #   line_breaks <int>, format <fct>, re_subj <fct>, exclaim_subj <dbl>,
## #   urgent_subj <fct>, exclaim_mess <dbl>, number <fct>

EDA

Fit the model

Summarize the model fit

Draw conclusions

LS0tCnRpdGxlOiAiTGFiMDc6IENsYXNzaWZpY2F0aW9uIgpvdXRwdXQ6IAogIGh0bWxfZG9jdW1lbnQ6IAogICAgdG9jOiB0cnVlCiAgICB0b2NfZmxvYXQ6IHRydWUKICAgIGNvZGVfZG93bmxvYWQ6IHRydWUKICAgIGNzczogc3RhdDQxLWxhYi10aGVtZS5jc3MKLS0tCgpgYGB7ciBzZXR1cCwgaW5jbHVkZT1GQUxTRX0KbGlicmFyeSh0aWR5dmVyc2UpCmxpYnJhcnkodGlkeW1vZGVscykKa25pdHI6Om9wdHNfY2h1bmskc2V0KG1lc3NhZ2UgPSBGQUxTRSwgd2FybmluZyA9IEZBTFNFKQpgYGAKCkZvciB0b2RheSdzIGxhYiwgd2UnbGwgdGFrZSBhIGNsb3NlciBsb29rIGF0ICpsb2dpc3RpYyByZWdyZXNzaW9uKi4gVGhlcmUncyBhIGh1Z2UgcmFuZ2Ugb2YgYmFja2dyb3VuZHMgaW4gdGhpcyBjbGFzcywgc28gc29tZSBvZiB5b3UgYXJlIHByb2JhYmx5IHByZXR0eSBmYW1pbGlhciB3aXRoIGxvZ2lzdGljIHJlZ3Jlc3Npb24sIGFuZCBpdCBtYXkgYmUgdG90YWxseSBuZXcgdG8gb3RoZXJzISBTbyB0b2RheSwgeW91IGhhdmUgdHdvIG9wdGlvbnMuIFRoZSBmaXJzdCBpcyB0byB3b3JrIHRocm91Z2ggYW5vdGhlciBPcGVuSW50cm8gdHV0b3JpYWwsIGFuZCB0aGUgc2Vjb25kIGlzIHRvIGRvIGEgY2FzZSBzdHVkeSBvbiB5b3VyIG93bi4gWW91IG9ubHkgbmVlZCB0byBkbyAqKm9uZSoqIGZvciB0aGUgbGFiIHRvIGJlIGNvbnNpZGVyZWQgY29tcGxldGUgdG9kYXkuIAoKIyBPcHRpb24gMTogVHV0b3JpYWwgCgpJbiB0aGlzIFt0dXRvcmlhbCBvcHRpb25dKGh0dHBzOi8vb3BlbmludHJvLnNoaW55YXBwcy5pby9pbXMtMDQtbXVsdGl2YXJpYWJsZS1hbmQtbG9naXN0aWMtbW9kZWxzLTA0LyNzZWN0aW9uLWxvZ2lzdGljLXJlZ3Jlc3Npb24pLCB5b3UnbGwgZXhwbG9yZToKCiAgKyBIb3cgR1BBIGlzIHJlbGF0ZWQgdG8gbWVkaWNhbCBzY2hvb2wgYWNjZXB0YW5jZQogICsgRml2ZS15ZWFyIHN1cnZpdmFsIG9mIGhlYXJ0IHRyYW5zcGxhbnQgcGF0aWVudHMKCklmIHlvdSBjaG9vc2UgdGhpcyBvcHRpb24sIHlvdSBzaG91bGQgaW5jbHVkZSBzY3JlZW5zaG90cyBvZiB0aGUgdHdvIGdyYXBocyB0aGF0IGJlc3QgZGlzcGxheSB0aGUgcmVzdWx0cyBvZiB0aGUgbW9kZWwuIFRoaXMgaXMgYSBzdWJqZWN0aXZlIGRlY2lzaW9uISAKCiMgT3B0aW9uIDI6IENhc2UgU3R1ZHkgb24gWW91ciBPd24KCklmIHlvdSBmZWVsIHByZXR0eSBjb25maWRlbnQgd2l0aCBsb2dpc3RpYyByZWdyZXNzaW9uIChvciBhcmUganVzdCBzaWNrIG9mIHRoZSB0dXRvcmlhbHMpIHlvdSBjYW4gY2hvb3NlIHRvIGV4dGVuZCB0aGUgZGVtbyBmcm9tIGNsYXNzIG9uIHlvdXIgb3duISBZb3Ugc2hvdWxkIGluY2x1ZGUgYXQgbGVhc3QgdHdvIHZhcmlhYmxlcyBpbiB5b3VyIG1vZGVsLiAKCkluY2x1ZGVkIGhlcmUgaXMgdGhlIHJvdWdoIG91dGxpbmUgb2Ygc3RlcHMgeW91IHNob3VsZCBmb2xsb3cuIEluIHlvdXIgd3JpdGUtdXAsIHlvdSBzaG91bGQgaW5jbHVkZSAoMSkgYW4gRURBIGdyYXBoLCAoMikgQSBNb2RlbCBGaXQgU3VtbWFyeSBHcmFwaCwgYW5kICgzKSBhIHNob3J0ICgyLTMgc2VudGVuY2UpIGNvbmNsdXNpb24uIAoKIyMjIERvd25sb2FkIERhdGEKCmBgYHtyfQojIHJlbW90ZXM6Omluc3RhbGxfZ2l0aHViKCJPcGVuSW50cm9TdGF0L29wZW5pbnRybyIpIApsaWJyYXJ5KG9wZW5pbnRybykKZW1haWwKYGBgCgojIyMgRURBIAoKIyMjIEZpdCB0aGUgbW9kZWwgCgojIyMgU3VtbWFyaXplIHRoZSBtb2RlbCBmaXQKCiMjIyBEcmF3IGNvbmNsdXNpb25zCgo=