nCal (version 0.3-2)

get.abc: Compute a measure of distance between two curves.

Description

abc criterion is area between curves divivded by the width of t.range. S1 criterion between two curves is defined as the integrated squared distance divivded by the width of t.range. S2 is relative bias divivded by the width of t.range.

Usage

get.abc(p1, p2, t.range)
get.S1(p1, p2, t.range)
get.S2(p1, p2, t.range)
get.abs.dev(p1, p2, t.range, y.range)

Arguments

p1

a vector of coefficients specify the first curve

p2

a vector of coefficients specify the second curve

t.range

a vector two real numbers. The range of log standard concentrations. The first one is the minimum t and the second one is the maximum t

y.range

Value

• a real number

Details

p1 and p2 can be in either classical or gh parameterization.

Examples

library(RUnit)

checkEqualsNumeric(
    get.abc(p.eotaxin[1,], p.eotaxin[2,], t.range=log(c(0.51,1e4)))
    , 
    0.04151633
    , tolerance=1e-7)


checkEqualsNumeric(
    get.S1(p.eotaxin[1,], p.eotaxin[2,], t.range=log(c(0.51,1e4)))
    , 
    0.002697115
    , tolerance=1e-7)


checkEqualsNumeric(
    get.S2(p.eotaxin[1,], p.eotaxin[2,], t.range=log(c(0.51,1e4)))
    , 
    17.48958
    , tolerance=1e-7)


checkEqualsNumeric(
    get.abs.dev(p.eotaxin[1,], p.eotaxin[2,], t.range=log(c(0.51,1e4)), y.range=c(5,6))
    , 
    0.07903227
    , tolerance=1e-7)


checkException(
    get.abs.dev(p.eotaxin[1,], p.eotaxin[2,], t.range=log(c(0.51,1e4)), y.range=c(1,11))
    )