simulation_functions-thomas_process.md
In mobsim: Spatial Simulation and Scale-Dependent Analysis of Biodiversity Changes

sim_thomas_community() - default argument values are as expected

Code
  sim_thomas_community(s_pool = 3L, n_sim = 100L, seed = 42L)
Output
  $census
                x           y   species
  1   0.110913769 0.943050781 species_1
  2   0.982808321 0.677287570 species_1
  3   0.856746803 0.738675602 species_1
  4   0.730903256 0.409163714 species_1
  5   0.898514625 0.662671447 species_1
  6   0.794821457 0.715673905 species_1
  7   0.141834638 0.946063141 species_1
  8   0.280994147 0.422478515 species_1
  9   0.500112333 0.256151539 species_1
  10  0.148763344 0.960716593 species_1
  11  0.305077798 0.443650749 species_1
  12  0.921386233 0.688046632 species_1
  13  0.908459093 0.597920694 species_1
  14  0.929854267 0.720228049 species_1
  15  0.910257912 0.694875781 species_1
  16  0.898578180 0.685874316 species_1
  17  0.286894878 0.463807244 species_1
  18  0.921398634 0.736579334 species_1
  19  0.531549143 0.198677562 species_1
  20  0.750174091 0.464089480 species_1
  21  0.907404686 0.675992153 species_1
  22  0.834152237 0.730748726 species_1
  23  0.170293091 0.951939625 species_1
  24  0.648462481 0.955442369 species_1
  25  0.712448578 0.468592521 species_1
  26  0.913002316 0.669462654 species_1
  27  0.968989094 0.684256354 species_1
  28  0.301503110 0.467017128 species_1
  29  0.899187367 0.714427943 species_1
  30  0.646903948 0.936441052 species_1
  31  0.820928533 0.701265213 species_1
  32  0.826104829 0.715457118 species_1
  33  0.741257842 0.477006104 species_1
  34  0.727064836 0.475299794 species_1
  35  0.151017478 0.928815111 species_1
  36  0.713953541 0.433108543 species_1
  37  0.617458644 0.947008175 species_1
  38  0.662640541 0.914608074 species_1
  39  0.133893758 0.927128493 species_1
  40  0.510650831 0.252981821 species_1
  41  0.634219704 0.959494572 species_1
  42  0.939236868 0.695356079 species_1
  43  0.268244027 0.443978457 species_1
  44  0.926822408 0.759102604 species_1
  45  0.937327558 0.653024922 species_1
  46  0.612336804 0.937166295 species_1
  47  0.916310598 0.736255446 species_1
  48  0.924801981 0.664571227 species_1
  49  0.916804555 0.663046857 species_1
  50  0.904948496 0.656993548 species_1
  51  0.720113935 0.486481677 species_1
  52  0.734241923 0.486322791 species_1
  53  0.281268181 0.450993465 species_1
  54  0.897052463 0.648098610 species_1
  55  0.788814706 0.681396970 species_1
  56  0.509476092 0.246765444 species_1
  57  0.592879100 0.368309826 species_2
  58  0.741049959 0.049715743 species_2
  59  0.491715511 0.153457006 species_2
  60  0.606232635 0.320251254 species_2
  61  0.497304294 0.335716513 species_2
  62  0.154492979 0.880127008 species_2
  63  0.472829897 0.365476029 species_2
  64  0.110498855 0.912218226 species_2
  65  0.510590511 0.396029059 species_2
  66  0.769915777 0.070049467 species_2
  67  0.143919743 0.886271171 species_2
  68  0.451945663 0.164739862 species_2
  69  0.134210524 0.869188670 species_2
  70  0.768139535 0.043659458 species_2
  71  0.776117067 0.058884837 species_2
  72  0.608020731 0.301625973 species_2
  73  0.575663009 0.328749640 species_2
  74  0.614698248 0.363810591 species_2
  75  0.003811242 0.931268740 species_2
  76  0.168826475 0.919619832 species_2
  77  0.751255553 0.064503013 species_2
  78  0.462687026 0.379263583 species_2
  79  0.114312613 0.877564337 species_2
  80  0.485923636 0.390496295 species_2
  81  0.484478036 0.394465181 species_2
  82  0.786591351 0.056044986 species_2
  83  0.756199271 0.088667052 species_2
  84  0.138043359 0.383685013 species_3
  85  0.828085319 0.668727012 species_3
  86  0.831103341 0.656485310 species_3
  87  0.322695493 0.002309701 species_3
  88  0.807263700 0.701841268 species_3
  89  0.870979203 0.669694587 species_3
  90  0.880304623 0.690784906 species_3
  91  0.821971491 0.665013712 species_3
  92  0.465128076 0.444997724 species_3
  93  0.818314542 0.647280967 species_3
  94  0.783749076 0.700757131 species_3
  95  0.794417569 0.711638164 species_3
  96  0.301837430 0.041589058 species_3
  97  0.304865176 0.091952956 species_3
  98  0.862811923 0.628776519 species_3
  99  0.453195018 0.416092355 species_3
  100 0.138914335 0.342856485 species_3

  $x_min_max
  [1] 0 1

  $y_min_max
  [1] 0 1

  attr(,"class")
  [1] "community"

sim_thomas_coords() - default argument values are as expected

Code
  sim_thomas_coords(abund_vec = 1:4, seed = 42L)
Output
  $census
               x         y   species
  1  0.914806043 0.9370754 species_1
  2  0.808355637 0.5298764 species_2
  3  0.860678066 0.5172028 species_2
  4  0.095764710 0.5926002 species_3
  5  0.111911586 0.5576663 species_3
  6  0.135123186 0.5699768 species_3
  7  0.001201375 0.8709752 species_4
  8  0.867380696 0.2148270 species_4
  9  0.846330157 0.1886754 species_4
  10 0.016443149 0.9206982 species_4

  $x_min_max
  [1] 0 1

  $y_min_max
  [1] 0 1

  attr(,"class")
  [1] "community"

rThomas_rcpp() - default argument values are as expected

Code
  rThomas_rcpp(n_points = 100L, n_mother_points = 2L, xmother = NA, ymother = NA,
    sigma = 0.2)
Output
              x           y
  1   0.9466524 0.609527738
  2   0.9219323 0.549131305
  3   0.8757477 0.474185939
  4   0.8509816 0.778993750
  5   0.5714043 0.129247733
  6   0.9560058 0.213928075
  7   0.9691198 0.562272618
  8   0.9701162 0.421997298
  9   0.9558408 0.756434050
  10  0.8636285 0.867493739
  11  0.9773625 0.589050264
  12  0.7284516 0.812410349
  13  0.7536635 0.735256699
  14  0.6955612 0.332834810
  15  0.7426475 0.059791799
  16  0.5917583 0.043270789
  17  0.8857904 0.278411143
  18  0.9359088 0.201688358
  19  0.9320569 0.852062172
  20  0.8383393 0.183609483
  21  0.9280392 0.622799476
  22  0.8514236 0.707713305
  23  0.7166751 0.665703825
  24  0.7176527 0.806983714
  25  0.8527856 0.781734083
  26  0.8975846 0.108603731
  27  0.6463813 0.338347184
  28  0.6671556 0.284587128
  29  0.9019702 0.616091149
  30  0.9272314 0.735198762
  31  0.8688504 0.453463348
  32  0.5671628 0.221536168
  33  0.9069202 0.540646200
  34  0.9450972 0.169317741
  35  0.6475587 0.426289298
  36  0.9780671 0.761430030
  37  0.5146202 0.352894974
  38  0.6896146 0.704826047
  39  0.6846349 0.144975256
  40  0.7099392 0.439154335
  41  0.6745616 0.693533968
  42  0.6453982 0.747903496
  43  0.9445585 0.804030019
  44  0.5104063 0.680679139
  45  0.9919396 0.215836960
  46  0.8975704 0.460043816
  47  0.8800024 0.389273081
  48  0.8138862 0.375542685
  49  0.9864894 0.109636786
  50  0.8965584 0.549504054
  51  0.9120125 0.694632630
  52  0.7342248 0.184206654
  53  0.6955748 0.295949625
  54  0.7080313 0.138451384
  55  0.8516260 0.705681596
  56  0.6527710 0.301763519
  57  0.6501359 0.167820254
  58  0.8370989 0.152286760
  59  0.7995560 0.919655837
  60  0.9050854 0.778831629
  61  0.8416083 0.797195328
  62  0.9301769 0.748522286
  63  0.9138819 0.438187968
  64  0.9260159 0.634760013
  65  0.9254979 0.926921519
  66  0.9006136 0.134355404
  67  0.8315415 0.129161973
  68  0.8401563 0.868273388
  69  0.9274545 0.574398472
  70  0.9001144 0.008734224
  71  0.6926516 0.166059354
  72  0.9105206 0.420239149
  73  0.7673123 0.612743654
  74  0.9031987 0.376023626
  75  0.8917240 0.983719869
  76  0.7766940 0.230853318
  77  0.7845687 0.077189546
  78  0.8934405 0.769492035
  79  0.7094022 0.897995867
  80  0.8182665 0.619537486
  81  0.7384596 0.679033320
  82  0.7373277 0.829961070
  83  0.9965024 0.549941554
  84  0.7965875 0.768161582
  85  0.9907499 0.921282692
  86  0.7189106 0.754378003
  87  0.9254201 0.174734810
  88  0.7418242 0.571427393
  89  0.8200731 0.894639131
  90  0.6763112 0.780264733
  91  0.8510299 0.609808322
  92  0.9077231 0.547964407
  93  0.6716230 0.744399106
  94  0.9352641 0.538820929
  95  0.8786627 0.946089867
  96  0.8365616 0.187906718
  97  0.9755999 0.456193414
  98  0.8919547 0.445457540
  99  0.8660830 0.340074856
  100 0.7337550 0.222980551

Code
  rThomas_rcpp(n_points = 100L, n_mother_points = 2L, xmother = 0.5, ymother = 0.5,
    sigma = 0.2)
Output
               x          y
  1   0.46433478 0.47697280
  2   0.35640712 0.27085336
  3   0.63522529 0.67971992
  4   0.28158235 0.18740570
  5   0.70582814 0.68295497
  6   0.09804431 0.61650747
  7   0.46037513 0.29415824
  8   0.69491107 0.37626795
  9   0.72299436 0.41741103
  10  0.43565710 0.16426601
  11  0.90826263 0.15603993
  12  0.52231538 0.56931722
  13  0.81948265 0.43328292
  14  0.33019099 0.70128729
  15  0.68409051 0.25869219
  16  0.68550260 0.51399955
  17  0.13911801 0.27563548
  18  0.60228930 0.36147738
  19  0.53247614 0.48707862
  20  0.59250173 0.16214114
  21  0.45431330 0.43043438
  22  0.10737080 0.51680993
  23  0.77647463 0.65261941
  24  0.50247465 0.22794346
  25  0.18774684 0.56495413
  26  0.57251604 0.44504938
  27  0.56021068 0.79856222
  28  0.25277226 0.66342880
  29  0.61754325 0.51792846
  30  0.66151266 0.65357480
  31  0.09843536 0.60819470
  32  0.54869050 0.18075501
  33  0.36571699 0.46853191
  34  0.36000715 0.41133816
  35  0.70905517 0.87546591
  36  0.80153188 0.50550364
  37  0.51071420 0.41302032
  38  0.36089556 0.42997953
  39  0.20806921 0.70969147
  40  0.05498505 0.88619994
  41  0.49743405 0.53038946
  42  0.65854255 0.48448342
  43  0.53206548 0.41327161
  44  0.58854577 0.25797324
  45  0.45030341 0.58446408
  46  0.77752318 0.50153221
  47  0.81281390 0.17540481
  48  0.52406528 0.45052283
  49  0.12683723 0.54903581
  50  0.68998667 0.28223238
  51  0.66773387 0.36907694
  52  0.51189457 0.49393584
  53  0.70022265 0.64949040
  54  0.69027857 0.66726580
  55  0.62616368 0.41345896
  56  0.15345027 0.39489233
  57  0.59310252 0.61487125
  58  0.53584219 0.54791908
  59  0.36760175 0.34452616
  60  0.58992199 0.62361993
  61  0.66058654 0.38530463
  62  0.58047322 0.56697776
  63  0.27873494 0.42788773
  64  0.66895844 0.49160570
  65  0.74998399 0.72926268
  66  0.19995581 0.37860215
  67  0.41240908 0.54098200
  68  0.38016362 0.75138133
  69  0.64505804 0.69744163
  70  0.55312495 0.71534525
  71  0.58186256 0.36165025
  72  0.64362724 0.59789140
  73  0.71688793 0.40438873
  74  0.31670879 0.24963531
  75  0.16174293 0.45840693
  76  0.50055444 0.23435806
  77  0.46763848 0.18414169
  78  0.40499326 0.38498858
  79  0.35964877 0.20357679
  80  0.32446712 0.25742057
  81  0.37758732 0.54167620
  82  0.39843042 0.33981290
  83  0.55893847 0.57854825
  84  0.34001303 0.59002398
  85  0.37291370 0.25803186
  86  0.47900584 0.17827495
  87  0.44823177 0.84591164
  88  0.93177327 0.31794551
  89  0.40254843 0.77458156
  90  0.59107101 0.84498682
  91  0.56572807 0.79330213
  92  0.55500854 0.73600334
  93  0.78443865 0.63277532
  94  0.81985220 0.09811363
  95  0.52831457 0.49921056
  96  0.56486313 0.35599497
  97  0.64827201 0.48007865
  98  0.48010188 0.41216634
  99  0.31669307 0.84193771
  100 0.53174443 0.83714527

Code
  rThomas_rcpp(n_points = 100L, n_mother_points = 2L, xmother = c(0.2, 0.8),
  ymother = c(0.2, 0.8), sigma = 0.2)
Output
                x           y
  1   0.930225194 0.693433475
  2   0.280881462 0.321842209
  3   0.303956577 0.085806675
  4   0.258104896 0.284461271
  5   0.544823799 0.666001742
  6   0.920920725 0.796171846
  7   0.092662661 0.325616847
  8   0.713741879 0.700094210
  9   0.598011118 0.009503932
  10  0.396679028 0.285698543
  11  0.736523627 0.835800793
  12  0.929037989 0.752638608
  13  0.754331477 0.935071112
  14  0.331030278 0.770183787
  15  0.082380428 0.067940834
  16  0.889935464 0.590223340
  17  0.251906291 0.232312014
  18  0.711755469 0.624261184
  19  0.937761550 0.945016604
  20  0.738212478 0.709419423
  21  0.954296808 0.886177547
  22  0.939415440 0.697986768
  23  0.009082876 0.186285387
  24  0.962272500 0.627077284
  25  0.289235030 0.752813210
  26  0.015041067 0.472709206
  27  0.348479032 0.025141465
  28  0.315311754 0.144900108
  29  0.797379870 0.933593615
  30  0.574363901 0.869759914
  31  0.944780546 0.666233452
  32  0.554463608 0.647198818
  33  0.735409938 0.976035980
  34  0.843909602 0.874058884
  35  0.793152383 0.977557159
  36  0.181826999 0.473006273
  37  0.228959675 0.071878646
  38  0.498250245 0.826700767
  39  0.005321169 0.036034945
  40  0.636601670 0.700210819
  41  0.603787899 0.799600882
  42  0.161655553 0.491878433
  43  0.146626646 0.246637875
  44  0.234094694 0.440133632
  45  0.388384843 0.150277193
  46  0.866614930 0.546984328
  47  0.208244076 0.816177648
  48  0.825012739 0.493490537
  49  0.843242309 0.837532763
  50  0.503121516 0.479806850
  51  0.608499552 0.775118884
  52  0.895576628 0.889512601
  53  0.250002809 0.349530499
  54  0.988825414 0.610747840
  55  0.345252807 0.451160894
  56  0.869539346 0.982752517
  57  0.955085077 0.913363519
  58  0.418318116 0.900589252
  59  0.692121835 0.608450870
  60  0.924896926 0.681410654
  61  0.911139824 0.363370130
  62  0.112903118 0.013067452
  63  0.461087668 0.559064774
  64  0.488421878 0.228207343
  65  0.905745607 0.967015785
  66  0.987960205 0.832674939
  67  0.438532853 0.518434127
  68  0.536707268 0.541195965
  69  0.161094261 0.261121137
  70  0.587178057 0.477284506
  71  0.607696967 0.489814280
  72  0.306564115 0.343933151
  73  0.661943322 0.939751850
  74  0.352140649 0.121754630
  75  0.503456426 0.910864128
  76  0.449048530 0.236275489
  77  0.800927277 0.870142702
  78  0.952764987 0.679322525
  79  0.071904246 0.325515271
  80  0.916594432 0.602686513
  81  0.674952248 0.651522141
  82  0.080503417 0.233050285
  83  0.368983263 0.286188790
  84  0.740308802 0.743277725
  85  0.897049544 0.759595187
  86  0.515533108 0.554497473
  87  0.488374328 0.191643579
  88  0.952734841 0.602811614
  89  0.277687815 0.031021348
  90  0.882137422 0.887029341
  91  0.857758672 0.818162264
  92  0.198625403 0.547154422
  93  0.579637026 0.619981971
  94  0.083579611 0.049223299
  95  0.332008302 0.149879913
  96  0.343784019 0.048772237
  97  0.245312715 0.284328716
  98  0.911662474 0.899192256
  99  0.942229883 0.613260803
  100 0.133158269 0.187003851