Error in .subset2 about spthin HOT 7 CLOSED

@kgahm Do you have a species column in your data frame? That's required for this program to work correctly.

from spthin.

@GisKid Sorry for the long delay. You still having this error? It looks like it may be a case of having too few data points, but I'd need to know more about the dataset you're dealing with to be sure.

from spthin.

from spthin.

Oh interesting. I have to admit that I never tested the code on anything other than long/lat in decimal-degree format. Is this UTM? The distance between points is determined using the function rdist.earth in the fields package. Best I can tell, it doesn't return reasonable distances with the 6 points you have above. Is it possible to convert them to a decimal degree format? I'll have to think about adding conversions into future editions of spThin.

from spthin.

I'm getting the same error as the above, but my coordinates are in decimal degrees.
This is the code that produces the error:
thinned <- thin(loc.data = data, lat.col = "lat", long.col = "lon", thin.par = 1, reps = 10)

and this is the error message itself:
Error in .subset2(x, i, exact = exact) : attempt to select less than one element in get1index.

Here's a bit of my data:

head(data)
lon lat
-91.34398 -0.29502
-91.34030 -0.29467
-91.34172 -0.29463
-91.33602 -0.29402
-91.33930 -0.29394
-91.33914 -0.29369

Any idea what might be going on here? Are my coordinates in the wrong format, or could there be a different issue?

Thank you!

from spthin.

@mlammens thanks for getting back to me so quickly! That was indeed the problem, and it's working now.

I think I got confused about this because the documentation says that the loc.data object "can include several columns, but must include at minimum a column of latitude and a column of longitude values." So I didn't realize a species column was required.

If you get the chance, it might be helpful to include a note that the species column is also required. But anyway I probably should have known to try that...

Thanks again for your help!

from spthin.

Thanks @kgahm. I'll add a note to fix that documentation.

from spthin.