set.seed(1)#NB: the seed was not set for the table in the publication
table1_data<-lapply(c(.1,1,2),function(xi){
message(paste0("Run simulations for xi=",xi));Simulation_data(popmodelfunction = model.Pareto.bernstrat,
theta=4,
xi=xi,
conditionalto=list(N=10000,sampleparam=list(tauh=c(0.01,0.1))))})
table1<-lapply(table1_data,simulation.summary)
table1
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.