1
我正在與一位項目的醫生合作,監測抗生素使用劑量的合規情況。來跟蹤是不符合事件的比例,醫生喜歡用P chartsR多條線的控制圖
我想以產生P-圖表與3限制線(對應於1,2,和3個SDS)上方和下方的中線。我還沒有找到辦法做到這一點。我還希望情節有幾次休息,將數據分成幾個時間段,我可以在qicharts包中進行,但不能在其他包中進行。
R有幾個包用於生成P圖表。我最喜歡的是qicharts。來自qicharts的標準P-Chart以及我見過的所有其他軟件包都會生成一條中心線,上限控制限和下限控制限,分別位於中心線+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')
我懷疑你需要真正修改一個包的源代碼來達到這個目的。在** qic.R **中,可從https://cran.r-project.org/web/packages/qicharts/index.html的** qicharts_0.5.1.tar.gz **獲得,第776-780行是可能是一個很好的開始 - 這個包在這裏計算限制。 – tluh
謝謝,但我希望有一個更簡單的方法 - 可能與另一個包或解決方法。 – user3072084