Calculate Total
calctot.RdCalculate Total
Examples
ff <- ks::kde.boundary(AllPlantsby2020[,2:3],
xmin=c(0,0),xmax=c(35,15))
probs1 <- quadratprobs(ff, quadratstart=c(0,0),quadratsize=5)
results <- table(replicate(100,calctot(15,
probs=probs1,
incprobs=BonThomasAIncProbsMat[15,],
qx=QuadratChange5_2021)))
results
#>
#> 0 1.18189339321593 2.36378678643186 3.5456801796478
#> 4 76 19 1
Hmisc::wtd.mean(as.numeric(names(results)), results)
#> [1] 1.382815
sqrt(Hmisc::wtd.var(as.numeric(names(results)), results))
#> [1] 0.5832561