Dear list:

 I have a ppp (i.e. spatstat) object representing a mapped stand of
 ~4000 trees. The polygon enclosing the stand is somewhat irregular
and
 has 9 edges. I am using the 'siplab' package for spatially explicit
 individual-based vegetation models, to calculate a simple
competition
 index (Hegyi’s index, based on diameter of competitor tree divided
 by diameter of target tree, times the reciprocal of distance between
 them) using the example given in the package documentation, i.e., 

 pairwise(myStandData, maxR = 6, kernel=powersker, kerpar =

 list(pi=1, pj=1, pr=1, r0=0.3048, smark="diameter")) 

  However, I am getting the error “Error in imarks[[smark]] :
 subscript out of bounds”. 

 imarks[[smark]] refers to the diameter of the target tree, yet all
 trees have diameter records -- there are no 'NA's. Can you suggest
how
 I might troubleshoot this problem? I can post data to github if
 necessary. 
   Seth Bigelow

        [[alternative HTML version deleted]]

_______________________________________________
R-sig-ecology mailing list
[email protected]
https://stat.ethz.ch/mailman/listinfo/r-sig-ecology

Reply via email to