xgrsr.ad(spc)
xgrsr.ad()所属R语言包：spc
Compute steady-state ARLs of Shiryaev-Roberts schemes
计算稳态连串长度的Shiryaev - 罗伯茨计划
译者：生物统计家园网 机器人LoveR
描述----------Description----------
Computation of the steady-state Average Run Length (ARL)
计算的稳态平均运行长度（ARL）
用法----------Usage----------
xgrsr.ad(k, g, mu1, mu0 = 0, zr = 0, sided = "one", r = 30)
参数----------Arguments----------
参数：k
reference value of the Shiryaev-Roberts scheme.
参考值的Shiryaev - 罗伯茨方案。
参数：g
control limit (alarm threshold) of Shiryaev-Roberts scheme.
控制上限（报警阈值）Shiryaev - 罗伯茨计划。
参数：mu1
out-of-control mean.
控制的意思。
参数：mu0
in-control mean.
在控制的意思。
参数：zr
reflection border to enable the numerical algorithms used here.
反射边界，使这里的数值算法。
参数：sided
distinguish between one- and two-sided schemes by choosing "one" and"two", respectively. Currently only one-sided schemes are implemented.
区分和双面计划的选择"one"和"two"，分别。目前只是片面的计划实施。
参数：r
number of quadrature nodes, dimension of the resulting linear equation system is equal to r+1.
正交节点的数量，产生的线性方程组的维数等于r+1。
Details
详细信息----------Details----------
xgrsr.ad determines the steady-state Average Run Length (ARL) by numerically solving the related ARL integral equation by means of the Nystroem method based on Gauss-Legendre quadrature.
xgrsr.ad确定的稳态平均运行长度（ARL）相关的ARL的Nystroem方法基于高斯 - 勒让德正交积分方程的数值求解。
值----------Value----------
Returns a single value which resembles the steady-state ARL.
类似的稳态ARL返回单个值。
（作者）----------Author(s)----------
Sven Knoth
参考文献----------References----------
S. Knoth (2006), The art of evaluating monitoring schemes – how to measure the performance of control charts? S. Lenz, H. & Wilrich, P. (ed.), Frontiers in Statistical Quality Control 8, Physica Verlag, Heidelberg, Germany, 74-99.
G. Moustakides, A. Polunchenko, A. Tartakovsky (2009), Numerical comparison of CUSUM and Shiryaev-Roberts procedures for detectin changes in distributions, Communications in Statistics: Theory and Methods, to appear.
参见----------See Also----------
xewma.arl and xcusum-arl for zero-state ARL computation of EWMA and CUSUM control charts, respectively, and xgrsr.arl for the zero-state ARL.
xewma.arl和xcusum-arl零状态EWMA和CUSUM控制图的ARL计算，分别为，和xgrsr.arl零状态的ARL。
实例----------Examples----------
## Small study to identify appropriate reflection border to mimic[＃小的研究，以确定合适的反射边界模仿]
## unreflected schemes[＃未反射的计划]
k <- .5
g <- log(390)
zrs <- -(0:10)
ZRxgrsr.ad <- Vectorize(xgrsr.ad, "zr")
ads <- ZRxgrsr.ad(k, g, 0, zr=zrs)
data.frame(zrs, ads)
## Table 2 from Knoth (2006)[＃2 Knoth（2006年）]
## original values are[＃原始值]
# mu arl[亩ARL]
# 0 689[0 689]
# 0.5 30[0.5 30]
# 1 8.9[1 8.9]
# 1.5 5.1[1.5 5.1]
# 2 3.6[2 3.6]
# 2.5 2.8[2.5 2.8]
# 3 2.4[3 2.4]
#[]
k <- .5
g <- log(390)
zr <- -5 # see first example[第一个例子]
mus <- (0:6)/2
Mxgrsr.ad <- Vectorize(xgrsr.ad, "mu1")
ads <- round(Mxgrsr.ad(k, g, mus, zr=zr), digits=1)
data.frame(mus, ads)
## Table 4 from Moustakides et al. (2009)[＃表4从Moustakides等。 （2009年）]
## original values are[＃原始值]
# gamma A STADD/steady-state ARL[伽马一个STADD /稳态ARL]
# 50 28.02 4.37[50 28.02 4.37]
# 100 56.04 5.46[100 56.04 5.46]
# 500 280.19 8.33[500 280.19 8.33]
# 1000 560.37 9.64[1000 560.37 9.64]
# 5000 2801.75 12.79[5000 2801.75 12.79]
# 10000 5603.7 14.17[10000 5603.7 14.17]
Gxgrsr.ad <- Vectorize("xgrsr.ad", "g")
As <- c(28.02, 56.04, 280.19, 560.37, 2801.75, 5603.7)
gs <- log(As)
theta <- 1
zr <- -6
ads <- round(Gxgrsr.ad(theta/2, gs, theta, zr=zr, r=100), digits=2)
data.frame(As, ads)
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注：
注1：为了方便大家学习，本文档为生物统计家园网机器人LoveR翻译而成，仅供个人R语言学习参考使用，生物统计家园保留版权。
注2：由于是机器人自动翻译，难免有不准确之处，使用时仔细对照中、英文内容进行反复理解，可以帮助R语言的学习。
注3：如遇到不准确之处，请在本贴的后面进行回帖，我们会逐渐进行修订。
|
发表于 2012-9-30 14:33:45
