Nothing
start.time <- proc.time()
suppressMessages(library(cubfits, quietly = TRUE))
set.seed(1234)
.CF.AC$renew.iter <- 3
# .CF.CT$type.p <- "lognormal_bias"
# .CF.CONF$scale.phi.Obs <- FALSE
# .CF.CONF$estimate.bias.Phi <- TRUE
ex.train$phi.Obs <- ex.train$phi.Obs / mean(ex.train$phi.Obs)
ret.time <- system.time({
ret <- cubfits(ex.train$reu13.df, ex.train$phi.Obs, ex.train$y, ex.train$n,
nIter = 20,
verbose = TRUE, report = 5,
model = "roc", adaptive = "simple")
})
print(ret.time)
x <- rowMeans(do.call("cbind", ret$phi.Mat)[, 11:20])
y <- ex.train$phi.Obs
x <- log10(x / mean(x))
y <- log10(y / mean(y))
print(mean(x))
print(summary(lm(y ~ x))$r.squared)
# warning: iterations terminated because half-step sizes are very small
print(proc.time() - start.time)
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.