[Biomod-commits] Re : discrepancies
Bruno Lafourcade
brunolafourcade at aol.com
Thu Sep 16 15:47:23 CEST 2010
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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.r-forge.r-project.org/pipermail/biomod-commits/attachments/20100916/266e1ef7/attachment.htm>
More information about the Biomod-commits
mailing list