Nothing
cat("# boot chooser test\n")
library("qgcomp")
library("splines")
library("survival")
dgm <- function(N){
dat <- data.frame(id=1:N)
dat <- within(dat, {
time=(tmg <- pmin(.1,rweibull(N, 10, 0.1)))
d=1.0*(tmg<0.1)
u = 0
x1 = runif(N)*4 + u
x2 = runif(N)*4 + u
x3 = runif(N)*4 + u
x4 = runif(N)*4 + u
x5 = runif(N)*4 + u
x6 = runif(N)*4 + u
y = rnorm(N, x1+x2, 2)
})
dat[,c("y", "time", "d", paste0("x", 1:6))]
}
set.seed(11232)
dat = dgm(100)
Xnm = c(paste0("x", 1:6))
f = y~x1 + x2 + x3 + I(x3^2)
qgcomp(f, expnms = c("x1", "x2"), data = dat)
qgcomp.noboot(f, expnms = c("x1", "x2"), data = dat)
# should cause errors
f = y~. + .^2
res = try(qgcomp(f, data = dat), silent=TRUE)
stopifnot(class(res)=="try-error")
fs = Surv(time,d)~. + .^2
res = try(qgcomp(fs, data = dat), silent=TRUE)
stopifnot(class(res)=="try-error")
f = y~.
res = try(fit1 <- qgcomp(f, expnms=Xnm, data = dat, parallel=TRUE), silent=TRUE)
stopifnot(class(res)=="try-error")
fs = Surv(time,d)~.
res = try(fit1 <- qgcomp(fs, expnms=Xnm, data = dat, B=5, MCsize=100), silent=TRUE)
stopifnot(class(res)=="try-error")
fs = Surv(time,d)~.
res = try(fit1 <- qgcomp(fs, expnms=Xnm, family=cox(), data = dat), silent=TRUE)
stopifnot(class(res)=="try-error")
#' \dontrun{
#' f = y~. + .^2
#' res = try(fit1 <- qgcomp(f, expnms=Xnm, data = dat), silent=TRUE)
#' stopifnot(inherits(res,"qgcompfit"))
#'
#' fs = Surv(time,d)~. + .^2
#' res = try(fit1 <- qgcomp(fs, expnms=Xnm, data = dat, B=5, MCsize=100), silent=TRUE)
#' stopifnot(inherits(res,"qgcompfit"))
#'
#' fs = Surv(time,d)~. + .^2
#' res = try(fit1 <- qgcomp(fs, expnms=Xnm, data = dat, B=5, MCsize=100, parallel=TRUE), silent=TRUE)
#' stopifnot(inherits(res,"qgcompfit"))
#'
#'
#'
#' # splines splines do work
#' f = y ~ x2 + x3 + x4 + x5 + x6 + splines::ns(x1, df=2)
#' res = try(qgcomp(f, data = dat), silent=TRUE)
#' stopifnot(class(res)=="try-error") # should give error that expnms not defined
#' res = try(fit1 <- qgcomp(f, expnms=Xnm, q=8, data = dat, deg=2), silent=TRUE)
#' stopifnot(inherits(res,"qgcompfit"))
#'
#'
#'
#' # splines splines do work
#' f = y ~ x2 + x3 + x4 + x5 + x6 + splines::ns(x1, df=2)
#' res = try(qgcomp(f, data = dat), silent=TRUE)
#' stopifnot(class(res)=="try-error") # should give error that expnms not defined
#' res = try(fit1 <- qgcomp(f, expnms=Xnm, q=8, data = dat, deg=2), silent=TRUE)
#' stopifnot(inherits(res,"qgcompfit"))
#'
#' # splines splines + bayes
#' f = y ~ x2 + x3 + x4 + x5 + x6 + splines::ns(x1, df=2)
#' res = try(fit2 <- qgcomp(f, expnms=Xnm, q=8, data = dat, deg=2, bayes=TRUE), silent=TRUE)
#' stopifnot(inherits(res,"qgcompfit"))
#'
#'
#' # indicator functions
#' f = y ~ factor(x1) + x2 + x3 + x4 + x5 + x6
#' res = try(qgcomp(f, data = dat), silent=TRUE)
#' stopifnot(class(res)=="try-error") # should give error that expnms not defined
#' res = try(fit1 <- qgcomp(f, expnms=Xnm, q=8, data = dat, deg=3), silent=TRUE)
#' stopifnot(inherits(res,"qgcompfit"))
#' }
cat("done")
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.