Estoy trabajando con médicos en un proyecto para monitorear el cumplimiento de la dosis adecuada de antibióticos. Para rastrear la proporción de eventos que no cumplen, a los médicos les gusta usarGráficos P

Me gustaría generar un P-Chart con3 líneas de límite (correspondiente a 1, 2 y 3 SD) encima y debajo de la línea central. No he encontrado una manera de hacer esto. También me gustaría que la trama tenga varios descansos que separen los datos en varios períodos de tiempo, lo que puedo hacer en el paquete qicharts pero no en otros paquetes.

Hay varios paquetes para R para generar Gráficos P. El que más me gusta es qicharts. El P-Chart estándar de qicharts, y todos los demás paquetes que he visto, genera un gráfico con una Línea central y un Límite de control superior y un Límite de control inferior a +3 y -3 SD desde la línea central.

Me gustaría descubrir cómo generar líneas de control SD +1, +2 y -1, -2 adicionales en el mismo diagrama. Alguna opción como

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

Aquí está el código, modificado der-proyectos, para generar datos, crear el gráfico e incluir dos saltos:

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