Skip to contents

Simulate data from multivariate MixAR models under the assumptions of multivariate Gaussian innovarion

Usage

mixVAR_sim(model, n, init, nskip = 100, flag = FALSE)

Arguments

model

model from which to simulate, an object inheriting from class MixVAR.

n

size of simulated multivariate series.

init

initial values, a numeric matrix. If missing, a matrix of 0 values is generated.

nskip

number of burn-in values.

flag

if TRUE returns also the regimes.

Details

mixVAR_sim simulates a series of length nskip + n and returns the last n values. init provides initial values for the algorithm. Each row is considered as a time point. The number of rows must be at least equal to the maximal AR order.

Value

a numeric matrix with n rows.

Author

Davide Ravagli

See also

Examples

AR <- list()
AR[[1]] <- array(c(0.5,-0.3,-0.6,0,0,0.5,0.4,0.5,-0.3), dim = c(3,3,1))
AR[[2]] <- array(c(-0.5,0.3,0,1,0,-0.5,-0.4,-0.2, 0.5), dim = c(3,3,1))

prob <- c(0.75, 0.25)
shift <- cbind(c(0,0,0), c(0,0,0))

Sigma1 <- cbind(c(1, 0.5, -0.4), c(0.5, 2, 0.8), c(-0.4, 0.8, 4))
Sigma2 <- cbind(c(1,0.2, 0), c(0.2, 2, -0.15), c(0, -0.15, 4))
Sigma <- array(c(Sigma1, Sigma2), dim = c(3,3,2))

m <- new("MixVARGaussian", prob=prob, vcov=Sigma, arcoef=AR, shift=shift)
mixVAR_sim(m, n=500, init=matrix(rep(0,3), ncol=3), nskip=100, flag=FALSE)
#>          [,1]         [,2]        [,3]
#>  -2.336238570 -6.280878130  1.44687446
#>  -0.879145854  2.676059013 -0.95350422
#>  -0.763205473 -0.414873518  1.41222381
#>   0.014177101 -0.372360850  0.85024212
#>  -0.152367756 -0.772784953 -1.21912027
#>   0.884956942 -2.631206642 -1.34986079
#>  -0.055224648 -1.756466761 -2.40594307
#>  -0.466465877 -0.960724622 -1.46020788
#>  -0.337766196 -0.041510511  1.55175686
#>  -1.880045720 -1.597079769  0.44349826
#>   0.076058948 -1.274074367 -1.02053097
#>   0.669356450 -0.021795888 -2.80529677
#>   1.740644332  0.580027709 -2.26346814
#>   1.617856460 -3.504782954 -0.64014806
#>  -0.575999855 -1.673557076 -3.13737320
#>  -1.063179893 -3.552539023 -5.33489627
#>  -1.044257561 -0.322868730 -0.43697243
#>  -1.333417112  0.993685566  0.98372651
#>   1.785275285 -1.941277442 -0.38301796
#>   0.791626077 -0.357706498  1.26935145
#>  -3.154534077  0.886838348  0.80668615
#>   0.057509919  0.749000862 -3.49773259
#>  -0.863038192 -2.215639311 -0.68781774
#>  -2.155573305 -0.205265056  5.66359823
#>   0.845685292  2.651944034  2.98376673
#>   1.197778767 -0.313243368 -2.68262074
#>  -0.815542853  1.026386636  2.31059825
#>   1.256426920 -0.451142797  1.20265614
#>  -0.955055654  0.483118728  1.12815045
#>   0.811796344  0.906820482  0.79070622
#>   0.009165060  2.466375513  3.01926368
#>  -0.663830118  0.928558705  1.63791045
#>   1.031997126  2.062866943 -0.39578949
#>   0.527190573  0.781737000 -4.32161153
#>  -0.395691478 -1.682704171 -2.32287408
#>   0.360489197  0.136761540 -1.65636596
#>   1.171627875  2.214638065  2.17712177
#>  -0.474608272 -0.807758194  0.29173160
#>   0.204729098 -0.891704087 -4.96124586
#>  -2.427595544 -4.220385391  0.73374418
#>  -1.636554782  0.555244222 -2.23957397
#>  -2.400393791 -0.910328414  4.52863051
#>   1.436639091  2.525118626 -2.96882668
#>   0.389264371 -3.940282097 -1.56946765
#>  -1.262478295 -2.640948901 -3.64221848
#>  -0.764914082 -1.469439211 -1.11343646
#>   0.052944141 -0.205181972 -2.41070845
#>  -1.657014042 -0.131575719  4.66759093
#>   2.753387287  4.343058808 -0.53505355
#>  -0.630069071 -2.549130985  0.34028715
#>   0.296522736  0.568456865 -0.94477246
#>   2.661977124  0.990019256 -1.33239692
#>   0.535139478 -1.635800307 -2.22093002
#>   0.323922689  0.335937264 -1.33508903
#>  -0.795493625 -0.000237571 -1.78320634
#>   0.669508533  0.524798882 -1.12240254
#>  -0.110346092 -2.618220573 -1.02078602
#>  -2.473452868 -1.740444326  1.39015316
#>  -0.272757672  1.379673155 -1.66824761
#>  -1.581526943 -2.262758729  1.29301044
#>   0.626567730  1.554315214 -1.83667979
#>   0.067803959 -0.235791236  2.63100208
#>   0.403907577 -0.446668243 -1.34077561
#>   0.327360218  1.401225954  0.14053625
#>   1.602607244 -0.106810792  3.68559118
#>   3.526744213  1.320340528 -2.63898741
#>   0.599673728  1.768699931 -3.11459173
#>  -0.858480752 -2.947142457 -0.89578934
#>  -1.304059855  0.940035370 -0.13802671
#>   2.899314421 -0.517066172 -1.09619657
#>  -2.480808405  0.981459401 -0.56356308
#>  -1.439347796  1.315323405  4.04231664
#>   1.275033314  4.156893971  2.61284228
#>  -1.304382663 -0.108702863  2.02288707
#>  -1.182849888 -0.106262086  0.59763479
#>  -1.317273310 -1.599185063  1.90609651
#>  -0.678756553  1.787282109  0.46215415
#>   1.248292323  1.221922383  2.29114280
#>   0.782177440 -0.489543516 -2.75643301
#>  -0.049003317 -0.110537242 -2.36061881
#>  -0.026053112 -0.989582754  0.70650860
#>  -0.048132565  2.059031445  2.43928277
#>   2.118667744  0.458292176  1.34272810
#>   0.947289786 -1.303821922 -0.86030579
#>  -1.383259628  1.493062770  3.77081919
#>   1.629973445  1.619650551  0.38646276
#>   2.844906529  1.247482674 -3.34629851
#>   1.011647111 -2.931488396 -1.54876956
#>  -3.940264539  0.051469632  2.08071678
#>   0.284079336  2.726374894  1.43231268
#>   0.909281784  1.623804318 -1.10638808
#>   0.322801454  0.899086095  2.91479655
#>   1.587293817 -0.343584635 -2.59781487
#>  -1.277009580 -3.912753246 -0.51350647
#>  -0.352523921 -2.760703591 -0.33663618
#>  -2.714716639 -1.174470428  2.78528057
#>   0.544907207  2.259462041 -1.69024827
#>   1.828600988 -1.921878581  1.85633565
#>   1.669648130 -1.231187242 -7.05392564
#>   0.280243112  1.162845408 -5.69167824
#>   0.593814577 -1.564099190 -1.07063174
#>   1.107208509  1.129874187  0.43375940
#>   1.047476309  0.330328738  0.96935798
#>  -0.074205902 -1.209415689 -2.73755319
#>  -1.351322520 -1.033310557  1.51452369
#>   0.679136590  3.418424488  1.29968268
#>   2.411128031  2.095954454 -0.17416263
#>   1.943647128  0.015384367 -0.11454176
#>   0.184243511 -0.025955171 -0.71070318
#>  -0.474426688  1.021328878  2.33815323
#>   0.975772008 -4.182283974 -1.01980762
#>   0.744351273 -1.507419678 -1.89434838
#>   0.156485070  0.214730982 -1.78954032
#>  -0.811606833 -0.164088673  3.16302737
#>   0.965225879  1.847687087 -3.10042450
#>   3.254069677  1.225957821 -1.20408046
#>   1.161993527  0.605576659 -0.77046125
#>   0.232962692 -0.639797363  1.46571548
#>  -0.078213291  0.905387319  0.40717950
#>  -0.927132881  0.103930561 -0.54650974
#>   1.121528862  0.872152011 -2.68721477
#>   0.506478388 -1.308885462 -0.62995665
#>  -0.694197391 -2.199027932 -0.26732932
#>  -3.388952355 -0.508896049  1.88015299
#>  -0.671776240  3.168055390  1.51622531
#>   2.049999191  0.697032590  4.00304157
#>   0.967198952 -0.690081139 -3.88287063
#>  -2.348538227 -0.286917140  2.89168683
#>   1.145441579  1.652559963  0.20205138
#>   2.202462291  1.868006046  0.34943879
#>   0.688244766 -1.052384236 -0.08863272
#>  -1.591391725 -0.783587833  3.66201433
#>  -0.028799575  2.608575829  0.65510308
#>   1.034753305  1.575422895  4.95651549
#>   2.798724401  1.625651080 -1.30847713
#>   1.027600014 -4.370169404 -1.31357267
#>  -0.034851038 -1.270199306 -2.37194188
#>  -2.799687457 -1.159475265  0.49992657
#>  -2.270318421 -3.058556844 -0.50128457
#>  -1.143790562  0.208218390 -1.25702428
#>  -1.084960401 -1.286513215  2.27362783
#>   1.256771892  3.166753748  1.01624915
#>   0.497049060  0.326127448  3.03612952
#>  -0.497635653  0.025837176  3.89821977
#>   0.511595249  2.312992968  1.82067416
#>   1.091038187 -0.034462538 -0.54568098
#>   1.533148292 -0.112495100 -0.53349553
#>   0.046629459 -1.849930331 -3.91871548
#>  -0.746735514 -1.048567472 -2.27530636
#>   0.662459583  0.719950014 -0.23680878
#>   0.690915028  2.544410484  0.11560673
#>   1.766309872  0.914980870  1.68124974
#>   1.668060465  3.865085992 -0.46966248
#>   0.192979555 -0.223103960  2.05493840
#>   1.423627470 -0.153455915 -1.45615460
#>   1.585275709 -1.543361506 -1.56431808
#>  -1.151132268 -0.065946115  2.28541246
#>   0.021960117 -0.570904602 -2.01040272
#>  -0.612264468  1.629277441  0.30044730
#>   0.476661657 -1.181237286 -2.38755695
#>   0.539846321  0.431343749 -2.96444583
#>   0.030181503 -1.085842062  0.00748961
#>   1.484289117 -1.010415413 -1.31875266
#>  -0.647460329 -0.811551491  2.10568691
#>   0.696409705  2.902212834 -0.37068537
#>   1.856657685 -1.609171599 -4.32381624
#>  -1.164849328 -5.896495247 -4.18395692
#>  -0.623575955 -1.195340511  1.34145934
#>   0.233084771 -0.149903628  1.27239356
#>  -0.522180913  0.755681555  2.28333426
#>  -1.900105472 -0.740870405  1.98517683
#>  -0.406098295  2.055889317  2.02227841
#>   1.431204133  1.610623307  0.82038191
#>  -0.620405461 -0.503021614  1.02519529
#>   0.232626257 -1.617075430 -4.28137270
#>  -1.392605254 -1.459315253 -0.88225562
#>  -0.704117995  1.104247395  2.38471460
#>   1.466536673  1.485304991 -2.86907406
#>  -0.774524183 -0.700927247  1.91060113
#>  -0.397081556  1.394765936 -0.42541076
#>   0.898157485  0.054890823  0.08999199
#>  -0.169023967 -1.166598910  1.84062105
#>   0.609004537  0.280556982 -0.42173193
#>   0.965298332 -0.243489856 -0.94235729
#>  -0.437678095 -3.483721011 -0.60090622
#>  -0.319820854 -0.773786383 -3.39116664
#>  -0.619934058 -1.947993544  0.80996385
#>   2.010640667  0.099647934 -2.64461158
#>   0.281203683  0.763872649 -2.34453586
#>  -0.183346698 -1.115227223 -1.21590857
#>  -1.635256642 -0.260033673 -1.37927037
#>  -2.013275544  1.296252932  1.87240944
#>  -0.800349857  2.802822285  4.15161075
#>   2.067305254  1.003402976 -0.96853874
#>   0.694731713 -2.999973152 -1.59351475
#>  -0.353668589 -1.042541294 -1.01085745
#>  -0.204873664  1.370425622 -0.12925070
#>   0.687262551 -0.266121440 -1.16473901
#>   0.169553567  1.394979213  1.28829765
#>   2.296765779  2.983452789  3.75765888
#>  -1.087284172 -1.149659966 -0.69572694
#>  -1.278022949 -0.495303023  1.14778476
#>  -1.616700099 -1.437547136 -0.27014012
#>  -1.602964433 -0.266083926  2.61991851
#>   1.260325668  1.942613458 -2.81680316
#>   2.344909956  0.482423047  3.27750590
#>  -3.903439934 -0.457171828  2.94091458
#>  -0.354030719  3.714586310  2.49094956
#>   0.487769150  0.675871582 -0.50477428
#>   0.681112830 -2.886138864 -1.01272318
#>   0.441522713 -0.446354163  1.14378422
#>   1.295845936 -0.569595686 -1.29669059
#>   0.284718656 -0.796772992 -1.65279901
#>  -1.314078965 -1.186899237  1.76212116
#>   0.978952241  0.720105920 -0.54972206
#>  -0.610423362 -0.722461112 -1.82534300
#>  -1.165142309 -2.793898885  1.07093083
#>  -2.303281512 -1.745799840  3.88520137
#>   1.373449779 -0.161796525  0.70090351
#>   2.262967263  0.300306830 -2.30329962
#>  -0.951858978 -3.451434464 -1.82411867
#>  -3.531080608 -1.753844149  1.51148657
#>  -1.516690154  3.776863768  1.19728976
#>   0.083238332  1.497744427  0.68298867
#>   0.378531198 -0.851644616  3.15521802
#>   2.013081203  1.154325257 -3.89841387
#>  -0.489989048 -2.124449554 -0.06302918
#>  -1.307377871 -1.443647560 -1.28218666
#>  -0.405924713  0.828463441 -2.00594342
#>   0.865756326 -2.680180940 -1.10804179
#>   0.544766299  0.007013980  0.85556367
#>   0.814050383 -0.045988243  2.04454392
#>   1.927714953  1.443865649 -2.87328518
#>  -0.695818356 -3.767990459  0.01778958
#>   0.854269246 -0.078942613 -1.66954033
#>   0.358827452 -3.431118500 -3.27829322
#>  -1.965766169 -1.527953364 -0.56869506
#>  -0.229681367 -1.210797040 -0.78460048
#>  -2.740612319  2.149870453  2.30742270
#>  -1.100170357  1.447209041  1.79706428
#>   3.071309187 -0.504279924  1.00122554
#>  -1.687996050  1.826173239 -2.61806032
#>  -1.388955467 -1.908209391  2.30283486
#>  -1.418841876  1.166338404 -2.44615224
#>  -3.945639869 -3.951060417 -0.72017601
#>  -1.767132483  1.547798248  1.95487509
#>  -0.634585986  1.736956614  2.19013980
#>   0.045645818  0.028257161  2.56306846
#>   2.977376384  4.647443943  1.22746749
#>   2.470228312  0.465410107 -2.31510903
#>   0.351905278 -2.522099147 -4.42924850
#>  -0.896767756  1.366271724 -2.86419767
#>   1.175544536 -0.212748037 -8.34320366
#>  -0.831585126 -3.433686927 -0.32903009
#>  -2.886974773 -0.342815340 -1.75999353
#>  -2.588112040  0.820700366  0.17998809
#>  -2.306874214 -0.023809925  4.26778571
#>   0.378263034  3.735250219  1.78598265
#>   0.251371294 -0.815042629  1.21902677
#>   0.683158453  2.302373223 -1.53684543
#>   2.321884420  1.479109215 -3.27687118
#>   1.651796170  2.068972041 -3.20822068
#>  -1.505277039 -2.019979406  0.71893276
#>  -2.363245556  0.747504677  0.32001456
#>  -2.760129581 -0.284318788  1.47948345
#>  -0.949998517  1.803601127 -1.82342152
#>  -2.449635494 -1.104447264  4.91162759
#>   0.115184941  2.136085841  2.25685605
#>   0.992604103  0.293416094  2.27885701
#>   2.224641221  2.523850920 -2.43656342
#>   1.403930451 -4.262168773 -1.97510432
#>  -4.739732857  1.117243809  0.18588941
#>  -1.744184009  2.395480846  4.70363975
#>   1.952937911  1.502875856  2.16561198
#>   1.678997657 -1.110154410 -4.74445547
#>  -1.640746261 -3.967146066 -3.91994660
#>  -1.896780452  1.038776557 -2.59460139
#>  -2.293432816 -2.588418387  1.99885135
#>  -1.527802883 -4.190949152  7.00237935
#>  -7.311031143 -0.720323339  9.71104591
#>  -0.992656693  6.158496767  3.11270436
#>  -0.503656995  0.745899965  5.75743080
#>   0.439008606 -1.576564530  7.07902819
#>   4.990405214  3.550253097 -5.09899382
#>   0.048703345 -2.142761340  1.83347508
#>  -2.261284091 -0.324713592  2.35513242
#>  -0.344812955  3.220180021  0.20213534
#>   1.803090222  1.716579408  3.54772409
#>   4.002255767  1.068336725 -5.02601864
#>  -0.645165493 -2.739886381  2.26416525
#>  -3.044556664  0.513514571  4.35097507
#>   1.407975877  4.304753220  0.42590263
#>   1.687758354 -0.578005160  1.58398882
#>   1.020600309  0.948936679 -0.58725290
#>   0.571120862  0.434304788  0.09779588
#>  -0.567824886  1.084658873  0.80859358
#>   1.575414924 -0.442883210 -3.88587693
#>  -0.204298170 -2.445556504  1.04074027
#>  -0.813221547 -1.789498289  4.50173450
#>   1.110517007  2.197138205 -2.76751864
#>  -1.315804448 -3.379398717  1.30295760
#>   0.526508818  0.495419160 -4.61464214
#>   1.710077229  1.661647060 -4.01477973
#>  -1.679814033 -4.114739260  1.97936962
#>  -0.116448974  5.940695488  1.25517999
#>   0.396686417  0.448564126  1.27110296
#>   0.020638509 -0.025437859 -3.04726981
#>  -0.606820730  1.412829353 -6.17437691
#>  -2.497218655 -1.135683411  0.10538932
#>   0.495551092  1.446978461  1.57544999
#>   1.791041438 -0.234350876  0.76295050
#>   1.301015556  2.360314839  3.80819805
#>   2.599683267  2.943984767 -1.28408131
#>   0.524634485 -2.601656318 -2.46448168
#>  -0.774453731 -1.246132428  0.53063613
#>  -0.572840135 -0.840887065 -0.76661024
#>   1.235576138  0.900561366  0.52269058
#>   1.277412670  2.807356863  3.15429446
#>   1.524322427  0.894845149 -0.98020596
#>   0.356003855 -0.494053200 -0.97868155
#>  -1.613821947 -0.831644824  2.82301530
#>   0.660029278 -1.500179924  3.50087378
#>   1.683271685  1.187570270 -3.54046625
#>  -1.799074598  0.033887474  2.92613319
#>  -2.050064536 -1.832514332  4.56710382
#>   0.816304985  3.399285446  0.65508655
#>   1.701783250  1.276504539  1.83705542
#>   3.361523171  3.152268076 -0.28540208
#>   1.315066426 -1.095305639  0.37212164
#>  -2.525202075  1.226711353 -1.87362724
#>  -0.474547063 -3.031125356  2.06630126
#>  -3.634084999  0.230047572  1.70379962
#>   1.134851717  2.304820947  1.42053624
#>   2.348253686  1.632426166  0.40224296
#>   1.328960074 -1.202481567 -3.42591506
#>   0.165014982 -1.612622412 -2.42274821
#>  -0.804150087  2.278895004  2.08562261
#>  -0.575180224  2.134639441  1.63739298
#>   2.081718372 -1.689378658 -0.62383705
#>  -1.308576620 -2.502945834 -2.37699660
#>  -0.644736967 -1.872900554  0.37658383
#>  -0.444194098 -1.398616909 -0.45780868
#>   0.661326211  1.366327996  2.69327414
#>   1.793132540  1.822316602 -4.93946905
#>  -2.340158789 -3.510067039  1.68332358
#>   0.237997372  2.538971869 -1.66629604
#>  -1.003257376 -1.635383513  0.50312095
#>  -1.385791706  0.805945113  1.95694797
#>  -0.475742214  0.316170014  0.33641169
#>   0.310739176  1.593181364 -4.58356615
#>  -2.572693651 -1.591801928  3.60249017
#>  -0.798601552  2.584433042  3.27102866
#>   0.215034624  4.140784118  2.03165760
#>   1.981113364  0.191814260 -1.17892942
#>   1.064989037 -0.619535255 -1.53706241
#>  -1.304775140  2.447275639 -3.18997284
#>  -0.399870431 -2.159498241 -1.13186977
#>   0.696122846  0.869126031 -1.63830205
#>  -0.900187765 -2.127069716  0.99191774
#>  -0.251134125 -0.842028731 -0.64394529
#>  -1.077271085  0.864466101  1.27478827
#>  -0.170533705 -0.248707839 -0.42619075
#>  -1.246817071  1.086192962  3.15009641
#>   0.527199352 -0.515338539  0.59929017
#>  -0.644960206  2.688566964  0.99844561
#>   1.630225439  0.997840334 -1.43597122
#>  -1.112087378 -2.843668476 -0.30306353
#>  -0.250427184  1.380010838 -3.62478126
#>  -3.061199507 -4.239652361  2.67148008
#>  -1.187411015 -0.578162378 -3.09579581
#>  -2.805983518 -0.874331423  2.68396343
#>  -1.388268632  3.413265657  1.81832261
#>   0.251591580  1.850551061  2.87908023
#>   2.082788075  0.873798477 -0.28242869
#>   3.315229874 -1.540461143 -0.85021958
#>  -0.166137787 -1.216051369 -0.29849800
#>  -2.948757882 -0.883069348  4.26062377
#>  -0.913583233  0.996018873  0.72029622
#>  -0.688885446  2.773983037  3.45337679
#>   2.416082064 -2.857795974  1.12163074
#>  -4.306937162 -2.435391833  3.33723577
#>  -0.741805020  1.525165610 -0.35831394
#>   1.832672111  1.274105160 -1.33483937
#>   0.868299864  1.209892302  0.42368957
#>   0.453900089  1.316148143 -0.59336007
#>  -0.311296486  1.132530232  4.33276394
#>   0.742669577  1.737741049 -0.61193482
#>   1.077722611  2.214191026  4.19734456
#>   2.350676877  2.783259886 -2.25926723
#>  -0.349982151 -3.876662769 -1.73275347
#>  -0.164654231  0.436564222 -0.68771283
#>  -0.833247568 -0.363527554 -1.66764394
#>  -2.821889832 -1.782292312  0.15117305
#>  -2.659364691 -2.154307063 -0.43553406
#>  -1.059152287  2.180301290  3.91852182
#>   1.373806131  2.094395664  1.67415400
#>   0.381671372 -1.654418226  4.07160216
#>   3.146531413  3.396018018 -0.12668687
#>   2.577402045  0.648325437  1.76300546
#>   2.896435656 -0.660746670 -2.61997230
#>  -0.023228062 -2.783661553  2.78913778
#>  -0.011408144  1.092616365 -3.54493284
#>  -1.291420572 -1.528497333  1.54922349
#>  -3.686927150 -2.799623599  0.83058252
#>  -2.470254156  1.648089175  2.30896863
#>  -0.602895694  0.819378776  1.01319264
#>   1.705983238  3.685930401  4.45567071
#>   3.997694704 -0.587282926 -4.81919445
#>   0.322867739 -4.206458895 -4.22435344
#>  -3.415668735 -0.540390968  2.89656307
#>  -0.237282675  2.737866216  0.67017778
#>   1.979927155 -2.199447498 -1.38295115
#>  -0.636723427 -1.089308722 -0.26194582
#>  -0.830674766  1.857520841  3.69765374
#>   1.907719662 -2.810513021 -2.09902275
#>   0.119845336 -1.218718387 -0.56041811
#>  -1.693997543 -0.945535522  0.34295904
#>  -2.561916517  1.026481448  4.51091860
#>  -0.791728125  4.762084275  3.34588248
#>   2.283221759  4.124248069  2.08594311
#>   1.860023229 -0.824885994  0.81382394
#>  -0.793175876  0.720583980  3.16222199
#>   1.710647068  0.417559195  0.04050092
#>   0.402759120 -0.872002004 -1.12632797
#>  -0.968332349  0.120461081 -1.73682368
#>  -1.069269274 -3.022299867 -0.09278887
#>  -2.713656915 -1.244675003  1.14048952
#>   0.352898240  2.229213755 -0.92081913
#>   2.200795266  0.152658884 -2.32754021
#>  -0.357837002 -1.724098325  3.62992034
#>   1.768154252  2.596572050 -2.79024480
#>   0.005313559 -1.375977885  1.89187889
#>   1.699955826  3.217964869 -4.68349435
#>  -0.266891591 -2.176632789  0.63096434
#>  -2.313247279 -1.454000543 -0.04422331
#>  -2.919896209  1.768034658  0.25297625
#>  -1.156271897  1.270014962  6.40101385
#>  -0.286462391  1.516773581  2.79651814
#>   2.343114372  2.075734908  0.58617008
#>   1.084809382 -0.146432064  0.93646153
#>   0.818724469 -1.695596579  0.91500714
#>  -3.133830343  1.337028696  1.01594457
#>  -0.832451280  1.954412269  0.58276683
#>  -0.515144797 -0.261400688  2.48610668
#>   0.777941775 -0.556218908 -1.59262428
#>  -0.976811113 -2.077104333 -0.32811608
#>  -1.994525925  2.470468110  1.35084880
#>  -2.583537232 -0.335013343  2.83274843
#>  -2.469531254  0.995300416  0.73158510
#>   0.459012788  4.030850825  4.59801125
#>   3.200071227  2.105884453  2.25717080
#>   2.908827102  2.749510437 -0.43628606
#>   1.899885772 -0.201898763  1.36681432
#>   1.745245417 -0.386468999 -1.60915804
#>   2.724942010 -2.001996065 -0.48454434
#>   2.165454931 -0.157858223 -6.89316432
#>  -2.712017778 -4.442921771 -0.91691302
#>  -1.973777693  1.179942176 -0.14674539
#>  -2.132905416  0.371001070  5.29313568
#>   0.880474486 -2.268961068  5.57467350
#>   3.612986013  3.649508701 -5.11500301
#>  -1.295054630 -2.355035002  0.31166707
#>  -0.821316045  0.289262153 -3.60192883
#>   2.840983716  2.399194034 -2.23228811
#>  -1.181396777 -3.314922109 -0.81524863
#>   0.286711097  1.713172452 -0.01104652
#>  -0.369814662 -0.213904723  2.04438180
#>   2.007457742  3.427525968 -2.59501056
#>   0.656160902 -1.887364618 -2.11545619
#>  -0.667155807 -2.352925057 -3.11792674
#>  -2.060467159 -4.317789412  1.64525588
#>   1.210143964  3.731647324  0.47023626
#>   3.198534851  0.521294318 -1.86461290
#>  -1.071866549  1.529931097  0.29929416
#>  -1.677927643 -0.258960877  0.78909734
#>  -0.580631819  0.023575792  2.27751534
#>   0.511881727  2.813066975 -0.94005158
#>   0.862774236 -0.644493005  2.38940669
#>   1.399653393  1.247567283  1.33449229
#>   1.530886592  0.838568386 -0.82021341
#>   0.179736877  3.715168583 -1.34961415
#>   4.336196126  0.795687321 -2.19933417
#>   2.444899080 -3.761958290 -2.74585396
#>  -1.204887183 -3.196421294 -1.54360892
#>  -0.897211434  4.013878851  1.70665260
#>  -0.344013219 -0.223617241  2.63002689
#>   1.630526377  1.111601299 -0.41683202
#>   0.669406464 -0.996246086  1.10780611
#>  -0.136438559  2.814518580 -0.39758973
#>   0.413713337 -0.242515636 -1.05354968
#>   0.485324293 -0.926972652 -2.19605610
#>  -0.957358971 -3.448021637  2.39631048
#>  -4.833269421 -1.097004129  6.62992531
#>  -0.893722190  6.935217809  2.60329585
#>   0.814646111 -1.330517091  2.08575838
#>  -0.329420709 -0.386759897 -0.87171999
#>  -1.118203747  0.964385827  0.77149548
#>  -1.482777236  1.379098035  0.18933901
#>  -1.437761841  2.003760609  2.60783429
#>   0.980280819  2.109439389 -0.29598846