Thanks bert!!!it worked out perfectly well.thankyou onceagain,
Eliza

> Date: Tue, 29 Oct 2013 08:30:13 -0700
> Subject: Re: [R] maximum value replacement
> From: gunter.berton@gene.com
> To: smartpink111@yahoo.com; eliza_botto@hotmail.com
> CC: r-help@r-project.org
> 
> To Eliza: What if the max in a column is not unique?
> 
> Given the small size of A, the solution given by Arun seems completely
> adequate. However, I was wondering if it could be done without the
> R-level loop in sapply by taking advantage of pmax() . Of course it
> can. Here's code to illustrate how:
> 
> A<- matrix(sample(15),nr=3)
> A
> B <- matrix(0,nr=1,nc=5)
> m <- do.call(pmax,data.frame(t(A)))
> d <- dim(A)
> A[which(A==matrix(m,nr=d[1],nc=d[2],byrow=TRUE))] <- B
> A
> 
> Note that:
> 
> 1. This does not generalize to functions other than max or min, afaik.
> 2. I don't even know if it would be faster for large data, because the
> data.frame call may slow things down.
> 
> But it is fully vectorized (I think). So for illustration only, maybe...
> 
> Cheers,
> Bert
> 
> 
> 
> On Tue, Oct 29, 2013 at 7:27 AM, arun <smartpink111@yahoo.com> wrote:
> > Hi,
> > Try:
> >  sapply(seq_len(ncol(A)),function(i) {indx <- which(A[,i]%in% max(A[,i])); A[,i][indx] <- B[,i]; A[,i]})
> >
> > A.K.
> >
> >
> >
> >
> > On Tuesday, October 29, 2013 10:16 AM, eliza botto <eliza_botto@hotmail.com> wrote:
> > Dear Users,
> > I have two matrices, one with 12 rows and 124 columns(A) and the other with 1 row and 124 column(B). i want to replace the maximum value in all columns of A with each (single) column value of B.
> > How can i do it??
> > Thanks indeed in advance,
> >
> > Eliza
> >     [[alternative HTML version deleted]]
> >
> > ______________________________________________
> > 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.
> >
> >
> > ______________________________________________
> > 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.
> 
> 
> 
> -- 
> 
> Bert Gunter
> Genentech Nonclinical Biostatistics
> 
> (650) 467-7374
 		 	   		  
	[[alternative HTML version deleted]]

