I have a 3D ordination with two groups, and want to identify the points in one group that lie within the 3D convex hull of the other. I would appreciate any enlightenment about why the command at the end (below in a toy example) crashes R when I try to ID the points in the hull using pip3d(), from package ptinpoly. convhulln() is from package geometry. Have I not supplied the right ?Faces? argument to pip3d()?
The objects are not large and are reproduced at the bottom of this message.
Thanks
------------------------
sessionInfo()
R version 3.0.3 (2014-03-06)
Platform: x86_64-apple-darwin10.8.0 (64-bit)
locale:
[1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8
attached base packages:
[1] stats graphics grDevices utils datasets methods base
other attached packages:
[1] ptinpoly_2.0 misc3d_0.8-4 geometry_0.3-4 magic_1.5-6 abind_1.4-0
loaded via a namespace (and not attached):
[1] tools_3.0.3
??