?

Log in

No account? Create an account

Previous Entry | Next Entry

Вчера обновил на своих компьютерах R. Это такая очень удобная программа для разных статистических вычислений; в последнее время она стала одним из моих основных инструментов, вместе с пакетом knitr для автоматической генерации отчетов. Я использовал этот пакет для написания этой заметки.

Среди прочего в дистрибутив R входит больше ста разных наборов данных для тестирования, от количества авиапассажиров за каждый месяц с 1949 по 1960 год до роста и веса американских женщин в семидесятые. Один из этих наборов — статистика выживаемости на "Титанике". Это четырехмерный массив с четырьмя категориями: класс, пол, возраст, выжил человек или нет. Из него легко узнать, например, были ли больше шансы выжить у взрослых или детей:

apply(Titanic, c("Age","Survived"),sum)
##        Survived
## Age       No Yes
##   Child   52  57
##   Adult 1438 654

Мы видим, что спаслось больше половины детей, но только около трети взрослых:

prop.table(apply(Titanic, c("Age","Survived"),sum),1)
##        Survived
## Age            No       Yes
##   Child 0.4770642 0.5229358
##   Adult 0.6873805 0.3126195

Женщин тоже спасали в первую очередь:

apply(Titanic, c("Sex","Survived"),sum)
##         Survived
## Sex        No Yes
##   Male   1364 367
##   Female  126 344
prop.table(apply(Titanic, c("Sex","Survived"),sum),1)
##         Survived
## Sex             No       Yes
##   Male   0.7879838 0.2120162
##   Female 0.2680851 0.7319149

Однако самое интересное — это влияние на выживаемость социального статуса. На "Титанике" были пассажиры первого, второго и третьего класса, а также экипаж. Вот статистика по ним:

apply(Titanic, c("Class","Survived"),sum)
##       Survived
## Class   No Yes
##   1st  122 203
##   2nd  167 118
##   3rd  528 178
##   Crew 673 212
prop.table(apply(Titanic, c("Class","Survived"),sum),1)
##       Survived
## Class         No       Yes
##   1st  0.3753846 0.6246154
##   2nd  0.5859649 0.4140351
##   3rd  0.7478754 0.2521246
##   Crew 0.7604520 0.2395480

Видно, как по мере уменьшения статуса падают и шансы выжить. Это можно изобразить на графике:

mosaicplot(prop.table(apply(Titanic, c("Class","Survived"),sum),1),
color=TRUE, main="Survival and Class")
plot of chunksurvival.and.class

Интересно также посмотреть на выживание детей:

apply(Titanic[,,"Child",],
c("Class","Survived"),sum)
##       Survived
## Class  No Yes
##   1st   0   6
##   2nd   0  24
##   3rd  52  27
##   Crew  0   0
prop.table(apply(Titanic[,,"Child",],
c("Class","Survived"),sum),1)
##       Survived
## Class         No       Yes
##   1st  0.0000000 1.0000000
##   2nd  0.0000000 1.0000000
##   3rd  0.6582278 0.3417722
##   Crew       NaN       NaN

Спасли всех детей, плывших в первом и втором классе, но только треть детей из третьего (среди экипажа детей не было). Более того, выживаемость детей в третьем классе меньше общей выживаемости в первом и втором.

Посмотрим еще на выживаемость взрослых мужчин и женщин по классам:

prop.table(Titanic[,"Male","Adult",],1)
##       Survived
## Class         No        Yes
##   1st  0.6742857 0.32571429
##   2nd  0.9166667 0.08333333
##   3rd  0.8376623 0.16233766
##   Crew 0.7772622 0.22273782
prop.table(Titanic[,"Female","Adult",],1)
##       Survived
## Class          No       Yes
##   1st  0.02777778 0.9722222
##   2nd  0.13978495 0.8602151
##   3rd  0.53939394 0.4606061
##   Crew 0.13043478 0.8695652

Ну и общая картинка (в абсолютных числах):

mosaicplot(Titanic, color=TRUE)
plot of chunk titanic

Да, обещанный совет. Он, боюсь, тривиален, но все же.

Если вы плывете на "Титанике", постарайтесь быть вызывающим жалость пассажиром первого класса.

Comments

( 36 comments — Leave a comment )
evgknyaginin
Mar. 27th, 2015 05:38 pm (UTC)
Россия напоминает огромный "Титаник", на котором количество шлюпок ограничено.
sivka2006
Mar. 27th, 2015 06:06 pm (UTC)
Не Россия, а планета.

Проблема России в том, что она в 3-м классе.
(no subject) - evgknyaginin - Mar. 27th, 2015 06:08 pm (UTC) - Expand
(no subject) - sivka2006 - Mar. 27th, 2015 06:10 pm (UTC) - Expand
(no subject) - evgknyaginin - Mar. 27th, 2015 06:29 pm (UTC) - Expand
(no subject) - sivka2006 - Mar. 27th, 2015 06:52 pm (UTC) - Expand
(no subject) - abba65 - Mar. 27th, 2015 09:15 pm (UTC) - Expand
(no subject) - sivka2006 - Mar. 27th, 2015 09:28 pm (UTC) - Expand
(no subject) - ygam - Mar. 27th, 2015 06:41 pm (UTC) - Expand
suricat
Mar. 27th, 2015 05:42 pm (UTC)
Интересно, делал ли кто-нибудь подобные разборы для "Эстонии".
agasfer
Mar. 27th, 2015 05:50 pm (UTC)
Все эти цифры давно были в открытом доступе.
maina_woopka
Mar. 27th, 2015 05:58 pm (UTC)
По-моему, из этого следует, что не стоит путешествовать вторым классом - тут уж или всё, "гуляй, рванина, на все деньги", либо эконом.

Кстати, действительно, отчего такой провал по мужчинам во втором классе?
anna_frid
Mar. 27th, 2015 06:09 pm (UTC)
Ага, на kaggle.com тоже самый первый туториал - про "Титаник", только там еще данные более подробные. Правда, остальные переменные гораздо менее информативные.
kray_zemli
Mar. 27th, 2015 07:45 pm (UTC)
Статистика -- лженаука. Можно получить любой результат, который хочешь, если правильно считать.

Я когда считал для себя, получил совершенно другой результат.



1. Если смотреть выживаемость мужчин (я же мужчина, бабы пусть за себя сами думают), то минимум наблюдается именно во втором классе. Из 3-го класса спаслось вдвое больше мужчин.

2. Надо понимать, что спастись из 3-го класса было сложнее по фундаментальным причинам. Это видно как раз по числу спасшихся женщин и детей. Третий класс ехал в жопе трюма, туда вода потекла в первую очередь, плюс свет из-за потопа отрубился, где выход неясно, да ещё те самые решётки (какие-то требования иммиграционного законодательства того времени), преграждающие выход, открывать которые никто не торопился.

2а. Если отнормировать процент спасшихся мужчин по проценту спасшихся женщин-детей, считая, что шансы выбежать из трюма наверх были у всех примерно равны, а уже наверху шансы женщин были такие же как у остальных классов, то получается: из числа тех мужчин 3-го класса, кто уже сумел выбраться из трюма наверх,
спаслось столько же процентов, как и из 1-го класса.

Итого имеем результат как в экономике: именно средний класс самый уязвимый. Именно средний класс платит больше всего налогов -- богачи от них шифруются, с лимиты взять нечего. И при этом по ним неблагоприятные экономические факторы бьют в первую очередь.

Иными словами. Богачи купили места в шлюпках ха деньги. Лимита пробилась локтями, им не впервой. А вот средний класс, не имея ни того, ни другого оказался между молотом и наковальней.

Edited at 2015-03-27 07:47 pm (UTC)
grihanm
Mar. 27th, 2015 11:28 pm (UTC)
Хороший анализ
(no subject) - affidavid - Mar. 28th, 2015 11:19 am (UTC) - Expand
(no subject) - kray_zemli - Mar. 28th, 2015 11:43 am (UTC) - Expand
(no subject) - affidavid - Mar. 28th, 2015 11:54 am (UTC) - Expand
(no subject) - kray_zemli - Mar. 28th, 2015 12:05 pm (UTC) - Expand
(no subject) - affidavid - Mar. 28th, 2015 01:21 pm (UTC) - Expand
(no subject) - kray_zemli - Mar. 28th, 2015 01:34 pm (UTC) - Expand
selfmade
Mar. 27th, 2015 08:12 pm (UTC)
По-моему, более правильный совет - будь женщиной как минимум 2-го класса. Тем более сейчас это доступно.
thawing_wind
Mar. 27th, 2015 09:04 pm (UTC)
А женщин таки в команде выжило чуть больше, чем во втором классе, несмотря на то, что общая доля выживших в команде чуть меньше, чем в третьем классе.
p_govorun
Mar. 27th, 2015 09:18 pm (UTC)
Правильный совет такой: выучите, куда бежать в случае катастрофы. После этого можете брать билет третьего класса и быть там статистическим выбросом.

(Выше уже написали: каюты третьего класса были, естественно, в дебрях, и многие не добежали оттуда до шлюпок.)
ppk_ptichkin
Mar. 27th, 2015 09:26 pm (UTC)
>постарайтесь быть вызывающим жалость пассажиром первого класса.

Или хотя бы успешно притворяться таким в очереди на шлюпку. Из шлюпки уже не выкинут.
dmpogo
Mar. 27th, 2015 10:58 pm (UTC)
Ну с нижних палуб выбраться всегда сложнее
vladoff_revolt
Mar. 27th, 2015 11:03 pm (UTC)
у писателя Никонова в его книжке про феминизм: "А на самом деле бедные хуже богатых. Фарид Закария в своей книге «Будущее свободы» приводит удивительный пример. Когда затонул «Титаник», среди спасшихся пассажиров первого класса было практически 100% женщин и детей. Среди пассажиров второго класса процент женщин и детей достигал 80%. А среди сброда, ехавшего в третьем классе, выжившими оказались одни мужчины, то есть самые сильные. Они, отталкивая женщин и детей, заняли спасительные шлюпки...».

Я эту цитату обнаружил на фикшнбуке уже очень давно, там же было обсуждение в комментариях, уже выпиленное, в котором комментаторы Никонова опровергали.
egh0st
Mar. 27th, 2015 11:33 pm (UTC)
никонов в целом быдлО и гонит. но это всем известно кому интересно, конечно :)
(no subject) - scholar_vit - Mar. 28th, 2015 01:11 am (UTC) - Expand
irene221b
Mar. 27th, 2015 11:56 pm (UTC)
Это же классика инфографики:
eta_ta
Mar. 28th, 2015 02:40 pm (UTC)
ничего не поняла.
берём сиреневый 3й класс: ленточка идёт к точке МЕЖДУ м. и ж. , а оттуда - к точке МЕЖДУ да/нет.
Т.е. если ты гермафродит, путешествующий 3м классом, то...то ли спасёшься, то ли нет.
а что ещё за другая полоса вертикальная,шо оно значит?

какая белиберда

Edited at 2015-03-28 02:41 pm (UTC)
(no subject) - irene221b - Mar. 30th, 2015 10:02 pm (UTC) - Expand
(no subject) - eta_ta - Apr. 2nd, 2015 11:38 am (UTC) - Expand
robert_lee_kmd
Mar. 28th, 2015 04:03 am (UTC)
Лучше всего - надеть чистую одежду и запеть "Ближе, мой Бог, к Тебе, ближе к Тебе...", как завещали нам дирижер Уоллес Хартли и кочегар Джо Доусон.
(Deleted comment)
scholar_vit
Mar. 29th, 2015 01:40 am (UTC)
Пожалуйста
( 36 comments — Leave a comment )

Profile

knot
scholar_vit
scholar_vit

Latest Month

December 2017
S M T W T F S
     12
3456789
10111213141516
17181920212223
24252627282930
31      

Tags

Powered by LiveJournal.com
Designed by Paulina Bozek