Contingency analysis for the ‘Titanic’ data.
Let’s consider the titanic dataset which contains a complete list of
passengers and crew members on the RMS Titanic. It includes a variable
indicating whether a person did survive the sinking of the RMS Titanic
on April 15, 1912. A data frame contains 2456 observations on 14
variables.
The website http://www.encyclopedia-titanica.org/ offers detailed
information about passengers and crew members on the RMS Titanic.
According to the website 1317 passengers and 890 crew member were
aboard.
8 musicians and 9 employees of the shipyard company are listed as
passengers, but travelled with a free ticket, which is why they have NA
values in fare. In addition to that, fare is truely missing for a few
regular passengers.
Data Cleansing
## Status Disembarked.at Home.Country Age Year.of.Birth
## DE GRASSE, Mr J. <NA> Cherbourg NA NA
## EVANS, Miss <NA> Cherbourg NA NA
## MULLEN, <NA> Cherbourg NA NA
## WOTTON, Mr Henry Swaffin <NA> Cherbourg 54 1858
## BRAND, Mr <NA> Cherbourg NA NA
## FLETCHER, Miss N. <NA> Cherbourg NA NA
## Crew.or.Passenger. Gender Class...Department
## DE GRASSE, Mr J. Passenger Male 2nd Class
## EVANS, Miss Passenger Female 2nd Class
## MULLEN, Passenger Female 2nd Class
## WOTTON, Mr Henry Swaffin Passenger Male 1st Class
## BRAND, Mr Passenger Male 1st Class
## FLETCHER, Miss N. Passenger Female 1st Class
## Embarked Job Job.details
## DE GRASSE, Mr J. Southampton
## EVANS, Miss Southampton
## MULLEN, Southampton
## WOTTON, Mr Henry Swaffin Southampton Butcher Butcher's Shop Proprietor
## BRAND, Mr Southampton
## FLETCHER, Miss N. Southampton
## Ticket.Number Fare.Price Fare_GBP Fare_today
## DE GRASSE, Mr J. 761 P1 1.0 82.110
## EVANS, Miss 88 P1 1.0 82.110
## MULLEN, 404 P1 1.0 82.110
## WOTTON, Mr Henry Swaffin 86 P1 10s 1.5 123.165
## BRAND, Mr 8 P1 10s 1.5 123.165
## FLETCHER, Miss N. 405 P1 10s 1.5 123.165
## Profile.on.Encyclopedia.Titanica
## DE GRASSE, Mr J. http://www.encyclopedia-titanica.org/titanic-biography/j-de-grasse.html
## EVANS, Miss http://www.encyclopedia-titanica.org/titanic-biography/evans.html
## MULLEN, http://www.encyclopedia-titanica.org/titanic-biography/mullen.html
## WOTTON, Mr Henry Swaffin http://www.encyclopedia-titanica.org/titanic-cross-channel-passenger/henry-swaffin-wotton.html
## BRAND, Mr http://www.encyclopedia-titanica.org/titanic-biography/brand.html
## FLETCHER, Miss N. http://www.encyclopedia-titanica.org/titanic-biography/n-fletcher.html
Survival by Passenger Class
We will create a contingency table to analyze the relationship
between survival and passenger class.
##
## 1st Class 2nd Class 3rd Class Deck Crew Engineering Crew
## Victim 123 166 528 23 253
## Survivor 201 119 180 43 71
##
## Restaurant Staff Victualling Crew
## Victim 66 337
## Survivor 3 94
##
## Female Male
## Victim 130 1366
## Survivor 359 352
##
## Female Male
## Victim 130 1366
## Survivor 359 352

Statistics number
## [1] "Phi value"
## [1] 0.3387276
## [1] "Contingency Coeficient"
## [1] 0.3208222
## [1] "Cramer Value"
## [1] 0.3387276
## [1] "Tschuprow"
## [1] 0.2164277
Association between survival status and passenger class is moderate
based on coefficients.
## [1] "Phi value"
## [1] 0.4703662
## [1] "Contingency Coeficient"
## [1] 0.4256325
## [1] "Cramer Value"
## [1] 0.4703662
## [1] "Tschuprow"
## [1] 0.4703662
Association between survival status and gender is significantly
higher.
Chi-Square
##
## Pearson's Chi-squared test
##
## data: contingency_table
## X-squared = 253.22, df = 6, p-value < 2.2e-16
##
## Pearson's Chi-squared test with Yates' continuity correction
##
## data: contingency_table_gender
## X-squared = 485.87, df = 1, p-value < 2.2e-16
P values are small, so it also sugest significant correlation between
survival status and passenger class or gender.
The results indicate that there is a statistically significant
relationship between survival status and the passenger class or gender.
Moreover we can conclude that gender correlation to survival status is
slightly stronger that correlation between passenger class and survival
status.
Plots


LS0tDQp0aXRsZTogJ0Rlc2NyaXB0aXZlIFN0YXRpc3RpY3MnDQpzdWJ0aXRsZTogJ1RpdGFuaWMgQml2YXJpYXRlIEFuYWx5c2lzJw0KZGF0ZTogImByIFN5cy5EYXRlKClgIg0KYXV0aG9yOiAiVG9tYXN6IEtydWN6YWxhaywgU3RhbmlzxYJhdyBLcnp5c2ljYSwgSGF0aWNlIFRhdGxpIg0Kb3V0cHV0Og0KICBodG1sX2RvY3VtZW50OiANCiAgICB0aGVtZTogY2VydWxlYW4NCiAgICBoaWdobGlnaHQ6IHRleHRtYXRlDQogICAgZm9udHNpemU6IDEwcHQNCiAgICB0b2M6IHllcw0KICAgIGNvZGVfZG93bmxvYWQ6IHllcw0KICAgIHRvY19mbG9hdDoNCiAgICAgIGNvbGxhcHNlZDogbm8NCiAgICBkZl9wcmludDogZGVmYXVsdA0KICAgIHRvY19kZXB0aDogNQ0KZWRpdG9yX29wdGlvbnM6IA0KICBtYXJrZG93bjogDQogICAgd3JhcDogNzINCi0tLQ0KDQpgYGB7ciBzZXR1cCwJbWVzc2FnZSA9IEZBTFNFLAl3YXJuaW5nID0gRkFMU0UsCWluY2x1ZGUgPSBGQUxTRX0NCmxpYnJhcnkoZHBseXIpDQpsaWJyYXJ5KHRpZHl2ZXJzZSkNCmxpYnJhcnkoSFNBVVIzKQ0KbGlicmFyeShoYXZlbikNCmxpYnJhcnkoZ2dwbG90MikNCmxpYnJhcnkoZ3JpZEV4dHJhKQ0KbGlicmFyeShwcGNvcikgIyB0aGlzIHBhY2thZ2UgY29tcHV0ZXMgcGFydGlhbCBhbmQgc2VtaXBhcnRpYWwgY29ycmVsYXRpb25zLg0KbGlicmFyeShsdG0pICMgdGhpcyBwYWNrYWdlIGNvbXB1dGVzIHBvaW50LWJpc2VyaWFsIGNvcnJlbGF0aW9ucy4NCmxpYnJhcnkoZGV2dG9vbHMpIA0KI2luc3RhbGxfZ2l0aHViKCJtYXJraGVja21hbm4vcnlvdXJlYWR5IikgIyBwbGVhc2UgaW5zdGFsbCBwYWNrYWdlICJyeW91cmVhZHkiIGZyb20gZ2l0aHViISAodGhlbiAjIGl0KQ0KbGlicmFyeShyeW91cmVhZHkpICMgdGhpcyBwYWNrYWdlIGNvbXB1dGVzIG5vbmxpbmVhciAiZXRhIiBjb3JyZWxhdGlvbnMuDQpsaWJyYXJ5KEdHYWxseSkgIyB0aGlzIHBhY2thZ2UgY29tcHV0ZXMgY29ycmVsYXRpb24gbWF0cml4Lg0KbGlicmFyeShwc3ljaCkgIyB0aGlzIHBhY2thZ2UgY29tcHV0ZXMgcXVhbGl0YXRpdmUgY29ycmVsYXRpb25zLg0KbGlicmFyeShEZXNjVG9vbHMpICMgdGhpcyBwYWNrYWdlIGNvbXB1dGVzIHF1YWxpdGF0aXZlIGNvcnJlbGF0aW9ucy4NCmBgYA0KDQojIyBDb250aW5nZW5jeSBhbmFseXNpcyBmb3IgdGhlICdUaXRhbmljJyBkYXRhLg0KDQpMZXQncyBjb25zaWRlciB0aGUgdGl0YW5pYyBkYXRhc2V0IHdoaWNoIGNvbnRhaW5zIGEgY29tcGxldGUgbGlzdCBvZiBwYXNzZW5nZXJzIGFuZCBjcmV3IG1lbWJlcnMgb24gdGhlIFJNUyBUaXRhbmljLiBJdCBpbmNsdWRlcyBhIHZhcmlhYmxlIGluZGljYXRpbmcgd2hldGhlciBhIHBlcnNvbiBkaWQgc3Vydml2ZSB0aGUgc2lua2luZyBvZiB0aGUgUk1TIFRpdGFuaWMgb24gQXByaWwgMTUsIDE5MTIuDQpBIGRhdGEgZnJhbWUgY29udGFpbnMgMjQ1NiBvYnNlcnZhdGlvbnMgb24gMTQgdmFyaWFibGVzLg0KDQpgYGB7ciBsb2FkLWRhdGEyLCB3YXJuaW5nPVRSVUUsIGluY2x1ZGU9RkFMU0V9DQpkb3dubG9hZC5maWxlKCJodHRwczovL2dpdGh1Yi5jb20va2ZsaXNpa293c2tpL2RzL2Jsb2IvbWFzdGVyL3RpdGFuaWMuY3N2P3Jhdz10cnVlIiwgZGVzdGZpbGUgPSJ0aXRhbmljLmNzdiIsbW9kZT0id2IiKQ0KdGl0YW5pYyA8LSByZWFkLmNzdigidGl0YW5pYy5jc3YiLHJvdy5uYW1lcz0xLHNlcD0iOyIpDQpgYGANCg0KVGhlIHdlYnNpdGUgaHR0cDovL3d3dy5lbmN5Y2xvcGVkaWEtdGl0YW5pY2Eub3JnLyBvZmZlcnMgZGV0YWlsZWQgaW5mb3JtYXRpb24gYWJvdXQgcGFzc2VuZ2VycyBhbmQgY3JldyBtZW1iZXJzIG9uIHRoZSBSTVMgVGl0YW5pYy4gQWNjb3JkaW5nIHRvIHRoZSB3ZWJzaXRlIDEzMTcgcGFzc2VuZ2VycyBhbmQgODkwIGNyZXcgbWVtYmVyIHdlcmUgYWJvYXJkLg0KDQo4IG11c2ljaWFucyBhbmQgOSBlbXBsb3llZXMgb2YgdGhlIHNoaXB5YXJkIGNvbXBhbnkgYXJlIGxpc3RlZCBhcyBwYXNzZW5nZXJzLCBidXQgdHJhdmVsbGVkIHdpdGggYSBmcmVlIHRpY2tldCwgd2hpY2ggaXMgd2h5IHRoZXkgaGF2ZSBOQSB2YWx1ZXMgaW4gZmFyZS4gSW4gYWRkaXRpb24gdG8gdGhhdCwgZmFyZSBpcyB0cnVlbHkgbWlzc2luZyBmb3IgYSBmZXcgcmVndWxhciBwYXNzZW5nZXJzLiANCg0KIyMgRGF0YSBDbGVhbnNpbmcNCmBgYHtyIGVjaG89RkFMU0UsIHdhcm5pbmc9VFJVRX0NCg0KIyBSZW1vdmUgcm93cyB3aXRoIE5BIHZhbHVlcyBpbiB0aGUgJ1N0YXR1cycgY29sdW1uDQp0aXRhbmljIDwtIHRpdGFuaWMgJT4lDQogIGZpbHRlcighaXMubmEoU3RhdHVzKSkNCg0KIyBDb252ZXJ0ICdTdGF0dXMnIHRvIGZhY3RvciB3aXRoIGxldmVscyAnVmljdGltJyBhbmQgJ1N1cnZpdm9yJw0KdGl0YW5pYyRTdGF0dXMgPC0gZmFjdG9yKHRpdGFuaWMkU3RhdHVzLCBsZXZlbHMgPSBjKCJWaWN0aW0iLCAiU3Vydml2b3IiKSkNCg0KIyBDb252ZXJ0ICdHZW5kZXInIHRvIGZhY3Rvcg0KdGl0YW5pYyRHZW5kZXIgPC0gZmFjdG9yKHRpdGFuaWMkR2VuZGVyKQ0KDQpoZWFkKHRpdGFuaWMpDQoNCmBgYA0KDQojIyBTdXJ2aXZhbCBieSBQYXNzZW5nZXIgQ2xhc3MNCldlIHdpbGwgY3JlYXRlIGEgY29udGluZ2VuY3kgdGFibGUgdG8gYW5hbHl6ZSB0aGUgcmVsYXRpb25zaGlwIGJldHdlZW4gc3Vydml2YWwgYW5kIHBhc3NlbmdlciBjbGFzcy4NCg0KYGBgIHtyIGVjaG89RkFMU0UsIHdhcm5pbmc9VFJVRX0NCg0KIyBDcmVhdGUgYSBjb250aW5nZW5jeSB0YWJsZQ0KY29udGluZ2VuY3lfdGFibGUgPC0gdGFibGUodGl0YW5pYyRTdGF0dXMsIHRpdGFuaWMkQ2xhc3MuLi5EZXBhcnRtZW50KQ0KY29udGluZ2VuY3lfdGFibGUNCg0KIyBDcmVhdGUgYSBjb250aW5nZW5jeSB0YWJsZSBmb3IgU3RhdHVzIGFuZCBHZW5kZXINCmNvbnRpbmdlbmN5X3RhYmxlX2dlbmRlciA8LSB0YWJsZSh0aXRhbmljJFN0YXR1cywgdGl0YW5pYyRHZW5kZXIpDQpjb250aW5nZW5jeV90YWJsZV9nZW5kZXINCg0KIyBDcmVhdGUgYSB2aXN1YWwgdGFibGUNCnk9YygxMzAsMzU5LDEzNjYsMzUyKQ0KZGltKHkpPWMoMiwyKQ0KY29udGluZ2VuY3lfdGFibGVfZ2VuZGVyX3ZpczwtYXMudGFibGUoeSkNCmRpbW5hbWVzKGNvbnRpbmdlbmN5X3RhYmxlX2dlbmRlcl92aXMpPWxpc3QoR2VuZGVyPWMoJ0ZlbWFsZScsJ01hbGUnKSxTdGF0dXM9YygnVmljdGltJywnU3Vydml2b3InKSkNCmNvbnRpbmdlbmN5X3RhYmxlX2dlbmRlcg0KZm91cmZvbGRwbG90KGNvbnRpbmdlbmN5X3RhYmxlX2dlbmRlcl92aXMpDQpgYGANCg0KIyMgU3RhdGlzdGljcyBudW1iZXINCmBgYHtyIGVjaG8gPSBGQUxTRSwgd2FybmluZz1UUlVFfQ0KcHJpbnQoIlBoaSB2YWx1ZSIpDQpQaGkoY29udGluZ2VuY3lfdGFibGUpDQoNCnByaW50KCJDb250aW5nZW5jeSBDb2VmaWNpZW50IikNCkNvbnRDb2VmKGNvbnRpbmdlbmN5X3RhYmxlKQ0KDQpwcmludCgiQ3JhbWVyIFZhbHVlIikNCkNyYW1lclYoY29udGluZ2VuY3lfdGFibGUpDQoNCnByaW50KCJUc2NodXByb3ciKQ0KVHNjaHVwcm93VChjb250aW5nZW5jeV90YWJsZSkNCmBgYA0KQXNzb2NpYXRpb24gYmV0d2VlbiBzdXJ2aXZhbCBzdGF0dXMgYW5kIHBhc3NlbmdlciBjbGFzcyBpcyBtb2RlcmF0ZSBiYXNlZCBvbiBjb2VmZmljaWVudHMuDQoNCmBgYHtyIGVjaG8gPSBGQUxTRSwgd2FybmluZz1UUlVFfQ0KcHJpbnQoIlBoaSB2YWx1ZSIpDQpQaGkoY29udGluZ2VuY3lfdGFibGVfZ2VuZGVyKQ0KDQpwcmludCgiQ29udGluZ2VuY3kgQ29lZmljaWVudCIpDQpDb250Q29lZihjb250aW5nZW5jeV90YWJsZV9nZW5kZXIpDQoNCnByaW50KCJDcmFtZXIgVmFsdWUiKQ0KQ3JhbWVyVihjb250aW5nZW5jeV90YWJsZV9nZW5kZXIpDQoNCnByaW50KCJUc2NodXByb3ciKQ0KVHNjaHVwcm93VChjb250aW5nZW5jeV90YWJsZV9nZW5kZXIpDQpgYGANCkFzc29jaWF0aW9uIGJldHdlZW4gc3Vydml2YWwgc3RhdHVzIGFuZCBnZW5kZXIgaXMgc2lnbmlmaWNhbnRseSBoaWdoZXIuDQoNCiMjIENoaS1TcXVhcmUgDQoNCmBgYCB7ciAgZWNobz1GQUxTRSwgd2FybmluZz1UUlVFfQ0KY2hpX3NxdWFyZV90ZXN0IDwtIGNoaXNxLnRlc3QoY29udGluZ2VuY3lfdGFibGUpDQpjaGlfc3F1YXJlX3Rlc3QNCg0KDQojIFBlcmZvcm0gYSBDaGktU3F1YXJlIFRlc3Qgb2YgSW5kZXBlbmRlbmNlIGZvciBTdGF0dXMgYW5kIEdlbmRlcg0KY2hpX3NxdWFyZV90ZXN0X2dlbmRlciA8LSBjaGlzcS50ZXN0KGNvbnRpbmdlbmN5X3RhYmxlX2dlbmRlcikNCmNoaV9zcXVhcmVfdGVzdF9nZW5kZXINCmBgYA0KUCB2YWx1ZXMgYXJlIHNtYWxsLCBzbyBpdCBhbHNvIHN1Z2VzdCBzaWduaWZpY2FudCBjb3JyZWxhdGlvbiBiZXR3ZWVuIHN1cnZpdmFsIHN0YXR1cyBhbmQgcGFzc2VuZ2VyIGNsYXNzIG9yIGdlbmRlci4NCg0KVGhlIHJlc3VsdHMgaW5kaWNhdGUgdGhhdCB0aGVyZSBpcyBhIHN0YXRpc3RpY2FsbHkgc2lnbmlmaWNhbnQgcmVsYXRpb25zaGlwIGJldHdlZW4gc3Vydml2YWwgc3RhdHVzIGFuZCB0aGUgcGFzc2VuZ2VyIGNsYXNzIG9yIGdlbmRlci4gTW9yZW92ZXIgd2UgY2FuIGNvbmNsdWRlIHRoYXQgZ2VuZGVyIGNvcnJlbGF0aW9uIHRvIHN1cnZpdmFsIHN0YXR1cyBpcyBzbGlnaHRseSBzdHJvbmdlciB0aGF0IGNvcnJlbGF0aW9uIGJldHdlZW4gcGFzc2VuZ2VyIGNsYXNzIGFuZCBzdXJ2aXZhbCBzdGF0dXMuDQoNCiMjIFBsb3RzDQoNCmBgYHtyIGVjaG8gPSBGQUxTRSwgd2FybmluZz1GQUxTRX0NCiMgUGxvdCB0aGUgY29udGluZ2VuY3kgdGFibGUNCmNvbnRpbmdlbmN5X2RmIDwtIGFzLmRhdGEuZnJhbWUoY29udGluZ2VuY3lfdGFibGUpDQoNCmdncGxvdChjb250aW5nZW5jeV9kZiwgYWVzKHggPSBWYXIxLCB5ID0gRnJlcSwgZmlsbCA9IFZhcjIpKSArDQogIGdlb21fYmFyKHN0YXQgPSAiaWRlbnRpdHkiLCBwb3NpdGlvbiA9ICJkb2RnZSIpICsNCiAgbGFicyh0aXRsZSA9ICJTdXJ2aXZhbCBDb3VudHMgYnkgUGFzc2VuZ2VyIENsYXNzIiwNCiAgICAgICB4ID0gIlN1cnZpdmFsIFN0YXR1cyIsDQogICAgICAgeSA9ICJDb3VudCIsDQogICAgICAgZmlsbCA9ICJQYXNzZW5nZXIgQ2xhc3MiKSArDQogIHRoZW1lX21pbmltYWwoKQ0KDQoNCiMgUGxvdCB0aGUgY29udGluZ2VuY3kgdGFibGUgZm9yIFN0YXR1cyBhbmQgR2VuZGVyDQpjb250aW5nZW5jeV9kZl9nZW5kZXIgPC0gYXMuZGF0YS5mcmFtZShjb250aW5nZW5jeV90YWJsZV9nZW5kZXIpDQoNCmdncGxvdChjb250aW5nZW5jeV9kZl9nZW5kZXIsIGFlcyh4ID0gVmFyMiwgeSA9IEZyZXEsIGZpbGwgPSBWYXIxKSkgKw0KICBnZW9tX2JhcihzdGF0ID0gImlkZW50aXR5IiwgcG9zaXRpb24gPSAiZG9kZ2UiKSArDQogIGxhYnModGl0bGUgPSAiU3Vydml2YWwgQ291bnRzIGJ5IEdlbmRlciIsDQogICAgICAgeCA9ICJHZW5kZXIiLA0KICAgICAgIHkgPSAiQ291bnQiLA0KICAgICAgIGZpbGwgPSAiU3RhdHVzIikgKw0KICB0aGVtZV9taW5pbWFsKCkNCmBgYA==