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