8 messages
Open this post in threaded view
|

 Dear all, I am working on a multitype point pattern, and I'd like to estimate how many of each type of point occurs into each quadrant. I know it is possible to use the quandracount on split marks as follows using spatstats: quadratcount(split(marks)). But the result produces as many windows as they are marks. I am wondering is there is a way to know many occurrence of each type there is per quadrant and to plot it in a single grid. Thanks, Guy This email and any files transmitted with it are confide...{{dropped:10}} _______________________________________________ R-sig-Geo mailing list [hidden email] https://stat.ethz.ch/mailman/listinfo/r-sig-geo
Open this post in threaded view
|

## Re: quadracount on multitype points

 Hi Guy, You can know how many points of each type are in each cell with sapply, i.e.:  >data (lansing)  >qs<- quadratcount(split(lansing)) #trees in the fist row, first column:  > sapply(qs, function(x) x[1,1]) blackoak  hickory    maple     misc   redoak whiteoak        13       66        0        0        4       23 #trees in the fist row, second column:  > sapply(qs, function(x) x[1,2]) blackoak  hickory    maple     misc   redoak whiteoak        16       38        3        0        9       16 # ETC. On the other hand I don't know if it is a good idea to have all this numbers printed in just one cell or how to print them without the result being a mess... Cheers, Marcelino. El 12/09/2017 a las 16:11, Guy Bayegnak escribió: > Dear all, > > > > I am working on a multitype point pattern, and I'd like to estimate how many of each type of point occurs into each quadrant. I know it is possible to use the quandracount on split marks as follows using spatstats: quadratcount(split(marks)). But the result produces as many windows as they are marks. I am wondering is there is a way to know many occurrence of each type there is per quadrant and to plot it in a single grid. > > Thanks, > > Guy > > > > This email and any files transmitted with it are confide...{{dropped:10}} > > _______________________________________________ > R-sig-Geo mailing list > [hidden email] > https://stat.ethz.ch/mailman/listinfo/r-sig-geo> . > -- Marcelino de la Cruz Rot Depto. de Biología y Geología Física y Química Inorgánica Universidad Rey Juan Carlos Móstoles España _______________________________________________ R-sig-Geo mailing list [hidden email] https://stat.ethz.ch/mailman/listinfo/r-sig-geo
Open this post in threaded view
|

## Re: quadracount on multitype points

 In reply to this post by Guy Bayegnak On 13/09/17 02:11, Guy Bayegnak wrote: > Dear all, > > I am working on a multitype point pattern, and I'd like to estimate > how many of each type of point occurs into each quadrant. I know it > is possible to use the quandracount on split marks as follows using > spatstats: quadratcount(split(marks)). But the result produces as > many windows as they are marks. I am wondering is there is a way to > know many occurrence of each type there is per quadrant and to plot > it in a single grid. > > Thanks. You really should start by mentioning that you are dealing with the spatstat package. It's not clear to me what you are after.  A minimal reproducible example would be helpful.  I presume that by "quadrant" you mean one of the four equal sub-windows formed by bisecting your (rectangular) window vertically and horizontally. If my presumption is correct then perhaps      lapply(split(X),quadratcount,nx=2) (where "X" is your point pattern) does what you want.  E.g.: > lapply(split(amacrine),quadratcount,nx=2) > \$off >          x > y         [0,0.801) [0.801,1.6] >   [0.5,1]        36          36 >   [0,0.5)        34          36 > > \$on >          x > y         [0,0.801) [0.801,1.6] >   [0.5,1]        35          39 >   [0,0.5)        42          36 Is this something like what you wish to achieve? cheers, Rolf Turner -- Technical Editor ANZJS Department of Statistics University of Auckland Phone: +64-9-373-7599 ext. 88276 _______________________________________________ R-sig-Geo mailing list [hidden email] https://stat.ethz.ch/mailman/listinfo/r-sig-geo
Open this post in threaded view
|

## Re: quadracount on multitype points

Open this post in threaded view
|

## Re: quadracount on multitype points

Open this post in threaded view
|