(аналогично графику).
https://www.nejm.org/na101/home/literatum/publisher/mms/journals/content/nejm/2013/nejm_2013.368.issue-13/nejmoa1301228/20130322/images/img_xlarge/nejmoa1301228_f2.jpeg
Я разобрался, как сделать похожий forest plot с текстом и т.п. Вопросы скорее по статистике:
1. Правильно ли я понимаю, что для каждой подгруппы (давайте для примера возьмем возраст и пол) нужно сделать сабсет и сравнить терапии?
У меня Overall модель для всех пациентов такого типа:
coxph(Surv(TIME, EVENT) ~ group, data = df)
Cоответственно, для каждой категории каждой подгруппы модель остается та же, но датасет сабсетится по-разному:
coxph(Surv(TIME, EVENT) ~ group, data = df[df["Sex"] == "M" ])
coxph(Surv(TIME, EVENT) ~ group, data = df[df["Sex"] == "F" ])и т.д.
Так можно получить все нужные результаты, кроме interaction p-value
Что за либа так отрисовывает красиво?
Обсуждают сегодня