"quadrats", "variance/population" rpoispp(25,win=owin(xrange=c(0,10),yrange=c(0,10)))   15% of quadrats sampled
225, 5.26547936692055
289, 6.35541388991654
361, 6.04057634141989
441, 6.39677776583233
529, 6.14111986823566
625, 5.73627896313843
729, 5.26696141331783
841, 6.74603229323467
961, 5.52194293512602
1089, 6.68607704593296
1225, 6.62173313199676
"quadrats", "variance/population" rpoispp(50,win=owin(xrange=c(0,10),yrange=c(0,10)))  15% of quadrats sampled
225, 6.90426797262287
289, 5.93712574474672
361, 5.69054454275825
441, 5.94236267771851
529, 5.79793223996722
625, 6.93873936093873
729, 6.14211049912364
841, 7.40233045854867
961, 5.6113432368888
1089, 6.38721680659351
1225, 6.15218655587803
"quadrats", "variance/population" rpoispp(function(x,y) {65*exp(-.01*x*y)},win=owin(xrange=c(0,10),yrange=c(0,10))) 15% of quadrats sampled
225, 12.626805860584
289, 10.7731801737102
361, 10.4681830189
441, 8.55636141882028
529, 7.73301091866841
625, 6.73681056410954
729, 8.39707888519799
841, 7.61243074579075
961, 6.02651863375618
1089, 7.0433516252464
1225, 6.72608969807015
"quadrats", "variance/population" houses = rpoispp(function(x,y) {30*exp(-.01*x*y)},win=owin(xrange=c(0,10),yrange=c(0,10))) 15% of quadrats sampled
225, 9.65385053000827
289, 6.23488370211011
361, 6.08291646127987
441, 8.52779561583964
529, 5.56813072248531
625, 5.67598465186832
729, 6.72241085359692
841, 6.49258174217286
961, 6.28399334209182
1089, 7.26188632291648
1225, 5.97714596082854
"quadrats", "variance/population" houses = rpoispp(function(x,y) {30*exp(-.01*x*y)},win=owin(xrange=c(0,10),yrange=c(0,10))) 15% of quadrats sampled
25, 19.453259800255
49, 18.5376218587603
81, 10.1344991703878
121, 11.8908929642853
169, 9.10359708015107
225, 8.83392916819609
289, 7.23049891811094
361, 7.42406479693757
"quadrats", "variance/population" rpoispp(25,win=owin(xrange=c(0,10),yrange=c(0,10))) 15% of quadrats sampled
25, 3.88841830134814
49, 5.95930215330414
81, 6.17631553903617
121, 7.02894352923651
169, 6.11820078343732
225, 5.05349108749279
289, 5.2453600268258
361, 5.15309355944328
"quadrats", "variance/population" rpoispp(50,win=owin(xrange=c(0,10),yrange=c(0,10)))  15% of quadrats sampled
16, 6.82767665097139
25, 3.80306254915077
36, 8.60579178095201
49, 7.24407188996162
64, 6.21445495192498
81, 7.11065050121424
100, 5.41684450247739
"quadrats", "variance/population" rpoispp(function(x,y) {30*exp(-.01*x*y)},win=owin(xrange=c(0,10),yrange=c(0,10))) 15% of quadrats sampled
16, 47.4054916397221
25, 22.4752773789684
36, 19.896522736861
49, 13.8761028705114
64, 13.14643764039
81, 12.0512572224728
100, 10.05325520091
"quadrats", "variance/population" rpoispp(function(x,y) {65*exp(-.01*x*y)},win=owin(xrange=c(0,10),yrange=c(0,10))) 15% of quadrats sampled
16, 84.902585790251
25, 41.2899653380262
36, 37.1720166660655
49, 30.310047185491
64, 15.8994522377162
81, 16.0521300622901
100, 13.2678460737787
"coverage", "variance/population" rpoispp(50,win=owin(xrange=c(0,10),yrange=c(0,10))) 400 quadrats
0.05, 15.04186468162
0.075, 13.5090129465037
0.1, 9.85178085598618
0.125, 8.97382189334219
0.15, 6.57741890231218
0.175, 4.99974816965592
0.2, 4.25989844954166
0.225, 3.49692120352781
0.25, 3.43727947259747
0.275, 2.9068252344945
0.3, 2.79104005804761
"coverage", "variance/population" rpoispp(25,win=owin(xrange=c(0,10),yrange=c(0,10))) 400 quadrats
0.05, 20.6237506648567
0.075, 13.0457006817482
0.1, 10.2194185778815
0.125, 7.54879547818079
0.15, 7.51504759182961
0.175, 4.7324347824819
0.2, 3.87006095842562
0.225, 3.5365780200909
0.25, 3.39775274476649
0.275, 2.76588758947559
0.3, 2.56702758570807
"coverage", "variance/population" rpoispp(function(x,y) {65*exp(-.01*x*y)},win=owin(xrange=c(0,10),yrange=c(0,10))) 400 quadrats
0.05, 28.7296610120297
0.075, 21.8892700786266
0.1, 14.3465188863068
0.125, 11.769560161222
0.15, 8.66751418496667
0.175, 8.26518848397094
0.2, 6.19811115163805
0.225, 5.25558811327276
0.25, 4.71280576931635
0.275, 3.52492286372912
0.3, 2.85793786241781
"coverage", "variance/population" rpoispp(function(x,y) {30*exp(-.01*x*y)},win=owin(xrange=c(0,10),yrange=c(0,10))) 400 quadrats
0.05, 20.1409590324504
0.075, 14.9291412672726
0.1, 10.2381114612611
0.125, 9.70688081245615
0.15, 6.75431502413872
0.175, 6.66698748804887
0.2, 5.25620409351033
0.225, 4.28986250916997
0.25, 3.52790003285191
0.275, 2.83387749652819
0.3, 2.8669397340933
(unless otherwise noted, quadrat 10 times more expensive than individual)
"quadrats", "cost" rpoispp(50,win=owin(xrange=c(0,10),yrange=c(0,10))) .15 coverage
225, 1079.81
289, 1158.13
361, 1275.695
441, 1397.345
529, 1523.165
625, 1680.93
729, 1827.805
841, 1995.265
961, 2176.765
1089, 2365.665
1225, 2578.185
"coverage", "cost" rpoispp(50,win=owin(xrange=c(0,10),yrange=c(0,10))) 400 quadrats
0.05, 445.4
0.075, 667.045
0.1, 892.54
0.125, 1117.525
0.15, 1336.505
0.175, 1559.46
0.2, 1778.51
0.225, 2006.67
0.25, 2228.155
0.275, 2449.515
0.3, 2677.215
"quadrats", "cost" rpoispp(25,win=owin(xrange=c(0,10),yrange=c(0,10))) .15 coverage
225, 725.185
289, 810.66
361, 919.59
441, 1042.135
529, 1170.15
625, 1325.16
729, 1470.805
841, 1639.83
961, 1821.195
1089, 2013.66
1225, 2222.135
"coverage", "cost" rpoispp(25,win=owin(xrange=c(0,10),yrange=c(0,10))) 400 quadrats
0.05, 327.26
0.075, 490.385
0.1, 654.285
0.125, 816.74
0.15, 981.52
0.175, 1145.675
0.2, 1309.62
0.225, 1472.66
0.25, 1636.865
0.275, 1798.235
0.3, 1963.37
"quadrats", "cost" rpoispp(function(x,y) {65*exp(-.01*x*y)},win=owin(xrange=c(0,10),yrange=c(0,10))) .15 coverage
225, 1126.565
289, 1203.74
361, 1316.48
441, 1440.425
529, 1571.94
625, 1722.965
729, 1866.85
841, 2043.46
961, 2220.615
1089, 2412.175
1225, 2624.485
"coverage", "cost" rpoispp(function(x,y) {65*exp(-.01*x*y)},win=owin(xrange=c(0,10),yrange=c(0,10))) 400 quadrats
0.05, 459.635
0.075, 689.145
0.1, 919.575
0.125, 1151.59
0.15, 1383.725
0.175, 1612.06
0.2, 1850.765
0.225, 2073.475
0.25, 2306.835
0.275, 2532.71
0.3, 2766.83
"quadrats", "cost" rpoispp(function(x,y) {30*exp(-.01*x*y)},win=owin(xrange=c(0,10),yrange=c(0,10))) .15 coverage
225, 704.85
289, 788.11
361, 899.98
441, 1021.87
529, 1150.68
625, 1304.955
729, 1452.23
841, 1622.085
961, 1802.415
1089, 1992.83
1225, 2203.69
"coverage", "cost" rpoispp(function(x,y) {30*exp(-.01*x*y)},win=owin(xrange=c(0,10),yrange=c(0,10))) 400 quadrats
0.05, 320.18
0.075, 482.59
0.1, 642.025
0.125, 800.695
0.15, 963.145
0.175, 1126.195
0.2, 1285.365
0.225, 1444.985
0.25, 1603.775
0.275, 1765.03
0.3, 1927.81