Estou trabalhando com médicos em um projeto para monitorar a conformidade com a dose adequada de antibióticos. Para rastrear a proporção de eventos que não são compatíveis, os médicos gostam de usarGráficos P

Gostaria de gerar um gráfico P com3 linhas de limite (correspondente a 1, 2 e 3 SDs) acima e abaixo da linha central. Não encontrei uma maneira de fazer isso. Também gostaria que o gráfico tivesse várias quebras que separam os dados em vários períodos, o que posso fazer no pacote qicharts, mas não em outros pacotes.

Existem vários pacotes para R para gerar gráficos P. O que eu mais gosto é o qicharts. O P-Chart padrão do qicharts, e todos os outros pacotes que eu já vi, gera um gráfico com uma Linha Central e um Limite Superior de Controle e um Limite Inferior de Controle em +3 e -3 SD da linha central.

Gostaria de descobrir como gerar linhas de controle adicionais +1, +2 e -1, -2 SD no mesmo gráfico. Alguma opção como

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

Aqui está o código, modificado der-projetos, para gerar dados, crie o gráfico e inclua duas quebras:

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