每年的药品剂量的可视化的方法

99 阅读2分钟

在过去的几年里,我经常服用泮托拉唑--一种质子泵抑制剂(PPI)--来控制我的胃酸反流。然而,众所周知,定期使用PPI会导致各种健康问题。我不会在这里列举这些问题,但相信我,当我说这些问题很多的时候。去年,我决定慢慢戒掉PPI,在任何一天都只限制在20毫克的剂量,并交替使用或完全放弃一些日子。我在手机上的谷歌日历中记录了我的剂量信息,以便我以后可以处理。现在我有了一年多的数据,我想把它可视化,以了解我的用药习惯。

把谷歌日历的数据做成可用的格式是第一个挑战。值得庆幸的是,我找到了一个在线服务--calendarlabs.com--它允许我轻松地下载CSV格式的数据。

我心目中的那种可视化就像下面的热图。虽然我的可视化只显示我拍摄PPI的那一天的填充方块。

由于我想用R来实现数据的可视化,我发现最接近的软件包是Paul Bleicher的日历热图代码。

生成图表的第一步是读取包含每日剂量数据的CSV文件。

> pan.data <- read.csv("gcal-p.csv")

# A few sample rows from the data file
> head(pan.data)
  dosage      date
1      1  2018-1-1
2      1  2018-1-4
3      1  2018-1-7
4      1  2018-1-8
5      1  2018-1-9
6      1 2018-1-10

我们需要 "日历 "代码来生成图表。我们可以在当前目录下下载并使用它,或者直接使用在线版本。

> source("https://raw.githubusercontent.com/iascchen/VisHealth/master/R/calendarHeat.R")

该代码将加载额外的库--lattice、grid和chron,所以要确保这些库已经安装。

现在加载了数据后,我们只需运行日历代码即可生成如下所示的年度图表。

> calendarHeat(pan.data$date, pan.data$dosage, varname="Pan 20 dosage",color = "w2b")

calendarHeat()函数需要以下选项。

date: 需要绘制数据的日期。
values:与这些日期相关的值。
color: 调色板。默认是r2g(红到绿)。其他预定义的选项是rbb(红到蓝)和w2b(白到蓝)。你可以通过定义一个矢量创建你自己的调色板,如下图所示。
ncolors。热图的颜色数量
varname: 图表的标题

 calendarHeat(pan.data$date, pan.data$dosage, varname="Pan 20 dosage",color = "r2g")

网格的颜色可以从给定的三种风格中改变。

r2b #红变蓝
r2g #红变绿
w2b #白变蓝
g2r #绿变红