Uses wtd.quantile() from Hmisc.
weightcdf(df1 = my.samp, ss = ss, SampID = "Sample.ID", xvar = "cond", nt = 25, log = TRUE, np = 61)
df1 | ; default = my.samp |
---|---|
ss | Species crosstabed data; default = "ss" from global environment. |
SampID | Site/sample id column; default = "Sample.ID" |
xvar | variable on which to base calculations; default = "cond" |
nt | Minimum number of occurence; default = 25. |
log | A boolean for if the values should be log (base 10) transformed; default = TRUE |
np | Number of bins; default = 61. |
data frame: taxaname, XC95.all, N
### Step 1A. use samples to calcualte HC05 ### basic models for season, and ecoregion allyear <- weightcdf(df1 = bio.sample[!is.na(bio.sample$lgSO4HCO3),], ss = ss, SampID = "Sample.ID", xvar = "lgSO4HCO3", nt = 25)