102 <strong>np</strong>regbwsubsetna.actionxdatydatbwsenvironment from which the function is called.an optional vector specifying a subset of observations to be used in the fittingprocess.a function which indicates what should happen when the data contain NAs. <strong>The</strong>default is set by the na.action setting of options, and is na.fail if that isunset. <strong>The</strong> (recommended) default is na.omit.a p-variate data frame of regressors which bandwidth selection will be performed.<strong>The</strong> datatypes may be continuous, discrete (unordered and orderedfactors), or some combination thereof.a one (1) dimensional numeric or integer vector of dependent data, each elementi corresponding to each observation (row) i of xdat.a bandwidth specification. This can be set as a rbandwidth object returnedfrom a previous invocation, or as a vector of bandwidths, with each elementi corresponding to the bandwidth for column i in xdat. In either case, thebandwidth supplied will serve as a starting point in the numerical search foroptimal bandwidths. If specified as a vector, then additional arguments willneed to be supplied as necessary to specify the bandwidth type, kernel types,selection methods, and so on. This can be left unset.... additional arguments supplied to specify the bandwidth type, kernel types, selectionmethods, and so on, detailed below.regtypebwmethodbwscalinga character string specifying which type of kernel regression estimator to use.lc specifies a local-constant estimator (Nadaraya-Watson) and ll specifies alocal-linear estimator. Defaults to lc.which method to use to select bandwidths. cv.aic specifies expected Kullback-Leibler cross-validation (Hurvich, Simonoff, and Tsai (1998)), and cv.ls specifiesleast-squares cross-validation. Defaults to cv.ls.a logical value that when set to TRUE the supplied bandwidths are interpretedas ‘scale factors’ (c j ), otherwise when the value is FALSE they are interpretedas ‘raw bandwidths’ (h j for continuous datatypes, λ j for discrete datatypes). Forcontinuous datatypes, c j and h j are related by the formula h j = c j σ j n −1/(2P +l) ,where σ j is the standard deviation of continuous variable j, n the number of observations,P the order of the kernel, and l the number of continuous variables.For discrete datatypes, c j and h j are related by the formula h j = c j n −2/(2P +l) ,where here j denotes discrete variable j. Defaults to FALSE.bwtype character string used for the continuous variable bandwidth type, specifying thetype of bandwidth to compute and return in the bandwidth object. Defaultsto fixed. Option summary:fixed: compute fixed bandwidthsgeneralized_nn: compute generalized nearest neighborsadaptive_nn: compute adaptive nearest neighborsbandwidth.computea logical value which specifies whether to do a numerical search for bandwidthsor not. If set to FALSE, a rbandwidth object will be returned with bandwidthsset to those specified in bws. Defaults to TRUE.ckertypecharacter string used to specify the continuous kernel type. Can be set as gaussian,epanechnikov, or uniform. Defaults to gaussian.
<strong>np</strong>regbw 103ckerorderukertypeokertypenmultireminitmaxftoltolsmallnumeric value specifying kernel order (one of (2,4,6,8)). Kernel order specifiedalong with a uniform continuous kernel type will be ignored. Defaults to2.character string used to specify the unordered categorical kernel type. Can beset as aitchisonaitken or liracine. Defaults to aitchisonaitken.character string used to specify the ordered categorical kernel type. Can be setas wangvanryzin or liracine. Defaults to wangvanryzin.integer number of times to restart the process of finding extrema of the crossvalidationfunction from different (random) initial points. Defaults to min(5,ncol(xdat)).a logical value which when set as TRUE the search routine restarts from locatedminima for a minor gain in accuracy. Defaults to TRUE.integer number of iterations before failure in the numerical optimization routine.Defaults to 10000.tolerance on the value of the cross-validation function evaluated at located minima.Defaults to 1.19e-07 (FLT_EPSILON).tolerance on the position of located minima of the cross-validation function.Defaults to 1.49e-08 (sqrt(DBL_EPSILON)).a small number, at about the precision of the data type used. Defaults to 2.22e-16 (DBL_EPSILON).Details<strong>np</strong>regbw implements a variety of methods for choosing bandwidths for multivariate (p-variate)regression data defined over a set of possibly continuous and/or discrete (unordered, ordered) data.<strong>The</strong> approach is based on Li and Racine (2003) who employ ‘generalized product kernels’ thatadmit a mix of continuous and discrete datatypes.<strong>The</strong> cross-validation methods employ multivariate numerical search algorithms (direction set (Powell’s)methods in multidimensions).Bandwidths can (and will) differ for each variable which is, of course, desirable.Three classes of kernel estimators for the continuous datatypes are available: fixed, adaptive nearestneighbor,and generalized nearest-neighbor. Adaptive nearest-neighbor bandwidths change witheach sample realization in the set, x i , when estimating the density at the point x. Generalizednearest-neighbor bandwidths change with the point at which the density is estimated, x. Fixedbandwidths are constant over the support of x.<strong>np</strong>regbw may be invoked either with a formula-like symbolic description of variables on whichbandwidth selection is to be performed or through a simpler interface whereby data is passed directlyto the function via the xdat and ydat parameters. Use of these two interfaces is mutuallyexclusive.Data contained in the data frame xdat may be a mix of continuous (default), unordered discrete(to be specified in the data frame xdat using factor), and ordered discrete (to be specified in thedata frame xdat using ordered). Data can be entered in an arbitrary order and data types will bedetected automatically by the routine (see <strong>np</strong> for details).Data for which bandwidths are to be estimated may be specified symbolically. A typical descriptionhas the form dependent data ~ explanatory data, where dependent data is
- Page 1 and 2:
The np PackageFebruary 16, 2008Vers
- Page 3 and 4:
Italy 3Examplesdata("cps71")attach(
- Page 5 and 6:
wage1 5Examplesdata("oecdpanel")att
- Page 7 and 8:
gradients 7## S3 method for class '
- Page 9 and 10:
np 9A variety of bandwidth methods
- Page 11 and 12:
npcmstest 11npcmstestKernel Consist
- Page 13 and 14:
npcmstest 13ReferencesAitchison, J.
- Page 15 and 16:
npcdens 15npcdensKernel Conditional
- Page 17 and 18:
npcdens 17Valuenpcdens returns a co
- Page 19 and 20:
npcdens 19# Gaussian kernel (defaul
- Page 21 and 22:
npcdens 21# (1993) (see their descr
- Page 23 and 24:
npcdensbw 23fit
- Page 25 and 26:
npcdensbw 25na.actionxdatydatbwspro
- Page 27 and 28:
npcdensbw 27data. The approach is b
- Page 29 and 30:
npconmode 29# depending on the spee
- Page 31 and 32:
npconmode 31tydatexdateydata one (1
- Page 33 and 34:
npconmode 33lwt,family=binomial(lin
- Page 35 and 36:
npudens 35npudensKernel Density and
- Page 37 and 38:
npudens 37Author(s)Tristen Hayfield
- Page 39 and 40:
npudens 39# EXAMPLE 1 (INTERFACE=DA
- Page 41 and 42:
npudensbw 41library("datasets")data
- Page 43 and 44:
npudensbw 43bwsa bandwidth specific
- Page 45 and 46:
npudensbw 45fvalobjective function
- Page 47 and 48:
npudensbw 47# previous examples.bw
- Page 49 and 50:
npudensbw 49# previous examples.bw
- Page 51 and 52: npksum 51Argumentsformuladatanewdat
- Page 53 and 54: npksum 53Usage IssuesIf you are usi
- Page 55 and 56: npksum 55# the bandwidth object its
- Page 57 and 58: npksum 57ss
- Page 59 and 60: npplot 59plot.behavior = c("plot","
- Page 61 and 62: npplot 61xtrim = 0.0,neval = 50,com
- Page 63 and 64: npplot 63xdatydatzdatxqyqzqxtrimytr
- Page 65 and 66: npplot 65DetailsValuenpplot is a ge
- Page 67 and 68: npplot 67year.seq
- Page 69 and 70: npplot 69# npplot(). When npplot()
- Page 71 and 72: npplreg 71## S3 method for class 'c
- Page 73 and 74: npplreg 73residR2MSEMAEMAPECORRSIGN
- Page 75 and 76: npplreg 75# Plot the regression sur
- Page 77 and 78: npplregbw 77and dependent data), an
- Page 79 and 80: npplregbw 79Detailsnpplregbw implem
- Page 81 and 82: npplregbw 81x2
- Page 83 and 84: npqcmstest 83npqcmstestKernel Consi
- Page 85 and 86: npqcmstest 85Author(s)Tristen Hayfi
- Page 87 and 88: npqreg 87ftol = 1.19209e-07,tol = 1
- Page 89 and 90: npqreg 89Li, Q. and J.S. Racine (20
- Page 91 and 92: npreg 91Usagenpreg(bws, ...)## S3 m
- Page 93 and 94: npreg 93residR2MSEMAEMAPECORRSIGNif
- Page 95 and 96: npreg 95summary(model)# Use npplot(
- Page 97 and 98: npreg 97# - this may take a few min
- Page 99 and 100: npreg 99# then a noisy samplen
- Page 101: npregbw 101## S3 method for class '
- Page 105 and 106: npregbw 105ReferencesAitchison, J.
- Page 107 and 108: npsigtest 107bw
- Page 109 and 110: npsigtest 109Author(s)Tristen Hayfi
- Page 111 and 112: npindex 111Usagenpindex(bws, ...)##
- Page 113 and 114: npindex 113MAEMAPECORRSIGNif method
- Page 115 and 116: npindex 115x2
- Page 117 and 118: npindex 117# x1 is chi-squared havi
- Page 119 and 120: npindexbw 119# plotting via persp()
- Page 121 and 122: npindexbw 121methodnmultithe single
- Page 123 and 124: npindexbw 123allows one to deploy t
- Page 125 and 126: npscoef 125x1
- Page 127 and 128: npscoef 127Valueeydatezdaterrorsres
- Page 129 and 130: npscoef 129# We could manually plot
- Page 131 and 132: npscoefbw 131optim.abstol,optim.max
- Page 133 and 134: npscoefbw 133optim.maxattemptsmaxim
- Page 135 and 136: npscoefbw 135ReferencesAitchison, J
- Page 137 and 138: uocquantile 137uocquantileCompute Q
- Page 139: INDEX 139npconmode, 29npindex, 110n