Entering edit mode
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]]