similar to: Memory problem?

Hello list, I am trying to apply the paired t.test between diseased and not diseased patients to identify genes that are more expressed in the one situation under the other. In order to retrieve the genes that are more expressed in the positive disease state I do: p.values<-c() for(i in 1:length(Significant[,1])){ p.values[i]<-try(t.test(positive[i,],negative[i,],alternative

Hello R-community, It's been a week now that I am struggling with the implementation of a cox model in R. I have 80 cancer patients, so 80 time measurements and 80 relapse or no measurements (respective to censor, 1 if relapsed over the examined period, 0 if not). My microarray data contain around 18000 genes. So I have the expressions of 18000 genes in each of the 80 tumors (matrix

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

Hi all, I am trying to draw a Kaplan-Meier curve and I found online that Kaplan - Meier estimates are computed with a function called km in the event package. Is there an update for that because when I choose to download packages in R,. there is no package called event, even though I have selected all the repositories. Thanks in advance, Eleni [[alternative HTML version deleted]]

Hi list! I have recently tried to take my office work home, meaning that I tried to transfer my "... .RData" workspace from my PC on my laptop. The office PC runs on Windows XP and my laptop runs on Windows Vista. I have saved the workspace at the office PC and kept it in a usb drive. When I tried to open it on my laptop I got an error: "Fatal Error: Unable to restore saved data in

Hi all! I am new to R and I would like to ask you the following question:How can I substitute the NA values with 0 in a data frame? I cannot find a command to check if a value is NA... Thank you very much! Eleni

Hi all, I am trying to install the package "GEOquery" in unix. I have downloaded the standard version of R and this package is not contained in the default. I know that I can select repositories under windows but I don't know how to do it in unix. Does anyone have any idea on this? Thank you in advance, Eleni [[alternative HTML version deleted]]

Hello all, I have the following problem. I want to insert blank rows in an data frame so as to make it comparable with another data frame with more rows. How can I insert these rows? If this is not solvable, then, alternatively, is it easy to delete rows from a data frame? Thanks for helping, Eleni [[alternative HTML version deleted]]

Hello R community, I am processing raw Affymetrix CEL files and I am using the Michigan custom CDF library hgu133plus2hsentrezgprobe. I have been looking for documentation on the function that it contains...I am specifically interested in converting probe names to gene symbols. Does anybody know where I can find it? Thank a lot! Eleni [[alternative HTML version deleted]]

Hello everyone, I wanted to know if it was possible to change the color of certain leaves in a hclust object in order to make my graph more readable. I know I can color certain groups but I would like to enter a vector telling the plot function which leaves to color in which color. Thanks in advance, Eleni [[alternative HTML version deleted]]

Dear list, I am trying to build a cox model and then perform ROC analysis in order to retrieve some genes that are correlated with breast cancer. When I calculate the hazard score taking into account different numbers of genes and their coefficients ( I am trying to find the pest predictor number of genes), I retrieve from around 1 values (for few genes included ) to size of e+80 values (for many

Dear list, I am having a set of human oligo ids (H200006022 H200002025 H300004703 H200000442 H200005719 H300018350) which I want to map to Ensembl or RefSeq. I am sure R has a function to do that. I downloaded the {oligo} package and tried to use the probeNames function. Although the factor of ologo ids is an object (as the argument to probeNames should be) I retrieve the following error:

Hello list! I have a proble trying to perform a SAM analysis using the function samr from the samr package. I have put the option *center.arrays=TRUE *in order to scale all the experiments to median=0. I would like to retrieved the scaled data but it seems that samr does not return it...Does anyone have any idea on this? Thanks a lot!!! E. [[alternative HTML version deleted]]

Hello all! I am performingsome clustering analysis on microarray data using agnes{cluster} and I have created my own dissimilarity matrix according to a distance measure different from "euclidean" or "manhattan" etc. My question is, if I choose for example method="complete", how are the distances between the elements calculated? Are they taken form the dissimilarity

For any given pre-specified gene or short list of genes, yes the Cox model works fine. Two important caveats: 1. Remeber the rule of thumb for a Cox model of 20 events per variable (not n=20). Many microarray studies will have very marginal sample size. 2. If you are looking at many genes then a completely different strategy is required. There is a large and growing literature; I like Newton

Dear list, I have a couple of questions concerning ridge regression. I am using the lm.ridge(...) function in order to fit a model to my microarray data. Thus *model=lm.ridge(...)* I retrieve some coefficients and some scales for each gene. First of all, I would like to ask: the real coefficients of the model are not included in the first argument of the output but in the result of coef(model),

Dear list, I have been searching for a week to fit a simple linear model to my data. I have looked into the previous posts but I haven't found anything relevant to my problem. I guess it is something simple...I just cannot see it. I have the following data frame, named "data", which is a subset of a microarray experiment. The columns are the samples and the rows are the probes. I

Dear R users, I d like to assess the effect of "treatment" covariate on a disease relapse risk with the package cmprsk. However, the effect of this covariate on survival is time-dependent (assessed with cox.zph): no significant effect during the first year of follow-up, then after 1 year a favorable effect is observed on survival (step function might be the correct way to say that ?).

Hello list, I am trying to perform ROC analysis and count the AUC in order to validate my results. I use package ROCR. I would like to count the AUC not under the cutoff found by "performance" but to use another cutoff that I calculate. How could I change the following command in order to get what I want? perform=performance(pred,measure="auc",x.measure="cutoff"),

Dear friends: I'm interested to make a stacked plot of cumulative incidence. that's, the cuminc model is fitted [fit=cuminc(time, relapse)] and cumulative incidence is in place. I'd like to stack the cuminc plots (relapse of luekemia and death free from leukemia, for example) , then the constituent ratio of leukemia relapse and treatment related mortality is very clear. Can