This calculation divides by values centered around zero.  The only context that 
I can think of that would require such silliness is a homework problem, and 
this list has a no-homework policy. If not, then mentioning the theory you are 
applying might help someone point you at an existing function that achieves 
your goals while avoiding divide-by-zero errors. 

Since you also posted in HTML I gather that you have not read the Posting Guide 
mentioned below.  Avoiding HTML on this list is to your benefit,  since using 
it inevitably leads to others seeing a garbled version of what you sent. Please 
read the PG for more important guidance. 
-- 
Sent from my phone. Please excuse my brevity.

On December 16, 2015 4:18:56 PM PST, Matteo Richiardi 
<matteo.richia...@gmail.com> wrote:
>I have to evolve each element of a matrix W
>
>W <- matrix(0,2,3)
>
>according to some function which uses the indices of the matrix [i,j]
>as
>arguments:
>w.fun = function(i,j) {
>  return A[i]*B[j]/(C[i,j])
>}
>
>where
>A<-c(100,100)
>B<-c(200,200,200)
>C <- matrix( rnorm(6,mean=0,sd=1), 2, 3)
>
>How can I do it, without recurring to a loop? Also, in my application I
>need to pass the function another argument.
>
>Thanks a lot for your suggestions.
>Matteo
>
>       [[alternative HTML version deleted]]
>
>______________________________________________
>R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
>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.

        [[alternative HTML version deleted]]

______________________________________________
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
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.

Reply via email to