Description Usage Arguments Examples
View source: R/qtlRaster_v2.1.R
Select top X percent of probability surface (defined by area percentage or cumulative probability)
1 | qtlRaster(pdR, threshold, thresholdType)
|
pdR |
probability density map which must be raster or RasterLayer |
threshold |
numeric from 0 to 1, top percentage to choose |
thresholdType |
numeric 1 or 2. 1 is to choose top percentage by cumulative probability; 2 is to choose top percentage by area |
genplot |
generate summary plot? (T or F) |
pdf |
save summary plot as PDF? (T or F) |
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 | # load North America SpatialPolygonsDataFrame data in assignR package
data("naMap")
data(d2h_world)
d = subOrigData(taxon = c("Homo sapiens"), mask = naMap)
# Exclude some outliers. This step is optional, which depends on your data quality
d <-as.data.frame(d)
dd = d[d$coords.x1<(-80),]
dd <- SpatialPointsDataFrame(dd[,2:3], as.data.frame(dd[,1]))
crs(dd) <- "+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0"
# Rescale from environmental isoscape to tissue isoscape (known must be SpatialPointsDataFrame with coordinate reference system )
r = calRaster(known = dd, isoscape = d2h_world, mask = naMap)
# four unknown-origin example
id = c("A", "B", "C", "D")
d2H = c(-110, -180, -130, -150)
un = data.frame(id,d2H)
# assignment
asn = pdRaster(r,unknown=un,mask=naMap)
# top 10% of probability surface (defined by % area)
qtlRaster(asn, threshold = 0.1, thresholdType = 2, pdf = T)
# top 10% of probability surface (defined by % cumulative probability)
qtlRaster(asn, threshold = 0.1, thresholdType = 1, pdf = F)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.