Significance Analysis of Microarrays (SAM)
0
0
Entering edit mode
@eleni-christodoulou-2653
Last seen 6.2 years ago
Singapore
Dear list, I am trying to perform a significance analysis of a microarray experiment with survival data using the {samr} package. I have a matrix containing my data which has 17816 rows corresponding to genes, and 286 columns corresponding to samples. The name of this matrix is data.matrix2. Some of the first values of this matrix are: data.matrix2[1:3,1:5] GSM36777 GSM36778 GSM36779 GSM36780 GSM36781 [1,] 1.009274 1.0740659 1.048540 1.015946 1.022650 [2,] 1.007992 0.8768410 0.962442 1.111742 1.121150 [3,] 0.981853 0.9606492 1.024987 1.053302 1.063408 I also have the time in which each patient-sample is examined for relapse. This information is in vector y, which has length 286, and is declared in months. Indicatively: y[1:5] [1] 101 118 9 106 37 Finally, I have a variable censored, which is 1 if the patient has relapsed when examined at the examined time and 0 if not. Indicatively: censored[1:5] [1] 0 0 1 0 1 I am trying to perform the following sam analysis: d=list(data.matrix2,y,censored) samr.obj=samr(d,resp.type="Survival", nperms=20) When I am running the above commands I get the error: Error in check.format(y, resp.type = resp.type, censoring.status = censoring.status) : Error in input response data: response type Survival specified; error in censoring indicator In addition: Warning message: In is.na(x) : is.na() applied to non-(list or vector) of type 'NULL' I really cannot understand what is wrong with my code. Could anyone please help me with this? Thank you all, Eleni [[alternative HTML version deleted]]
Survival Survival • 2.1k views
ADD COMMENT

Login before adding your answer.

Traffic: 811 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6