Implementation of the random number generators for multivariate extreme-value distributions and max-stable processes based on the two algorithms described in Dombry, Engelke and Oesting (2016).
Usage
rmev(
n,
d,
param,
asy,
sigma,
model = c("log", "alog", "neglog", "aneglog", "bilog", "negbilog", "hr", "br", "xstud",
"smith", "schlather", "ct", "sdir", "dirmix", "pairbeta", "pairexp", "wdirbs",
"wexpbs"),
alg = c("ef", "sm"),
weights = NULL,
vario = NULL,
coord = NULL,
grid = FALSE,
dist = NULL,
...
)Arguments
- n
number of observations
- d
dimension of sample
- param
parameter vector for the logistic, bilogistic, negative bilogistic and extremal Dirichlet (Coles and Tawn) model. Parameter matrix for the Dirichlet mixture. Degree of freedoms for extremal student model. See Details.
- asy
list of asymmetry parameters, as in function
rmvevdfrom pcakageevd, of \(2^d-1\) vectors of size corresponding to the power set ofd, with sum to one constraints.- sigma
covariance matrix for Brown-Resnick and extremal Student-t distributions. Symmetric matrix of squared coefficients \(\lambda^2\) for the Husler-Reiss model, with zero diagonal elements.
- model
for multivariate extreme value distributions, users can choose between 1-parameter logistic and negative logistic, asymmetric logistic and negative logistic, bilogistic, Husler-Reiss, extremal Dirichlet model (Coles and Tawn) or the Dirichlet mixture. Spatial models include the Brown-Resnick, Smith, Schlather and extremal Student max-stable processes.
- alg
algorithm, either simulation via extremal function (
'ef') or via the spectral measure ('sm'). Default toef.- weights
vector of length
mfor themmixture components. Must sum to one- vario
semivariogram function whose first argument must be distance. Used only if provided in conjunction with
coordand ifsigmais missing- coord
dbykmatrix of coordinates, used as input in the variogramvarioor as parameter for the Smith model. IfgridisTRUE, unique entries should be supplied.- grid
Logical.
TRUEif the coordinates are two-dimensional grid points (spatial models).- dist
symmetric matrix of pairwise distances. Default to
NULL.- ...
additional arguments for the
variofunction
Details
The vector param differs depending on the model
log: one dimensional parameter greater than 1alog: \(2^d-d-1\) dimensional parameter fordep. Values are recycled if needed.neglog: one dimensional positive parameteraneglog: \(2^d-d-1\) dimensional parameter fordep. Values are recycled if needed.bilog:d-dimensional vector of parameters in \([0,1]\)negbilog:d-dimensional vector of negative parametersct, dir, negdir, sdir:d-dimensional vector of positive (a)symmetry parameters. Fordirandnegdir, a \(d+1\) vector consisting of thedDirichlet parameters and the last entry is an index of regular variation in \((-\min(\alpha_1, \ldots, \alpha_d), 1]\) treated as shape parameterxstud: one dimensional parameter corresponding to degrees of freedomalphadirmix:dbym-dimensional matrix of positive (a)symmetry parameterspairbeta, pairexp:d(d-1)/2+1vector of parameters, containing the concentration parameter and the coefficients of the pairwise beta, in lexicographical order e.g., \(\beta_{12}, \beta_{13}, \ldots\)wdirbs, wexpbs:2dvector ofdconcentration parameters followed by thedDirichlet parameters
Stephenson points out that the multivariate asymmetric negative logistic model given in e.g. Coles and Tawn (1991) is not a valid distribution function in dimension \(d>3\) unless additional constraints are imposed on the parameter values.
The implementation in mev uses the same construction as the asymmetric logistic distribution (see the vignette). As such it does not match the bivariate implementation of rbvevd.
The dependence parameter of the evd package for the Husler-Reiss distribution can be recovered taking
for the Brown--Resnick model \(2/r=\sqrt(2\gamma(h))\) where \(h\) is the lag vector between sites and \(r=1/\lambda\) for the Husler--Reiss.
Warning
As of version 1.8 (August 16, 2016), there is a distinction between models hr and br. The latter is meant to be used in conjunction with variograms. The parametrization differs between the two models.
The family of scaled Dirichlet is now parametrized by a parameter in \(-\min(\alpha)\) appended to the the d vector param containing the parameter alpha
of the Dirichlet model. Arguments model='dir' and model='negdir' are still supported internally, but not listed in the options.
References
Dombry, Engelke and Oesting (2016). Exact simulation of max-stable processes, Biometrika, 103(2), 303--317.
Examples
set.seed(1)
rmev(n=100, d=3, param=2.5, model='log', alg='ef')
#> [,1] [,2] [,3]
#> [1,] 1.3241844 2.7993629 1.0646407
#> [2,] 0.8080132 1.6992939 0.7179083
#> [3,] 26.8322925 23.0714108 37.5944866
#> [4,] 4.2548221 3.3191385 3.0360654
#> [5,] 1.7926629 2.8146683 1.9201000
#> [6,] 0.3436559 0.2275397 0.1736779
#> [7,] 4.1524779 1.6554806 1.1395523
#> [8,] 0.6924634 1.6015045 0.4729466
#> [9,] 0.7722329 1.0479223 0.6663154
#> [10,] 0.9734914 2.2388627 1.2705100
#> [11,] 9.3166649 5.0719498 4.3993434
#> [12,] 1.6068357 2.0191002 2.5851826
#> [13,] 1.7208045 6.4176355 2.0864478
#> [14,] 0.7217394 3.6688121 1.2280653
#> [15,] 0.2550256 0.3148263 0.3432364
#> [16,] 1.4594709 2.6277023 1.2050499
#> [17,] 0.7952011 0.9130174 0.7408276
#> [18,] 5.6864751 10.8918433 9.1068160
#> [19,] 0.6507906 0.5649713 2.5489719
#> [20,] 0.2943543 0.3187133 0.3552160
#> [21,] 3.9921646 5.2547137 4.3501015
#> [22,] 45.4919358 44.1298669 53.0023479
#> [23,] 9.2109274 14.1265549 14.8903439
#> [24,] 6.5100218 2.3194623 5.4323191
#> [25,] 0.7788335 1.5107596 0.8286699
#> [26,] 0.4689831 0.3076351 0.2766227
#> [27,] 2.2297072 4.5375915 2.9498290
#> [28,] 3.5158139 8.9986099 2.1972200
#> [29,] 1.1075761 0.7893003 0.4072057
#> [30,] 4.4474789 3.5740845 1.3516939
#> [31,] 0.1880307 0.7157074 0.3655790
#> [32,] 4.1948912 3.6034486 2.4692792
#> [33,] 2.4452410 2.0918579 1.7233233
#> [34,] 2.1078167 1.9527438 2.3825863
#> [35,] 20.6483841 7.3692035 12.9039323
#> [36,] 5.7910208 0.9927368 1.6359326
#> [37,] 0.4405659 0.5622922 0.2563285
#> [38,] 0.2737696 0.4675466 0.8151963
#> [39,] 1.9077606 1.3977859 3.2930667
#> [40,] 518.7965587 654.0947432 559.9352772
#> [41,] 0.5044269 0.4816114 0.5386546
#> [42,] 0.6525537 0.6020635 0.6410463
#> [43,] 0.8297775 0.3931504 0.5103706
#> [44,] 5.5615656 9.4178517 2.5822721
#> [45,] 0.8782299 0.6859625 0.8512454
#> [46,] 0.7255494 0.5939896 0.8462600
#> [47,] 32.9434266 24.5289388 76.1289342
#> [48,] 0.8242095 0.7357299 0.9947758
#> [49,] 0.7569966 0.3463297 1.0311805
#> [50,] 4.4421658 6.2614536 3.0787288
#> [51,] 2.0389109 0.8444871 1.2737657
#> [52,] 6.7758576 4.8729472 2.8735788
#> [53,] 5.6882860 5.4424085 4.8457198
#> [54,] 0.6472684 0.8741934 0.8616940
#> [55,] 1.0202678 1.8892189 2.9672814
#> [56,] 1.5625309 1.5043384 0.8932051
#> [57,] 5.3784519 4.7888378 2.3316439
#> [58,] 2.7429866 6.5206491 5.7368401
#> [59,] 0.4239659 0.3392086 0.3693929
#> [60,] 1.4285849 1.1103833 0.8652395
#> [61,] 0.5262386 0.4187189 0.4797577
#> [62,] 1.4667154 0.7042429 2.1704179
#> [63,] 0.3941662 0.5197489 0.3099264
#> [64,] 0.5009843 1.8451471 0.6201354
#> [65,] 0.2289987 1.7555116 0.4913419
#> [66,] 0.7424231 1.3076994 1.1753006
#> [67,] 0.3944793 1.0649244 0.4124939
#> [68,] 0.3826278 0.5261374 0.5137356
#> [69,] 2.5422436 1.9500485 1.8923649
#> [70,] 5.2454130 10.4540583 7.5350216
#> [71,] 6.8215580 1.6971728 1.5727520
#> [72,] 2.4431299 0.9174740 1.1897401
#> [73,] 0.2947914 0.2014534 0.7053788
#> [74,] 9.8183223 7.4168505 28.9644807
#> [75,] 2.0977387 1.4587730 1.2596296
#> [76,] 2.1063441 0.6056433 0.7216470
#> [77,] 0.3791312 0.4946861 1.2087285
#> [78,] 1.4132108 1.6579087 1.3394142
#> [79,] 1.7290945 4.4431699 1.3579806
#> [80,] 6.3547801 5.0884003 15.6788521
#> [81,] 1.8221446 1.1368440 3.6184991
#> [82,] 3.4655900 3.2022965 4.4324373
#> [83,] 10.2291760 9.8751531 13.2110210
#> [84,] 2.5906485 6.9595485 5.9532892
#> [85,] 0.5516058 1.2833323 0.5715908
#> [86,] 0.7759360 0.6456412 0.7055688
#> [87,] 0.5620228 0.5964398 0.7658183
#> [88,] 1.5912064 3.0733047 5.5349735
#> [89,] 3.1445205 0.6125466 0.6581886
#> [90,] 1.4322970 2.4934474 0.9413238
#> [91,] 0.4132584 1.7433116 1.4005017
#> [92,] 4.2212600 6.5874352 3.4886966
#> [93,] 2.6039127 2.2772943 1.1978808
#> [94,] 0.6981758 0.7403026 5.6111802
#> [95,] 1.6999919 1.0957491 1.3305272
#> [96,] 4.2275592 4.1182303 2.7306551
#> [97,] 0.9989890 3.1715032 3.0418132
#> [98,] 0.2552138 0.2086764 0.3668845
#> [99,] 6.0318994 3.3560933 2.2229644
#> [100,] 0.3473329 0.4059501 0.5636032
rmev(n=100, d=4, param=c(0.2,0.1,0.9,0.5), model='bilog', alg='sm')
#> [,1] [,2] [,3] [,4]
#> [1,] 1.1884635 1.7856587 2.8667964 2.2733176
#> [2,] 0.9032010 0.6720215 1.2541370 0.8445604
#> [3,] 1.7075216 1.9608242 0.3366486 0.9901735
#> [4,] 0.4156491 0.3476536 0.1668837 0.2283162
#> [5,] 1.1298928 0.5550409 4.6648697 0.4910048
#> [6,] 0.3224053 0.4380449 9.6972984 0.4065408
#> [7,] 1.7367523 2.3840136 2.1111006 1.6840591
#> [8,] 1.9166787 2.4495006 1.6636824 3.1150988
#> [9,] 0.5067202 0.4726854 0.2515280 0.4411369
#> [10,] 11.7437015 13.4095504 1.6479289 10.8853319
#> [11,] 0.6995022 0.6393852 0.3757318 0.8326771
#> [12,] 0.3800904 0.3449422 1.1407585 0.2197184
#> [13,] 3.7710060 3.7483615 0.7849528 3.7461276
#> [14,] 1.7425539 2.3762805 0.5773511 1.4901389
#> [15,] 0.3480587 0.2756752 0.7426411 0.3763338
#> [16,] 4.4209590 4.4931624 1.4037907 4.7047912
#> [17,] 0.6159474 0.3805702 1.9014536 2.6146804
#> [18,] 7.0803539 8.0412040 1.1650497 6.5926116
#> [19,] 4.5906471 4.9678162 1.4667451 1.8034935
#> [20,] 2.2942669 1.5583731 1.4383188 0.9897646
#> [21,] 56.9384856 43.5959255 21.7835460 23.4911680
#> [22,] 3.1496360 3.5144286 2.2835179 1.5948546
#> [23,] 1.4426739 1.7089270 4.2473665 3.3691313
#> [24,] 2.6630914 2.2005614 0.7211449 3.3945207
#> [25,] 0.8094212 0.8864928 0.7710882 0.4694609
#> [26,] 1.5989515 1.7162134 0.3248020 1.2279848
#> [27,] 0.4937936 0.7036582 0.6242557 0.4632328
#> [28,] 0.6923108 1.1141193 0.3119199 1.6666024
#> [29,] 0.9973205 0.9246299 0.3460527 0.5627246
#> [30,] 0.4002020 0.3645767 169.3113196 0.3172569
#> [31,] 0.3530947 0.4394674 1.5978777 0.6109954
#> [32,] 4.7824812 3.5195027 0.7250942 2.5419976
#> [33,] 1.5728188 1.6341189 5.1859387 1.7909067
#> [34,] 0.3429683 0.4285018 3.7552176 0.3764642
#> [35,] 21.4590685 13.4515813 13.0075224 27.2801653
#> [36,] 1.5529179 1.8094334 9.2302150 1.0811398
#> [37,] 0.3387633 0.2861620 5.5813367 0.5004467
#> [38,] 1.3000031 1.6151976 1.2505467 1.8911989
#> [39,] 5.2531148 2.4354263 133.3924922 1.0389956
#> [40,] 0.7605952 0.7120925 0.8856417 0.4306784
#> [41,] 0.5705526 0.6126565 1.0061543 0.4820991
#> [42,] 2.2220782 2.8037560 0.9964778 3.2412322
#> [43,] 0.8132851 0.6505220 0.4941494 1.3546971
#> [44,] 1.4475011 1.4125120 0.6792525 0.4936465
#> [45,] 0.7393947 0.8615995 1.0265909 0.5293244
#> [46,] 0.9569868 0.9507842 1.0416360 0.9049183
#> [47,] 0.4854348 0.3900925 0.4753514 0.5542988
#> [48,] 1.3934476 1.7967266 1.3667497 1.5258232
#> [49,] 0.6879177 0.6372396 3.3790040 0.5296521
#> [50,] 1.6948870 1.4645749 0.8877760 1.4178395
#> [51,] 14.7184457 9.2388045 8.8106998 10.2894886
#> [52,] 0.5881312 0.6773375 1.7791944 0.9337532
#> [53,] 2.1488720 2.2719036 1.5975857 1.0084198
#> [54,] 9.5549346 8.8806253 1.3528289 6.1275327
#> [55,] 0.5989547 0.6636190 0.6192832 3.6334874
#> [56,] 1.7850294 0.9748040 0.6469666 0.7049369
#> [57,] 33.6368542 42.4667696 26.0254380 22.1432179
#> [58,] 0.9651103 0.9931536 1.0563475 0.6784568
#> [59,] 0.6457674 0.4758477 0.9404102 0.3763188
#> [60,] 2.5029137 2.2874233 1.0633913 2.2840566
#> [61,] 0.8395874 0.7187367 2.2695518 0.4581007
#> [62,] 0.7556948 0.7349402 0.6403012 0.7347460
#> [63,] 16.1254142 14.6595358 6.2462223 11.7277311
#> [64,] 8.5763334 11.7596210 4.1187242 13.1519686
#> [65,] 1.4075540 1.5583706 5.4806099 3.8322557
#> [66,] 0.7773871 0.8902242 0.7012725 0.8922112
#> [67,] 3.1252387 4.1993378 1.7147613 1.7754124
#> [68,] 28.4086749 18.3648855 2.0059032 30.5148516
#> [69,] 0.6330372 0.7132324 0.4107640 3.8127257
#> [70,] 3.2438869 3.4645575 2.5188149 6.5789334
#> [71,] 0.3310286 0.3428823 5.2519884 0.6023043
#> [72,] 8.5553953 6.6879366 2.8070044 3.6172843
#> [73,] 0.5212343 1.1856712 0.8254436 0.4810649
#> [74,] 1.6063783 1.7594667 0.4895264 1.6241634
#> [75,] 0.4783367 0.5442460 0.8207148 0.2072106
#> [76,] 0.7711092 0.7460709 0.8061218 3.7351572
#> [77,] 14.5657010 16.0249078 4.8870310 12.0103480
#> [78,] 2.5828553 2.6662212 0.6794608 3.0703792
#> [79,] 1.8446194 2.3708271 0.7091086 1.2919117
#> [80,] 3.8542471 3.5560984 1.3017582 8.5111474
#> [81,] 0.4397355 0.4484222 2.8764417 0.2363978
#> [82,] 0.4210346 0.5019033 0.5087373 2.5070201
#> [83,] 3.1001613 3.3345045 1.2804778 1.9047906
#> [84,] 0.3857180 0.4356819 12.4961168 1.2548797
#> [85,] 0.7289760 0.9559479 2.2172534 0.6009609
#> [86,] 0.2960657 0.3088411 1.2612274 0.6581833
#> [87,] 2.7655550 2.4564235 4.0441158 1.3555999
#> [88,] 3.6246254 3.3613546 0.7627075 4.4099907
#> [89,] 2.8719349 2.3708033 6.8367612 7.4823751
#> [90,] 0.6752532 0.5590593 1.3493152 0.5473113
#> [91,] 3.3226255 3.2965453 2.4067569 2.0675351
#> [92,] 0.6707315 0.8748905 1.3554857 1.3167451
#> [93,] 0.8514512 0.7275407 3.7170197 0.5119197
#> [94,] 6.9704161 9.0544583 1.5059288 10.7108686
#> [95,] 7.8529006 9.4703257 6.2684291 11.9846924
#> [96,] 0.9038928 0.7684667 0.6516356 1.3644579
#> [97,] 6.2625974 5.5626443 2.3676784 5.0974838
#> [98,] 1.4696896 1.9987126 1.4589666 1.5832782
#> [99,] 0.8136669 0.7285609 1.0887153 0.7871845
#> [100,] 0.3300757 0.3073172 1.8663563 0.2142376
## Spatial example using power variogram
#NEW: Semi-variogram must take distance as argument
semivario <- function(x, scale, alpha){ scale*x^alpha }
#grid specification
grid.coord <- as.matrix(expand.grid(runif(4), runif(4)))
rmev(n=100, vario=semivario, coord=grid.coord, model='br', scale = 0.5, alpha = 1)
#> [,1] [,2] [,3] [,4] [,5] [,6]
#> [1,] 0.8426083 0.7302149 0.8595143 1.2658793 0.8395195 0.7320301
#> [2,] 1.5199974 1.5449539 1.2240541 1.6035839 1.5196261 1.5532248
#> [3,] 0.6088198 0.5904609 0.5934316 0.5064741 0.6223486 0.5905136
#> [4,] 2.1171572 2.7276682 2.8080420 1.3993248 2.1080347 2.8131175
#> [5,] 2.0326237 2.0538142 1.8373249 2.3039332 2.0137120 2.0455075
#> [6,] 2.6826745 2.4342983 1.4418201 2.0512716 2.6473059 2.4379174
#> [7,] 19.5229087 14.1133826 14.5960848 21.5587371 19.1862665 14.0922850
#> [8,] 0.5620293 0.8013278 1.0972927 0.6738161 0.5705288 0.7904279
#> [9,] 1.6908310 1.4060703 1.3095035 2.1935243 1.6901523 1.4180811
#> [10,] 2.7573322 2.7079406 2.2670682 3.3608515 2.7396211 2.7001485
#> [11,] 0.9161147 0.7492908 0.9298025 0.7648425 0.9451168 0.7568701
#> [12,] 1.8346379 2.0902827 2.3185569 1.6776362 1.8478125 2.0977678
#> [13,] 2.0805400 2.0010957 1.8926709 2.8205763 2.1073916 1.9961118
#> [14,] 0.3870196 0.4317807 0.4121914 0.4409632 0.3874604 0.4317450
#> [15,] 0.5982203 0.7239396 0.9266797 0.6628174 0.6019651 0.7182503
#> [16,] 1.8498994 1.8085588 1.6540520 3.5483621 1.8667318 1.8568906
#> [17,] 0.4940947 0.5078744 0.5222667 0.4727587 0.4969251 0.5132593
#> [18,] 0.2438105 0.2691679 0.1839277 0.3128320 0.2387531 0.2652175
#> [19,] 1.1340436 1.5392345 0.9191043 1.1896614 1.1146288 1.5482971
#> [20,] 0.4872884 0.5445631 0.6121230 0.4413535 0.4880010 0.5423699
#> [21,] 0.9508563 0.8336151 1.0771294 1.0778311 0.9264179 0.8364528
#> [22,] 130.8864583 136.7697063 113.7737251 97.8790596 133.0114150 137.0302104
#> [23,] 1.8242597 1.9913570 1.4777887 2.7793295 1.8298311 1.9851785
#> [24,] 1.0448629 1.1824350 0.9807948 0.7729364 1.0438824 1.1782034
#> [25,] 2.3113371 2.4990377 2.2407845 2.5106443 2.2995544 2.4809869
#> [26,] 52.9133772 52.2316092 35.4347863 54.8314604 53.0388026 52.1259826
#> [27,] 9.7925055 7.8564244 7.0333323 8.9547918 9.7243409 8.0425023
#> [28,] 0.9950331 0.9847861 1.4047819 1.0536189 1.0212370 0.9830035
#> [29,] 1.6293715 1.5954972 1.8194734 1.6194653 1.5921860 1.5699304
#> [30,] 3.6325551 4.2271676 3.2845134 6.2561412 3.6918101 4.2325798
#> [31,] 0.9669727 0.8654438 0.5394738 0.8127238 0.9461243 0.8859033
#> [32,] 1.0770786 1.2132190 1.0808176 1.0150982 1.0571392 1.1984023
#> [33,] 24.2306240 26.1922981 22.3818810 21.9422163 23.8967535 26.3232389
#> [34,] 2.7137924 2.7495218 2.1766820 2.0877566 2.7351239 2.6990720
#> [35,] 6.7412122 6.6351219 5.8668637 6.0652172 6.7594036 6.6464301
#> [36,] 1.1436936 1.0146938 1.2783765 0.8522676 1.1436905 1.0184643
#> [37,] 0.5547986 0.4976018 0.4252762 0.4646366 0.5634325 0.4909869
#> [38,] 0.6392770 0.6621808 0.6362456 0.7011622 0.6524610 0.6658638
#> [39,] 0.8020017 0.8058205 0.9201135 0.7305441 0.8027111 0.8118322
#> [40,] 0.9469137 0.8870885 1.2062306 1.1449414 0.9671403 0.8719274
#> [41,] 2.5244447 2.3877954 1.9135019 1.9839181 2.5161425 2.3750274
#> [42,] 3.3951860 3.1827546 2.8483555 2.5694390 3.3800744 3.2198799
#> [43,] 3.3207250 3.6517007 2.4388614 5.5951737 3.3794053 3.6148381
#> [44,] 3.4354300 3.3695747 2.7007622 3.0230622 3.4444952 3.4264089
#> [45,] 1.9083845 1.9504984 3.6798106 2.5657510 1.8903617 1.9274457
#> [46,] 0.3619218 0.4224601 0.5614226 0.3615642 0.3694996 0.4346221
#> [47,] 0.9279558 0.9060527 0.6575254 0.5939395 0.9116666 0.8964040
#> [48,] 9.3723008 11.3793127 10.5238224 10.9632575 9.5025198 11.3044104
#> [49,] 60.4465918 61.0368380 48.2552933 49.1500247 60.4124785 60.3287150
#> [50,] 0.2027938 0.2123719 0.1771356 0.2006902 0.2021040 0.2155369
#> [51,] 0.6673273 0.6259050 0.6395109 0.5783991 0.6500526 0.6254674
#> [52,] 7.6047825 7.1460732 5.6763485 8.3686956 7.6901919 7.2161146
#> [53,] 0.7110485 0.7592389 0.8810381 0.7153269 0.7277134 0.7595429
#> [54,] 2.9053299 3.1881021 2.3806705 2.8193385 2.9620920 3.1545509
#> [55,] 0.9894241 0.8644508 0.8582722 0.8439871 0.9678272 0.8712305
#> [56,] 0.8936386 0.6743083 0.7225482 0.9610249 0.8851675 0.6699073
#> [57,] 0.2334174 0.2589149 0.2802148 0.2842939 0.2323557 0.2588801
#> [58,] 0.7053754 0.8544064 1.0782276 0.7455042 0.7183791 0.8674258
#> [59,] 0.9536414 1.0399452 0.9691168 0.5984984 0.9328715 1.0207021
#> [60,] 0.6578081 0.6765497 0.6006873 0.5209725 0.6603897 0.6819939
#> [61,] 2.0268900 2.1296205 1.8123543 2.0767627 2.0164571 2.1067793
#> [62,] 1.3631737 1.1090441 1.1214431 1.4623760 1.3968597 1.1027889
#> [63,] 5.4865635 5.8465401 3.7809679 6.1960004 5.4968646 5.6695705
#> [64,] 3.4741536 3.5898800 5.8624341 3.3267169 3.5015546 3.5699625
#> [65,] 4.1610579 4.7938905 4.9197444 3.0310767 4.2066585 4.8056127
#> [66,] 2.2923848 2.5126159 2.4763333 2.7003003 2.3534844 2.5400992
#> [67,] 1.7056776 1.5722805 1.3510829 2.1718512 1.6984085 1.5843918
#> [68,] 9.7478480 9.5475024 8.6160351 11.8807480 9.5809984 9.6542213
#> [69,] 3.7622693 4.1142359 4.5073091 4.0568295 3.6870502 4.0746306
#> [70,] 0.6117022 0.6163562 0.7345999 0.6943126 0.6136670 0.6150636
#> [71,] 11.0083810 12.4604510 12.2096756 16.4447586 11.2629770 12.3793965
#> [72,] 0.7113457 0.7078233 0.8331269 1.1494892 0.7087351 0.7065447
#> [73,] 4.4283967 3.7867946 4.9288401 3.9094045 4.4164570 3.8611842
#> [74,] 10.8412252 11.0949766 14.2416297 9.9874671 10.8642363 11.1399998
#> [75,] 1.1867807 1.1742420 0.9294485 1.2235953 1.1623313 1.1520978
#> [76,] 0.2896588 0.2639968 0.3010517 0.3320183 0.2814946 0.2626035
#> [77,] 0.3342088 0.3671696 0.3376078 0.3827003 0.3378268 0.3577576
#> [78,] 1.6845455 1.8398543 1.9321961 2.2001262 1.6637603 1.8740125
#> [79,] 0.5213473 0.6107365 0.7973775 0.4112081 0.5033820 0.6113133
#> [80,] 125.2067763 141.2579273 122.7770717 98.4393427 126.3668475 141.3580343
#> [81,] 0.7712729 0.8388798 0.6918421 0.7798409 0.7879326 0.8529047
#> [82,] 1.7338738 1.9761059 1.5595382 1.1576748 1.7527301 1.9881877
#> [83,] 40.9806601 34.0269067 31.6840801 40.8808869 40.2866002 34.5719366
#> [84,] 0.6001580 0.5075262 0.4315789 0.4538916 0.5862538 0.5017500
#> [85,] 3.5461905 2.5099762 3.0959663 2.4733719 3.4902819 2.5136255
#> [86,] 8.7639700 9.5896901 8.5261619 11.6671219 8.8075829 9.4878695
#> [87,] 1.6638204 1.6408884 1.1056086 1.6086369 1.6452492 1.6156497
#> [88,] 3.5119978 2.9417079 3.7140267 4.4969722 3.5091925 2.9588413
#> [89,] 0.2499934 0.2547207 0.3019126 0.2506668 0.2528157 0.2560502
#> [90,] 0.6317557 0.6955508 0.4605090 0.5108540 0.6184765 0.6868962
#> [91,] 0.2979791 0.2926200 0.4349974 0.2642916 0.2987700 0.2971711
#> [92,] 1.2292259 1.2110646 1.5899029 1.3237745 1.1936552 1.2225458
#> [93,] 0.7774696 0.6532188 0.8487412 0.7302286 0.7682122 0.6749382
#> [94,] 0.4617700 0.3747055 0.4055236 0.3530516 0.4628303 0.3759452
#> [95,] 18.9043440 20.4338412 21.3297595 13.6127109 19.1391302 20.4445650
#> [96,] 0.6474815 0.5879732 0.7262747 0.6810981 0.6508330 0.5804630
#> [97,] 0.8680973 0.7906039 0.6141117 0.7787126 0.8673459 0.7800437
#> [98,] 13.1695169 14.2515859 13.4842031 9.9088226 13.1520611 14.4397445
#> [99,] 0.6486717 0.6436363 0.5179893 0.5214322 0.6484908 0.6280337
#> [100,] 2.4553895 2.4330333 3.2985355 2.6598540 2.5063540 2.3893450
#> [,7] [,8] [,9] [,10] [,11] [,12]
#> [1,] 0.8327566 1.2560504 0.6557407 0.6675480 0.5583301 1.0678820
#> [2,] 1.2202750 1.5653417 3.2493438 3.2160608 2.9612941 4.0489951
#> [3,] 0.5869662 0.5054560 0.7830661 0.9665690 0.8018338 0.4846389
#> [4,] 2.7871801 1.4252057 3.1986589 3.7743797 3.1441579 2.3789840
#> [5,] 1.8481963 2.3640931 1.9369016 2.0535495 2.3211250 1.8389804
#> [6,] 1.4177496 2.0616654 1.0567668 1.2385469 0.9710156 0.8670988
#> [7,] 14.4080706 21.2407469 10.9858319 10.4823932 9.7766079 12.8002806
#> [8,] 1.1045148 0.6807221 0.7663265 0.8418036 1.0185099 0.6719683
#> [9,] 1.3088083 2.2427207 2.4618356 2.4590485 1.5693856 2.9494863
#> [10,] 2.2741667 3.3021584 2.3191006 2.0604314 2.0237426 2.0529067
#> [11,] 0.9275677 0.7599246 1.0248755 1.1109968 1.4135643 1.0862740
#> [12,] 2.2945343 1.6668284 1.8258300 1.6516551 1.7885314 1.6829577
#> [13,] 1.9065334 2.7646328 1.1968902 1.2932446 0.9318013 2.1083621
#> [14,] 0.4099023 0.4493639 0.3488305 0.3717157 0.2898711 0.3803369
#> [15,] 0.9201046 0.6541192 0.7514847 0.6997322 0.8000001 0.5863166
#> [16,] 1.6173936 3.5604545 2.0611896 2.1654188 2.2303093 3.1090555
#> [17,] 0.5238285 0.4532858 0.3625784 0.2999189 0.3502250 0.4792490
#> [18,] 0.1859557 0.3106789 0.2353987 0.2369068 0.2683563 0.2497575
#> [19,] 0.9208690 1.1608073 0.9256311 1.0294211 1.1044359 0.8451894
#> [20,] 0.6071159 0.4462044 0.2316054 0.3093494 0.4013808 0.2341877
#> [21,] 1.0523298 1.1080439 0.5875544 0.6090133 0.7897838 0.6230546
#> [22,] 113.7827493 97.8440825 149.8543846 125.7274797 167.2679099 168.9301848
#> [23,] 1.4843966 2.8264934 1.9800827 2.3371788 2.0275799 3.6113154
#> [24,] 0.9807941 0.7669417 0.9983796 1.1464428 1.5319643 1.0763389
#> [25,] 2.2412199 2.4975123 1.6518737 1.2660256 1.3354434 1.1341826
#> [26,] 35.4789318 55.3310328 21.4325301 19.0006845 19.9707794 20.9677844
#> [27,] 7.0976424 8.8856778 21.2594429 17.8272348 11.8366064 15.1908863
#> [28,] 1.4217661 1.0778858 1.6517310 1.7034084 1.4540444 1.2603681
#> [29,] 1.7969492 1.5945701 1.2677645 1.2023824 1.1064341 1.9380365
#> [30,] 3.2503767 6.2645805 8.8311522 12.6018059 7.5720479 9.2828105
#> [31,] 0.5430418 0.7985802 0.9482873 0.9750697 1.1642797 1.0030182
#> [32,] 1.0808478 1.0150503 0.7987009 0.8065194 0.9795856 0.9777417
#> [33,] 22.6600982 22.1500481 30.6205255 30.3229354 36.0416505 38.9860394
#> [34,] 2.1741224 2.0790004 2.4256817 2.6208281 2.0556643 3.9784049
#> [35,] 5.9160824 6.1224530 3.5901933 3.9771856 6.6480483 4.4689927
#> [36,] 1.2679657 0.8595163 0.9923684 1.1253013 0.9086796 0.9049682
#> [37,] 0.4278154 0.4616109 1.1412999 1.2086568 1.2311754 1.3396608
#> [38,] 0.6349455 0.6988468 0.5034284 0.5264484 0.4015423 0.6265682
#> [39,] 0.9155948 0.7344261 0.9032604 0.8983678 1.1621279 0.6766302
#> [40,] 1.1866301 1.1574085 1.8215388 1.8345997 1.5133395 1.7667248
#> [41,] 1.9128037 1.9797422 2.4108030 2.4415413 1.3692384 3.1095613
#> [42,] 2.8794621 2.5897007 1.8031866 2.1596559 3.4151957 1.4243539
#> [43,] 2.4472745 5.6335781 3.6072167 3.6468950 3.2112939 5.1592206
#> [44,] 2.6708724 2.9226659 3.3465654 3.4584780 2.8867777 4.1977383
#> [45,] 3.5876312 2.6104609 1.5946506 1.5458651 1.5433350 1.7933679
#> [46,] 0.5527721 0.3651031 0.7600383 0.6504084 0.7915350 0.8827272
#> [47,] 0.6619081 0.5915694 0.3672173 0.4580523 0.4821844 0.3909723
#> [48,] 10.5641362 10.9923336 5.6013218 6.4229534 9.7729485 6.6906957
#> [49,] 48.8907673 49.4698676 28.2112388 33.0195021 51.8609854 42.3732651
#> [50,] 0.1824876 0.1991927 0.2408143 0.2346960 0.2363074 0.2345646
#> [51,] 0.6377801 0.5744148 0.5151280 0.6127882 0.5833437 0.5050398
#> [52,] 5.5919518 8.3004662 2.6825565 3.9316850 4.2389475 2.6810031
#> [53,] 0.8843105 0.7008193 0.4877147 0.5975779 0.7527168 0.6136767
#> [54,] 2.3642281 2.7973771 3.4687460 2.9561199 2.3966877 3.2277530
#> [55,] 0.8599228 0.8499212 0.8602408 0.7773863 0.8492159 0.9385647
#> [56,] 0.7212844 0.9809029 0.7111322 0.6957430 0.6937791 0.5856716
#> [57,] 0.2843839 0.2869878 0.4434099 0.4856307 0.4164912 0.3786716
#> [58,] 1.0593126 0.7381490 0.4461725 0.5244459 0.5027276 0.4755898
#> [59,] 0.9673642 0.6018039 0.7596030 0.6789657 0.6619295 0.5084909
#> [60,] 0.6030086 0.5218243 0.6543594 0.6135507 0.9897508 0.7874383
#> [61,] 1.8093523 2.0420776 1.4162070 1.5906339 2.1664705 1.1476253
#> [62,] 1.1167809 1.4731507 1.0001074 1.2661611 0.8877320 0.9703858
#> [63,] 3.8007778 6.2242685 4.8593281 4.4234316 4.6970430 4.2551848
#> [64,] 5.8251062 3.3137192 2.7571586 2.2472189 2.4256863 3.0176948
#> [65,] 4.9461600 3.0395836 4.7238590 4.7333726 3.7579359 3.1874711
#> [66,] 2.4471799 2.7632394 3.1846279 2.6398936 3.5463713 4.0190203
#> [67,] 1.3777531 2.1547203 0.9511904 0.9790176 1.1928718 0.8727624
#> [68,] 8.7013827 11.5954873 9.4056954 9.2759870 8.6081831 8.0375052
#> [69,] 4.4263707 4.0405200 3.6061903 3.8910196 2.8846797 2.9892399
#> [70,] 0.7352764 0.6884231 1.0285528 1.0003258 0.8865015 1.2438979
#> [71,] 11.8465388 16.5956305 10.8126345 13.3387790 13.5666349 11.0207850
#> [72,] 0.8253410 1.1548827 1.4985010 1.6421192 1.9209603 1.3428917
#> [73,] 5.0060112 3.9998930 4.0369723 3.6947510 2.8955976 2.4415009
#> [74,] 14.7040987 10.0763190 12.1873180 11.6566269 13.3499720 13.4313137
#> [75,] 0.9553195 1.2353689 4.1511362 4.3775995 3.6904803 4.6095284
#> [76,] 0.3037800 0.3368242 0.6133827 0.6271274 0.7876305 0.5512530
#> [77,] 0.3356046 0.3817858 0.4230596 0.5237951 0.7129618 0.4068187
#> [78,] 1.9115090 2.1763018 4.4521429 5.0186420 8.0286080 3.1953208
#> [79,] 0.8024421 0.4069146 0.9545429 1.1890632 1.0889542 0.7424340
#> [80,] 122.2695722 96.9067190 30.4911330 38.3139672 41.0244214 38.2550405
#> [81,] 0.6913379 0.7993532 1.0541382 1.0308733 0.7707784 1.0785932
#> [82,] 1.5768766 1.1523198 1.2820985 1.3725818 1.3696760 1.8827896
#> [83,] 31.7733362 41.4117796 17.4712515 21.1544887 26.5437864 16.7090136
#> [84,] 0.4264609 0.4601483 0.4776772 0.4760081 0.5680203 0.4628562
#> [85,] 3.0984054 2.4220470 2.7731943 2.7727827 1.4586283 3.5133094
#> [86,] 8.4535361 11.6230306 14.0553513 14.2393191 13.5595556 14.7181156
#> [87,] 1.0982229 1.5438235 1.8445090 2.0464888 1.3225480 1.0049953
#> [88,] 3.6847632 4.4146195 2.0459273 2.6544048 3.4789540 2.4883238
#> [89,] 0.3063149 0.2520365 0.3252853 0.3279702 0.3312064 0.2514629
#> [90,] 0.4588755 0.5067020 0.7401046 0.6760834 0.4559503 0.8558664
#> [91,] 0.4409784 0.2690826 0.3560064 0.3855393 0.3470677 0.3348954
#> [92,] 1.5575316 1.3417864 1.5074456 1.5259563 2.4045851 1.8354152
#> [93,] 0.8356742 0.7505371 1.2305568 1.1766346 1.0851128 1.0470381
#> [94,] 0.4037395 0.3532286 1.8969368 2.0296393 1.9891769 2.8062485
#> [95,] 22.0099773 13.5578631 5.0261551 5.4669510 7.2131244 5.8545052
#> [96,] 0.7122054 0.6795955 0.5168275 0.6156701 0.4813555 0.5596008
#> [97,] 0.6025621 0.7854466 1.2314205 1.3323393 1.5963386 1.1234846
#> [98,] 13.6701928 9.9807631 29.0228757 27.1251151 20.7612078 33.2565877
#> [99,] 0.5184963 0.5288918 0.6668451 0.5721096 0.6394911 0.6180270
#> [100,] 3.2253971 2.6351400 1.2138865 1.4851694 1.6974659 1.4549563
#> [,13] [,14] [,15] [,16]
#> [1,] 0.7801504 0.9832116 0.6131255 0.7043107
#> [2,] 2.5098488 2.8210042 2.3332712 3.1209359
#> [3,] 0.7998675 0.7192381 0.8026200 0.4754337
#> [4,] 2.5997023 2.7094006 3.2020482 1.9226119
#> [5,] 2.5551043 2.2949440 2.3195108 2.4518573
#> [6,] 0.8835520 0.9557190 1.2113228 1.0015142
#> [7,] 14.4298977 14.7284785 7.8744849 15.4525114
#> [8,] 0.7707189 0.7274023 1.0132372 0.6926530
#> [9,] 2.7539897 2.1598659 2.0176613 1.6215997
#> [10,] 2.6186546 2.5771822 2.0333760 2.2731384
#> [11,] 1.2005121 1.3141498 1.4017004 0.9510328
#> [12,] 1.4936301 1.3981869 1.5049462 1.1457495
#> [13,] 1.9015717 1.8978951 1.0128161 2.0312602
#> [14,] 0.3157133 0.3380934 0.3666783 0.4200485
#> [15,] 0.7564417 0.7851358 0.5970236 0.6932737
#> [16,] 2.9495683 2.8108437 2.4759466 4.0847171
#> [17,] 0.4135180 0.3600391 0.3553818 0.4704095
#> [18,] 0.2615115 0.2491214 0.2099516 0.2955285
#> [19,] 1.6002938 1.4794120 1.2834312 1.1542792
#> [20,] 0.2910850 0.3037398 0.3234771 0.2727842
#> [21,] 1.0045998 1.0321444 0.9113530 0.8635390
#> [22,] 183.8625564 182.3387782 184.3207216 186.1988739
#> [23,] 1.4512299 1.4229659 1.5830127 2.5313631
#> [24,] 1.0290970 1.1882520 0.8249629 0.9855272
#> [25,] 1.4123187 1.3936328 1.9755000 1.4632853
#> [26,] 25.7891850 25.1742806 31.5793899 31.7082366
#> [27,] 15.1576563 19.1266084 12.4124627 14.0953633
#> [28,] 1.2616858 1.1447525 1.5743840 1.0577564
#> [29,] 1.8065344 1.5724174 1.3704209 1.5029258
#> [30,] 11.1549077 11.0604101 7.1132151 7.7364311
#> [31,] 0.9418333 0.9563673 0.9540322 0.6177266
#> [32,] 0.9027846 0.9147217 0.9260617 1.1384308
#> [33,] 29.1193015 26.7752318 34.5913402 30.9770195
#> [34,] 2.9165617 2.3779797 2.6348713 2.9073611
#> [35,] 4.8046534 4.7019221 5.5760688 3.9255718
#> [36,] 0.8392400 0.9180816 1.0852719 0.7508497
#> [37,] 1.0471167 1.1132491 0.9117734 1.2885507
#> [38,] 0.4877064 0.4375998 0.3435985 0.9638696
#> [39,] 1.1908028 1.1794494 0.8288894 0.8175128
#> [40,] 1.9585483 1.9428099 1.0957512 1.8058861
#> [41,] 2.3859204 2.0625963 1.2115996 1.8237428
#> [42,] 2.1178400 2.1035119 2.9518985 1.8955397
#> [43,] 3.7096835 3.5013176 3.9185220 4.2970495
#> [44,] 3.1688152 3.4875674 2.9968712 3.0827858
#> [45,] 1.5807567 1.6929497 1.9227273 2.0631873
#> [46,] 0.9432532 0.8114366 0.7285444 1.0372871
#> [47,] 0.5046325 0.4825983 0.4297452 0.6667551
#> [48,] 7.0272715 5.6691089 8.2209736 8.7560829
#> [49,] 40.1541406 41.4680829 53.6648304 39.8154577
#> [50,] 0.2065137 0.1617573 0.2190104 0.2510260
#> [51,] 0.4943799 0.5310160 0.5921210 0.5105304
#> [52,] 4.2377400 4.7786346 4.5154413 4.6633187
#> [53,] 0.6917022 0.6552052 0.7167207 0.5273342
#> [54,] 2.2386584 2.4595884 2.4745560 2.2066433
#> [55,] 1.0594522 1.0610739 1.0663532 0.9552143
#> [56,] 0.7206595 0.6674797 0.7547826 0.6702581
#> [57,] 0.3187216 0.3076663 0.4186603 0.3826275
#> [58,] 0.4591697 0.5643836 0.6398599 0.4888548
#> [59,] 0.7370848 0.6866777 0.7187212 0.6543839
#> [60,] 0.7401341 0.8528639 1.0519319 0.7922992
#> [61,] 1.8942436 1.5628761 2.0968906 1.5954533
#> [62,] 1.3136025 1.3524492 0.7239225 1.2124028
#> [63,] 6.3306740 6.1738812 6.2427431 7.0208927
#> [64,] 2.2434220 2.3920850 3.1704581 2.1090812
#> [65,] 3.2482261 3.4781089 5.3739642 3.1164402
#> [66,] 3.3936345 3.0954939 3.1634538 3.5837885
#> [67,] 0.9283848 0.8732206 1.1044690 0.9840491
#> [68,] 9.1210320 9.3205592 8.6280703 8.9304411
#> [69,] 5.9267402 4.9900483 3.1062650 3.3277248
#> [70,] 0.8468930 0.9558838 0.7943335 0.8102525
#> [71,] 13.1281793 11.7497177 11.7964443 14.5092992
#> [72,] 1.7608203 1.5669338 1.7306061 1.8410296
#> [73,] 3.7580723 4.2145312 4.8607861 2.5782951
#> [74,] 13.8068947 15.4717478 12.6089940 11.1670747
#> [75,] 4.3958758 4.5581285 3.7774977 3.6412822
#> [76,] 0.5605331 0.5225130 0.5693564 0.5501085
#> [77,] 0.5180642 0.5746807 0.7829309 0.2895740
#> [78,] 4.6389317 5.7068845 7.3988768 4.4336274
#> [79,] 0.7550198 0.7649316 0.8729182 0.9406590
#> [80,] 44.6010456 48.4424396 51.4979355 48.0715458
#> [81,] 0.9256853 0.8999535 0.6692380 1.5870822
#> [82,] 1.8707728 1.8134628 1.5598892 1.7986883
#> [83,] 19.7155314 20.4566052 22.5681671 25.0347147
#> [84,] 0.5001006 0.4981018 0.4678546 0.5384461
#> [85,] 2.7054062 2.3717674 1.3635604 3.6438142
#> [86,] 9.4995985 10.1281932 9.3419386 11.0683392
#> [87,] 1.1481881 1.1938974 1.1125243 0.8592134
#> [88,] 1.8562812 1.7141613 1.7523162 1.9282786
#> [89,] 0.2676033 0.2673805 0.3319031 0.2839174
#> [90,] 0.5358542 0.5127611 0.4035959 0.5586551
#> [91,] 0.2681198 0.2497549 0.1737191 0.3115668
#> [92,] 1.5659524 1.4676368 2.1215335 1.8706740
#> [93,] 0.9762998 1.1405319 0.9227644 0.9817692
#> [94,] 1.5073507 1.5298211 1.7016845 1.7503515
#> [95,] 12.0326962 8.9856883 13.6522288 8.5910626
#> [96,] 0.7745168 0.9398637 0.6134516 0.8809995
#> [97,] 1.5501478 1.5405510 1.2023118 1.3905242
#> [98,] 30.6058503 29.0812542 25.5560371 39.4427496
#> [99,] 0.5110001 0.5336657 0.6126535 0.8810912
#> [100,] 1.2814733 1.3567249 2.1080053 1.0063125
#using the Brown-Resnick model with a covariance matrix
vario2cov <- function(coord, semivario,...){
sapply(1:nrow(coord), function(i) sapply(1:nrow(coord), function(j)
semivario(sqrt(sum((coord[i,])^2)), ...) +
semivario(sqrt(sum((coord[j,])^2)), ...) -
semivario(sqrt(sum((coord[i,]-coord[j,])^2)), ...)))
}
rmev(n=100, sigma=vario2cov(grid.coord, semivario = semivario, scale = 0.5, alpha = 1), model='br')
#> [,1] [,2] [,3] [,4] [,5] [,6]
#> [1,] 0.9552671 0.8239410 0.9239665 0.8649631 0.9609080 0.8051049
#> [2,] 1.6792763 1.3804349 1.3864842 2.3955415 1.6680526 1.3954869
#> [3,] 0.7116208 0.7689730 0.6784106 0.9869712 0.7237560 0.7804772
#> [4,] 1.6688561 1.3469676 1.3446439 1.1621214 1.6663137 1.3528515
#> [5,] 7.0573601 7.4404407 10.4258531 7.8196776 7.0196561 7.3775143
#> [6,] 115.4834847 117.0859303 90.4620688 148.4135705 112.9691566 118.3874161
#> [7,] 1.4389955 1.4308989 1.1889382 1.1107546 1.4199826 1.4376228
#> [8,] 0.3630235 0.4139229 0.5229194 0.3849355 0.3661393 0.4107830
#> [9,] 1.2761116 1.1927764 1.4723699 1.4611034 1.2321557 1.2036320
#> [10,] 0.8273277 0.7778192 0.8770229 0.6088137 0.8214141 0.7747612
#> [11,] 7.1027805 7.0123653 6.9435105 8.2220703 6.9579361 7.1550690
#> [12,] 2.3249471 2.4815695 2.8614459 2.3491276 2.3211762 2.5007465
#> [13,] 0.8589982 0.9476493 0.6766544 0.9736369 0.8809654 0.9551996
#> [14,] 1.8344014 1.9846538 3.0691413 1.4354778 1.8411711 1.9681971
#> [15,] 2.0682482 2.1238726 2.3642905 1.9654595 2.0696341 2.1258852
#> [16,] 1.6084239 1.7525447 1.5926353 1.2562319 1.5850635 1.7505342
#> [17,] 21.0153473 21.6972917 22.5015236 30.8069597 20.8894315 22.0674527
#> [18,] 6.6483499 6.6584946 4.6310059 4.9190483 6.6634374 6.6674027
#> [19,] 2.0500993 1.8337225 1.7074603 3.2181394 2.0681284 1.8638484
#> [20,] 2.4697758 2.5739317 3.2444240 2.3597201 2.4455217 2.5907415
#> [21,] 3.1763100 2.7748754 2.9260783 4.8793949 3.1248016 2.7398827
#> [22,] 0.9212928 0.9247121 0.5101077 0.8937884 0.9055090 0.9241193
#> [23,] 1.4940513 1.3021833 1.0555241 1.9797422 1.5119413 1.3268184
#> [24,] 2.5860023 2.6128797 2.8603012 3.0345403 2.5819016 2.5978561
#> [25,] 33.5565776 33.6563161 37.2847330 44.9412860 33.8740006 33.8440050
#> [26,] 0.8080512 0.7215496 1.3465935 0.7950086 0.7912732 0.7292398
#> [27,] 2.1665648 1.8743660 1.2855674 1.9897507 2.1341732 1.8993848
#> [28,] 3.8545865 3.5799341 3.1329661 4.0114111 3.8429322 3.5705577
#> [29,] 1.2737153 1.1161459 1.0737414 1.0340607 1.2969921 1.1312603
#> [30,] 2.5494535 2.7784962 2.6359623 2.7267642 2.5596730 2.7621680
#> [31,] 17.1731383 14.1935116 10.3934871 14.5750721 17.1380111 14.1542393
#> [32,] 1.8122162 1.7744005 1.8492381 1.6689874 1.8144018 1.7816962
#> [33,] 1.5554239 1.8131881 1.8142226 1.2860842 1.5444537 1.7985796
#> [34,] 1.5873055 1.7840292 1.3577625 1.5197614 1.5598532 1.7680469
#> [35,] 5.8178994 5.5172212 4.2411801 4.5442846 6.0268445 5.5250249
#> [36,] 0.2591402 0.2331119 0.2445813 0.2810436 0.2566708 0.2314786
#> [37,] 8.2774336 8.3824816 12.3802806 11.8082186 8.1190488 8.3737953
#> [38,] 0.6248818 0.5137158 0.7581339 0.6262986 0.6193309 0.5178714
#> [39,] 0.8970550 0.8451497 0.5559656 1.2085795 0.8760455 0.8505935
#> [40,] 3.4546984 4.1215569 6.0375148 2.3627259 3.4963190 4.1102780
#> [41,] 0.3240369 0.3527772 0.4121134 0.3675735 0.3273126 0.3372511
#> [42,] 0.5221244 0.5632913 0.5382160 0.7595857 0.5149304 0.5584240
#> [43,] 2.1836178 1.9122341 2.1285857 2.0728696 2.1636915 1.9177345
#> [44,] 0.8230207 0.7663601 0.8317321 1.0395485 0.8170397 0.7657896
#> [45,] 0.3154054 0.2901851 0.2450141 0.3100774 0.3118271 0.2943616
#> [46,] 2.8744926 2.4709649 1.9912518 1.8389364 2.8985735 2.4273954
#> [47,] 31.8678901 32.4209808 22.4523650 33.4681593 31.6518221 32.2647413
#> [48,] 0.6439568 0.5793491 0.6138764 0.7017163 0.6261408 0.5632026
#> [49,] 5.8783228 6.4924426 6.9585010 8.1920795 5.8494155 6.6017423
#> [50,] 1.7843720 1.5721478 1.0419077 1.2559903 1.7973786 1.5861460
#> [51,] 12.7300491 13.7788792 14.8977064 18.1199375 12.6893712 14.0583808
#> [52,] 0.5584420 0.5112161 0.3753817 0.6360252 0.5495245 0.5079303
#> [53,] 0.8373551 0.8190413 0.9629574 0.8050972 0.8256573 0.8213187
#> [54,] 5.7597803 5.6392310 5.7869796 5.3858897 5.6945335 5.6347436
#> [55,] 0.5314800 0.6443570 0.4848834 0.4582246 0.5402581 0.6507913
#> [56,] 0.6603934 0.7351723 0.8560847 0.5365490 0.6594406 0.7264439
#> [57,] 0.7252982 0.8192870 1.1138514 0.6212145 0.7336798 0.8376423
#> [58,] 18.5970408 21.7386329 15.1835215 18.7080167 18.6651685 22.6443881
#> [59,] 0.4190977 0.3952901 0.3830081 0.2996396 0.4253309 0.3957508
#> [60,] 0.3372667 0.3469880 0.3088262 0.3380081 0.3345892 0.3458392
#> [61,] 0.7216769 0.8310813 0.7469832 0.7146127 0.6998135 0.8180594
#> [62,] 0.4121905 0.4174860 0.4769077 0.3825429 0.4070892 0.4045093
#> [63,] 1.0851038 1.0336014 1.0169100 1.0584668 1.0440040 1.0259809
#> [64,] 2.7020666 2.9775249 3.0370208 3.6475725 2.7974089 2.9493141
#> [65,] 8.8908645 8.2899240 13.3151995 10.7407399 8.9277782 8.1782719
#> [66,] 1.8241323 1.7909194 1.8934064 2.3251136 1.8224495 1.7616442
#> [67,] 0.5808203 0.4687600 0.4994300 0.5146759 0.5673811 0.4778570
#> [68,] 0.7483445 0.9032570 1.2724710 0.8462357 0.7476898 0.8982078
#> [69,] 0.4315063 0.5270225 0.6784097 0.5528307 0.4373871 0.5297048
#> [70,] 1.6993956 1.9233384 1.5146843 2.0819263 1.7566680 1.9485753
#> [71,] 34.8962918 35.2877905 38.4252481 44.4426910 34.1063138 34.7279516
#> [72,] 0.9444279 0.9942802 0.8101688 1.0717454 0.9442053 0.9819765
#> [73,] 1.9831821 2.1254927 1.4115970 1.5468113 1.9625122 2.1780648
#> [74,] 0.7281167 0.6350117 0.5761523 0.9107843 0.7247970 0.6326282
#> [75,] 0.2886131 0.2580392 0.2279038 0.2678929 0.2890978 0.2590903
#> [76,] 0.9284595 0.7279607 0.5176531 1.0332164 0.9273088 0.7151318
#> [77,] 3.6249238 3.1109648 2.4524139 3.5930707 3.6347981 3.0685890
#> [78,] 0.9699325 1.0275578 0.8966013 0.9504462 0.9387382 1.0283654
#> [79,] 1.0290044 0.9430781 1.1863030 0.9905569 1.0270085 0.9237549
#> [80,] 5.6852514 5.1857073 5.0438681 4.4729294 5.7842636 5.1545939
#> [81,] 101.8945097 129.7589880 177.7059723 106.7405758 103.8513021 128.2821782
#> [82,] 14.6178138 13.8444660 7.9614868 11.3790029 14.6581980 13.6263697
#> [83,] 1.9981264 2.2838800 1.6491452 1.3539511 1.9657855 2.3034218
#> [84,] 1.0976254 1.0507418 0.8251977 1.0347192 1.0909318 1.0568619
#> [85,] 2.8897934 2.8679755 3.3266434 2.7548547 2.9334071 2.8851997
#> [86,] 35.5798447 37.0927609 31.9375694 37.5498948 34.9762997 38.0407685
#> [87,] 1.1930785 1.0820073 0.9290892 0.8840933 1.2078100 1.0414333
#> [88,] 6.9087908 6.8596673 5.5112696 5.8315646 6.9470310 6.8602417
#> [89,] 3.1372227 3.1130273 3.3919333 2.0355384 3.1868734 3.1209416
#> [90,] 0.9190769 0.9368618 0.9159166 0.8452622 0.9266443 0.9394305
#> [91,] 0.3931971 0.3005444 0.3292441 0.3164939 0.3889489 0.3067694
#> [92,] 1.4187806 1.4739575 0.9939337 1.6530507 1.3997986 1.4936667
#> [93,] 1.8102510 2.0078190 1.7470584 1.7780243 1.7796812 2.0172892
#> [94,] 0.3908485 0.3374217 0.3076896 0.6170973 0.3928350 0.3344726
#> [95,] 1.0361291 1.1417540 1.4929033 1.6099914 1.0262216 1.1634291
#> [96,] 3.6698356 3.6923245 4.5641011 4.1305256 3.6360781 3.7460189
#> [97,] 0.8234268 0.7996377 0.5703694 0.5130850 0.8287182 0.7985064
#> [98,] 34.6791075 31.9246750 27.3320101 33.1318508 34.9379381 31.6991628
#> [99,] 1.4470835 1.4568453 0.7986627 1.4603915 1.4398484 1.4551076
#> [100,] 27.7630657 27.7760193 23.6317602 27.1142706 28.1228115 28.1276211
#> [,7] [,8] [,9] [,10] [,11] [,12]
#> [1,] 0.9276581 0.8711558 0.8345461 0.9144693 0.9612704 0.5618196
#> [2,] 1.3959686 2.3603998 1.3522902 1.2188319 1.1039952 1.4546472
#> [3,] 0.6760000 0.9732182 0.7405774 0.6049569 0.7019892 0.7681583
#> [4,] 1.3779540 1.1769832 1.8140069 1.7827704 1.8535070 1.4179226
#> [5,] 10.2612620 7.8237650 8.3017115 8.5026379 6.2058924 12.5943861
#> [6,] 90.2688772 150.1072163 158.1675799 184.2461092 127.6222267 143.2065977
#> [7,] 1.1679030 1.0913959 1.0159751 1.2033727 1.7989657 0.9059334
#> [8,] 0.5185797 0.3930768 0.4334739 0.4394312 0.4553650 0.5457015
#> [9,] 1.4765642 1.4721726 1.6265212 1.4334376 1.4826593 2.1182089
#> [10,] 0.8783827 0.6145625 1.1636958 1.1612188 1.2478997 1.2112001
#> [11,] 6.7569330 8.1553524 7.3561438 6.3899213 5.5266641 8.3069006
#> [12,] 2.8714873 2.3794120 2.8483668 3.5116056 3.5110587 2.6383367
#> [13,] 0.6716830 0.9652916 0.4527130 0.4913808 0.5106477 0.5674330
#> [14,] 3.0568437 1.4314222 1.2533103 1.6109107 1.0882844 1.1365843
#> [15,] 2.3395838 1.9805091 1.7113392 1.4903562 1.2225371 1.9987503
#> [16,] 1.6279924 1.2520409 1.3154838 0.9594803 1.2862647 1.0112481
#> [17,] 22.2217510 30.8280916 26.3631892 23.7987031 23.0624894 31.0555021
#> [18,] 4.6679371 4.8651645 3.0121249 3.3138945 2.3370107 3.4699653
#> [19,] 1.6817344 3.2478374 1.9622523 1.5280881 1.7339669 2.4798573
#> [20,] 3.2488562 2.3696062 2.9193725 3.0562303 2.7332286 3.0677359
#> [21,] 2.8663577 4.8920149 3.5807738 3.0449386 3.7859618 3.6052103
#> [22,] 0.5109555 0.8851277 0.4820385 0.5695639 0.4762065 0.5378858
#> [23,] 1.0441980 1.9651511 0.8418195 0.9751856 0.8450702 1.1882387
#> [24,] 2.8805716 2.9840201 1.6697138 1.8146611 2.2820942 1.4303349
#> [25,] 37.5881136 45.0914151 12.6716261 13.0122706 10.7937968 21.0418855
#> [26,] 1.3377815 0.7949111 0.8845952 0.7953465 0.6880555 0.9986477
#> [27,] 1.2683425 1.9627909 2.3804744 2.9704886 2.9705295 1.9001269
#> [28,] 3.0831522 4.0844702 2.0924588 1.9319866 1.8109059 2.0500149
#> [29,] 1.0905893 1.0188699 0.7826432 0.7040166 0.6322806 0.7988269
#> [30,] 2.6340930 2.7257823 4.2511838 4.2873744 2.1194863 4.3233710
#> [31,] 10.5856945 14.5354414 9.9664714 9.4073547 12.2086214 15.9714683
#> [32,] 1.8413088 1.6850419 1.4850649 1.2458285 1.1373714 1.6346149
#> [33,] 1.8493517 1.2578037 1.5369816 1.7787223 1.9127577 1.2121638
#> [34,] 1.3460723 1.5091209 0.8530474 0.7841969 0.9744479 0.9020648
#> [35,] 4.2360260 4.5685756 5.6312317 6.0868696 5.1569211 3.6887826
#> [36,] 0.2467479 0.2778073 0.5018238 0.5114197 0.4968671 0.3501075
#> [37,] 12.2567064 11.7942424 10.7378411 9.5762080 12.7493093 11.4549484
#> [38,] 0.7343191 0.6127161 0.7723970 0.8131345 0.6770285 0.9196005
#> [39,] 0.5613365 1.2097913 1.0455354 1.2022589 1.1300707 1.3515585
#> [40,] 5.9854614 2.4329355 3.8041464 4.3449054 5.0808820 3.4399767
#> [41,] 0.4190005 0.3622250 0.3221913 0.3724543 0.3522903 0.2616515
#> [42,] 0.5301623 0.7546655 0.5576066 0.6562208 0.5167878 0.7961032
#> [43,] 2.0865806 2.0760142 1.9687346 1.8638016 1.9534835 1.6474749
#> [44,] 0.8463389 1.0231345 0.9138032 0.9212682 1.0720119 0.7755999
#> [45,] 0.2416482 0.3106846 0.3894938 0.4958418 0.4934578 0.4508711
#> [46,] 1.9795674 1.8777563 3.1498484 2.9695849 2.4290216 2.1534417
#> [47,] 22.6754507 33.0872556 58.3605039 41.9442816 61.6212960 74.6327271
#> [48,] 0.6081016 0.7156690 1.8239852 1.8152447 1.3930599 1.2790941
#> [49,] 6.8342425 8.2072452 14.9787879 15.8583958 10.4650076 11.9159859
#> [50,] 1.0370975 1.2475686 1.5209585 1.4715544 1.6087200 1.3811728
#> [51,] 14.9099118 18.2819672 15.4956336 16.9454145 14.4395648 28.1972956
#> [52,] 0.3745153 0.6352147 0.5624924 0.6664396 0.5033078 0.5931516
#> [53,] 0.9790109 0.8052395 0.9863786 0.9502209 1.0509421 1.2920417
#> [54,] 5.7711708 5.4835680 10.4035898 11.8692503 9.6268634 11.3934139
#> [55,] 0.4885832 0.4611290 0.5153245 0.4927255 0.4027637 0.4683790
#> [56,] 0.8511615 0.5289000 0.4834076 0.4000947 0.3555422 0.4599993
#> [57,] 1.1393045 0.6253108 1.4090169 1.3206787 1.3203572 1.1719475
#> [58,] 15.3393034 19.0653351 10.3146740 11.0290343 7.8208585 6.3892621
#> [59,] 0.3816798 0.2970917 0.4313216 0.4481727 0.4803382 0.4499174
#> [60,] 0.3080139 0.3357060 1.0220383 0.8257369 0.4567392 0.7543959
#> [61,] 0.7335348 0.7129037 0.2993866 0.3599542 0.4725668 0.3254475
#> [62,] 0.4795676 0.3798617 0.5830165 0.5903291 0.8150045 0.5417720
#> [63,] 1.0033880 1.0492971 0.6240121 0.5953918 0.5513940 0.5829681
#> [64,] 3.0575732 3.6789760 2.2358514 1.9908756 1.7342959 2.7287913
#> [65,] 13.3116025 10.8387423 6.0848315 7.2283451 9.8236661 5.8722041
#> [66,] 1.8820880 2.2808933 1.1847087 1.2922557 2.3984328 1.0282010
#> [67,] 0.4897675 0.5218234 1.0222660 1.1460242 0.8742729 0.8430843
#> [68,] 1.2669302 0.8511037 1.0672264 1.1849762 1.1137209 0.7303788
#> [69,] 0.6787380 0.5437174 0.5226750 0.5123628 0.6709751 0.5342334
#> [70,] 1.5096215 2.0458298 0.7668750 0.8334828 0.6945066 0.8358163
#> [71,] 36.8879335 44.0203586 20.9311735 19.2729453 25.6102451 22.9135836
#> [72,] 0.8119489 1.0889944 1.3099489 1.3813836 1.1086771 0.9140541
#> [73,] 1.4404999 1.5551567 2.9380874 2.7830043 2.1918715 3.5771396
#> [74,] 0.5776620 0.9106247 0.9527030 0.9016141 0.7414328 0.9793334
#> [75,] 0.2310827 0.2624110 0.4040236 0.3819930 0.2931648 0.3485965
#> [76,] 0.5209974 1.0401147 0.7778763 0.8799144 1.0130236 0.7365147
#> [77,] 2.4865165 3.5774189 2.3362514 2.0917854 2.3646037 2.1783240
#> [78,] 0.9144644 0.9727734 1.2524457 1.2671617 1.0613560 1.2475258
#> [79,] 1.2184995 0.9901187 1.2818801 1.4487827 1.1371578 1.7494384
#> [80,] 5.0333956 4.4519609 3.4265297 4.5306773 6.0688340 3.1652451
#> [81,] 177.5649579 107.2553365 86.6246157 86.8824169 75.8795800 56.7527169
#> [82,] 8.0347938 11.4904584 10.3914213 12.2930546 12.3860621 9.3002868
#> [83,] 1.6485516 1.3726146 1.0750327 1.0044464 1.0009631 1.1604727
#> [84,] 0.8252038 1.0227531 1.3535502 1.3009545 0.8727275 1.7824100
#> [85,] 3.3847235 2.7081252 3.0819249 2.7065553 3.2272662 3.0770470
#> [86,] 32.2465960 37.7221570 24.0166775 24.1140936 12.8449122 14.2877147
#> [87,] 0.9361515 0.8756168 2.7559580 2.5004395 2.2857260 2.4839009
#> [88,] 5.5037760 5.9204111 7.0420391 8.0166585 4.6564044 6.7122190
#> [89,] 3.3753804 2.0721142 3.7009022 3.8495715 3.9157499 4.3691371
#> [90,] 0.8900034 0.8371123 1.3622845 1.4124928 1.1964872 1.5419117
#> [91,] 0.3314075 0.3131258 0.4426594 0.4583077 0.3701031 0.3375262
#> [92,] 0.9971145 1.6612244 1.6605159 2.0311503 1.8781906 1.6287944
#> [93,] 1.7181269 1.8025459 1.2242579 1.5303369 1.3252719 1.3312726
#> [94,] 0.3052065 0.6081318 0.4595480 0.3783969 0.4639368 0.4615956
#> [95,] 1.4965619 1.5887283 1.4392368 1.2698642 1.0884747 1.3802794
#> [96,] 4.6098473 4.1490417 3.6470306 3.3314620 3.4395327 4.5351950
#> [97,] 0.5750426 0.5197097 0.4942522 0.5412504 0.6785530 0.4847582
#> [98,] 26.9075586 34.0431721 20.1769295 19.9736033 16.4247242 24.7779439
#> [99,] 0.8039579 1.4559722 1.0469096 0.9626254 0.7674626 1.2958720
#> [100,] 22.7065169 27.0708847 43.0179408 41.9064635 28.5676799 41.8112895
#> [,13] [,14] [,15] [,16]
#> [1,] 0.7004544 0.6954706 1.1159438 0.5500656
#> [2,] 1.6116646 1.6718327 1.6872022 1.5680935
#> [3,] 0.5371208 0.5800821 0.5072287 0.5508722
#> [4,] 1.4603837 1.5996908 1.4544103 1.6986315
#> [5,] 11.4573961 9.8042757 10.2363401 17.3113129
#> [6,] 124.2929378 137.3900497 79.2390159 123.1950886
#> [7,] 0.8769198 0.9780093 1.2638389 0.9987225
#> [8,] 0.4165916 0.4523917 0.4936134 0.4129824
#> [9,] 1.5072663 1.3590433 1.7270184 1.7812140
#> [10,] 0.9291882 1.1269511 1.2369357 0.6882308
#> [11,] 8.2027324 8.8744262 6.4212383 9.7304057
#> [12,] 3.2607947 3.4880601 3.3613172 2.3356870
#> [13,] 0.5319334 0.4553311 0.4527653 0.6066661
#> [14,] 1.1713595 1.3043197 1.0372446 1.2606793
#> [15,] 1.6836696 1.6987869 1.0745561 2.2160328
#> [16,] 1.3111090 1.2232295 1.5480935 1.4699050
#> [17,] 18.8180597 19.3997762 24.6318549 33.6030140
#> [18,] 3.4388680 3.4043437 2.8389870 5.1153883
#> [19,] 2.4242289 2.3412059 1.9826303 2.9397149
#> [20,] 2.0070791 1.8204382 3.1031261 2.0700823
#> [21,] 4.1375800 5.1308056 4.2398299 4.2121138
#> [22,] 0.5788159 0.5044782 0.6883838 0.4490326
#> [23,] 1.3452730 1.4598261 1.1164470 1.4301322
#> [24,] 1.5371703 1.8767509 2.4334540 1.4267214
#> [25,] 21.1205443 22.7709707 18.9228362 30.5963171
#> [26,] 0.6235878 0.6499179 0.6625614 0.8531072
#> [27,] 1.9589203 2.3927045 2.4842570 1.5035062
#> [28,] 2.3716510 2.5013376 1.6589273 2.4177917
#> [29,] 0.7451405 0.7964479 0.6989973 0.8648907
#> [30,] 3.7310894 3.2519835 2.2864218 4.2777013
#> [31,] 13.0779219 13.7070181 12.0505800 13.0797156
#> [32,] 1.5982789 1.5659633 1.3322754 2.3944701
#> [33,] 2.5344433 2.5305880 2.1012637 1.5451380
#> [34,] 1.2711157 0.9493422 1.1482516 1.2345941
#> [35,] 5.6924121 6.2731608 5.8641830 5.6886859
#> [36,] 0.4065413 0.3508270 0.2739387 0.4242328
#> [37,] 11.5402120 12.3731706 11.0086799 16.5916151
#> [38,] 0.9740256 0.9074632 0.7871957 0.9196522
#> [39,] 1.5280900 1.4696137 1.3704188 1.8120499
#> [40,] 4.6104028 5.9946893 6.0945876 4.0868296
#> [41,] 0.3704825 0.3476781 0.2967906 0.3139205
#> [42,] 0.7484104 0.8361242 0.7079949 0.7684774
#> [43,] 1.3987694 1.2565395 1.2647602 1.4634188
#> [44,] 1.0007435 0.8031443 0.7113674 0.9603913
#> [45,] 0.3839303 0.4073644 0.3010277 0.3704116
#> [46,] 3.4757015 3.9313996 4.3327927 3.5700351
#> [47,] 61.3610036 50.9390590 54.6404077 69.4829415
#> [48,] 2.0841323 2.3874638 1.5869389 1.6470562
#> [49,] 13.4118677 12.8401452 8.2676744 12.2870584
#> [50,] 1.3746573 1.5099291 1.6220857 1.5498900
#> [51,] 14.1994687 11.6535610 14.4158453 15.6194654
#> [52,] 0.8032660 0.6183875 0.7310849 0.6623546
#> [53,] 0.9973230 0.9211543 0.8656659 1.2380795
#> [54,] 9.1754554 8.3166759 6.4347060 9.7506856
#> [55,] 0.4045753 0.3585052 0.3985622 0.5217842
#> [56,] 0.6596404 0.6726087 0.6995014 0.4336552
#> [57,] 1.2185731 1.1679828 1.3155526 1.6581463
#> [58,] 11.3075856 15.7102422 12.8783730 10.1345000
#> [59,] 0.5590052 0.5463224 0.5345983 0.5411223
#> [60,] 0.5417659 0.4651287 0.4873044 0.6063740
#> [61,] 0.4557064 0.5363349 0.6431590 0.4493267
#> [62,] 0.6041061 0.5894839 0.8566687 0.4812772
#> [63,] 0.6737876 0.6986896 0.4613802 0.8785750
#> [64,] 1.9400777 2.0408771 1.8913685 2.5719659
#> [65,] 7.2998380 7.9303347 9.6117613 5.5389799
#> [66,] 1.3805284 1.3211816 2.0789673 1.8107521
#> [67,] 0.7286837 0.7523700 0.6421802 0.6756338
#> [68,] 1.1163030 1.0076389 1.1200342 1.0459828
#> [69,] 0.4923139 0.4469592 0.4974626 0.4891862
#> [70,] 0.8443584 0.8949859 0.8401181 0.6627140
#> [71,] 28.3252616 27.8744369 25.4549759 30.5620998
#> [72,] 0.9683329 0.8364375 1.0375586 0.8604148
#> [73,] 2.5200796 2.4036721 2.0842551 2.9471482
#> [74,] 0.6839707 0.7458511 0.7568274 0.9662435
#> [75,] 0.2927912 0.3178749 0.3066878 0.3620511
#> [76,] 0.9586500 0.9315728 0.8928675 0.5502877
#> [77,] 2.2018388 1.9607483 2.2413014 3.5249455
#> [78,] 1.6716094 1.5419292 1.1408904 1.4829938
#> [79,] 1.6784467 1.4674304 0.9483450 1.7436843
#> [80,] 5.7309965 4.6193502 6.8709310 4.3185504
#> [81,] 77.3443236 71.1597443 70.1592595 79.6690656
#> [82,] 10.2913827 12.8048983 10.0482813 7.7306954
#> [83,] 0.9575612 1.0331514 1.1564878 1.2831273
#> [84,] 0.9560426 0.9777432 0.8249844 1.2370079
#> [85,] 3.0376983 3.3024882 3.5827953 2.7707433
#> [86,] 14.1033549 13.8492178 17.8139690 16.9979909
#> [87,] 1.9460200 1.8419539 1.6216618 1.5237177
#> [88,] 6.2792919 7.9577765 6.5726359 5.5107908
#> [89,] 3.9459534 3.5204438 3.7709779 3.1013408
#> [90,] 1.0488221 0.9877006 0.9788676 0.9916548
#> [91,] 0.5674201 0.5972529 0.4061772 0.3575116
#> [92,] 2.1224583 2.1531078 1.4289631 1.6433819
#> [93,] 1.2416561 1.3341480 1.6493077 1.0808190
#> [94,] 0.3463568 0.3895353 0.4095485 0.3723621
#> [95,] 1.0211914 1.0060158 1.0888198 1.9207936
#> [96,] 5.1161636 5.1840377 3.8852212 5.1319459
#> [97,] 0.6096065 0.5679481 0.6687911 0.4358321
#> [98,] 20.5834326 18.5607625 13.9492318 21.3083370
#> [99,] 0.9888304 1.1753384 0.7454012 1.1111812
#> [100,] 36.4122254 35.7510762 24.7819894 43.7696347
# asymmetric logistic model - see function 'rmvevd' from package 'evd '
asy <- list(0, 0, 0, 0, c(0,0), c(0,0), c(0,0), c(0,0), c(0,0), c(0,0),
c(.2,.1,.2), c(.1,.1,.2), c(.3,.4,.1), c(.2,.2,.2), c(.4,.6,.2,.5))
rmev(n=1, d=4, param=0.3, asy=asy, model="alog")
#> [,1] [,2] [,3] [,4]
#> [1,] 2.954491 0.4630662 3.958977 0.7101828
#Example with a grid (generating an array)
rmev(n=10, sigma=cbind(c(2,1), c(1,3)), coord=cbind(runif(4), runif(4)), model='smith', grid=TRUE)
#> , , 1
#>
#> [,1] [,2] [,3] [,4]
#> [1,] 12.98884 15.42073 14.37246 10.870865
#> [2,] 12.51703 14.61044 13.72464 10.611359
#> [3,] 10.61326 11.93471 11.40633 9.254866
#> [4,] 13.03755 15.52108 14.44763 10.889014
#>
#> , , 2
#>
#> [,1] [,2] [,3] [,4]
#> [1,] 0.4416108 0.4208311 0.4342221 0.4364541
#> [2,] 0.4582573 0.4293436 0.4464993 0.4587588
#> [3,] 0.4571634 0.5088429 0.4365965 0.4707586
#> [4,] 0.4379975 0.4185344 0.4313044 0.4319857
#>
#> , , 3
#>
#> [,1] [,2] [,3] [,4]
#> [1,] 2.355028 2.196908 2.289272 2.365336
#> [2,] 2.355478 2.160341 2.268922 2.396360
#> [3,] 2.167283 2.760912 2.293635 2.267986
#> [4,] 2.349694 2.197954 2.287457 2.355088
#>
#> , , 4
#>
#> [,1] [,2] [,3] [,4]
#> [1,] 15.17152 24.74092 19.90900 9.987556
#> [2,] 19.68619 31.56283 25.59890 13.127075
#> [3,] 32.09137 49.56827 40.90207 22.011330
#> [4,] 14.51361 23.73305 19.07374 9.534645
#>
#> , , 5
#>
#> [,1] [,2] [,3] [,4]
#> [1,] 21.95604 18.09137 19.96637 24.22236
#> [2,] 22.34710 18.10365 20.13750 24.97237
#> [3,] 21.36592 16.67508 18.87137 24.55907
#> [4,] 21.84458 18.04898 19.89431 24.04944
#>
#> , , 6
#>
#> [,1] [,2] [,3] [,4]
#> [1,] 2.173076 1.776555 1.661140 3.061310
#> [2,] 2.099504 1.502484 1.590329 2.995885
#> [3,] 1.790214 1.011296 1.329146 2.627635
#> [4,] 2.180324 1.821240 1.669138 3.065154
#>
#> , , 7
#>
#> [,1] [,2] [,3] [,4]
#> [1,] 2.349687 3.551346 2.960053 2.028936
#> [2,] 2.592107 3.851796 3.235799 2.056771
#> [3,] 2.958021 4.234604 3.619322 2.521738
#> [4,] 2.307532 3.497212 2.911234 2.054575
#>
#> , , 8
#>
#> [,1] [,2] [,3] [,4]
#> [1,] 2.639554 2.470003 2.570146 2.644878
#> [2,] 2.294501 2.110972 2.213883 2.328839
#> [3,] 1.551176 1.374855 1.466983 1.619437
#> [4,] 2.693963 2.527841 2.626993 2.693801
#>
#> , , 9
#>
#> [,1] [,2] [,3] [,4]
#> [1,] 1.153961 0.9933452 1.0743339 1.231656
#> [2,] 1.020632 0.8637854 0.9415786 1.103427
#> [3,] 1.095722 0.8444621 0.9612698 1.271527
#> [4,] 1.174456 1.0137649 1.0950274 1.250933
#>
#> , , 10
#>
#> [,1] [,2] [,3] [,4]
#> [1,] 2.514571 2.813753 2.695322 2.200923
#> [2,] 2.454153 2.699921 2.606679 2.175798
#> [3,] 2.139679 2.267773 2.227574 1.951271
#> [4,] 2.518834 2.826266 2.703874 2.200085
#>
## Example with Dirichlet mixture
alpha.mat <- cbind(c(2,1,1),c(1,2,1),c(1,1,2))
rmev(n=100, param=alpha.mat, weights=rep(1/3,3), model='dirmix')
#> [,1] [,2] [,3]
#> [1,] 1.7256102 0.5846787 1.1427186
#> [2,] 0.3499582 1.0519567 0.9258928
#> [3,] 0.3256478 0.2464800 0.1828500
#> [4,] 0.5038384 0.2373393 0.7204479
#> [5,] 0.9115554 1.9472892 0.8383267
#> [6,] 2.5796791 5.8154164 3.0827486
#> [7,] 0.5328958 0.6627461 0.5054585
#> [8,] 3.3110091 2.5891914 1.2088149
#> [9,] 7.1328600 1.1367244 1.3666279
#> [10,] 1.5899568 3.0591357 0.5509226
#> [11,] 1.1638077 0.7385675 1.5980405
#> [12,] 251.2600971 29.7094401 262.1621587
#> [13,] 9.8034612 0.9468627 2.1962425
#> [14,] 1.8339534 0.2095615 1.3366477
#> [15,] 1.2540402 3.4624393 0.9967765
#> [16,] 0.5138806 0.8744057 0.3119824
#> [17,] 0.7501497 3.5970049 1.0371341
#> [18,] 2.2099216 2.0022036 2.3293946
#> [19,] 0.7798161 1.7906425 1.5114915
#> [20,] 14.1354388 2.8254560 3.4329760
#> [21,] 19.8999725 9.5319927 63.0690240
#> [22,] 0.2544788 0.3816362 0.4610124
#> [23,] 1.0432453 0.5826977 1.5833175
#> [24,] 2.0639982 1.5465400 0.4877852
#> [25,] 2.1993896 1.1702199 2.3341204
#> [26,] 1.9557143 1.1727705 1.4510191
#> [27,] 0.7937444 0.8399898 0.7272106
#> [28,] 3.4189593 2.0882206 3.7517532
#> [29,] 34.9650995 31.4647833 25.6985467
#> [30,] 2.7376336 15.5567680 2.6707741
#> [31,] 0.8313366 0.6446365 0.5471616
#> [32,] 2.5945287 8.9313678 2.0714382
#> [33,] 5.4058132 1.0511733 2.3075330
#> [34,] 13.4604881 5.9007424 7.1857961
#> [35,] 0.7878278 3.2986657 1.4374290
#> [36,] 2.8630059 4.9078085 1.2888347
#> [37,] 10.9379506 8.5744541 0.8883530
#> [38,] 0.4652805 0.8296353 0.3147421
#> [39,] 0.3947759 0.3343116 0.3178591
#> [40,] 0.3928706 0.5818001 0.3184558
#> [41,] 1.0826516 4.3914798 2.9067917
#> [42,] 1.2107099 0.7773342 0.8352794
#> [43,] 0.9732368 0.5738284 0.6662570
#> [44,] 2.2176350 0.6195366 0.9186807
#> [45,] 0.3819071 0.2975143 0.9137178
#> [46,] 9.1440676 14.3515874 2.7175076
#> [47,] 12.5697306 8.1269787 4.2305131
#> [48,] 1.0568433 1.6460488 0.7995212
#> [49,] 0.5382643 1.2040129 1.4423600
#> [50,] 5.7909121 3.7751498 3.9181270
#> [51,] 6.6164372 2.3820202 4.2380598
#> [52,] 0.4993421 0.2467188 0.3669028
#> [53,] 0.6208297 0.8145474 83.1811189
#> [54,] 1.4175540 0.6253529 1.6085974
#> [55,] 0.7355190 2.6328965 0.6516679
#> [56,] 192.9421926 3.2329951 68.8271512
#> [57,] 1.9927580 1.1923993 0.5699933
#> [58,] 1.9557912 1.8414261 0.6659310
#> [59,] 4.4314033 4.7621051 0.7821811
#> [60,] 0.9340977 0.6243714 2.3328791
#> [61,] 4.6175969 4.9114225 2.8341926
#> [62,] 0.5922554 0.6155032 0.2303388
#> [63,] 1.5850214 0.8623912 1.2573554
#> [64,] 0.9588624 0.4067635 0.2899709
#> [65,] 1.2885044 0.5575733 1.6546550
#> [66,] 0.6909401 2.7401999 0.7745740
#> [67,] 1.5438068 1.4281182 0.9236493
#> [68,] 2.9168079 1.1019767 0.9019602
#> [69,] 0.8865898 8.0058187 2.0059042
#> [70,] 30.4420835 2.3326463 6.6015322
#> [71,] 1.9341816 0.8653862 2.7544850
#> [72,] 1.1807932 0.7538528 1.5536932
#> [73,] 33.6632400 85.8494061 58.1541682
#> [74,] 0.3245085 0.5107813 0.5167024
#> [75,] 5.2374610 2.9525616 3.2391658
#> [76,] 6.1480735 1.4185436 0.9164292
#> [77,] 2.2149889 8.2491304 1.4824927
#> [78,] 1.5097219 22.0040100 4.0067921
#> [79,] 1.4278904 3.0039171 2.1123461
#> [80,] 1.7909433 0.7928260 3.4485898
#> [81,] 0.7437094 0.8385937 0.8332425
#> [82,] 4.9235325 1.9314674 2.2036865
#> [83,] 0.8351811 0.7264963 0.3674909
#> [84,] 2.7487421 1.1917754 1.5711182
#> [85,] 9.0056659 49.6340675 1.5135840
#> [86,] 3.8024370 3.8457471 2.1923028
#> [87,] 0.5237128 0.2880730 1.2051598
#> [88,] 1.4188974 1.3842895 3.3272805
#> [89,] 7.6188992 7.8856157 13.2300381
#> [90,] 0.9718113 1.4439526 1.8759952
#> [91,] 0.5947192 2.1617586 1.0385477
#> [92,] 0.5783418 1.3617116 0.6488635
#> [93,] 1.0534374 1.0104602 2.7699398
#> [94,] 5.9457199 3.9297503 6.4678743
#> [95,] 0.7635072 3.0211491 4.1473391
#> [96,] 31.4177375 18.2079488 2.1950209
#> [97,] 0.9431628 4.8408951 8.6891505
#> [98,] 2.1741161 0.9885503 3.4733264
#> [99,] 0.6741489 0.5926507 2.8430918
#> [100,] 41.2323814 40.7045096 17.0570588
rmev(n=10, param=c(0.1,1,2,3), d=3, model='pairbeta')
#> [,1] [,2] [,3]
#> [1,] 0.3227291 1.8610507 0.6340582
#> [2,] 3.2009836 3.1587561 5.2080782
#> [3,] 0.8553240 0.5748510 1.6925412
#> [4,] 0.6101622 0.3649300 0.8713223
#> [5,] 1.4052551 0.7505828 1.1049207
#> [6,] 1.5066514 3.6511266 0.3876806
#> [7,] 6.6746807 1.7482743 1.0857901
#> [8,] 1.7978592 0.9683812 2.1479411
#> [9,] 1.2401136 0.5715295 0.3390979
#> [10,] 5.6892034 7.0140057 2.0015396