plotGeoRatesCon: plotGeoRatesCon

In michellelawing/ppgm: PaleoPhyloGeographic Modeling of Climate Niches and Species Distributions

plotGeoRatesCon

Description

plotGeoRatesCon

Usage

plotGeoRatesCon(geo_center, geo_size, time_int, trees, path="")

Arguments

geo_center	change in geographic center of suitable climate envelope
geo_size	change in geographic size of suitable climate envelope
time_int	time intervals to plot
trees	distribution of phylogenies
path	path to the directory where the results to be saved

Details

Creates plot with gray background of all pairwise comparisons of change in geo center and area through time. Blue points on top show the sequential change in geo center and expansion/contraction for all lineages

Value

plots of geo rate

Author(s)

A. Michelle Lawing, Alexandra F. C. Howard

See Also

getGeoRates

Examples

data(sampletrees)
data(occurrences)
tree <- sampletrees[[25]]
occurrences <- getBioclimVars(occurrences, which.biovars=1)
sp_data_min<- tapply(occurrences[,4],occurrences$Species,min)
sp_data_max<- tapply(occurrences[,4],occurrences$Species,max)
treedata_min <- geiger::treedata(tree,sp_data_min,sort=TRUE,warnings=F)
treedata_max <- geiger::treedata(tree,sp_data_max,sort=TRUE,warnings=F)
full_est <- nodeEstimateEnvelopes(treedata_min,treedata_max)
node_est <- full_est$est
example_getEnvelopes <- getEnvelopes(treedata_min, treedata_max, node_est)
example_getGeoRate <- getGeoRate(example_getEnvelopes, tree,which.biovars=1)

## Not run: plotGeoRatesCon(example_getGeoRate$geo_center,example_getGeoRate$geo_size,
example_getGeoRate$time_int, trees = trees[[1]], path=tempdir())
## End(Not run)

michellelawing/ppgm documentation built on April 24, 2024, 10:17 p.m.