search for: epxressed

m0<-epxression((4*theta1*theta2-theta3^2)/(2*x*theta3^2)-0.5*theta1*x) params<-all.vars(m0) this reads all the params from m0 so theta1,2 and 3 correct? params<-params[-which(params=="x")] checks which params are multiplied by x? np<-length(params) for(i in 1:6){ esp<-get(sprintf("m%d",i-1))

...2.1 1 2 2.2 1.1 0.7 2.7 Gene2 1.5 1.4 1.7 2.2 1.3 1.2 Gene3 1.4 1.7 1.8 2.7 1.6 1.5 Gene4 2.2 2.4 2.1 2.3 2.1 1.9 Gene5 2.6 3.4 2.1 1.3 2.6 2.9 Objective: find genes who are differentially epxressed. -- View this message in context: http://www.nabble.com/Statistical-Questions%3A-finding-differentially-expressed-genes-tp15873163p15873163.html Sent from the R help mailing list archive at Nabble.com.

...2 1.1 0.7 2.7 >Gene2 1.5 1.4 1.7 2.2 1.3 1.2 >Gene3 1.4 1.7 1.8 2.7 1.6 1.5 >Gene4 2.2 2.4 2.1 2.3 2.1 1.9 >Gene5 2.6 3.4 2.1 1.3 2.6 2.9 >Objective: find genes who are differentially epxressed. I'm not sure what you are asking, but to find whether one of your genes is significantly expressed is relatively straightforward in R or excel, and you have already outlined the procedure yourself. Have you tried to perform a paired t test or log transform in either software yet, and if so...