inst/unitTests/runittwLogSumExp.R

test.one <- function(){
    x <- rnorm(100)
    r0 <- log(sum(exp(x)))
    #mtrace(twLogSumExp)
    res <- twLogSumExp(x)
    checkEqualsNumeric(r0,res)
    
    x2 <- -1000+x
    exp(x2) #zero
    log(mean(exp(x2)))
    res <- twLogMeanExp(x2)
    hist(x2)
    abline(v=res,col="red")
    abline(v=mean(x2),col="blue")
    checkTrue( res < max(x2))
    checkTrue( res > mean(x2))
}

Try the twMisc package in your browser

Any scripts or data that you put into this service are public.

twMisc documentation built on May 2, 2019, 6:11 p.m.