[R] outer

Kaspar Pflugshaupt pflugshaupt at cirsium.ethz.ch
Thu Mar 15 18:50:30 CET 2001


On Thursday 15 March 2001 18:11, Olivier Martin wrote:

> Dear r-plus users,

r-plus? That's new! :-)

>
> i would like to use outer in the following case  outer(x,y,FUN="fun")
> i suppose that my function fun is of the following form:
> fun<-function(x,y)
> {
> if(y>x) return(x+y)
> if(y<=x) return(0)
> }
> My problem is that the command outer(x,y,FUN="fun") return me a
> null matrix instead of an upper triangular matrix.
> Is somebody have a solution ?

you would have to change the if() clause: x and y are your data, not indices.

Try this:

m <- outer(x,y,FUN="+")
m[lower.tri(m, diag=T)] <- 0

Cheers

Kaspar



-- 


Kaspar Pflugshaupt
Geobotanical Institute
ETH Zurich, Switzerland

http://www.geobot.umnw.ethz.ch
mailto:pflugshaupt at geobot.umnw.ethz.ch
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._



More information about the R-help mailing list