for line 29 of gstudio/R/partition.R, I believe the line should be changed from:
ret[[lvl]] <- droplevels( x[ x[[stratum]]==lvl, ] )
to:
ret[[lvl]] <- droplevels( x[ which(x[[stratum]]==lvl), ] )
in order to avoid empty rows of NA within each population corresponding to every line of the input data frame that doesn't match the stratum.
for line 29 of gstudio/R/partition.R, I believe the line should be changed from:
ret[[lvl]] <- droplevels( x[ x[[stratum]]==lvl, ] )
to:
ret[[lvl]] <- droplevels( x[ which(x[[stratum]]==lvl), ] )
in order to avoid empty rows of NA within each population corresponding to every line of the input data frame that doesn't match the stratum.