<div dir="ltr"><div>Hey Roberto!</div><div><br></div>Never done it in Grass but its a typical operation. <div><br></div><div>For me it seems a bit confusing with for me since you are not showing me a glm fitting but instead just using the occurrence map in the logistic link function.</div><div><br></div><div><br></div><div>What you need to do is to fit a GLM to the binary data you have, something you can do in R for example check here: </div><div><br></div><blockquote style="margin:0 0 0 40px;border:none;padding:0px"><blockquote style="margin:0 0 0 40px;border:none;padding:0px"><div><a href="http://www.statmethods.net/advstats/glm.html">http://www.statmethods.net/advstats/glm.html</a></div></blockquote></blockquote><div><br></div><div>Then you will build your output in grass using the r.mapcalc: F(x) = 1 / ( 1+ e^(-GLM)  (That will be your probability map). </div><div><br></div><div>But, you can actually do everything using R only, no need to come into grass but you can if you want. </div><div><br></div><div>Good luck! </div><div>Ciao!</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On 5 September 2014 10:46, Roberto Marzocchi <span dir="ltr"><<a href="mailto:roberto.marzocchi@gmail.com" target="_blank">roberto.marzocchi@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div></div><div> Dear all, <br></div>
<br>
I have dichotomous data (occurrence 1 and non-occurence 0) and I want to
 use the  logistic regression method in order to obtain susceptibility 
map of landslide occurrence. <br>
<br>
I try to use r.regression.multi, but I am not sure to correctly define the mapy variable..<br>
<br>
I try to use r.mapcalc<br>
<br>
<i>r.mapcalc expression="logit_occurrence=log(occurrence_map+0.00001)-log(1-occurrence_map+0.00001)"</i><br>
<br>
obtaining a map that go from -11.5 (not occurrence) to +11.5 (occurrence) <br>
<br>
when I try to use r.regression.multi I obtain a lot of null values in residuals and estimates map.<br>
<br>
i do not understand why.. <br>
<br>
Thanks in advanced for the attention<span class="HOEnZb"><font color="#888888"><br>
Roberto</font></span></div>
<br>_______________________________________________<br>
grass-user mailing list<br>
<a href="mailto:grass-user@lists.osgeo.org">grass-user@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/mailman/listinfo/grass-user" target="_blank">http://lists.osgeo.org/mailman/listinfo/grass-user</a><br></blockquote></div><br><br clear="all"><div><br></div>-- <br><div dir="ltr"><br>Nuno César de Sá<br>+351 91 961 90 37<br><br></div>
</div>