[R] Test of Presence Matrix HOWTO?

ONKELINX, Thierry Thierry.ONKELINX at inbo.be
Mon Feb 26 16:38:25 CET 2007


> a <- c("A","B","C","D")
> b <- c("A","B","E","F")
> c <- c("A","C","E","G")
> Df <- cbind(a, b, c)
> apply(Df, 2, function(x)(LETTERS[1:7] %in% x))
         a     b     c
[1,]  TRUE  TRUE  TRUE
[2,]  TRUE  TRUE FALSE
[3,]  TRUE FALSE  TRUE
[4,]  TRUE FALSE FALSE
[5,] FALSE  TRUE  TRUE
[6,] FALSE  TRUE FALSE
[7,] FALSE FALSE  TRUE
> 
> apply(Df, 2, function(x)(as.numeric(LETTERS[1:7] %in% x)))
     a b c
[1,] 1 1 1
[2,] 1 1 0
[3,] 1 0 1
[4,] 1 0 0
[5,] 0 1 1
[6,] 0 1 0
[7,] 0 0 1


Cheers,

Thierry
------------------------------------------------------------------------
----

ir. Thierry Onkelinx

Instituut voor natuur- en bosonderzoek / Reseach Institute for Nature
and Forest

Cel biometrie, methodologie en kwaliteitszorg / Section biometrics,
methodology and quality assurance

Gaverstraat 4

9500 Geraardsbergen

Belgium

tel. + 32 54/436 185

Thierry.Onkelinx op inbo.be

www.inbo.be 

 

Do not put your faith in what statistics say until you have carefully
considered what they do not say.  ~William W. Watt

A statistical analysis, properly conducted, is a delicate dissection of
uncertainties, a surgery of suppositions. ~M.J.Moroney


> -----Oorspronkelijk bericht-----
> Van: r-help-bounces op stat.math.ethz.ch [mailto:r-help-
> bounces op stat.math.ethz.ch] Namens Johannes Graumann
> Verzonden: maandag 26 februari 2007 16:25
> Aan: r-help op stat.math.ethz.ch
> Onderwerp: [R] Test of Presence Matrix HOWTO?
> 
> Hello,
> 
> Imagine 3 lists like so:
> 
> > a <- list("A","B","C","D")
> > b <- list("A","B","E","F")
> > c <- list("A","C","E","G")
> 
> What I need (vennDiagram) is a matrix characterizing with 1 or 0
whether
> any
> given member is present or not like so:
>      x1 x2 x3
> [1,]  1  1  1
> [2,]  1  1  0
> [3,]  1  0  1
> [4,]  1  0  0
> [5,]  0  1  1
> [6,]  0  1  0
> [7,]  0  0  1
> 
> (where the rows represent "A"-"G" and the columns a-c, respectively).
> 
> > table(c(a,b,c))
> will give me a quick answer for the "1 1 1" case, but how to deal with
the
> other cases efficiently without looping over each string and looking
for
> membership %in% each list?
> 
> Thanks for enlightening the learning,
> 
> Joh
> 
> ______________________________________________
> R-help op stat.math.ethz.ch 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.



More information about the R-help mailing list