Gives natural statistics (number of observations, average, standard deviations) for each variable and each zones.
ANOVA of the model "variable~zone" is run, and difference tests are computed to find groups.
Usage
summaryrshp(rshp, test = "hsd", p.adj = "none", alpha = 0.05)
"lsd" (corresponding to Least Significant Difference) or "hsd" (corresponding to the Tukey Test)
p.adj
LSD.test parameter to adjust the p-value for the multiple comparison
alpha
Threshold of significance for determining the gorups
Value
Returns a list containing, for each variable: the results of the ANOVA; the summary of the variable; and the usual statistics by zones (n,average, sd,group)