英国石油公司

Travis构建状态 AppVeyor生成状态

CRAN_状态_标签

bp:血压分析对于R

心血管疾病(CVD)是世界范围内的主要死亡原因高血压,特别是影响11亿多人的高血压每年。该包的目标是提供一个全面的工具箱使用各种统计数据分析血压(BP)数据为CVD带来更清晰的指标和可视化。

安装

您可以从安装发布版本的bpCRAN(起重机)具有:

安装.包(“bp”)

您可以从安装开发版本github具有:

#install.packages(“devtools”)
开发工具::安装github(“johnschwenck/bp”)

对于带有渐晕图的安装:

开发工具::安装github(“johnschwenck/bp”构建渐晕图= 真的)

示例数据

该软件包包括六个样本数据集:

预期功能

这个英国石油公司包的设计允许用户通过处理_数据通过指定以下变量的任意组合出现在用户提供的数据集中(最低要求是SBP公司数据库管理员包括在内)。该包裹将然后利用处理过的数据帧计算来自医学和统计文献,并提供可视化。也许该软件包最有用的用户友好功能是能够生成可视化报告以识别关系并评估受试者的血压阶段进展。

该软件包能够利用以下功能生理变量(以整数表示):

还有动脉压(AP)的处理功能包含以下输入的数据:*对应于AP数据(AP公司)*经过的时间(时间重叠(_E))

此外,在处理_数据功能,如:

识别并处理所有可用变量后处理后的结果数据帧用于所有其他功能。

唯一的英国石油公司包是创建最初可能不在提供的数据集。目前,以下附加列将是创建:

有关详细信息,请参见下面的示例。

可用指标

然后,软件包将利用上述变量进行计算医学和统计文献中的各种指标量化读数的可变性并将其分类为高血压的不同类别(正常、升高或高血压)。

以下指标目前通过英国石油公司包裹:

功能 指标名称 来源
反向传播 平均实际变量 梅纳等人al(2005)
bp_中心 平均值和中位数 阿马罗·利贾西奥等人(2006年)
业务流程控制 变异系数 Munter等人(2011)
bp_标记 血压幅度(峰值和低谷) Munter等人(2011)
bp_范围 血压范围 Levitan等人(2013)
业务流程视图 连续变化 Munter等人(2011)
bp_睡眠_指标 血压睡眠指标 (多个-请参阅文档)
bp_阶段 血压分期分类 美国心脏协会
bp_状态 统计汇总汇总 不适用
dip_calc(倾斜_倾斜) 夜间浸渍百分比和分类 Okhubo等人al(1997)

以下可视化功能目前通过提供这个英国石油公司包裹:

功能 可视化类型
bp_散射 BP阶段散点图
bp_ts_plots(bp_ts_图) 时间序列图
业务流程历史记录 BP阶段直方图
dip类plot 倾斜%类别图
bp_报告 业务伙伴汇总的可导出报告

示例-HYPNOS数据

有两个主要步骤涉及英国石油公司包裹:数据处理步骤和功能/分析步骤。

  1. 将数据加载并处理到一个新的可用数据帧中,以便进一步使用分析处理_数据功能
#开发工具::install_github(“johnschwenck/bp”)
图书馆(bp)

##加载bp_催眠
数据(bp_催眠)

##流程bp_催眠
催眠_程序<- 处理_数据(bp_催眠,
                     丁苯橡胶= “系统”
                     数据库处理程序= “离散”
                     日期_时间= '日期时间'
                     小时= “小时”
                     聚丙烯= “PP”
                     地图= “MaP”
                     快速定位程序= “卢比”
                     标识= “id”
                     访问= '访问'
                     尾流,尾迹= “唤醒”)
#>2个值超过DUL或DLL阈值,被强制为NA

注:这个处理_数据函数是对提供的数据列名的大小写不敏感。对于此示例中,即使原始列名“SYST”存在于bp_催眠,“syst”仍然是可以接受的名称到如图所示的功能。为了强调,上面所有的列名故意使用错误的大写字母输入。

SBP公司数据库管理员必须为任何其他功能正常工作。

  1. 使用新处理的催眠_程序,我们现在可以计算各种指标。现在包括生的bp_催眠数据集已处理为催眠_程序,我们现在可以依赖这个新的数据帧了计算各种指标和可视化。计算夜间浸泡分类如下所示,仅使用两名受试者(受试者70417和70435)进行比较:
dip_calc(倾斜_倾斜)(催眠_程序,子j= c(c)(7041770435))
#> [[1]]
#>#A tibble:8 x 6
#>#组:ID,VISIT[4]
#>ID VISIT WAKE avg_SBP avg_DBP N(ID访客唤醒平均值_SBP平均值_DBP N)
#><fct><fct>
#> 1 70417 1     0        123.    六十点五一零
#> 2 70417 1     1        128     66.6    20
#> 3 70417 2     0        136.    六十点五八
#> 4 70417 2     1        136.    六十五点六一七
#> 5 70435 1     0        106.    63       6
#> 6 70435 1     1        129.    八十二点一二三
#> 7 70435 2     0        136.    七十九点二九
#> 8 70435 2     1        123.    七十二点五二零
#>
#> [[2]]
#>#A tibble:4 x 6
#>#组:ID[2]
#>ID访问dip_sys class_sys dip_dias class_dias
#><fct><fct><dbl><chr><dbl><chr
#>1 70417 1 0.0359非铲斗0.0916非铲斗
#>2 70417 2-0.00450反向0.0784非铲斗
#>3 70435 1 0.179铲斗0.233极限
#>4 70435 2-0.104反向-0.0927反向

就统计指标而言bp_状态功能将许多可变性和中心指标聚合到一个表中这使得将不同的度量值进行比较方便。让我们假设对于这个示例,我们想进一步通过他们的BP_分类:我们会包括add_groups=“BP_CLASS”作为附加参数(请注意,资本化并不重要)。

bp_状态(催眠_程序,主题= c(c)(7041770435),添加组(_G)= “bp_class”bp_类型= 1)
#>#A tibble:35 x 16
#>#组:ID、VISIT、WAKE[8]
#>ID N VISIT WAKE BP_CLASS SBP_平均SBP_med SD ARV SV CV
#><fct><int><fct>
#>1 70417 5 1 0正常114。114  3.51  3.5   4.74 3.07
#>2 70417 3 1 0高架125。125  1.53  1.5   1.58 1.22
#>3 70417 2 1 0 ISH-S2 143 143 2.83 4 4 1.98
#>4 70417 3 1 1正常116。118  4.04  3.5   4.95 3.49
#>5 70417 9 1 1高架125。124  2.54  3.25  3.71 2.03
#>6 70417 6 1 1 ISH-S1 134。133  2.42  2.4   2.83 1.81
#>7 70417 2 1 1 ISH-S2 144 144 1.41 2 2 0.982
#>8 70417 3 2 0高架124 122 4.36 4 5 3.52
#>9 70417 2 0 ISH-S1 133 1.41 2 2 1.06
#>10 70417 3 2 0 ISH-S2 151 151 0 0 0 0
#> # ... 还有25行和5个变量:SBP_max<dbl>、SBP_min<dbl>,
#>#SBP_range<dbl>,峰值<dbl=,波谷<dbl:

下面是一个示例bp_散射的函数主题70417:

bp_散射(催眠_程序,主题= 70417)