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)
}
相关推荐
wjs20242 小时前
状态模式(State Pattern)
开发语言
我命由我123452 小时前
Kotlin 数据容器 - List(List 概述、创建 List、List 核心特性、List 元素访问、List 遍历)
java·开发语言·jvm·windows·java-ee·kotlin·list
liulilittle2 小时前
C++ TAP(基于任务的异步编程模式)
服务器·开发语言·网络·c++·分布式·任务·tap
励志要当大牛的小白菜4 小时前
ART配对软件使用
开发语言·c++·qt·算法
Chef_Chen4 小时前
从0开始学习R语言--Day55--弹性网络
r语言
爱装代码的小瓶子6 小时前
数据结构之队列(C语言)
c语言·开发语言·数据结构
Maybe_ch7 小时前
.NET-键控服务依赖注入
开发语言·c#·.net
超浪的晨7 小时前
Java UDP 通信详解:从基础到实战,彻底掌握无连接网络编程
java·开发语言·后端·学习·个人开发
终焉暴龙王7 小时前
CTFHub web进阶 php Bypass disable_function通关攻略
开发语言·安全·web安全·php
Edingbrugh.南空8 小时前
Aerospike与Redis深度对比:从架构到性能的全方位解析
java·开发语言·spring