Hi, Try: example1<- as.matrix(read.table("example1.txt",header=TRUE,stringsAsFactors=FALSE,sep="\t",row.names=1,check.names=FALSE)) example2 <- example1 example1New <- cbind(example1, `MAX without Restriction`=apply(example1[,1:4],1,max,na.rm=TRUE)) #or library(matrixStats) `MAX without Restriction`<- rowMaxs(example2[ ,colnames(example2)[!grepl("Restriction", colnames(example2))]],na.rm=TRUE) example2New <- cbind(example2, `MAX without Restriction`) identical(example1New,example2New) #[1] TRUE
A.K. Hello I've got a matrix of the following form BARN BCGE BCVN BEAN Restriction 1 Restriction 2 Restriction 3 Restriction 4 Restriction 5 alpha.1 0.000172449 7.87E-05 -0.003271044 0.000921609 9.28E-19 2.00E-05 -0.000608211 NA NA alpha.2 0.000896744 0.000579453 -0.000623357 0.001260358 -1.36E-19 -5.22E-05 NA NA NA alpha.3 0.000832748 0.00076229 0.002170229 0.001159895 3.09E-19 -7.86E-05 NA NA NA alpha.4 0.000920545 NA 0.001680433 0.000459149 -3.08E-19 -3.59E-05 NA NA NA alpha.5 0.001385238 0.000527484 0.000593311 0.000549358 7.72E-19 -6.99E-05 NA NA NA alpha.6 0.000644262 0.000305714 -0.00044514 0.000407448 -9.68E-20 -5.56E-05 NA NA NA alpha.7 -0.00022808 -0.00017047 0.000109545 0.000601197 0 3.50E-05 NA NA NA alpha.8 -1.16E-05 -0.000105657 0.001403036 0.00058719 3.88E-19 8.64E-06 NA NA NA alpha.9 0.000633559 -4.33E-05 0.000724611 0.000841646 -4.82E-20 -3.29E-05 NA NA NA (see also file) Now I'd like to calculate the maximal value of each row, but exclude Restriction 1-5. An add a new column to the given matrix with the name 'MAX without Restriction' and the max of each row. I tried max<-apply(example1,1,function(x)max(x)) new<-cbind(example,max) but it gave a strange output, also I couldn't manage to exclude the columns Restriction 1-5. Thank you for your help!! ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.