我正在尝试绘制医院中如何开出单一药物的图解。在这个虚拟数据库中,我在 2017 年 1 月 1 日之后遇到了 1000 名患者。
绘图的目的是了解该药物的给药模式:在接近入院、出院或患者住院期间是否更频繁/高剂量给药。
#Get_random_dates that we will use multiple times
gen_random_dates <- function(N, st, et) {
st <- as.POSIXct(as.Date(st))
et <- as.POSIXct(as.Date(et))
dt <- as.numeric(difftime(et,st,unit="sec"))
ev <- runif(N, 0, dt)
rt <- st + ev
return(rt)
}
#Generate admission and discharge dates
admission <- gen_random_dates(1000, "2017/01/01", "2017/01/10")
discharge <- gen_random_dates(1000, "2017/01/11", "2017/01/20")
patient <- sort(sample(1:1000, 1000))
patient_data <- data.frame(patient_ID = patient, admission_date = admission, discharge_date = discharge)
#Grow the database
patient_data <- patient_data[sort(sample(1000, 100000, replace=TRUE)), ]
#Medication admin date and dose
patient_data$admin_date <- gen_random_dates(100000, patient_data$admission_date, patient_data$discharge_date)
patient_data$admin_dose <- abs(as.integer(rnorm(100000, 50, 100)))
我尝试了这个 ggplot 函数,但它没有帮助我可视化该模式。
ggplot(patient_data, aes(x = admin_date, y = admin_dose)) +
xlab("Use of Drug in Patient Encounters") + ylab("Dose (mg)") +
geom_jitter()
ggplot