Nothing
cat("# numeric test\n")
library("qgcomp")
# are results at a given seed numerically stable across versions?
set.seed(50)
N=50
dat <- data.frame(time=(tmg <- pmin(.1,rweibull(N, 10, 0.1))),
d=1.0*(tmg<0.1), x1=runif(N), x2=runif(N), z=runif(N))
expnms=paste0("x", 1:2)
##### binomial
set.seed(123123)
f0 = d ~ x1 + x2
obj0a <- qgcomp(f0, expnms = expnms, data = dat, B=2, family=binomial())
obj0b <- qgcomp(f0, expnms = expnms, data = dat, family=binomial())
stopifnot(all.equal(
coef(obj0a),
c(-0.63497247,0.06082349 ),
check.names=FALSE, tolerance = 1e-4))
stopifnot(all.equal(
coef(obj0b),
c(-0.63497247,0.06082349 ),
check.names=FALSE, tolerance = 1e-4))
#' \dontrun{
#' ##### survival
#' set.seed(123123)
#' f1 = survival::Surv(tmg, d) ~ x1 + x2
#' obj1a <- qgcomp(f1, expnms = expnms, data = dat)
#' obj1b <- qgcomp.cox.boot(f1, expnms = expnms, B=2, data = dat)
#'
#' stopifnot(all.equal(
#' coef(obj1a),
#' c(0.02461938),
#' check.names=FALSE, tolerance = 1e-4))
#' stopifnot(all.equal(
#' coef(obj1b),
#' c(0.02654175),
#' check.names=FALSE, tolerance = 1e-4))
#'
#'
#' ##### zi
#' set.seed(123123)
#' f2 = d ~ x1 + x2 | x1 + x2
#' tm = pscl::zeroinfl(f2 , data=dat, dist="poisson")
#'
#' obj2a <- qgcomp.zi.noboot(f2, expnms = expnms, data = dat, dist="poisson")
#' obj2a <- qgcomp(f2, expnms = expnms, data = dat, family=poisson())
#' obj2b <- qgcomp.zi.boot(f2, expnms = expnms, data = dat, dist="poisson", B = 2)
#'
#' stopifnot(all.equal(
#' coef(obj2a),
#' list(count=c( -0.42823824,-0.05519007 ),
#' zero=c(9.057221, -40.356026)),
#' check.names=FALSE, tolerance = 1e-4))
#' stopifnot(all.equal(
#' coef(obj2b),
#' list(count=c( -0.43554489,-0.05279048 ),
#' zero=c(7.369426 , -18.411273)),
#' check.names=FALSE, tolerance = 1e-4))
#' }
cat("\n# done\n")
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.