Skip to content

Multitype point pattern analysis

2 messages · Ervan Rutishauser, Marcelino de la Cruz

#
Maybe solutionset can help you. Try the followin example:

data(lansing)
lansing.spl <- split(lansing)
  blackoak <- with (lansing.spl, density.ppp(blackoak))
  hickory <- with (lansing.spl, density.ppp(hickory))
maple <- with (lansing.spl, density.ppp(maple))

# where density of maple is > 500 and density of hickory is > 500?
plot(solutionset(maple>500 & hickory>500))

# Three types case: where density of maple is > 500 and density of 
hickory is > 500 and density of blackoak > 100?
  plot(solutionset(maple>500 & hickory>500 & blackoak>100))

HTH

Marcelino
At 13:57 15/11/2009, Ervan Rutishauser wrote:
***************************************
Marcelino de la Cruz Rot
Depto. Biologia Vegetal
EUIT Agricola
Universidad Politecnica de Madrid
28040-Madrid
SPAIN