R lines.saddle.distn 将鞍点近似添加到绘图中

说明

此函数将对应于鞍点密度或分布函数近似值的线添加到当前绘图中。

用法

## S3 method for class 'saddle.distn'
lines(x, dens = TRUE, h = function(u) u, J = function(u) 1, 
      npts = 50, lty = 1, ...)

参数

x	"saddle.distn" 类的对象(请参阅 saddle.distn.object 表示分布的鞍点近似。
dens	一个逻辑变量，指示是否应绘制鞍点密度(TRUE；默认值)或鞍点分布函数(FALSE)。
h	所需的变量的任何转换。它的第一个参数必须是执行近似的值，并且该函数必须进行矢量化。
J	当dens=TRUE 时，此函数指定可能需要的任何转换的雅可比行列式。 J 的第一个参数必须是执行近似的值，并且函数必须进行矢量化。如果提供了h，则还必须提供J，并且两者必须具有相同的参数列表。
npts	用于绘图的点数。这些点将在用于查找鞍点近似的点范围内均匀分布。
lty	要使用的线类型。
...	h 和 J 的任何其他参数。

细节

该函数使用smooth.spline 生成鞍点曲线。当dens=TRUE 时，样条曲线采用对数标度，而当dens=FALSE 时，样条曲线采用概率标度。

值

sad.d 隐形返回。

副作用

一条线将添加到当前绘图中。

例子

# In this example we show how a plot such as that in Figure 9.9 of
# Davison and Hinkley (1997) may be produced.  Note the large number of
# bootstrap replicates required in this example.
expdata <- rexp(12)
vfun <- function(d, i) {
     n <- length(d)
     (n-1)/n*var(d[i])
}
exp.boot <- boot(expdata,vfun, R = 9999)
exp.L <- (expdata - mean(expdata))^2 - exp.boot$t0
exp.tL <- linear.approx(exp.boot, L = exp.L)
hist(exp.tL, nclass = 50, probability = TRUE)
exp.t0 <- c(0, sqrt(var(exp.boot$t)))
exp.sp <- saddle.distn(A = exp.L/12,wdist = "m", t0 = exp.t0)

# The saddlepoint approximation in this case is to the density of
# t-t0 and so t0 must be added for the plot.
lines(exp.sp, h = function(u, t0) u+t0, J = function(u, t0) 1,
      t0 = exp.boot$t0)

参考

Davison, A.C. and Hinkley, D.V. (1997) Bootstrap Methods and Their Application. Cambridge University Press.

也可以看看

saddle.distn