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

R 复制代码
pbdb_temporal_resolution (data, do.plot=TRUE)

Arguments

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

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


Value

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


Example

R 复制代码
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

R 复制代码
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)
}
相关推荐
Kiri霧33 分钟前
Linux下的Rust 与 C 的互操作性解析
c语言·开发语言·rust
雪芽蓝域zzs1 小时前
uniapp AES 加密解密
开发语言·uni-app·c#
雨夜的星光1 小时前
Python JSON处理:load/loads/dump/dumps全解析
开发语言·python·json
fen_fen2 小时前
Java打包时,不将本地Jar打包到项目的最终 JAR 中
开发语言·python·pycharm
wjs20244 小时前
SQL AND & OR 操作符详解
开发语言
~光~~4 小时前
【环境配置 】WSL2 +ubuntu20.04 +Qt配置+Kits配置
开发语言·qt·ubuntu
·心猿意码·5 小时前
C++右值语义解析
开发语言·c++
小龙报5 小时前
《彻底理解C语言指针全攻略(2)》
c语言·开发语言·c++·visualstudio·github·学习方法
zzzsde6 小时前
【c++】深入理解string类(4)
开发语言·c++
郝学胜-神的一滴6 小时前
Effective Python 第44条:用纯属性与修饰器取代旧式的 setter 与 getter 方法
开发语言·python·程序人生·软件工程