## variation on example from help("example_reftable")
if (Infusion.getOption("example_maxtime")>2) {
blurred <- function(mu,s2,sample.size) {
s <- rnorm(n=sample.size,mean=mu,sd=sqrt(s2))
s <- exp(s/4)
return(c(mean=mean(s),var=var(s)))
}
set.seed(123)
dSobs <- blurred(mu=4,s2=1,sample.size=40)
parsp_j <- init_reftable(lower=c(mu=2.5, s2=0.25, sample.size=40),
upper=c(mu=5.2, s2=2.4, sample.size=40))
dsimuls <- add_reftable(,Simulate="blurred", parsTable=parsp_j,verbose=FALSE)
## Then instead of
#
# mufit <- project("mu",stats=c("mean","var"),data=dsimuls,verbose=FALSE)
# s2fit <- project("s2",stats=c("mean","var"),data=dsimuls,verbose=FALSE)
# dprojectors <- list(MEAN=mufit,VAR=s2fit)
#
## one can use
#
dprojectors <- def_projectors(dsimuls, pars=c("mu","s2"),
stats=c("mean","var"), npp = 1L,
projNames = c("MEAN", "VAR"))
## More speculatively, for non-default 'npp':
dprojectors <- def_projectors(dsimuls, pars=c("mu","s2"),
stats=c("mean","var"), npp = 2L)
}
Run the code above in your browser using DataLab