if (FALSE) {
data(pollution)
Hpar.hr <- list(mean.lambda=0, sd.lambda=3)
PNS.hr <- fExtDep(x=PNS, method="BayesianPPP", model = "HR",
Nsim=5e+4, Nbin=3e+4, Hpar=Hpar.hr, MCpar=0.35, seed=14342)
bic(PNS.hr)
}
Run the code above in your browser using DataLab