Background
Background: Histologic grade in breast cancer provides clinically important prognostic information. Researchers examined whether histologic grade was associated with gene expression profiles of breast cancers and whether such profiles could be used to improve histologic grading. In this tutorial we will assess the association between histologic grade and the expression of the KPNA2 gene that is known to be associated with poor BC prognosis. The patients, however, do not only differ in the histologic grade, but also on their lymph node status. The lymph nodes were not affected (0) or chirugically removed (1).
Data analysis
Import KPNA2 data in R
kpna2 <- read_tsv("https://raw.githubusercontent.com/GTPB/PSLS20/master/data/kpna2.txt")
kpna2
# A tibble: 24 x 3
grade node gene
<dbl> <dbl> <dbl>
1 3 1 368.
2 3 1 590.
3 1 1 347.
4 1 1 258.
5 1 0 154.
6 3 0 644.
7 3 1 818.
8 1 1 329.
9 3 0 746.
10 3 0 380.
# ... with 14 more rows
Transform the variable grade and node to a factor
Data exploration
Histologic grade and lymph node status can be associated with the kpna2 gene expression. Moreover, it is also possible that the differential expression associated with histological grade is different in patients that have unaffected lymph nodes and patients for which the lymph nodes had to be removed.
What does the plot suggests
Model
- Specify the model, interpret each model parameter in the model.
- Formulate the relevant research questions
- Translate them in a null and alternative hyptheses.
- Perform the hypothesis tests and calculate the confidence intervals.
Formulate your conclusions in terms of the research question
LS0tCnRpdGxlOiAiMy4gQnJlYXN0Y2FuY2VyIEdlbmUgRXhwcmVzc2lvbiBTdHVkeTogS1BOQTIgZ2VuZSIKYXV0aG9yOiAiTGlldmVuIENsZW1lbnQiCmRhdGU6ICJzdGF0T21pY3MsIEdoZW50IFVuaXZlcnNpdHkgKGh0dHBzOi8vc3RhdG9taWNzLmdpdGh1Yi5pbykiCm91dHB1dDoKICAgIGh0bWxfZG9jdW1lbnQ6CiAgICAgIGNvZGVfZG93bmxvYWQ6IHRydWUgICAgCiAgICAgIHRoZW1lOiBjb3NtbwogICAgICB0b2M6IHRydWUKICAgICAgdG9jX2Zsb2F0OiB0cnVlCiAgICAgIGhpZ2hsaWdodDogdGFuZ28KICAgICAgbnVtYmVyX3NlY3Rpb25zOiB0cnVlCi0tLQoKYGBge3Igc2V0dXAsIGluY2x1ZGU9RkFMU0V9CmtuaXRyOjpvcHRzX2NodW5rJHNldChpbmNsdWRlID0gVFJVRSwgY29tbWVudCA9IE5BLCBlY2hvID0gVFJVRSwKICAgICAgICAgICAgICAgICAgICAgIG1lc3NhZ2UgPSBGQUxTRSwgd2FybmluZyA9IEZBTFNFKQpsaWJyYXJ5KHRpZHl2ZXJzZSkKYGBgCgojQmFja2dyb3VuZApCYWNrZ3JvdW5kOiBIaXN0b2xvZ2ljIGdyYWRlIGluIGJyZWFzdCBjYW5jZXIgcHJvdmlkZXMgY2xpbmljYWxseSBpbXBvcnRhbnQgcHJvZ25vc3RpYyBpbmZvcm1hdGlvbi4gUmVzZWFyY2hlcnMgZXhhbWluZWQgd2hldGhlciBoaXN0b2xvZ2ljIGdyYWRlIHdhcyBhc3NvY2lhdGVkIHdpdGggZ2VuZSBleHByZXNzaW9uIHByb2ZpbGVzIG9mIGJyZWFzdCBjYW5jZXJzIGFuZCB3aGV0aGVyIHN1Y2ggcHJvZmlsZXMgY291bGQgYmUgdXNlZCB0byBpbXByb3ZlIGhpc3RvbG9naWMgZ3JhZGluZy4gSW4gdGhpcyB0dXRvcmlhbCB3ZSB3aWxsIGFzc2VzcyB0aGUgYXNzb2NpYXRpb24gYmV0d2VlbiBoaXN0b2xvZ2ljIGdyYWRlIGFuZCB0aGUgZXhwcmVzc2lvbiBvZiB0aGUgS1BOQTIgZ2VuZSB0aGF0IGlzIGtub3duIHRvIGJlIGFzc29jaWF0ZWQgd2l0aCBwb29yIEJDIHByb2dub3Npcy4KVGhlIHBhdGllbnRzLCBob3dldmVyLCBkbyBub3Qgb25seSBkaWZmZXIgaW4gdGhlIGhpc3RvbG9naWMgZ3JhZGUsIGJ1dCBhbHNvIG9uIHRoZWlyIGx5bXBoIG5vZGUgc3RhdHVzLiAKVGhlIGx5bXBoIG5vZGVzIHdlcmUgbm90IGFmZmVjdGVkICgwKSBvciBjaGlydWdpY2FsbHkgcmVtb3ZlZCAoMSkuCgoKI0RhdGEgYW5hbHlzaXMKIyNJbXBvcnQgS1BOQTIgZGF0YSBpbiBSCmBgYHtyfQprcG5hMiA8LSByZWFkX3RzdigiaHR0cHM6Ly9yYXcuZ2l0aHVidXNlcmNvbnRlbnQuY29tL0dUUEIvUFNMUzIwL21hc3Rlci9kYXRhL2twbmEyLnR4dCIpCmtwbmEyCmBgYAoKIyNUcmFuc2Zvcm0gdGhlIHZhcmlhYmxlIGdyYWRlIGFuZCBub2RlIHRvIGEgZmFjdG9yCmBgYHtyfQoKYGBgCgojI0RhdGEgZXhwbG9yYXRpb24KSGlzdG9sb2dpYyBncmFkZSBhbmQgbHltcGggbm9kZSBzdGF0dXMgY2FuIGJlIGFzc29jaWF0ZWQgd2l0aCB0aGUga3BuYTIgZ2VuZSBleHByZXNzaW9uLiBNb3Jlb3ZlciwgaXQgaXMgYWxzbyBwb3NzaWJsZSB0aGF0IHRoZSBkaWZmZXJlbnRpYWwgZXhwcmVzc2lvbiBhc3NvY2lhdGVkIHdpdGggaGlzdG9sb2dpY2FsIGdyYWRlIGlzIGRpZmZlcmVudCBpbiBwYXRpZW50cyB0aGF0IGhhdmUgdW5hZmZlY3RlZCBseW1waCBub2RlcyBhbmQgcGF0aWVudHMgZm9yIHdoaWNoIHRoZSBseW1waCBub2RlcyBoYWQgdG8gYmUgcmVtb3ZlZC4gCgpgYGB7cn0KCmBgYAoKV2hhdCBkb2VzIHRoZSBwbG90IHN1Z2dlc3RzCgojIyBNb2RlbAoKMS4gU3BlY2lmeSB0aGUgbW9kZWwsIGludGVycHJldCBlYWNoIG1vZGVsIHBhcmFtZXRlciBpbiB0aGUgbW9kZWwuIAoyLiBGb3JtdWxhdGUgdGhlIHJlbGV2YW50IHJlc2VhcmNoIHF1ZXN0aW9ucyAKMy4gVHJhbnNsYXRlIHRoZW0gaW4gYSBudWxsIGFuZCBhbHRlcm5hdGl2ZSBoeXB0aGVzZXMuICAKNC4gUGVyZm9ybSB0aGUgaHlwb3RoZXNpcyB0ZXN0cyAgYW5kIGNhbGN1bGF0ZSB0aGUgY29uZmlkZW5jZSBpbnRlcnZhbHMuIAoKIyBGb3JtdWxhdGUgeW91ciBjb25jbHVzaW9ucyBpbiB0ZXJtcyBvZiB0aGUgcmVzZWFyY2ggcXVlc3Rpb24KCg==