Skip to main content

Browse the archive

Revision 63c8e8a453ea587001e2438a8ce51cf0e1e1675c authored by Charles J. Geyer on 23 March 2009, 00:00:00 UTC, committed by Gabor Csardi on 23 March 2009, 00:00:00 UTC 
version 0.7-7
1 parent b524c08
Raw File
fampoi.R 

 library(aster)

 ifam <- fam.poisson()

 mu <- seq(0.1, 3.0, 0.2)
 theta <- log(mu)

 zeroth <- double(length(theta))
 first <- double(length(theta))
 second <- double(length(theta))

 for (i in seq(along = theta)) {
    zeroth[i] <- famfun(ifam, 0, theta[i])
    first[i] <- famfun(ifam, 1, theta[i])
    second[i] <- famfun(ifam, 2, theta[i])
 }

 all.equal(zeroth, mu)
 all.equal(first, mu)
 all.equal(second, mu)
back to top