R语言【paleobioDB】——pbdb_temporal_resolution()：展示化石数据的时空分辨率

Package?paleobioDB?version 0.7.0

paleobioDB 包在2020年已经停止更新，该包依赖PBDB v1 API。

可以选择在Index of /src/contrib/Archive/paleobioDB (r-project.org)下载安装包后，执行本地安装。

Usage

pbdb_temporal_resolution (data, do.plot=TRUE)

Arguments

参数【data】：输入的数据，数据帧格式。可以通过 pbdb_occurrences() 函数获得数据。

参数【do.plot】：TRUE/FALSE。默认为 TRUE。

Value

一个图和一个列表，其中包含数据的时间分辨率（最小值、最大值、第 1 和第 3 个四分位数、中位数和平均值）以及每个化石记录的时间分辨率 （马） 的摘要。

Example

library(paleobioDB)
library(RCurl)

options(RCurlOptions = list(cainfo = system.file("CurlSSL", "cacert.pem", package = "RCurl")))


data<- pbdb_occurrences (taxon_name= "Canidae", interval= "Quaternary")

pbdb_temporal_resolution (data)





$summary
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
 0.0117  0.3846  1.5000  2.2046  2.8440 23.0183 

$temporal_resolution
 [1]  3.1000  4.6000  3.1000  3.1000  3.1000  3.1000  3.1000  3.1000
 [9]  4.6000  3.1000  1.5000  1.5000  1.5000  1.5000  1.5000  1.5000
[17]  1.5000  1.0250  0.7693  1.8070  2.5763  2.5763  2.5763  2.5763
[25]  2.5763  5.3213  2.5763  0.7693 23.0183  0.6550  2.5763  0.6550
[33]  0.6550  8.5000  0.6550  3.1000  0.1143  0.1143  2.5880  0.1143
[41]  0.1143  1.8070  0.1143  0.1143  0.1143  2.5880  0.1143  2.5763
[49]  0.1143  1.8070  0.0117  0.6550  0.1143  4.6000  0.0117  0.0117
[57]  0.1143  2.5763  0.1143

Page

function (data, do.plot = TRUE) 
{
    if ("eag" %in% colnames(data)) {
        diff <- .numeric_age(data$eag) - .numeric_age(data$lag)
        tr <- list(summary = summary(diff), temporal_resolution = diff)
    }
    if ("early_age" %in% colnames(data)) {
        diff <- .numeric_age(data$early_age) - .numeric_age(data$late_age)
        tr <- list(summary = summary(diff), temporal_resolution = diff)
    }
    if (do.plot == TRUE) {
        hist(unlist(tr[[2]]), freq = T, col = "#0000FF", border = F, 
            xlim = c(max(unlist(tr[[2]]), na.rm = TRUE), 0), 
            breaks = 50, xlab = "Temporal resolution of the data (Ma)", 
            main = "", col.lab = "grey30", col.axis = "grey30", 
            cex.axis = 0.8)
    }
    return(tr)
}