[R] comparing row by row in matrix

[Huntsinger, Reid]

I don't know the Russell-Rao coefficient but maybe this will help:

You can compute the number of times y and x are both 1 (that's what your
function f does) by counting 1s in y*x, i.e., sum(y*x), aka t(y)%*%x.
Ordinary matrix multiplication does this row-by-column, so if M and N are
two binary matrices, M%*%t(N) will have (i,j) entry equal to the count of 1s
in common to row i of M and row j of N.

Hope I've understood the problem correctly.

Reid Huntsinger

-----Original Message-----
From: Petra Steiner [mailto:steinep at uni-muenster.de]
Sent: Tuesday, March 26, 2002 2:08 PM
To: r-help at stat.math.ethz.ch
Subject: [R] comparing row by row in matrix


Hello,

and thanks for the two responses to my questions on binary matrixes, which
showed me that the functions I needed do not exist.

To get a distance matrix with the Russell-Rao-coefficient, I first have to
compare each row of a binary matrix with each row and count how many
elements are
a. equal
and
b. 1.

by
f <- function(x,y) (sum(x == 1& y == x))

Now how can I iterate this over a matrix without a loop? I think apply etc.
will not work in this case.

Thanks for any help.

Regards,
Petra

-
---------------------------------------------------
Petra Steiner
Arbeitsbereich Linguistik
Universitaet Muenster
Huefferstrasse 27
48149  Muenster

Tel: 0251 / 83 39442
petra at marley.uni-muenster.de
http://santana.uni-muenster.de/~petra/





-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.
-.-
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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._.
_._


------------------------------------------------------------------------------
Notice:  This e-mail message, together with any attachments, contains information of Merck & Co., Inc. (Whitehouse Station, New Jersey, USA) that may be confidential, proprietary copyrighted and/or legally privileged, and is intended solely for the use of the individual or entity named in this message.  If you are not the intended recipient, and have received this message in error, please immediately return this by e-mail and then delete it.

==============================================================================

-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._