modelfree (version 1.0-2)

bandwidth_bootstrap: Bootstrap estimate of bandwidth

Description

Finds bootstrap estimate of the optimal bandwidth h for binomial data in local polynomial estimation of psychometric function (PF) with guessing and lapsing rates specified in lims.

Usage

bandwidth_bootstrap( r, m, x, H, N, h0 = NULL, link = c( "logit" ), guessing = 0, lapsing = 0, K = 2, p = 1, ker = c( "dnorm" ), maxiter = 50, tol = 1e-6, method = c( "all" ) )

Arguments

r

number of successes in points x

m

number of trials in points x

x

stimulus levels

H

minimum and maximum values of bandwidth to be considered

N

number of bootstrap replications

h0

pilot bandwidth; if not specified, then the scaled plug-in bandwidth is used

link

name of the link function to be used; default is "logit"

guessing

guessing rate; default is 0

lapsing

lapsing rate; default is 0

K

power parameter for Weibull and reverse Weibull link; default is 2

p

order of the polynomial; default is 1

ker

kernel function for weights; default "dnorm"

maxiter

maximum number of iterations in Fisher scoring; default is 50

tol

tolerance level at which to stop Fisher scoring; default is 1e-6

method

loss function to be used in bandwidth: choose from: 'ISEeta', 'ISE', 'deviance'; default is "all"

Value

• h
  bootstrap bandwidth for the chosen method; if no method was specified, then it a list of three elements with entries corresponding to the estimated bandwidths on p-scale (h$pscale), on eta-scale (h$etascale) and for mean likelihood (h$deviance)

Examples

data("01_Miranda")
h<- bandwidth_bootstrap( example01$r, example01$m, example01$x, c( 0.1, 10 ), 10 )