xgrsr.crit(spc)
xgrsr.crit()所属R语言包：spc
Compute alarm thresholds for Shiryaev-Roberts schemes
计算的报警阈值Shiryaev - 罗伯茨计划
译者：生物统计家园网 机器人LoveR
描述----------Description----------
Computation of the alarm thresholds (alarm limits)
计算的报警阈值（报警限值）
用法----------Usage----------
xgrsr.crit(k, L0, mu0 = 0, zr = 0, hs = NULL, sided = "one", r = 30)
参数----------Arguments----------
参数：k
reference value of the Shiryaev-Roberts scheme.
参考值的Shiryaev - 罗伯茨方案。
参数：L0
in-control ARL.
在控制ARL。
参数：mu0
in-control mean.
在控制的意思。
参数：zr
reflection border to enable the numerical algorithms used here.
反射边界，使这里的数值算法。
参数：hs
so-called headstart (give fast initial response). If hs=NULL, then the classical headstart -Inf is used (corresponds to 0 for the non-log scheme).
所谓的领先地位提供快速初始响应。如果hs=NULL，然后古典领先地位的Inf使用（对应于0的非对数方案）。
参数：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.crit determines the alarm threshold (alarm limit) for given in-control ARL L0 by applying secant rule and using xgrsr.arl().
xgrsr.crit决定报警阈值（报警下限）对于给定的控制ARL L0，应用割线规则，使用xgrsr.arl()。
值----------Value----------
Returns a single value which resembles the alarm limit
返回一个值，类似于报警极限
（作者）----------Author(s)----------
Sven Knoth
参考文献----------References----------
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----------
xgrsr.arl for zero-state ARL computation.
xgrsr.arl零状态ARL计算。
实例----------Examples----------
## Table 4 from Moustakides et al. (2009)[＃表4从Moustakides等。 （2009年）]
## original values are[＃原始值]
# gamma/L0 A/exp(g)[gamma/L0 A / EXP（G）]
# 50 28.02[50 28.02]
# 100 56.04[100 56.04]
# 500 280.19[500 280.19]
# 1000 560.37[1000 560.37]
# 5000 2801.75[5000 2801.75]
# 10000 5603.7[10000 5603.7]
theta <- 1
zr <- -6
r <- 100
Lxgrsr.crit <- Vectorize("xgrsr.crit", "L0")
L0s <- c(50, 100, 500, 1000, 5000, 10000)
gs <- Lxgrsr.crit(theta/2, L0s, zr=zr, r=r)
data.frame(L0s, gs, A=round(exp(gs), digits=2))
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注：
注1：为了方便大家学习，本文档为生物统计家园网机器人LoveR翻译而成，仅供个人R语言学习参考使用，生物统计家园保留版权。
注2：由于是机器人自动翻译，难免有不准确之处，使用时仔细对照中、英文内容进行反复理解，可以帮助R语言的学习。
注3：如遇到不准确之处，请在本贴的后面进行回帖，我们会逐渐进行修订。
|
发表于 2012-9-30 14:33:59
