context("Typical session")
test_that("Typical session can be run.", {
set.seed(0)
neg=100/365
off.r=5
w.shape=10
w.scale=0.1
pi=0.25
dateT=2008
expect_silent(simu <- simulateOutbreak(neg=neg,pi=pi,off.r=off.r,w.shape=w.shape,
w.scale=w.scale,dateStartOutbreak=2005,dateT=dateT))
expect_is(simu,'ctree')
expect_silent(plotCTree(simu))
expect_is(capture_output(print(simu)),'character')
expect_silent(plot(simu))
expect_silent(ttree<-extractTTree(simu))
expect_is(ttree,'ttree')
expect_is(capture_output(print(ttree)),'character')
expect_silent(plot(ttree,w.shape,w.scale))
expect_silent(ptree<-extractPTree(simu))
expect_silent(plot(ptree))
expect_is(ptree,'ptree')
expect_is(capture_output(print(ptree)),'character')
expect_silent(p<-phyloFromPTree(ptree))
expect_is(p,'phylo')
expect_silent(ptree<-ptreeFromPhylo(p,dateLastSample=max(simu$ctree[,1])))
expect_equal(dateLastSample(simu),dateLastSample(ttree))
expect_equal(dateLastSample(ptree),dateLastSample(ttree))
expect_is(ptree,'ptree')
expect_is(capture_output(res<-inferTTree(ptree,mcmcIterations=100,w.shape=w.shape,w.scale=w.scale,dateT=dateT,updateOff.p = T,verbose=T)),'character')
expect_is(res,'resTransPhylo')
expect_is(capture_output(res<-inferTTree(ptree,mcmcIterations=100,w.shape=w.shape,w.scale=w.scale,dateT=dateT,updateOff.p = T,optiStart=1,verbose=T)),'character')
expect_is(res,'resTransPhylo')
expect_is(capture_output(print(res)),'character')
expect_is(capture_output(summary(res)),'character')
expect_silent(last<-extractCTree(res,100))
expect_is(last,'ctree')
expect_silent(plot(last))
expect_silent(plot(res))
expect_silent(plotTraces(res,burnin=0.5,extend=T))
expect_silent(mcmc<-as.mcmc.resTransPhylo(res))
expect_is(mcmc,'mcmc')
expect_is(capture_output(cons<-consTTree(res,debug=T)),'character')
expect_is(cons,'ttree')
expect_silent(sel<-selectTTree(res))
expect_is(sel,'numeric')
expect_silent(med<-medTTree(res))
expect_is(med,'ctree')
expect_silent(plot(med))
expect_silent(plotTTree(extractTTree(med),w.shape,w.scale))
expect_silent(plotTTree2(extractTTree(med),w.shape,w.scale))
expect_silent(mat<-computeMatWIW(res))
expect_is(mat,'matrix')
expect_silent(mat<-computeMatTDist(res))
expect_is(mat,'matrix')
expect_silent(a<-getIncidentCases(res,show.plot = T))
expect_is(a,'list')
expect_silent(a<-getGenerationTimeDist(res,show.plot = T))
expect_is(a,'matrix')
expect_silent(a<-getSamplingTimeDist(res,show.plot = T))
expect_is(a,'matrix')
expect_silent(a<-getInfectionTimeDist(res,k='1',show.plot = T))
expect_is(a,'numeric')
expect_silent(a<-getOffspringDist(res,k='1',show.plot = T))
expect_is(a,'numeric')
})
test_that("Inference with multiple trees runs", {
set.seed(0)
neg=100/365
off.r=5
w.shape=10
w.scale=0.1
pi=0.25
dateT=2008
simu=simulateOutbreak(neg=neg,pi=pi,off.r=off.r,w.shape=w.shape,
w.scale=w.scale,dateStartOutbreak=2005,dateT=dateT)
ptree=extractPTree(simu)
expect_is(capture.output(res<-infer_multittree_share_param(list(ptree),mcmcIterations=100,updateOff.p=T,w.shape=w.shape,w.scale=w.scale,dateT=dateT)),'character')
expect_is(res,'resTransPhylo')
expect_is(capture.output(res<-infer_multittree_share_param(list(ptree),mcmcIterations=100,updateOff.p=T,w.shape=w.shape,w.scale=w.scale,dateT=dateT,share=c('off.r','off.p','neg','pi'))),'character')
expect_is(res,'resTransPhylo')
expect_is(capture.output(res<-infer_multittree_share_param(list(ptree),mcmcIterations=100,updateOff.p=T,w.shape=w.shape,w.scale=w.scale,dateT=dateT,share=c('off.r','off.p','neg','pi')),verbose=T),'character')
expect_is(res,'resTransPhylo')
})
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.