plot.h.ucv: Plot for Unbiased Cross-Validation

In kedd: Kernel Estimator and Bandwidth Selection for Density and Its Derivatives

Description

The plot.h.ucv function loops through calls to the h.ucv function. Plot for unbiased cross-validation function for 1-dimensional data.

Usage

## S3 method for class 'h.ucv'
plot(x, seq.bws=NULL, ...)
## S3 method for class 'h.ucv'
lines(x,seq.bws=NULL, ...)

Arguments

x	object of class h.ucv (output from h.ucv).
seq.bws	the sequence of bandwidths in which to compute the unbiased cross-validation function. By default, the procedure defines a sequence of 50 points, from 0.15*hos to 2*hos (Over-smoothing).
...	other graphics parameters, see par in package "graphics".

Value

Plot of 1-d unbiased cross-validation function are sent to graphics window.

kernel	name of kernel to use.
deriv.order	the derivative order to use.
seq.bws	the sequence of bandwidths.
ucv	the values of the unbiased cross-validation function in the bandwidths grid.

Author(s)

Arsalane Chouaib Guidoum [e-mail]

See Also

h.ucv.

Examples

par(mfrow=c(2,1))
plot(h.ucv(trimodal,deriv.order=0),seq.bws=seq(0.06,0.2,length=50))
plot(h.ucv(trimodal,deriv.order=1),seq.bws=seq(0.06,0.2,length=50))

Example output

$kernel
[1] "gaussian"

$deriv.order
[1] 0

$seq.bws
 [1] 0.06000000 0.06285714 0.06571429 0.06857143 0.07142857 0.07428571
 [7] 0.07714286 0.08000000 0.08285714 0.08571429 0.08857143 0.09142857
[13] 0.09428571 0.09714286 0.10000000 0.10285714 0.10571429 0.10857143
[19] 0.11142857 0.11428571 0.11714286 0.12000000 0.12285714 0.12571429
[25] 0.12857143 0.13142857 0.13428571 0.13714286 0.14000000 0.14285714
[31] 0.14571429 0.14857143 0.15142857 0.15428571 0.15714286 0.16000000
[37] 0.16285714 0.16571429 0.16857143 0.17142857 0.17428571 0.17714286
[43] 0.18000000 0.18285714 0.18571429 0.18857143 0.19142857 0.19428571
[49] 0.19714286 0.20000000

$ucv
 [1] -0.2302222 -0.2310167 -0.2316960 -0.2322752 -0.2327675 -0.2331840
 [7] -0.2335340 -0.2338254 -0.2340649 -0.2342581 -0.2344093 -0.2345227
[13] -0.2346012 -0.2346477 -0.2346646 -0.2346538 -0.2346173 -0.2345568
[19] -0.2344737 -0.2343698 -0.2342463 -0.2341047 -0.2339464 -0.2337728
[25] -0.2335852 -0.2333849 -0.2331733 -0.2329515 -0.2327208 -0.2324824
[31] -0.2322375 -0.2319871 -0.2317323 -0.2314742 -0.2312136 -0.2309514
[37] -0.2306886 -0.2304258 -0.2301637 -0.2299031 -0.2296446 -0.2293886
[43] -0.2291357 -0.2288863 -0.2286409 -0.2283998 -0.2281632 -0.2279316
[49] -0.2277050 -0.2274837

$kernel
[1] "gaussian"

$deriv.order
[1] 1

$seq.bws
 [1] 0.06000000 0.06285714 0.06571429 0.06857143 0.07142857 0.07428571
 [7] 0.07714286 0.08000000 0.08285714 0.08571429 0.08857143 0.09142857
[13] 0.09428571 0.09714286 0.10000000 0.10285714 0.10571429 0.10857143
[19] 0.11142857 0.11428571 0.11714286 0.12000000 0.12285714 0.12571429
[25] 0.12857143 0.13142857 0.13428571 0.13714286 0.14000000 0.14285714
[31] 0.14571429 0.14857143 0.15142857 0.15428571 0.15714286 0.16000000
[37] 0.16285714 0.16571429 0.16857143 0.17142857 0.17428571 0.17714286
[43] 0.18000000 0.18285714 0.18571429 0.18857143 0.19142857 0.19428571
[49] 0.19714286 0.20000000

$ucv
 [1] -0.3463442 -0.7239392 -1.0122201 -1.2297570 -1.3914786 -1.5094289
 [7] -1.5933020 -1.6508499 -1.6882123 -1.7101966 -1.7205197 -1.7220171
[13] -1.7168231 -1.7065222 -1.6922758 -1.6749246 -1.6550712 -1.6331452
[19] -1.6094533 -1.5842177 -1.5576045 -1.5297449 -1.5007497 -1.4707199
[25] -1.4397528 -1.4079466 -1.3754018 -1.3422224 -1.3085149 -1.2743884
[31] -1.2399522 -1.2053155 -1.1705852 -1.1358653 -1.1012552 -1.0668490
[37] -1.0327350 -0.9989947 -0.9657027 -0.9329264 -0.9007261 -0.8691550
[43] -0.8382591 -0.8080781 -0.7786451 -0.7499871 -0.7221257 -0.6950771
[49] -0.6688529 -0.6434602

kedd documentation built on May 2, 2019, 7:32 a.m.