Dear all,
i'm using gdistance to model animal movement across landscape.
I have imported 11 rasters with roads, freeways, slope, use-of-land, lakes
(...) after recoding them with GRASS with a HSI value ranging from 1 to 4.
I've assigned zero to the NAs and then transformed all the rasters in
TransitionLayers (function=mean,directions=4) and later summed all of them
into a new transition:
> MY_transitionclass : TransitionLayer
dimensions : 2181, 1648, 3594288 (nrow, ncol, ncell)
resolution : 100.0049, 99.98945 (x, y)
extent : 1460708, 1625516, 4947383, 5165460 (xmin, xmax, ymin, ymax)
coord. ref. : +proj=tmerc +lat_0=0 +lon_0=9 +k=0.9996 +x_0=1500000
+y_0=0 +ellps=intl +towgs84=-225,-65,9,0,0,0,0 +units=m +no_defs
values : conductance
matrix class: dsCMatrix
I applied two different geocorrections to this same transition:
> MY_correction<-geoCorrection(MY_transition,"c",F,scl=T)
> MY_Rcorrection<-geoCorrection(MY_transition,"r",F,scl=T)
My coords are 44 points where samples were taken (some further genetic
analyses are programmed), but i've also tried with only two points:
> coords<-read.csv("coords_md.csv",header=F)
> mycoords<-as.matrix(coords)
> my_spatialpoints<-SpatialPoints(mycoords)
Finally I calculated three geographical distance matrices:
> geodist <- pointDistance(my_spatialpoints,longlat=FALSE)
> summary(geodist) V1 V2 V3 V4
> V5
Min. : 0 Min. : 0 Min. : 0 Min. : 0
Min. : 0
1st Qu.:133788 1st Qu.:133788 1st Qu.:133788 1st Qu.:133788
1st Qu.:133788
Median :133788 Median :133788 Median :133788 Median :133788
Median :135168
Mean :114704 Mean :117371 Mean :120166 Mean :123097
Mean :126174
3rd Qu.:137971 3rd Qu.:137971 3rd Qu.:137971 3rd Qu.:137971
3rd Qu.:137971
Max. :137971 Max. :137971 Max. :137971 Max. :137971
Max. :137971
NA's : 1 NA's : 2 NA's : 3
NA's : 4
(....)
> costdist <- costDistance(MY_correction, my_spatialpoints)
> summary(costdist) Min. 1st Qu. Median Mean 3rd Qu. Max.
0.00 0.00 21.92 68.53 210.90 244.70
> resdist <- resistanceDistance(MY_Rcorrection, my_spatialpoints)
> summary(resdist) Min. 1st Qu. Median Mean 3rd Qu. Max. 0 0 10220000
17070000 42860000 48300000
Now, I'd like to represent this measures but I couldn't find any examples
of how to do it; plotting or imaging them doesn't return a meaningful
graph, I guess I should transform them somehow but I don't know how.
ANY help would be really appreciated, along with comments about the rest of
the work done.
Thanks in advance,
Lisa
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
