MGBT (version 1.0.4)

critK10: Single Grubbs--Beck Critical Values for 10-percent Test as used in Bulletin 17B

Description

Return the critical values at the 10-percent (\(\alpha_\mathrm{17B} = 0.10\)) significance level for the single Grubbs--Beck test as in Bulletin 17B (IACWD, 1982).

Usage

critK10(n)

Arguments

n

The number of observations.

Value

The critical value for sample size \(n\) unless it is outside the range \(10 \le n \le 149\) for which the critical value is NA.

References

Cohn, T.A., 2013--2016, Personal communication of original R source code: U.S. Geological Survey, Reston, Va.

Interagency Advisory Committee on Water Data (IACWD), 1982, Guidelines for determining flood flow frequency: Bulletin 17B of the Hydrology Subcommittee, Office of Water Data Coordination, U.S. Geological Survey, Reston, Va., 183 p.

See Also

CritK

Examples

# NOT RUN {
critK10(58)
#[1] 2.824

# }
# NOT RUN {
 # Modified slightly from TAC sources (Original has the # Not run:)
# KJRS() is the ".kngb()" function in TAC sources
n <- 10:149; KJRS <- function(n) -0.9043+3.345*sqrt(log10(n))-0.4046*log10(n)
result <- data.frame(n=n, Ktrue=sapply(n, critK10), # 17B single Grubbs--Beck
                          KJRS= sapply(n, KJRS   )) # name mimic of TAC sources 
# }
# NOT RUN {
# }
# NOT RUN {
 # Near verbatim from TAC sources, GGBK() does not work, issues a stop().
# KJRS() is the ".kngb()" function in TAC sources
n <- 10:149; KJRS <- function(n) -0.9043+3.345*sqrt(log10(n))-0.4046*log10(n)
result <- data.frame(n=n, Ktrue=sapply(n, critK10), # 17B single Grubbs--Beck
                          KJRS= sapply(n, KJRS   ), # name mimic of TAC sources
                          KTAC= sapply(n, GGBK   )) # name mimic of TAC sources
# }