search

AndrewLJackson / SIBER

ellipse and convex hull fitting package for stable isotope data
36 stars 14 forks source link

siberDensityPlot #69

Closed elenagoro closed 3 years ago

elenagoro commented 4 years ago

Hello, I am running into an error when trying to create siberDensityPlot for data with 2 communities, 2 groups in each: Error in hdrcde::hdr(dat[, j], probs, h = stats::bw.nrd0(dat[, j])) : Insufficient data The data set is below. Any help is much appreciated! /Elena

iso1 iso2 group community
-21.6 6.58 1 1
-21.96 7.12 1 1
-20.97 7.26 1 1
-23.06 7.26 1 1
-20.86 7.42 1 1
-20.55 7.58 1 1
-21.07 7.73 1 1
-21.57 6.18 1 1
-21.41 6.25 1 1
-22.7 6.31 1 1
-23.39 6.66 1 1
-20.6 6.77 1 1
-21.52 6.87 1 1
-21.15 7.09 1 1
-21.37 7.11 1 1
-21.38 7.35 1 1
-21.14 7.43 1 1
-21.36 7.19 1 1
-21.33 6.73 1 1
-21.16 7.24 1 1
-22.59 6.38 1 1
-22.76 7.42 2 1
-21.94 7.65 2 1
-21.39 7.91 2 1
-21.28 8.18 2 1
-21.47 8.25 2 1
-20.72 8.82 2 1
-19.81 11.89 2 1
-26.99 7.09 2 1
-24.07 5.68 2 1
-21.94 7.05 2 1
-21.88 7.06 2 1
-20.06 7.64 2 1
-21.95 6.88 1 2
-21.89 7.3 1 2
-21.18 7.54 1 2
-22.22 7.55 1 2
-21.4 7.55 1 2
-21.38 7.67 1 2
-21.27 7.89 1 2
-20.84 8.08 1 2
-20.62 9.59 1 2
-21.07 6.35 1 2
-21.86 6.83 1 2
-21.32 6.95 1 2
-22.61 6.96 1 2
-21.49 7.36 1 2
-20.95 8.14 1 2
-21.29 7.56 1 2
-21.16 7.36 1 2
-21.08 7.36 1 2
-24.33 5.94 2 2
-22.7 8 2 2
-22.43 8.14 2 2
-21.55 8.7 2 2
-23.19 8.83 2 2
-22 10.1 2 2
-19.11 12.42 2 2
-24.01 6.08 2 2
-23.58 6.38 2 2
-23.07 6.45 2 2
-23.6 6.85 2 2
-20.44 8.7 2 2
-23.59 7.29 2 2
elenagoro commented 4 years ago

Sorry, I have figured it out. Stupid me, there are only 2 groups in each community, which means that community hull metrics are screwed. No problem anymore...