Dear Diogo,
the following function should give you the p-values (physo.obs.p) you
are looking for:
# obs : observed phylosor-values (distance matrix)
# rand : list of distance matrices from the phylosor.rnd-function
ses.physo - function(obs, rand){
rand - t(as.data.frame(lapply(rand,
I have a doubt with using the picante::phylosor.rnd for testing hypothesis
about the phylogenetic dissimilarities between sites. This function gives
the randomized matrices obtained with the null model chosen. But how can I
do to explicitly test for a significant difference between the observed
