Ich arbeite mit Ärzten an einem Projekt, um die Einhaltung der richtigen Dosierung von Antibiotika zu überwachen. Um den Anteil von Ereignissen zu erfassen, die nicht konform sind, verwenden Ärzte gerneP Charts

Ich möchte ein P-Chart mit @ generier 3 Grenzlinien (entspricht 1, 2 und 3 SDs) über und unter der Mittellinie. Ich habe keinen Weg gefunden, dies zu tun. Ich möchte auch, dass die Handlung mehrere Pausen hat, die die Daten in mehrere Zeiträume aufteilen, was ich im Paket qicharts aber nicht in anderen Paketen tun kann.

Es gibt mehrere Pakete für R zum Generieren von P-Diagrammen. Das, was mir am besten gefällt, sind Qicharts. Das Standard-P-Chart von qicharts und alle anderen Pakete, die ich gesehen habe, erzeugen einen Plot mit einer Mittellinie und einer oberen Kontrollgrenze und einer unteren Kontrollgrenze bei +3 und -3 SD von der Mittellinie.

Ich möchte herausfinden, wie zusätzliche +1, +2 und -1, -2 SD-Steuerlinien auf demselben Plot generiert werden. Einige Optionen wie

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

Hier ist der Code, geändert von r-projects, um Daten zu generieren, erstellen Sie das Diagramm und fügen Sie zwei Pausen ein:

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