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