Я работаю с врачами над проектом по контролю за соблюдением правильных дозировок антибиотиков. Чтобы отслеживать долю несоответствующих событий, врачи предпочитают использовать P-диаграммы.
Я хочу создать P-диаграмму с тремя ограничительными линиями (соответствующими 1, 2 и 3 SD) выше и ниже центральной линии. Я не нашел способ сделать это. Я также хотел бы, чтобы на графике было несколько разрывов, которые разделяют данные на несколько периодов времени, что я могу сделать в пакете qicharts, но не в других пакетах.
Существует несколько пакетов R для создания P-диаграмм. Мне больше всего нравится qicharts. Стандартный P-Chart из qicharts и все другие пакеты, которые я видел, генерируют график с центральной линией и верхним контрольным пределом и нижним контрольным пределом на +3 и -3 SD от центральной линии.
Я хотел бы выяснить, как создать дополнительные контрольные линии +1, +2 и -1, -2 SD на одном графике. Какой-то вариант, например
LimitLines = c(1, 2, 3) where the default is LimitlLines = 3
Вот код, измененный с r-projects, чтобы сгенерировать данные, создать диаграмму и включить два разрыва:
# Setup parameters
m.beds <- 300
m.stay <- 4
m.days <- m.beds * 7
m.discharges <- m.days / m.stay
p.pu <- 0.08
# Simulate data
discharges <- rpois(24, lambda = m.discharges)
patientdays <- round(rnorm(24, mean = m.days, sd = 100))
n.pu <- rpois(24, lambda = m.discharges * p.pu * 1.5)
n.pat.pu <- rbinom(24, size = discharges, prob = p.pu)
week <- seq(as.Date('2014-1-1'),
length.out = 24,
by = 'week')
# Combine data into a data frame
d <- data.frame(week, discharges, patientdays,n.pu, n.pat.pu)
# Create a P-chart to measure the number of patients with pressure ulcers (n.pat.pu) each week (week) as a proportion of all discharges (discharges) with breaks one third (8) and two thirds (16) of the way through the data
qic(n.pat.pu,
n = discharges,
x = week,
data = d,
chart = 'p',
multiply = 100,
breaks = c(8,16),
main = 'Hospital acquired pressure ulcers (P chart)',
ylab = 'Percent patients',
xlab = 'Week')