This function computes the gap statistic which is used to test for an outlier using Daniels method

Share:

Description

This function computes the gap statistic which is used to test for an outlier using Daniels method

Usage

1
gapstat(beta, pse)

Arguments

beta

input - vector of coefficients from saturated model fit to the data

pse

input - Lenth's PSE statistic calculated from the elements of beta

Value

returned gap statistic

Author(s)

John Lawson

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
## The function is currently defined as
function (beta, pse) 
{
    p <- length(beta)
    psehe <- pse
    sel <- beta >= 0
    betap <- beta[sel]
    betap <- sort(betap)
    betas <- betap[1]
    sel <- beta < 0
    betan <- beta[sel]
    nn <- length(betan)
    betan <- sort(betan)
    betal <- betan[nn]
    zl <- qnorm((nn - 0.375)/(p + 0.25))
    zs <- qnorm((nn + 1 - 0.375)/(p + 0.25))
    gap <- ((betas - betal)/psehe)/(zs - zl)
    return(gap)
  }

Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker.