mapast: mapast

In macroecology/pasta: Paleogeography Combined with Spatial Paleobiodiversity

Description

Plots your fossil occurrence data onto the corresponding MAp of the PAST.

Usage

mapast(model = "SETON2012", data, map = NULL, do.plot = TRUE, save.as = NULL,
                         colland = "#66666660",
                         colsea = "#00509010", 
                         colpoints = "#65432190", 
                         pch = 16, cex = 1, ...)

Arguments

model	character. The model the map should be created with. "SETON2012" (default), "MULLER2016", "GOLONKA", "PALEOMAP" or "MATTHEWS2016".
data	data.frame. Fossil occurrences data.
map	(list of) SpatialPolygonDataFrames. Containing map(s) which can be created by getmap.
do.plot	logical. If a plot is created or not. By default do.plot = TRUE.
save.as	character. The format the plots should be saved. "tiff", "pdf", "jpeg" or "png".
colland	character. The color of the land masses. By default colland = "#66666660".
colsea	character. The color of the sea. By default colsea = "#00509010".
colpoints	character. The color of the occurrence-points. By default colpoints = "#65432190".
pch	numeric. Point symbol for plotting the occurences. By default pch = 16 (filled circle).
cex	numeric. Size of the fossil occurrences points. By default cex = 1.
...	Graphical parameters. Any argument that can be passed to image.plot and to plot, such as main="my own title" or main.col="red".

Value

Plot

Examples

## Not run: 

library(mapast)

#get data and preprocess it
data  <-  base::data.frame(paleobioDB::pbdb_occurrences(base_name = "Canis", 
                                                        min_ma = 0, max_ma = 10, 
                                                        show = c("coords", "phylo"), 
                                                        vocab = "pbdb", limit = 100))
df <- formatdata(data)
df_auto <- paleocoords(df, time = "automatic")
                                                        
#create a plot with fossils on the paleogeographical map
mapast(model = "SETON2012", data = df_auto)

#save the maps before so the function does not need to load them
maps <- getmap(ma = 2.5, model = "SETON2012", do.plot = FALSE)
mapast(model = "SETON2012", data = df_auto, map = maps)

#save maps as pdf
mapast(model = "SETON2012", data = df_auto, map = maps, save.as = "pdf")



## End(Not run)

macroecology/pasta documentation built on March 1, 2020, 11:22 p.m.