[Biomod-commits] Re : discrepancies

[Wilfried Thuiller]

You can also visualize the data to check if your My.Data is correct


plot(My.Data[,1], My.Data[,2], pch=".")

if you have everything here, it means our function level.plot has a  
problem somewhere.



Le 16 sept. 10 à 15:47, Bruno Lafourcade a écrit :

>
> Hi again,
>
> You have a point on the bottom right corner which totally stands out  
> of the rest (or maybe several points)
> 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 :
>
>
> level.plot(My.Data[-156,5], Coors, title='var 3')  that would be if  
> the line you wish to take off is line 156
> or create a new object
> My.Data.NEW <- My.Data[-156, ]
>
>
> Also, do consider that Arcgis has a way to represent the roundness  
> of the world, especially true for a large
> scale like yours. R doesn't do this and flattens the whole area  
> which produces a differences.
>
> 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
> you are trying to visualise which possibly leads to a saturation for  
> R's plotting device. How many points do you have ?
>
> Bruno
>
>
>
> -------
> Bruno Lafourcade
> Statistical tools engineer
>
> Laboratoire d'Ecologie Alpine, bureau 308
> CNRS - UMR 5553, 2233 rue de la piscine
> 38400 Saint Martin d'Hères
> -------
>
>
> -----E-mail d'origine-----
> De : Ashley Brooks <ashleycbrooks84 at gmail.com>
> A : biomod-commits at lists.r-forge.r-project.org
> Envoyé le : Jeudi, 16 Septembre 2010 15:31
> Sujet : [Biomod-commits] discrepancies
>
> Hello,
>
> 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!
>
> Thank you so much,
> Ashley
>
>
>  My.Data <- read.dbf("C:/Documents and Settings/New User/My  
> Documents/curoccdata.dbf")
> > Coors <- read.dbf("C:/Documents and Settings/New User/My Documents/ 
> Coors.dbf")
> > head(My.Data)
>           X        Y curbio1cli curbio2cli curbio7cli curbio12cl  
> curbio15cl
> 1 -95.14583 49.35417         22        112        491         
> 568         54
> 2 -95.10416 49.35417         22        112        491         
> 570         54
> 3 -95.06250 49.35417         22        112        490         
> 572         54
> 4 -95.02083 49.35417         22        111        491         
> 573         54
> 5 -95.14583 49.31251         21        113        492         
> 567         55
> 6 -95.10416 49.31251         22        113        492         
> 567         55
>   curbio17cl Occurance
> 1         60         0
> 2         60         0
> 3         60         0
> 4         60         0
> 5         57         0
> 6         57         0
> > head(Coors)
>           X        Y
> 1 -95.14583 49.35417
> 2 -95.10416 49.35417
> 3 -95.06250 49.35417
> 4 -95.02083 49.35417
> 5 -95.14583 49.31251
> 6 -95.10416 49.31251
> > level.plot(My.Data[,5], Coors, title='var 3')
> _______________________________________________
>
> Biomod-commits mailing list
>
> Biomod-commits at lists.r-forge.r-project.org
>
> https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/biomod-commits
>
>
> _______________________________________________
> Biomod-commits mailing list
> Biomod-commits at lists.r-forge.r-project.org
> https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/biomod-commits

-------------------------
Dr. Wilfried Thuiller
Laboratoire d'Ecologie Alpine, UMR CNRS 5553
Université Joseph Fourier
BP53, 38041 Grenoble cedex 9, France
tel: +33 (0)4 76 51 44 97
fax: +33 (0)4 76 51 42 79

Email: wilfried.thuiller at ujf-grenoble.fr
Home page: http://www.will.chez-alice.fr
Website: http://www-leca.ujf-grenoble.fr/equipes/tde.htm

FP6 European MACIS project: http://www.macis-project.net
FP6 European EcoChange project: http://www.ecochange-project.eu






-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.r-forge.r-project.org/pipermail/biomod-commits/attachments/20100916/b047a8ac/attachment-0001.htm>