I am new at ‘R’ and not even sure whether I am contacting the appropriate forum!! However, I have carried out several mantel tests (using ‘R’) on genetic data and would like to plot the results in the form of a histogam of the distribution – and preferably use a procedure that shows the critical ‘r’ value.
The tail end of my analysis is shown below. I just need help to dig out an instruction to follow on from this result.
> mantel.rtest(sample.dists, cat.dists, nrepet = 9999)
Monte-Carlo test
Observation: 0.4431225
Call: mantel.rtest(m1 = sample.dists, m2 = cat.dists, nrepet = 9999)
Based on 9999 replicates
Simulated p-value: 1e-04