demo/table1.R

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
DanielBonnery/pubBonneryBreidtCoquet2016 documentation built on May 6, 2019, 1:35 p.m.