On the third day of the “Practical Statistics for the Life Sciences (2021)” course, we will have exercises on linear regression, based on different datasets:
The puromycin dataset
Data on the velocity of an enzymatic reaction were obtained by Treloar (1974).
The number of counts per minute of radioactive product from the reaction was measured as a function of substrate concentration in parts per million (ppm) and from these counts the initial rate (or velocity) of the reaction was calculated (counts/min/min). The experiment was conducted once with the enzyme treated with Puromycin, and once with the enzyme untreated.
The goal is to assess if there is an association between the substrate concentration and rate for the treated enzyme.
The poison dataset:
In this experiment, 96 fish (dojofish, goldfish and zebrafish) were placed separately in a tank with two litres of water and a certain dose (in mg) of a certain poison EI-43,064. The resistance of the fish a against the poison was measured as the amount of minutes the fish survived upon adding the poison (Surv_time, in minutes). Additionally, the weight of each fish was measured.
In this tutorial, we will study the association between dose and survival time by using a linear regression model.
Armpit example
Smelly armpits are not caused by sweat, itself. The smell is caused by specific micro-organisms belonging to the group of Corynebacterium spp. that metabolise sweat. Another group of abundant bacteria are the Staphylococcus spp., which do not metabolise sweat in smelly compounds.
The CMET-groep at Ghent University does research on transplanting the armpit microbiome to help people with smelly armpits. To test the effect of transplanting the microbiome, they conducted an experiment on two groups of volunteers: one group was treated with a placebo, while the other had a microbiome transplant. The goal was to assess whether a microbial transplant can affect the relative abundance Staphylococcus spp. in the microbiome of the armpit.
In the previous tutorial, we analysed the armpit data using a two sample t-test. Add a section to the script file for the armpit example where you model the data using a linear model. Compare the output of the linear model to the results you optained with the two sample t-test and formulate your conclusion based on the output of the linear model.
LS0tCnRpdGxlOiAiRXhlcmNpc2VzIG9uIGNoYXB0ZXIgNjogU2ltcGxlIGxpbmVhciByZWdyZXNzaW9uIgphdXRob3I6ICJMaWV2ZW4gQ2xlbWVudCBhbmQgSmVyb2VuIEdpbGlzIgpkYXRlOiAic3RhdE9taWNzLCBHaGVudCBVbml2ZXJzaXR5IChodHRwczovL3N0YXRvbWljcy5naXRodWIuaW8pIgpvdXRwdXQ6CiAgICBodG1sX2RvY3VtZW50OgogICAgICBjb2RlX2Rvd25sb2FkOiB0cnVlICAgIAogICAgICB0aGVtZTogY29zbW8KICAgICAgdG9jOiB0cnVlCiAgICAgIHRvY19mbG9hdDogdHJ1ZQogICAgICBoaWdobGlnaHQ6IHRhbmdvCiAgICAgIG51bWJlcl9zZWN0aW9uczogdHJ1ZQotLS0KCk9uIHRoZSB0aGlyZCBkYXkgb2YgdGhlICJQcmFjdGljYWwgU3RhdGlzdGljcyBmb3IgdGhlIExpZmUgU2NpZW5jZXMgKDIwMjEpIiAKY291cnNlLCB3ZSB3aWxsIGhhdmUgZXhlcmNpc2VzIG9uIGxpbmVhciByZWdyZXNzaW9uLCBiYXNlZCBvbiBkaWZmZXJlbnQgCmRhdGFzZXRzOgoKIyBUaGUgcHVyb215Y2luIGRhdGFzZXQKCkRhdGEgb24gdGhlIHZlbG9jaXR5IG9mIGFuIGVuenltYXRpYyByZWFjdGlvbiB3ZXJlIG9idGFpbmVkIGJ5IFRyZWxvYXIgKDE5NzQpLiAgClRoZSBudW1iZXIgb2YgY291bnRzIHBlciBtaW51dGUgb2YgcmFkaW9hY3RpdmUgcHJvZHVjdCBmcm9tIHRoZSByZWFjdGlvbiB3YXMgCm1lYXN1cmVkIGFzIGEgZnVuY3Rpb24gb2Ygc3Vic3RyYXRlIGNvbmNlbnRyYXRpb24gaW4gcGFydHMgcGVyIG1pbGxpb24gKHBwbSkgCmFuZCBmcm9tIHRoZXNlIGNvdW50cyB0aGUgaW5pdGlhbCByYXRlIChvciB2ZWxvY2l0eSkgb2YgdGhlIHJlYWN0aW9uIHdhcyAKY2FsY3VsYXRlZCAoY291bnRzL21pbi9taW4pLiBUaGUgZXhwZXJpbWVudCB3YXMgY29uZHVjdGVkIG9uY2Ugd2l0aCB0aGUgZW56eW1lIAp0cmVhdGVkIHdpdGggUHVyb215Y2luLCBhbmQgb25jZSB3aXRoIHRoZSBlbnp5bWUgdW50cmVhdGVkLgoKVGhlIGdvYWwgaXMgdG8gYXNzZXNzIGlmIHRoZXJlIGlzIGFuIGFzc29jaWF0aW9uIGJldHdlZW4gdGhlIHN1YnN0cmF0ZSAKY29uY2VudHJhdGlvbiBhbmQgcmF0ZSBmb3IgdGhlICoqdHJlYXRlZCoqIGVuenltZS4gCgotIEV4ZXJjaXNlOiBbRXhlcmNpc2UxXSguLzA2XzNfcHVyb215Y2luLmh0bWwpCi0gRGF0YSBwYXRoOiBOb3QgcmVxdWlyZWQKLSBTb2x1dGlvbjogW1NvbHV0aW9uMV0oLi8wNl8zX3B1cm9teWNpbl9zb2wuaHRtbCkKCgojIFRoZSBwb2lzb24gZGF0YXNldDoKCkluIHRoaXMgZXhwZXJpbWVudCwgOTYgZmlzaCAoZG9qb2Zpc2gsIGdvbGRmaXNoIGFuZCB6ZWJyYWZpc2gpCndlcmUgcGxhY2VkIHNlcGFyYXRlbHkgaW4gYSB0YW5rIHdpdGggdHdvIGxpdHJlcyBvZiB3YXRlciBhbmQKYSBjZXJ0YWluIGRvc2UgKGluIG1nKSBvZiBhIGNlcnRhaW4gcG9pc29uIEVJLTQzLDA2NC4gVGhlIHJlc2lzdGFuY2UKb2YgdGhlIGZpc2ggYSBhZ2FpbnN0IHRoZSBwb2lzb24gd2FzIG1lYXN1cmVkIGFzIHRoZSBhbW91bnQgb2YKbWludXRlcyB0aGUgZmlzaCBzdXJ2aXZlZCB1cG9uIGFkZGluZyB0aGUgcG9pc29uIChTdXJ2X3RpbWUsIGluCm1pbnV0ZXMpLiBBZGRpdGlvbmFsbHksIHRoZSB3ZWlnaHQgb2YgZWFjaCBmaXNoIHdhcyBtZWFzdXJlZC4KCkluIHRoaXMgdHV0b3JpYWwsIHdlIHdpbGwgc3R1ZHkgdGhlIGFzc29jaWF0aW9uIGJldHdlZW4gZG9zZSBhbmQgc3Vydml2YWwgdGltZQpieSB1c2luZyBhIGxpbmVhciByZWdyZXNzaW9uIG1vZGVsLgoKLSBFeGVyY2lzZTogW0V4ZXJjaXNlMl0oLi8wNl8xX3BvaXNvbi5odG1sKQotIERhdGEgcGF0aDogaHR0cHM6Ly9yYXcuZ2l0aHVidXNlcmNvbnRlbnQuY29tL3N0YXRPbWljcy9QU0xTMjEvZGF0YS9wb2lzb24uY3N2Ci0gU29sdXRpb246IFtTb2x1dGlvbjJdKC4vMDZfMV9wb2lzb25fc29sLmh0bWwpCgotLS0KCiMgQXJtcGl0IGV4YW1wbGUgCgpTbWVsbHkgYXJtcGl0cyBhcmUgbm90IGNhdXNlZCBieSBzd2VhdCwgaXRzZWxmLiBUaGUgc21lbGwgaXMgY2F1c2VkIGJ5IHNwZWNpZmljIAptaWNyby1vcmdhbmlzbXMgYmVsb25naW5nIHRvIHRoZSBncm91cCBvZiBDb3J5bmViYWN0ZXJpdW0gc3BwLiB0aGF0IG1ldGFib2xpc2UgCnN3ZWF0LiBBbm90aGVyIGdyb3VwIG9mIGFidW5kYW50IGJhY3RlcmlhIGFyZSB0aGUgU3RhcGh5bG9jb2NjdXMgc3BwLiwgd2hpY2ggZG8Kbm90IG1ldGFib2xpc2Ugc3dlYXQgaW4gc21lbGx5IGNvbXBvdW5kcy4KClRoZSBDTUVULWdyb2VwIGF0IEdoZW50IFVuaXZlcnNpdHkgZG9lcyByZXNlYXJjaCBvbiB0cmFuc3BsYW50aW5nIHRoZSBhcm1waXQgCm1pY3JvYmlvbWUgdG8gaGVscCBwZW9wbGUgd2l0aCBzbWVsbHkgYXJtcGl0cy4gVG8gdGVzdCB0aGUgZWZmZWN0IG9mIAp0cmFuc3BsYW50aW5nIHRoZSBtaWNyb2Jpb21lLCB0aGV5IGNvbmR1Y3RlZCBhbiBleHBlcmltZW50IG9uIHR3byBncm91cHMgCm9mIHZvbHVudGVlcnM6IG9uZSBncm91cCB3YXMgdHJlYXRlZCB3aXRoIGEgcGxhY2Vibywgd2hpbGUgdGhlIG90aGVyIGhhZCBhIAptaWNyb2Jpb21lIHRyYW5zcGxhbnQuIFRoZSBnb2FsIHdhcyB0byBhc3Nlc3Mgd2hldGhlciBhIG1pY3JvYmlhbCB0cmFuc3BsYW50IGNhbiBhZmZlY3QgdGhlIHJlbGF0aXZlIGFidW5kYW5jZSBTdGFwaHlsb2NvY2N1cyBzcHAuIGluIHRoZSBtaWNyb2Jpb21lIG9mIHRoZSBhcm1waXQuCgpJbiB0aGUgcHJldmlvdXMgdHV0b3JpYWwsIHdlIGFuYWx5c2VkIHRoZSBhcm1waXQgZGF0YSB1c2luZyBhIHR3byBzYW1wbGUgdC10ZXN0LgpBZGQgYSBzZWN0aW9uIHRvIHRoZSBzY3JpcHQgZmlsZSBmb3IgdGhlIGFybXBpdCBleGFtcGxlIHdoZXJlIHlvdSBtb2RlbCB0aGUgZGF0YSB1c2luZyBhIGxpbmVhciBtb2RlbC4gQ29tcGFyZSB0aGUgb3V0cHV0IG9mIHRoZSBsaW5lYXIgbW9kZWwgdG8gdGhlIHJlc3VsdHMgeW91IG9wdGFpbmVkIHdpdGggdGhlIHR3byBzYW1wbGUgdC10ZXN0IGFuZCBmb3JtdWxhdGUgeW91ciBjb25jbHVzaW9uIGJhc2VkIG9uIHRoZSBvdXRwdXQgb2YgdGhlIGxpbmVhciBtb2RlbC4gIAoKLSBEYXRhIHBhdGg6IGh0dHBzOi8vcmF3LmdpdGh1YnVzZXJjb250ZW50LmNvbS9zdGF0T21pY3MvUFNMUzIxL2RhdGEvYXJtcGl0LmNzdiIKLSBFeGVyY2lzZTogW0V4ZXJjaXNlM10oLi8wNl8zX2FybXBpdC5odG1sKQoK