Вопрос, который я публикую здесь, тесно связан с другим вопросом, который я опубликовал два дня назад, об анализе старения gompertz.
Я пытаюсь построить объект выживания, см. Surv, в R. Надеюсь, он будет использован для выполнения анализа Гомперца, чтобы получить на выходе два значения (см. Исходный вопрос для получения дополнительных сведений).
У меня есть данные о выживаемости из эксперимента на мухах, в котором изучается скорость старения у разных генотипов. Мне доступны данные в нескольких макетах, поэтому выбор остается за вами, в зависимости от того, что лучше всего подходит для ответа.
Один фрейм данных (wide.df) выглядит так, где каждый генотип (Exp, из которых ~ 640) имеет строку, а дни идут последовательно по горизонтали от дня 4 до дня 98 с подсчетом новых смертей каждые два дня.
Exp Day4 Day6 Day8 Day10 Day12 Day14 ...
A 0 0 0 2 3 1 ...
Я делаю пример, используя это:
wide.df2<-data.frame("A",0,0,0,2,3,1,3,4,5,3,4,7,8,2,10,1,2)
colnames(wide.df2)<-c("Exp","Day4","Day6","Day8","Day10","Day12","Day14","Day16","Day18","Day20","Day22","Day24","Day26","Day28","Day30","Day32","Day34","Day36")
Другая версия похожа на эту, где каждый день имеет строку для каждого «опыта» и записывается количество смертей в этот день.
Exp Deaths Day
A 0 4
A 0 6
A 0 8
A 2 10
A 3 12
.. .. ..
Чтобы сделать этот пример:
df2<-data.frame(c("A","A","A","A","A","A","A","A","A","A","A","A","A","A","A","A","A"),c(0,0,0,2,3,1,3,4,5,3,4,7,8,2,10,1,2),c(4,6,8,10,12,14,16,18,20,22,24,26,28,30,32,34,36))
colnames(df2)<-c("Exp","Deaths","Day")
В каждом генотипе около 50 мух. В чем мне сейчас нужна помощь, так это в том, как перейти от одного из вышеупомянутых фреймов данных к работающему объекту выживания. Как выглядит этот объект? И как мне плавно добраться от вышеуказанного до объекта выживания?