<div dir="ltr">Thank you for the answer<br>I want to test whether space (lat+lon) has significant effect on the genetic structure. Therefore I would like to use spatial variables in the right side of the model. Can you suggest a better representation of the spatial structures than lat-lon?<br>

Thank you<br>Hanan    <br><br><div class="gmail_quote">On Thu, Jul 12, 2012 at 1:58 PM, Jombart, Thibaut <span dir="ltr"><<a href="mailto:t.jombart@imperial.ac.uk" target="_blank">t.jombart@imperial.ac.uk</a>></span> wrote:<br>

<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear Hanan,<br>
<br>
this is a tricky question, and I don't think there is a single universal answer. Technically speaking, the only requirement is that your residuals are independent, so you need to make sure there is no spatial autocorrelation left there. Otherwise minimizing the sum of squared residuals is no longer the ML solution.<br>


<br>
The real problem relates to the interpretation, and the assumption implicitly made by the model. There is several reasons why spatial genetic patterns can occur. Your model has the form:<br>
genetic pattern = lat+lon + environment + residuals<br>
<br>
Which means that beyond linear trends, genetic patterns are due to the environment. It makes sense to treat spatial autocorrelation as a confounding factor first removed from the analysis. But lat+lon is often not enough to capture all spatial structures. In this respect, using PCs from PCA on the left side is probably better than sPCA (no need to seek spatial structures to remove them afterwards).<br>


<br>
Cheers<br>
<br>
Thibaut<br>
<br>
________________________________________<br>
From: <a href="mailto:adegenet-forum-bounces@lists.r-forge.r-project.org">adegenet-forum-bounces@lists.r-forge.r-project.org</a> [<a href="mailto:adegenet-forum-bounces@lists.r-forge.r-project.org">adegenet-forum-bounces@lists.r-forge.r-project.org</a>] on behalf of Hanan Sela [<a href="mailto:dooshra@gmail.com">dooshra@gmail.com</a>]<br>


Sent: 12 July 2012 07:34<br>
To: <a href="mailto:adegenet-forum@lists.r-forge.r-project.org">adegenet-forum@lists.r-forge.r-project.org</a><br>
Subject: [adegenet-forum] Using PCA of SPCA in linear models with       environmental data.<br>
<br>
Hello all<br>
I am trying to estimate the major factors affecting the spatial distribution of wild wheat genotypes.  I am using a linear model where the PCA or the SPCA   first and second axis are the dependent variables and the environmental variables are the predictors. Additionally I am using the longitude and the latitude as predictors.   Since there is a spatial reference on the left side of the formula, I was wondering if using SPCA on the right side will not be a problem.<br>


Thank you<br>
Hanan</blockquote></div><br><br clear="all"><br>-- <br><div dir="ltr"><div>
<div>Hanan Sela Ph.D.</div>
<div>Curator of the Lieberman Cereal Germplasm Bank</div>
<div>The Institute for Cereal Crops Improvement<br>Tel-Aviv University</div>
<div>P.O. Box 39040</div>
<div>Tel Aviv 69978 </div>
<div>Israel </div> </div>
<div><a href="mailto:hans@tauex.tau.ac.il" target="_blank">hans@tauex.tau.ac.il</a> <br></div>
<div>Phone: 972-3-6405773</div>
<div>Cell: 972-50-5727458<br>Fax: 972-3-6407857</div></div><br>
</div>