Hi! I have managed to sort this problem out by adapting script from a published 'instruction'. This gives: < m1 <- mantelkdist(kdist(data gene, data geo),999) < plot (m1).
The script produced two sub-plots: one showing an xy plot of gene(x) end geo(y) data and another showing a histogram of gene-geo 'r', together with a marker for the 'observed value' - just what I wanted. I hope this helps other neophytes! -- View this message in context: http://n4.nabble.com/Mantel-results-shown-in-a-graphical-histogram-form-tp1009216p1120041.html Sent from the R help mailing list archive at Nabble.com. ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.