Я провожу анализ выживания, используя значения экспрессии гена (в пересчете на z-показатели) из определенного гена и времени выживания без событий для 249 пациентов. Я новичок в этом.
Вот мой код:
finaldata $ SurvObj <- with (finaldata, Surv (efst, status == 2)) </p>
finaldata $ SurvObj
res.cox1 <- coxph (SurvObj ~ Zscore, data = finaldata) </p>
сводка (res.cox1)
А вот и выходные данные:
coxph(formula = SurvObj ~ Zscore, data = finaldata)
n= 247, number of events= 140
(2 observations deleted due to missingness)
coef exp(coef) se(coef) z Pr(>|z|)
Zscore 0.5878 1.8001 0.2118 2.776 0.0055
---
Signif. codes: 0 0.001 0.01 0.05 ‘.’ 0.1 ‘ ’ 1
exp(coef) exp(-coef) lower .95 upper .95
Zscore 1.8 0.5555 1.189 2.726
Concordance= 0.561 (se = 0.025 )
Likelihood ratio test= 7.92 on 1 df, p=0.005
Wald test = 7.71 on 1 df, p=0.006
Score (logrank) test = 7.74 on 1 df, p=0.005
efst - это время выживания пациентов без событий, состояние - мертв ли они или живы (1 жив, 2 мертв), а Zscore является значением Z-Score для экспрессия рассматриваемого гена для каждого из пациентов (чтобы понять, является ли экспрессия относительно высокой или низкой).
Однако я нахожу крайне маловероятным, что я получаю такие p-значения, и это приводит меня думать, что я определенно сделал что-то неправильно, так как я впервые пробую это.
Как я могу найти любые ошибки, которые я сделал, как я могу исправить их, и как я должен интерпретировать эти данные?