Skip to content

Mantel results shown in a graphical (histogram) form.

2 messages · colin1

colin1
# 
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.
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

Colin.
14 days later
colin1
# 
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!