modelfree (version 1.0-2)

bandwidth_cross_validation: Cross-validation estimate of bandwidth

Description

Finds the cross-validation bandwidth for the local polynomial estimator of the psychometric function (PF) with guessing and lapsing rates specified in lims.

Usage

bandwidth_cross_validation( r, m, x, H, 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

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 cross-validation: choose from: 'ISEeta', 'ISE', 'deviance'; default is "all"

Value

• h
  cross-validation 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_cross_validation( example01$r, example01$m, example01$x, c( 0.1, 10 ) )