[BioC] Survival analysis (hazard ratio calculation)

[Natasha Sahgal]

Hi Luffy,

I recently had to generate HR's and KM plots.

I used the survival package - which is what is used for KM plots but it can also get you Hazard ratios, look at the function: summary.coxph() or coxph then summary of that.


HTH,
Natasha




-----Original Message-----
From: bioconductor-bounces at r-project.org [mailto:bioconductor-bounces at r-project.org] On Behalf Of Farscape2012
Sent: 09 July 2012 12:08
To: bioconductor at r-project.org
Subject: [BioC] Survival analysis (hazard ratio calculation)

Hi,

I have been trying to compute hazard ratio between groups. I failed to find similar examples.

I have data like this:
      status time group
1    LIVING   14     1
2    LIVING   44     5
3  DECEASED   32     4
4    LIVING    8     7
5    LIVING   19     4
6    LIVING   22     4
7    LIVING   63     6
8    LIVING   44     3
9    LIVING   25     3
10   LIVING   60     3
11 DECEASED   25     4
12   LIVING   62     2
13   LIVING   43     6
14   LIVING   42     6
15   LIVING   44     2
16   LIVING   21     4
17   LIVING   26     4
18   LIVING   25     4
19   LIVING   24     3
20   LIVING   16     1

I have successfully plotted kaplan meier curve for each group. Now I want to compute hazard ratio between the groups.  The result could be like a symmetric matrix.  Each entry Eij shows the hazard ratio between group i and group j.

How can I do this? Please help me.

Br,
Luffy

_______________________________________________
Bioconductor mailing list
Bioconductor at r-project.org
https://stat.ethz.ch/mailman/listinfo/bioconductor
Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor