<html><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div>You can also visualize the data to check if your My.Data is correct</div><div><br></div><div><br></div><div>plot(My.Data[,1], My.Data[,2], pch=".")</div><div><br></div><div>if you have everything here, it means our function level.plot has a problem somewhere. </div><div><br></div><div><br></div><br><div><div>Le 16 sept. 10 à 15:47, Bruno Lafourcade a écrit :</div><br class="Apple-interchange-newline"><blockquote type="cite"><font color="black" size="2" face="arial"> <div> <font size="2"><font face="Arial, Helvetica, sans-serif"><br> Hi again,<br> <br> You have a point on the bottom right corner which totally stands out of the rest (or maybe several points)<br> Try to find it in your data, take it off and plot again to see what you get. You can take the line out by doing this :<br> <br> <br> </font></font>level.plot(My.Data[<font color="blue"><font size="2"><font face="Arial, Helvetica, sans-serif">-156</font></font></font>,5], Coors, title='var 3')<font size="2"><font face="Arial, Helvetica, sans-serif"> <font color="blue">that would be if the line you wish to take off is line 156</font></font></font><br> </div> <div> <font size="2"><font face="Arial, Helvetica, sans-serif"><font color="blue">or create a new object</font><br> My.Data.NEW <- </font></font>My.Data[<font color="blue"><font size="2"><font face="Arial, Helvetica, sans-serif">-156</font></font></font>,<font size="2"><font face="Arial, Helvetica, sans-serif"> </font></font>]<br> <font size="2"><font face="Arial, Helvetica, sans-serif"><br> <br> Also, do consider that Arcgis has a way to represent the roundness of the world, especially true for a large<br> scale like yours. R doesn't do this and flattens the whole area which produces a differences.<br> <br> But it is true that the bottom part of your data is missing from the plot. Maybe it is due to the number of points<br> you are trying to visualise which possibly leads to a saturation for R's plotting device. How many points do you have ?<br> <br> Bruno<br> <br> <br> <br> </font></font></div> <div style="clear: both;">-------<br> Bruno Lafourcade<br> Statistical tools engineer<br> <br> Laboratoire d'Ecologie Alpine, bureau 308<br> CNRS - UMR 5553, 2233 rue de la piscine<br> 38400 Saint Martin d'Hères<br> -------</div> <div> <br> </div> <div> <br> </div> <div style="font-family: arial,helvetica; font-size: 10pt; color: black;">-----E-mail d'origine-----<br> De : Ashley Brooks <<a href="mailto:ashleycbrooks84@gmail.com">ashleycbrooks84@gmail.com</a>><br> A : <a href="mailto:biomod-commits@lists.r-forge.r-project.org">biomod-commits@lists.r-forge.r-project.org</a><br> Envoyé le : Jeudi, 16 Septembre 2010 15:31<br> Sujet : [Biomod-commits] discrepancies<br> <br> <div id="AOLMsgPart_3_de638867-6e6f-4d69-bae1-18a7b2c54de4"> Hello,<br> <br> Please forgive me for being slow at this and for having so many questions. I sampled my climate data in arcgis and then imported the .dbf file into R/Biomod to run. Shown below are the commands I used to retrieve this data. Then I did a level.plot to see my data in R. Attached are the images of what I got from R and what I see in arcgis. From what I can tell everything is the same (the number of coulmns and rows in my dbf are the same when viewed in R and in arcgis). I guess I don't understand why they look so different. Again, sorry for all the questions, but any help or advice or ideas would be most appreciated!<br> <br> Thank you so much,<br> Ashley<br> <br> <br> My.Data <- read.dbf("C:/Documents and Settings/New User/My Documents/curoccdata.dbf")<br> > Coors <- read.dbf("C:/Documents and Settings/New User/My Documents/Coors.dbf")<br> > head(My.Data)<br> X Y curbio1cli curbio2cli curbio7cli curbio12cl curbio15cl<br> 1 -95.14583 49.35417 22 112 491 568 54<br> 2 -95.10416 49.35417 22 112 491 570 54<br> 3 -95.06250 49.35417 22 112 490 572 54<br> 4 -95.02083 49.35417 22 111 491 573 54<br> 5 -95.14583 49.31251 21 113 492 567 55<br> 6 -95.10416 49.31251 22 113 492 567 55<br> curbio17cl Occurance<br> 1 60 0<br> 2 60 0<br> 3 60 0<br> 4 60 0<br> 5 57 0<br> 6 57 0<br> > head(Coors)<br> X Y<br> 1 -95.14583 49.35417<br> 2 -95.10416 49.35417<br> 3 -95.06250 49.35417<br> 4 -95.02083 49.35417<br> 5 -95.14583 49.31251<br> 6 -95.10416 49.31251<br> > level.plot(My.Data[,5], Coors, title='var 3')<br> </div> <!-- end of AOLMsgPart_3_de638867-6e6f-4d69-bae1-18a7b2c54de4 --> <div id="AOLMsgPart_6_de638867-6e6f-4d69-bae1-18a7b2c54de4" style="margin: 0px; font-family: Tahoma,Verdana,Arial,Sans-Serif; font-size: 12px; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);"> <pre style="font-size: 9pt;"><tt>_______________________________________________<br>
Biomod-commits mailing list<br>
<a href="mailto:Biomod-commits@lists.r-forge.r-project.org">Biomod-commits@lists.r-forge.r-project.org</a><br>
<a href="https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/biomod-commits" target="_blank">https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/biomod-commits</a><br>
</tt></pre> </div> <!-- end of AOLMsgPart_6_de638867-6e6f-4d69-bae1-18a7b2c54de4 --> <br> </div> </font> _______________________________________________<br>Biomod-commits mailing list<br><a href="mailto:Biomod-commits@lists.r-forge.r-project.org">Biomod-commits@lists.r-forge.r-project.org</a><br>https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/biomod-commits<br></blockquote></div><br><div> <span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div>-------------------------</div><div>Dr. Wilfried Thuiller</div><div>Laboratoire d'Ecologie Alpine, UMR CNRS 5553</div><div>Université Joseph Fourier</div><div>BP53, 38041 Grenoble cedex 9, France</div><div>tel: +33 (<span class="Apple-style-span" style="border-collapse: collapse; font-family: Arial, Helvetica, sans-serif; -webkit-border-horizontal-spacing: 10px; -webkit-border-vertical-spacing: 10px; ">0)4 76 51 44 97</span></div><div>fax: +33 (0)4 76 51 42 79</div><div><br></div><div>Email: <a href="mailto:wilfried.thuiller@ujf-grenoble.fr">wilfried.thuiller@ujf-grenoble.fr</a><br>Home page: <a href="http://www.will.chez-alice.fr/">http://www.will.chez-alice.fr</a><br>Website: <a href="http://www-leca.ujf-grenoble.fr/equipes/tde.htm">http://www-leca.ujf-grenoble.fr/equipes/tde.htm</a><br><br>FP6 European MACIS project: <a href="http://www.macis-project.net/">http://www.macis-project.net</a><br>FP6 European EcoChange project: <a href="http://www.ecochange-project.eu/">http://www.ecochange-project.eu</a></div><div><br></div><div><br></div><div><br></div></div></span><br class="Apple-interchange-newline"></div></span><br class="Apple-interchange-newline"> </div><br></body></html>