Probability density function of multiple structured Normal inverse Wishart
Source:R/mmsNiWlogpdf.R
mmsNiWlogpdf.RdProbability density function of structured Normal inverse Wishart (sNiW) for multiple inputs, on the log scale.
Arguments
- U_xi
a list of length n of observed mean vectors, each of dimension p
- U_psi
a list of length n of observed skew vectors of dimension p
- U_Sigma
a list of length n of observed covariance matrices, each of dimension p x p
- U_xi0
a list of length K of mean vector parameters for sNiW, each of dimension p
- U_psi0
a list of length K of mean vector parameters for sNiW, each of dimension p
- U_B0
a list of length K of structuring matrix parameters for sNiW, each of dimension 2 x 2
- U_Sigma0
a list of length K of covariance matrix parameters for sNiW, each of dimension p x p
- U_df0
a list of length K of degrees of freedom parameters for sNiW, each of dimension p x p
Examples
hyperG0 <- list()
hyperG0$b_xi <- c(-1.6983129, -0.4819131)
hyperG0$b_psi <- c(-0.0641866, -0.7606068)
hyperG0$kappa <- 0.001
hyperG0$D_xi <- 16.951313
hyperG0$D_psi <- 1.255192
hyperG0$nu <- 27.67656
hyperG0$lambda <- matrix(c(2.3397761, -0.3975259,-0.3975259, 1.9601773), ncol=2)
xi_list <- list()
psi_list <- list()
S_list <- list()
for(k in 1:1000){
NNiW <- rNNiW(hyperG0, diagVar=FALSE)
xi_list[[k]] <- NNiW[["xi"]]
psi_list[[k]] <- NNiW[["psi"]]
S_list[[k]] <- NNiW[["S"]]
}
mmsNiWlogpdf(U_xi=xi_list, U_psi=psi_list, U_Sigma=S_list,
U_xi0=list(hyperG0$b_xi), U_psi0=list(hyperG0$b_psi) ,
U_B0=list(diag(c(hyperG0$D_xi, hyperG0$D_psi))) ,
U_Sigma0=list(hyperG0$lambda), U_df0=list(hyperG0$nu))
#> [1] 6.24819325 7.09519233 2.38417823 7.98159771 9.24676782 8.57167819
#> [7] 9.02768584 8.27644594 7.23655526 7.36875081 2.22609620 2.27975006
#> [13] 3.54675072 5.87840349 5.16517348 2.16201055 6.34551187 5.06266144
#> [19] 6.89241595 4.24524902 2.34118165 3.67413378 2.84431127 6.70480856
#> [25] 1.50598683 7.24727396 7.49534584 5.97688757 6.75287492 6.26316248
#> [31] 7.98136017 6.72307861 7.11605033 -0.00844199 7.68952531 5.73086864
#> [37] 0.85760353 6.10315856 6.10475908 6.67449629 7.24543290 5.07016017
#> [43] 5.38342728 7.05341297 1.07731198 8.33117029 4.37472802 9.03847661
#> [49] 5.99911735 5.45105849 5.13788473 5.12899831 8.51482988 7.68324383
#> [55] 2.81392397 2.02725765 8.08984830 1.50476358 5.27034323 3.08779191
#> [61] -1.82980915 8.47602105 5.27429001 6.39281793 4.59746682 5.75316610
#> [67] 6.85962726 8.06167108 7.17607479 1.71838795 8.92930149 3.58182593
#> [73] 4.50215165 5.90762563 -0.59626800 4.00243306 5.94378882 6.87240830
#> [79] 6.92918810 2.37685468 7.54108106 1.40421693 7.16925197 6.45898903
#> [85] 7.62033478 2.35426464 6.20525148 3.91896051 7.52693741 8.15056870
#> [91] 7.54459704 3.12051126 6.20168888 4.88594128 0.07861976 6.41736329
#> [97] 3.04076652 4.15663054 2.24175495 5.96049055 5.91665166 8.13821483
#> [103] 7.15328418 8.31812283 6.55485838 6.88063121 6.31280459 5.76908227
#> [109] 8.09058348 4.61791495 4.91444454 8.12540879 9.41996194 3.34006293
#> [115] 2.68170230 7.67499482 8.03093925 6.54701669 8.06097064 7.92488168
#> [121] 8.06615805 6.77279083 8.46686316 8.81466189 7.56202036 1.05097703
#> [127] 6.14400272 6.80986389 4.82721934 7.47687569 6.86548211 5.33329385
#> [133] 7.26285298 1.12391491 5.08956516 7.92288112 -1.23459959 7.13352253
#> [139] 7.84743699 7.30627228 7.26176106 7.78325930 5.51997980 7.92917119
#> [145] 7.66249311 7.87439072 6.20939656 8.06686379 8.34890806 3.64634897
#> [151] 6.65339136 5.56928208 6.75687746 6.95185896 3.60055600 6.82061034
#> [157] 3.09066343 5.53843388 3.09899250 9.01952174 6.37750972 2.53630115
#> [163] 8.48605689 7.53380109 5.55679969 -1.06160945 8.62028322 7.81272128
#> [169] 4.65755821 4.06877417 8.00409349 3.86621478 5.19421347 5.95330529
#> [175] 1.98569370 9.14719031 7.76615337 8.78959740 7.59757386 2.17886054
#> [181] 8.87838571 8.62195289 6.72207279 7.78284208 7.56084922 3.71220678
#> [187] 5.87629888 7.61753536 6.42742214 9.36973440 7.05588242 4.59545325
#> [193] 4.02099059 7.82572125 2.97804159 1.31020870 3.48378368 5.37788323
#> [199] 1.57047610 5.65114587 1.59248726 5.22197652 6.68947122 2.81457892
#> [205] 6.25449736 6.58646181 2.75834163 8.67122520 7.52268635 5.79047234
#> [211] 2.89492950 2.55480107 8.84555014 7.00315831 5.30510691 3.99115224
#> [217] 3.44950476 8.01807879 2.93328071 2.74344038 3.50859696 7.75490120
#> [223] 4.78099583 3.60373866 6.63882546 4.26443080 8.51246437 -0.55530190
#> [229] 3.50833057 5.24763548 -1.88760287 5.72298460 9.00563784 5.89580215
#> [235] 2.36533330 4.93313929 8.41894729 7.85246798 4.89032109 5.65867688
#> [241] 6.66233098 3.51986221 8.29364511 5.03178848 6.66679230 6.12904567
#> [247] 8.24700539 5.85968440 6.48308040 2.42597926 6.22783673 6.16969537
#> [253] 2.81217891 7.85669593 3.99677272 5.82053254 8.08265244 8.03427826
#> [259] 7.73166708 6.87763265 6.87390466 3.21647591 5.28350866 0.89109271
#> [265] 8.75357821 7.76625353 7.12859197 5.91047621 5.15806218 2.87970690
#> [271] 5.58713092 4.19285305 7.41786987 2.70211527 5.87057871 6.25821215
#> [277] 7.10221477 4.70697703 5.65234655 4.17575410 5.98453297 8.49243525
#> [283] 6.86633674 6.64756074 5.72863066 9.02761154 7.10405911 3.40068587
#> [289] 6.35655646 8.49488780 9.43733902 6.82462600 8.41328053 4.25018615
#> [295] 6.73765472 7.61088305 4.51592699 0.34226172 0.65371226 2.71093998
#> [301] 4.27472572 7.95789523 5.87515761 5.72424573 7.42145790 6.79362394
#> [307] 9.70225920 6.27186513 8.46459071 6.38312993 7.13483575 1.41726032
#> [313] 7.78896018 7.59972680 2.85028667 6.48142660 9.19454327 7.33152120
#> [319] 6.83719795 7.66601777 0.39162955 4.09831717 3.20496781 2.15709717
#> [325] 5.31549040 5.01048622 7.43182539 5.44325554 5.06865323 6.71806799
#> [331] 1.84656287 5.39791557 7.90487588 5.50681152 8.86450320 5.08927384
#> [337] 5.88599153 7.24662990 5.79237457 7.43498474 2.78487977 7.01872315
#> [343] 8.48914971 7.60172092 5.59978698 6.87943410 7.00040068 5.80740260
#> [349] 8.87498781 8.80893651 3.97049621 4.70875803 5.63037564 8.13174373
#> [355] 3.65170262 6.72010782 3.26838574 8.15287879 4.49603600 5.70917352
#> [361] 1.44471424 1.95725971 7.07086265 2.11007877 5.84713026 5.68520377
#> [367] 6.11993699 5.82131097 5.79725846 8.01364983 5.15912127 4.51000517
#> [373] 5.08941895 8.02590621 8.47599795 2.10637351 4.89348847 6.15414412
#> [379] 9.20929880 6.91892666 6.18604479 7.26262929 6.42156139 7.48604178
#> [385] 4.82805116 7.54573637 -3.74316807 4.87980850 0.52042873 4.66867663
#> [391] 5.37272283 7.17035032 7.06389272 5.67230539 6.75219337 5.54446260
#> [397] 5.27256138 6.34457303 7.72885001 3.93143044 8.03413544 6.42381160
#> [403] 2.89132194 3.43225163 7.19899104 6.60733729 6.51673876 6.78744451
#> [409] 0.06622454 3.32707630 3.55740561 6.87359228 7.93628190 7.73750081
#> [415] 3.23363078 6.91884323 7.51099256 7.76592396 0.44834590 6.45286247
#> [421] 6.89395775 6.05149369 6.37595471 6.91434561 -0.78169038 7.95926401
#> [427] 7.42852094 7.48956735 8.11138458 8.51207507 8.06930625 4.81665137
#> [433] 4.83869758 7.50121776 4.25309980 5.52872407 4.40373042 4.52378442
#> [439] 5.54085078 4.64546728 3.97058522 8.16905655 7.97506606 6.10716360
#> [445] 5.67132307 6.69233033 6.53618728 7.26349846 4.39334871 6.15138719
#> [451] 5.93669871 5.59931854 3.86417250 7.04086403 6.87813875 1.05048041
#> [457] 5.11803324 7.71887547 6.90723556 5.62853357 5.11036356 4.66463320
#> [463] 9.20862421 6.87689123 5.39104732 6.73353008 -0.90157087 8.49202162
#> [469] 6.15209782 6.90485521 4.50055829 6.40721671 4.92887247 6.17557475
#> [475] 3.93425281 6.31014295 3.65122920 8.29283629 9.32949914 7.76080306
#> [481] 5.33260239 7.99856000 4.79599502 9.29206633 5.25856276 4.03227477
#> [487] 8.87359496 5.10124845 2.48454181 5.17802792 9.24126735 7.21091633
#> [493] 4.09636698 1.23579623 7.94777646 6.00231927 2.91882195 7.87292833
#> [499] 7.83715064 7.42068182 6.82783861 5.35307259 5.20621960 7.19665769
#> [505] 3.81485377 6.51549058 5.71559849 7.39253570 3.65297621 7.05038526
#> [511] 6.41002410 5.90797258 5.96131821 6.94508894 3.29546762 4.21132832
#> [517] 8.62853127 2.28365965 5.61041783 6.92673424 7.36331947 3.55452720
#> [523] 3.49717895 4.51009560 7.69153690 0.74212992 7.21377832 6.24481773
#> [529] 7.94816478 3.06444459 7.71398961 5.12545345 6.32953389 7.35457888
#> [535] 4.98704798 8.19103972 4.51110572 3.10943122 8.44241426 8.45149017
#> [541] 6.60286637 5.82648248 7.88006545 5.37352362 4.60079213 5.34657601
#> [547] 6.22174464 7.96631437 4.45251408 4.97354257 1.66611516 8.28643241
#> [553] -0.47612168 2.84598112 5.09863200 2.71964335 5.62854922 7.12243701
#> [559] 5.98976313 5.33721563 4.49901701 4.27063529 8.86443750 6.85352331
#> [565] 5.27060709 -1.71957636 5.88841526 4.24689251 8.05575372 5.39740480
#> [571] 6.84635769 6.52009759 6.93103349 2.93931989 3.96168988 7.07060297
#> [577] 5.37730319 4.28255335 9.29624535 4.73831748 0.53464292 4.32256465
#> [583] 7.00100554 6.88830161 5.37253773 6.17072702 -0.40705320 6.56198525
#> [589] 8.30311655 5.84202021 3.72098250 5.99144010 5.29906852 8.49332484
#> [595] 5.13508566 6.31861778 5.83653213 6.44967387 7.54281643 6.21092533
#> [601] 6.54663559 6.89140545 8.21130761 5.36814764 4.75443510 7.11033243
#> [607] 4.15229681 3.31112118 7.62053480 5.22505715 6.23381904 5.45991600
#> [613] 4.54221124 6.79242821 6.32966135 2.30139173 6.42790571 7.00268204
#> [619] 5.57092622 6.27406567 -2.79294980 7.22681804 8.04296846 4.24622081
#> [625] 7.66038669 7.74947104 7.43923898 6.00096319 5.37950133 5.25778183
#> [631] 6.02288138 6.36185001 9.11578645 7.89654236 6.04608328 5.70394820
#> [637] 6.82520770 8.35839201 3.39176529 7.79940129 7.37526500 2.73915976
#> [643] 7.15185236 5.38741431 5.44398507 4.23193384 4.36829068 1.00218307
#> [649] 7.84629687 6.75734763 8.41283102 6.36151717 6.21304500 4.90464858
#> [655] 8.85375224 1.99733216 3.80060866 7.64259683 8.20447792 3.84052260
#> [661] 8.99379274 6.17175328 6.38630283 5.84436326 5.07816026 6.08821535
#> [667] 7.69172169 3.33719360 4.73313574 7.72239544 8.03988784 8.35464035
#> [673] 6.28017920 7.10283666 7.45877351 4.64773962 6.57630426 3.77718274
#> [679] 7.05293260 8.20505097 8.18450930 5.05152546 0.17806917 4.17455394
#> [685] 7.87413911 3.75257819 5.18880737 5.47610747 4.17597293 1.59543124
#> [691] 6.93487589 5.88644381 3.91511463 5.28653980 6.46350334 8.29584229
#> [697] 2.63143553 3.66571388 7.28318483 5.19807480 8.62731060 6.34383556
#> [703] 2.82271701 8.41122922 6.75785971 4.26069223 6.12615408 6.41823825
#> [709] 2.88701885 3.89332098 5.42321982 5.98966542 8.14571852 3.44358760
#> [715] 7.25462375 8.74510386 4.52938530 5.38340574 5.16968160 3.35887591
#> [721] 4.90262238 7.40111057 6.61651282 4.47096965 8.10717794 7.34766815
#> [727] 7.57491522 7.62857858 5.81718406 4.27882815 6.37293271 4.87392947
#> [733] 6.66803299 3.32905362 6.82411343 4.51183259 6.57823016 7.31026833
#> [739] 7.13062148 6.69408170 4.38423629 6.81623925 1.05549484 4.60174277
#> [745] 4.52329228 7.13214554 5.39262385 7.44292156 3.78550878 7.46260507
#> [751] 2.72413996 -0.27464443 7.41307671 7.21765814 3.08709055 7.82015186
#> [757] 0.14038087 7.02893662 2.32107741 8.39344295 1.88674737 6.32405521
#> [763] 5.61897290 5.77390747 8.85992241 7.05372198 9.13421888 6.48162886
#> [769] 3.26692017 7.62934623 -2.78162166 4.84539209 5.36970697 -1.11006030
#> [775] 8.48055877 2.29309526 4.53512950 3.93366586 3.41591438 2.20267110
#> [781] 3.31201142 6.14017692 8.71228424 7.32656636 9.25037586 7.97694773
#> [787] -0.78186476 7.08658693 7.68341035 2.57441859 6.49033037 7.26273686
#> [793] 4.72072778 4.23604614 4.18841461 7.81901528 8.72174533 5.55546377
#> [799] 7.21029405 3.74971486 4.77881733 5.74315680 1.59000063 5.87755938
#> [805] 8.85652022 5.87134359 1.02573129 7.24212440 6.83455025 6.82654309
#> [811] 6.49815669 5.28974868 5.27438081 7.05277507 9.04230338 3.69039769
#> [817] 4.41435314 7.73357470 5.83418052 7.59554602 6.23713694 4.89942312
#> [823] 8.02471080 5.22807514 5.36493471 6.46287389 1.59018300 7.85493513
#> [829] 4.40431228 7.37937234 3.17959055 6.79171102 7.94752796 3.26217539
#> [835] 2.24181227 5.60168367 3.91732075 8.36005923 3.73520406 2.27518841
#> [841] 0.82757478 8.23489748 2.94360781 8.37146995 8.15435559 4.53881151
#> [847] 4.04365503 3.55114037 5.98116547 8.48933641 7.40196623 7.22117521
#> [853] 3.92864440 4.66662779 6.66005042 3.49377716 6.62243842 8.35999330
#> [859] 7.52375094 6.79518472 4.10581920 0.71086767 7.91663941 7.58785800
#> [865] 4.75497668 6.71051748 2.57397311 1.57216368 5.29290239 3.67652528
#> [871] -3.89538086 5.16395501 6.66257378 8.21053358 0.81990460 5.72453452
#> [877] 7.06398663 7.35055322 3.41852426 6.12839340 4.64308568 7.12513810
#> [883] 3.99876565 4.84665202 4.33637765 5.85158759 7.62954757 7.36209268
#> [889] 5.51254530 8.40069471 6.13236212 5.65985848 1.53539306 7.75252392
#> [895] 6.84918408 8.13869021 8.21332484 7.81891148 4.75385160 8.27796845
#> [901] 5.50020049 4.02171287 1.53384754 6.76867521 4.98831769 5.95919435
#> [907] 8.11283490 0.61545753 7.40278079 8.36613679 6.92625133 7.42561564
#> [913] 7.83543677 2.72475359 5.48903611 5.50254850 6.10971362 8.63953761
#> [919] 8.26466598 6.59791763 4.26246681 5.68337564 4.85170261 6.73689034
#> [925] 8.61804250 5.21720414 8.02952071 8.05693190 8.81396888 5.28094763
#> [931] 4.96949264 9.52720450 2.65927756 4.76236171 5.89647748 8.18494264
#> [937] 2.67998557 6.85111083 7.43683125 3.28125772 8.03232923 7.40850638
#> [943] 6.08364937 5.19064068 0.55622763 7.33720965 2.63109571 0.81523273
#> [949] 5.99735431 1.98223768 7.22103118 5.72436927 3.88764269 8.42236044
#> [955] 7.61258015 3.04020340 5.16476512 7.09860755 6.80837786 5.95471620
#> [961] 6.35393881 3.43527186 6.98171918 3.65805725 7.69641849 4.42598310
#> [967] 6.09361730 6.73641620 8.25896428 0.27618991 6.17303879 5.16337952
#> [973] 3.17318742 -5.12513450 5.59099344 1.22595434 8.93741145 2.54383798
#> [979] 4.78009125 8.37058231 5.96201998 5.98932175 6.00762637 6.34722210
#> [985] 8.33304948 8.42182945 3.90810103 8.55038336 4.64045379 4.59793564
#> [991] 7.62815911 5.52066757 5.49683910 5.54847696 5.95486210 7.37842650
#> [997] 4.92517655 6.66911096 8.32985962 2.79611352
hyperG0 <- list()
hyperG0$b_xi <- c(-1.6983129)
hyperG0$b_psi <- c(-0.0641866)
hyperG0$kappa <- 0.001
hyperG0$D_xi <- 16.951313
hyperG0$D_psi <- 1.255192
hyperG0$nu <- 27.67656
hyperG0$lambda <- matrix(c(2.3397761), ncol=1)
#'xi_list <- list()
psi_list <- list()
S_list <- list()
for(k in 1:1000){
NNiW <- rNNiW(hyperG0, diagVar=FALSE)
xi_list[[k]] <- NNiW[["xi"]]
psi_list[[k]] <- NNiW[["psi"]]
S_list[[k]] <- NNiW[["S"]]
}
mmsNiWlogpdf(U_xi=xi_list, U_psi=psi_list, U_Sigma=S_list,
U_xi0=list(hyperG0$b_xi), U_psi0=list(hyperG0$b_psi) ,
U_B0=list(diag(c(hyperG0$D_xi, hyperG0$D_psi))) ,
U_Sigma0=list(hyperG0$lambda), U_df0=list(hyperG0$nu))
#> [1] 2.906303875 2.671372120 1.325353660 1.057552200 1.049203088
#> [6] 2.935065612 1.040110921 -2.057000714 2.730023305 2.479453006
#> [11] 2.061061855 0.849743560 1.296964637 1.164313348 0.889914710
#> [16] 1.621657255 2.838385568 -0.714630166 1.179854769 2.538936571
#> [21] 2.120324857 2.000557936 1.355440499 1.926799302 1.958887585
#> [26] -1.141021178 2.841875595 1.674703897 1.253710328 2.015727825
#> [31] 2.834835080 2.518846052 1.743120856 2.538204395 -1.567358399
#> [36] 2.734103676 0.955616995 2.059975858 0.293845923 2.476201494
#> [41] 0.494220822 1.953708089 0.387236657 1.670179538 0.968235794
#> [46] 2.736937935 2.424902703 1.836295571 2.713275202 2.373466041
#> [51] 1.454679140 2.448166406 1.284789420 1.950122238 1.671794267
#> [56] 2.174297962 -0.978210318 2.340425949 -2.225397556 0.819309873
#> [61] 2.598870451 2.848271959 2.105866190 1.214921785 2.215728480
#> [66] 2.563401378 0.618143408 2.046114810 1.437890208 0.977646659
#> [71] 2.759866381 1.561814390 2.883125770 1.907314959 1.990645426
#> [76] 2.464989017 0.847170415 -0.076004650 2.612693674 1.536689268
#> [81] 2.894619803 1.092910254 -2.133967275 -2.244729497 2.085185164
#> [86] 2.841725483 2.799454354 0.920213210 2.633582613 1.656639838
#> [91] -0.036541812 -4.619866592 2.539746205 1.356097880 2.826578008
#> [96] 1.591853272 0.238789590 2.120467231 2.263254294 -0.880352237
#> [101] -0.819075811 2.466356151 2.405308344 1.759222657 1.903382393
#> [106] -0.420686801 -3.629912465 2.873066216 0.508039477 2.345859567
#> [111] 1.331695238 -2.621865041 -1.118953380 1.636120781 1.539335716
#> [116] 2.210144255 2.863053805 1.951278966 2.681360645 0.390965031
#> [121] 2.004792905 2.427252638 2.699661873 0.553152012 1.277232217
#> [126] 2.140451538 1.154965965 1.509696945 1.312776508 2.095594262
#> [131] 2.732206899 -0.240482646 1.312253480 0.994373868 2.744285675
#> [136] -1.345945658 2.475882342 1.200618431 2.814946815 2.212781506
#> [141] 0.483079259 2.251641279 2.226026911 2.367347626 2.654407699
#> [146] 1.982407771 2.702290063 1.299005178 2.349803185 1.988088178
#> [151] 1.887150555 0.367404788 2.338153077 0.621188580 -0.563457848
#> [156] 2.319097776 1.980876933 1.847542370 -0.175046438 2.716551478
#> [161] 2.053164466 -1.007623963 2.314757375 2.012733821 2.932420766
#> [166] 2.274674875 2.757788431 2.737979407 2.751100060 2.372315217
#> [171] 1.692607279 2.745649259 -4.015948493 2.279497797 1.832308556
#> [176] 1.107666859 2.001912781 0.634418666 1.578439435 -0.297697098
#> [181] 0.857011221 1.904876526 -1.624446482 1.490688898 1.863392530
#> [186] -1.354491498 1.216180106 -0.081244501 2.291310764 -0.133402697
#> [191] -0.197357515 1.532636877 2.660503920 2.766708662 2.573694932
#> [196] 2.951082970 2.953695701 2.573527008 0.778053145 2.192580758
#> [201] 1.651387474 0.148026793 0.055713608 2.299425624 1.648875932
#> [206] 1.671885699 1.640469068 2.238133363 1.610719069 0.482877304
#> [211] 0.919146968 2.344799477 1.805486814 1.712484102 0.276091209
#> [216] 0.421096313 2.406139064 2.361656125 1.465082483 2.368897744
#> [221] 1.680341227 1.464829970 -0.868506511 0.976896094 2.716863110
#> [226] 2.164490568 -2.273071709 -2.869760047 2.334985225 2.563360796
#> [231] 1.623631471 0.958448707 2.823719970 -0.695726489 0.537843859
#> [236] 2.031910381 2.433013223 1.288192049 0.568997754 2.646631712
#> [241] 2.882265236 2.736676048 0.759766891 0.822441748 2.581879979
#> [246] 1.156301193 2.958338605 1.846756670 0.259362633 0.517738710
#> [251] 1.125793077 2.261962556 2.016294761 1.375212804 -2.799953418
#> [256] 2.955307383 2.127756791 2.198923546 0.715620033 0.193571387
#> [261] 2.372041071 -0.379097550 -2.296526084 0.055996184 2.569123119
#> [266] 0.923021875 1.624389326 2.327999705 1.106662241 0.810848645
#> [271] 1.144685080 0.593623563 2.128120752 2.448801187 -4.595156540
#> [276] 3.000290202 0.058090905 2.843525948 0.799037364 1.934248001
#> [281] 0.679751456 2.103883576 2.111176353 1.925236724 0.907464156
#> [286] 1.896188475 2.090633415 2.310932643 2.660877649 1.109212809
#> [291] 2.461258901 0.500037376 1.820242586 0.051151229 1.790316620
#> [296] 0.119284700 2.375561988 -0.213729339 2.039203543 -2.967246511
#> [301] 1.459262465 -3.412822574 -0.421434149 0.332701152 2.643320933
#> [306] -0.787517334 1.529878924 2.596224144 -0.320517125 -2.933991891
#> [311] 0.521072366 1.157747131 1.336262431 1.096386060 2.229050948
#> [316] -3.783778954 1.859996732 1.896434025 -3.222866629 1.328196939
#> [321] 2.318876398 -0.347849059 2.830698311 1.580895382 -0.601392347
#> [326] -0.590399228 0.381109125 2.760678161 2.314842012 2.266137096
#> [331] 2.168528668 1.893126682 0.800948320 2.222284767 1.497053926
#> [336] 0.374091548 0.900546870 1.746269279 2.083827131 1.325099891
#> [341] 2.335969390 -0.006902181 1.000770997 -0.211341119 2.276059679
#> [346] -1.988568069 -2.433190435 1.503315007 -0.727938546 1.959191393
#> [351] 2.162121777 2.066802395 1.297229872 2.742794897 0.061099245
#> [356] 2.463724925 2.394546437 1.523302948 2.711215988 0.725807833
#> [361] 2.332935268 0.680154232 0.903773832 0.530608252 2.228806540
#> [366] 1.975941620 -4.097564867 1.213792163 0.812672134 2.022187864
#> [371] -2.533557887 1.126487945 1.764136698 1.307980204 2.525177167
#> [376] 0.014695297 1.007272172 -1.210730099 2.858342109 2.736185416
#> [381] 2.419264160 0.947132578 1.963150779 0.975255540 2.414165159
#> [386] 2.623632445 2.287022755 2.399876093 -0.059527140 1.923881029
#> [391] 0.878526230 1.207812080 1.961503347 2.727138678 -2.300345508
#> [396] 2.696176359 2.217046227 2.752208191 -1.054455869 0.789725253
#> [401] 2.871668558 1.730576495 -0.415284511 2.638859088 1.467645261
#> [406] 2.462198257 0.735679783 0.083637841 2.616223174 2.933983660
#> [411] 0.672827349 -1.217349697 0.405706366 2.504458287 -0.568433687
#> [416] 1.226040040 2.264251530 0.805872762 1.788135111 1.331611804
#> [421] 1.979352244 1.396860435 2.238256873 -0.557588989 2.470281325
#> [426] 2.932274346 -2.843481093 2.424195808 0.292642877 -0.667431733
#> [431] 2.830371010 -1.165087887 1.049402466 1.882470244 2.668816609
#> [436] 2.057281547 0.743801326 1.946680603 2.185067226 2.848530392
#> [441] 2.453631800 1.373836368 2.920677642 2.589064788 0.609689346
#> [446] -0.514841189 1.383460852 2.592997476 -0.580040770 2.685925313
#> [451] 0.150572684 2.759629729 -0.247948553 1.268093470 -0.654485614
#> [456] 1.759844941 -2.887629809 0.760475511 2.497767470 2.084477848
#> [461] 2.149954646 2.971759715 0.776251930 -1.694051228 -0.306913882
#> [466] 2.112335933 1.728326044 0.430239049 2.442007245 0.589402517
#> [471] 2.437210329 2.472631403 -0.240628357 -0.191945190 -1.860595114
#> [476] 0.095844670 0.306170118 0.502435720 1.826728663 2.030446116
#> [481] -3.780485354 1.282188805 1.896366144 2.751644902 0.944955948
#> [486] 1.943059630 1.848023271 2.346522975 2.522611407 0.499024114
#> [491] 1.985199086 2.850741710 1.412459680 2.780036753 0.956871460
#> [496] 2.621766641 2.704375356 1.357257050 1.359500769 2.023871981
#> [501] 0.769490206 2.636104012 2.554653425 2.228339009 1.123870636
#> [506] 1.782086066 0.215977867 -1.978038719 2.248979427 0.981700833
#> [511] 1.021497991 1.521566688 1.037458133 2.263308616 2.497896342
#> [516] 2.676866191 1.593063187 0.732983895 -1.904284098 2.440279661
#> [521] -1.388629949 1.892941129 2.380343254 2.266268108 1.265732686
#> [526] 2.030319660 2.213894619 2.703599093 1.083233317 -1.636567413
#> [531] 2.371373774 2.436889067 0.330478503 0.470699159 1.722981881
#> [536] 2.019827114 2.877680138 -0.763393655 0.253659383 -0.918610859
#> [541] -1.678670999 -2.700557441 1.384019615 1.888676289 2.729490245
#> [546] 2.310715358 0.351760716 0.245889749 0.095219677 2.197463093
#> [551] 0.360181289 2.039063639 2.720819787 0.935481887 0.569023345
#> [556] 2.810406465 2.109911836 2.129800942 2.597316740 2.856240770
#> [561] 2.961826212 2.328233825 -3.896307161 2.575655516 -0.008986700
#> [566] 2.206780167 1.470782634 1.690665844 2.893178993 1.748477613
#> [571] 1.283140833 2.604326073 1.766146195 0.393727783 2.522026830
#> [576] -0.858398652 2.686996831 0.132161107 0.397116879 1.781390818
#> [581] 1.647457842 2.331165542 2.866280281 1.819443397 2.248880399
#> [586] 1.543876390 2.301726988 0.012496766 0.654055373 1.865325579
#> [591] 1.566691058 1.483072323 1.508649803 -0.321383471 2.494003381
#> [596] 0.076502744 2.820809043 0.929689248 -2.256556820 2.965678032
#> [601] 2.437709675 1.293681618 0.218010992 1.425140018 1.460894655
#> [606] 1.649604379 -0.344593524 1.263839928 1.648113377 2.661647111
#> [611] 2.212810897 1.685708123 1.573532271 2.339955779 2.208440413
#> [616] 0.292194164 2.380808523 1.340682077 2.937516863 1.827329221
#> [621] 2.243284321 2.530250875 1.833566773 2.603273183 2.164894986
#> [626] 2.155853229 1.712356908 -0.078882527 1.007417480 2.201746338
#> [631] 1.623164879 1.487589270 1.412748849 2.639379494 1.980431500
#> [636] 2.733961375 1.940356032 -3.185845807 2.587175670 2.524467162
#> [641] 2.893095255 1.830619584 1.388773936 -0.463940728 2.515680944
#> [646] 0.181699103 -2.520188353 2.762628316 0.333987524 -1.728584339
#> [651] 0.834914981 2.610907875 2.418174306 2.781792217 2.365646901
#> [656] 0.651484574 2.471175025 0.885832639 2.833014570 2.690539621
#> [661] 1.249904095 1.032035685 -1.214794651 0.114471726 2.077043249
#> [666] 0.560288945 1.310910032 0.316818750 2.925813452 1.692150665
#> [671] 2.317544658 1.106411023 1.947087653 2.970608710 0.838316123
#> [676] 0.838552962 0.110511737 1.500907965 1.534947115 2.210081472
#> [681] 2.015055908 2.002308274 2.448722234 2.123979879 2.305616902
#> [686] 2.285163832 2.252250394 1.201324002 0.861825456 1.530228400
#> [691] 0.795434297 1.205509767 1.219559200 0.015175546 1.895455509
#> [696] 2.296679238 2.694528180 -0.613159218 0.921685038 1.845192673
#> [701] 2.134551508 0.988706167 1.835939757 1.356145940 1.775032842
#> [706] 2.211494751 2.249806066 2.299161811 1.298646335 2.678636407
#> [711] 2.009354048 1.657994621 0.288941029 0.662762449 0.831851834
#> [716] -2.497058535 1.298631676 2.843743217 1.499942680 2.756346283
#> [721] 2.604813206 1.821271123 1.795050401 2.624153101 1.636759418
#> [726] 0.158734197 2.757074119 0.891161569 -0.781016915 0.321690089
#> [731] 2.618101601 2.833417113 -3.758033164 -0.309997348 1.997649578
#> [736] -1.936939755 1.894764984 2.630714560 2.057073801 1.133311976
#> [741] 2.145321626 1.977599928 -0.629999632 2.531947584 2.482346051
#> [746] 1.977140249 0.366584897 1.924591160 1.208880263 2.778937182
#> [751] -0.252191329 -0.684709747 1.746030866 2.087889878 2.271750699
#> [756] 1.894644973 2.817616515 2.460626901 2.633341698 0.812319511
#> [761] 2.674133777 -1.442991203 0.969067023 2.868182076 1.222611225
#> [766] 2.085224938 2.311619088 2.952348288 2.685140725 -1.915890827
#> [771] 2.043832193 1.536625868 1.557234630 0.677816235 1.378410403
#> [776] 2.727214940 0.906891501 -0.345417627 2.913027285 2.273415816
#> [781] 1.917090347 1.041077819 1.378274734 2.719200692 2.126392976
#> [786] -1.791164535 1.582846377 2.629863563 1.318484144 2.186802800
#> [791] 2.789795030 1.943151191 -0.138956563 0.632532571 1.242901917
#> [796] 2.024535238 2.825292746 -3.983165637 0.571209721 2.383707241
#> [801] 2.423562951 -0.255578388 0.390177858 1.999304110 1.183117769
#> [806] 1.663409709 2.708230906 1.734487918 2.663693676 2.865759222
#> [811] 3.001655879 2.705857721 -1.101414260 -2.133066627 2.365829873
#> [816] -0.451889889 -0.281967452 2.597179019 1.142515827 2.356927884
#> [821] 2.562108022 -0.785378434 0.623957475 2.810794175 0.369308455
#> [826] 2.248577283 0.982159321 0.170850759 2.279862550 2.109226797
#> [831] 2.413815366 -2.755668654 -1.245523325 2.480211066 1.477074483
#> [836] 1.445494439 2.199440905 1.609155590 2.490348621 0.874756544
#> [841] 2.516024760 1.483687607 -0.137522193 2.740789357 -1.357528332
#> [846] 1.637720598 0.914897401 1.954979412 1.918016946 2.411599489
#> [851] -1.396228590 1.810661186 2.327797086 2.496293238 0.489902039
#> [856] 1.605879225 2.019544781 2.724969059 2.678046936 2.543942073
#> [861] 1.418252571 2.598694628 2.233045027 2.306492787 -1.280827070
#> [866] 2.370005163 1.887935229 0.749618763 -0.659954814 -0.287520121
#> [871] 2.709700423 0.628763225 1.440797532 2.634496567 2.896609347
#> [876] 0.696041789 2.076089363 1.725312590 2.259753749 -1.993360612
#> [881] 1.184763517 1.102545275 -5.029025670 -0.041833957 2.002593544
#> [886] 2.048273548 2.929366592 2.934960155 2.590673538 1.160287869
#> [891] 1.892855075 2.431871663 0.700515990 -0.856649814 2.623438451
#> [896] 2.752247634 0.068368197 -2.338217958 1.629621920 2.314554761
#> [901] 0.896532397 2.767141284 2.597680231 0.368074146 2.696799880
#> [906] 1.255178465 0.272788412 2.172450692 1.390352952 2.612112014
#> [911] 1.396605764 2.548910773 1.450846316 1.576092679 0.343212316
#> [916] 2.398780814 -4.494998432 1.965806229 1.359631038 0.974216014
#> [921] -1.012082935 2.715990721 1.557769924 1.619164763 2.824355334
#> [926] 1.258583337 2.993528398 1.714952816 2.081719991 -0.864684581
#> [931] 1.089331255 2.416319680 0.943398944 2.538179440 2.344896748
#> [936] 2.860280966 0.379979023 2.811519338 -0.427355887 -0.076355983
#> [941] 2.465485584 2.219989050 2.183186722 1.989535882 0.970154522
#> [946] 2.828108013 2.150477684 1.881006248 1.943442203 -0.262768339
#> [951] 2.417405406 -4.816296633 -0.643834140 2.041264446 0.858758454
#> [956] 2.700401969 -0.431047154 2.650879795 2.604148616 -0.246874148
#> [961] 1.473777806 1.937823522 -0.839535409 0.959031647 1.947023886
#> [966] 1.841923691 1.423778837 -1.423201460 0.966977450 2.762190190
#> [971] 2.019019497 2.531044088 2.500863318 2.649756960 2.384771184
#> [976] 1.811120304 1.852982830 2.014057469 1.178951539 0.974524735
#> [981] 2.077705438 2.312485168 1.042932401 2.880952313 2.001113044
#> [986] 2.666681896 2.639793164 1.813437750 1.864138382 2.119714650
#> [991] 2.019749538 1.752609440 2.097475553 1.477939643 2.654576109
#> [996] 2.088943642 2.546157394 2.492047445 2.052534080 2.326097848