Я работаю с врачами над проектом по контролю за соблюдением правильной дозировки антибиотиков. Чтобы отслеживать долю событий, которые не соответствуют, врачи любят использоватьP диаграммы

Я хотел бы создать P-Chart с3 предельные линии (соответствует 1, 2 и 3 SD) выше и ниже центральной линии. Я не нашел способ сделать это. Я также хотел бы, чтобы на графике было несколько разрывов, которые разделяют данные на несколько периодов времени, что я могу сделать в пакете qicharts, но не в других пакетах.

Есть несколько пакетов для R для генерации P диаграмм. Больше всего мне нравятся qicharts. Стандартная P-диаграмма от qicharts и все другие пакеты, которые я видел, генерируют график с центральной линией, верхним контрольным пределом и нижним контрольным пределом в +3 и -3 SD от центральной линии.

Я хотел бы выяснить, как создать дополнительные линии управления SD, +1, +2 и -1, -2 на одном графике. Некоторые варианты, такие как

LimitLines = c(1, 2, 3) where the default is LimitlLines = 3

Вот код, модифицированный изг-проекты, чтобы сгенерировать данные, создать диаграмму и включить два перерыва:

# 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')