[R] Mantel - like analysis

serrouya serrouya at ualberta.ca
Tue May 31 22:46:36 CEST 2011


Hello,

I've used mantel tests a fair bit for standard genetic/geographic/landscape
analyses.  I realize it involves permutations of rows and columns in the
matrix to obtain correlation coefficients.

However, for a certain question, I am interested in populations that are
nearest neighbours to each other, and identifying factors that affect
genetic distance between them.  So because of autocorrelation problems, I
can't use standard regression or correlation, because some populations
appear several times depending upon their spatial arrangement.  In essence I
am cherry picking a typical matrix, but only for populations that are
closest.

Is there code or a test that can handle this problem? It seems very similar
to a mantel, but only for a subset of possible pairs, not all possible
pairwise comparisons as a mantel would handle.  What about a bootstrapping
algorithm than handles populations as the sample unit instead of individual
data points?
I bet mixed models could handle this problem somehow, but I'd like to stick
to the theme of a mantel test.

The data look like this

pop1 pop2  genetic_distance(dependent)  (independent1)
popa popb    data    data
popb popd    data    data
popc pope    data    data
popf popg    data    data
popg  poph  data   data

So for example, popb appears twice, as does popg.  Clearly each row is not
independent.

Thanks for any insight

Rob.

--
View this message in context: http://r.789695.n4.nabble.com/Mantel-like-analysis-tp3564296p3564296.html
Sent from the R help mailing list archive at Nabble.com.



More information about the R-help mailing list