R : Copyright 2006, The R Foundation for Statistical Computing Version 2.3.1 (2006-06-01) ISBN 3-900051-07-0 R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > > library(aster) Loading required package: trust > > ifam <- fam.bernoulli() > > p <- seq(0.1, 0.9, 0.1) > theta <- log(p) - log(1 - p) > > zeroth <- double(length(p)) > first <- double(length(p)) > second <- double(length(p)) > > for (i in seq(along = p)) { + zeroth[i] <- famfun(ifam, 0, theta[i]) + first[i] <- famfun(ifam, 1, theta[i]) + second[i] <- famfun(ifam, 2, theta[i]) + } > > all.equal(zeroth, - log(1 - p)) [1] TRUE > all.equal(first, p) [1] TRUE > all.equal(second, p * (1 - p)) [1] TRUE > >