Main smoothLRC algorithm
Arguments
- x
sparse matrix; assay matrix
- u_init
matrix; u initialization matrix
- v_init
matrix; v initialization matrix
- w
matrix; distance matrix. If null, computed to specification with k nearest neighbors.
- index
w indicies.
- lambda
positive numeric; penalization parameter.
- epsilon
positive numeric; convergence criterion.
- maxiter
positive integer; maximum desired iterations
Examples
sce <- example_sce()
x <- SummarizedExperiment::assay(sce)
coords <- SummarizedExperiment::colData(sce)[, c("col", "row")]
k <- 10
init <- smooth_init(x, coords, k)
lambda <- 1
epsilon <- 1e-3
maxiter <- 5
smooth_model(x, init$u0, init$v0, init$w, init$index, lambda, epsilon, maxiter)
#> iteration: 1 | convergence: 0.106201 | 0.0371337 | 1
#> iteration: 2 | convergence: 0.0903484 | 0.0311731 | 0.00525504
#> iteration: 3 | convergence: 0.0870562 | 0.0367362 | 0.00560428
#> iteration: 4 | convergence: 0.0876035 | 0.0403189 | 0.00625686
#> iteration: 5 | convergence: 0.0890037 | 0.0420041 | 0.0071044
#> $u
#> [,1] [,2] [,3] [,4] [,5]
#> [1,] 0.2007066735 8.865972e-02 -0.0671240212 0.0221808448 -0.3089801298
#> [2,] 0.0068779493 4.225074e-02 0.0634413733 0.0059770036 -0.0051897911
#> [3,] 0.1474356856 3.878247e-02 0.3524605550 0.2631222898 0.0083108161
#> [4,] 0.1758078569 7.641218e-02 -0.0994096059 0.0392977980 -0.0755143651
#> [5,] 0.0523636971 5.986077e-02 0.0098076570 0.0304484059 -0.0348481628
#> [6,] 0.1442057406 1.416731e-01 0.3133375466 0.1137707557 0.1003648849
#> [7,] 0.0532971171 4.387710e-02 -0.0863071615 0.0754363119 -0.0814809179
#> [8,] 0.2205740278 3.011001e-01 0.0490613055 -0.3482692322 -0.1070543084
#> [9,] 0.0746570395 5.414608e-02 -0.1160163893 0.0591123971 0.0629343738
#> [10,] 0.0106591946 3.037733e-02 -0.0576274600 0.0576634995 -0.0239671931
#> [11,] 0.0734189036 1.813488e-02 -0.1137814347 0.0093114216 0.0376100823
#> [12,] 0.0789922001 1.507262e-02 -0.1302131670 0.0274423331 -0.0010450493
#> [13,] 0.0612681065 3.513416e-03 0.0265970617 -0.1501539770 -0.0597147888
#> [14,] 0.1937663701 2.341532e-02 0.1456436556 -0.1878269022 -0.0366633205
#> [15,] 0.1065638279 1.347485e-01 -0.1152510032 -0.0907983826 -0.0122281114
#> [16,] 0.1925211758 -2.790384e-02 -0.1243203230 0.1794781533 0.1976379866
#> [17,] 0.0333029908 1.154662e-01 0.0357108169 -0.1163251629 -0.0362815652
#> [18,] 0.0675422292 9.447625e-02 0.1667602378 0.0658661931 -0.0162974902
#> [19,] 0.0610069136 1.073274e-01 0.2555974145 0.1054721570 0.0167433454
#> [20,] -0.0256668678 3.990149e-02 -0.0409612167 0.0708914994 -0.0402666261
#> [21,] 0.1295508997 4.871966e-02 0.0757343049 -0.1185101262 -0.0313169427
#> [22,] 0.1229994476 1.001976e-02 0.1008258802 0.0895054918 -0.0310976643
#> [23,] 0.0045216409 3.818953e-02 -0.0948238753 0.0662675543 -0.0009356166
#> [24,] 0.1010423838 1.497169e-01 0.0601373379 -0.2851338137 -0.0511941482
#> [25,] 0.0605055328 1.475781e-01 0.0376472791 0.0259878629 -0.0164971186
#> [26,] -0.0151799690 9.034042e-02 -0.0617189230 0.0561316178 0.0042896253
#> [27,] 0.0867396688 1.207830e-01 -0.0333724719 0.0655977713 -0.0461601368
#> [28,] 0.1170477000 -5.921024e-02 -0.0570784354 0.0254073965 0.1987673309
#> [29,] 0.1098341952 1.357562e-01 -0.1061515406 -0.1612299855 0.0160733449
#> [30,] -0.0035738690 3.522324e-02 -0.0221909269 -0.0082103036 -0.0123589763
#> [31,] 0.2603915083 -8.218033e-02 -0.1273625625 0.2449432940 0.2583705736
#> [32,] 0.0168444599 5.903458e-02 -0.0148333523 0.0328390465 0.0081771004
#> [33,] 0.0291858334 4.754408e-02 0.0115379053 0.0256341993 0.0266797942
#> [34,] 0.0082246248 5.921923e-02 0.0723520957 0.0205583231 -0.0504449555
#> [35,] 0.0538666178 1.505109e-02 -0.0687611900 -0.0007870401 0.0310764628
#> [36,] -0.0210291773 6.749544e-02 -0.0227602559 0.0277247774 -0.0090546807
#> [37,] 0.1890950618 -1.549371e-01 0.1395826214 0.1159481212 -0.1623220098
#> [38,] 0.0792824202 -2.090137e-02 0.0052084003 0.0992171516 -0.2573076918
#> [39,] 0.0252860294 5.100394e-02 0.0306379392 0.0190346581 -0.0364963579
#> [40,] 0.0472026974 1.690404e-02 0.1980241203 0.1943563565 -0.0325927848
#> [41,] 0.1384828110 -3.427575e-02 0.0022426605 0.0356678836 -0.1759588872
#> [42,] -0.0031611172 9.505073e-02 -0.0297923511 0.0649329683 -0.0012840483
#> [43,] 0.2563036641 -4.689148e-01 -0.2009638709 0.0979694737 0.0073417811
#> [44,] 0.1162037849 1.833571e-01 -0.1439317651 0.0698613063 -0.0308369964
#> [45,] 0.0064284386 7.864658e-02 -0.1086020661 0.0548502460 -0.0496703852
#> [46,] 0.0725679327 7.660213e-02 -0.0158569175 0.1188521489 -0.0292352646
#> [47,] 0.1479485089 -3.704375e-02 -0.1071841346 0.0599274149 -0.1389458756
#> [48,] 0.0780417763 1.216165e-01 -0.1163015976 0.1888882950 0.0124732687
#> [49,] 0.1104965723 3.082087e-02 -0.0683676104 -0.0416542627 -0.0922662861
#> [50,] 0.0861975576 -2.096414e-02 0.0055342590 -0.0966823289 -0.1073734990
#> [51,] 0.0552893731 7.545413e-02 -0.0829414015 0.0290662754 0.0950316026
#> [52,] 0.0570641393 1.238584e-01 -0.0963074105 0.0467455776 -0.0509507345
#> [53,] -0.0047869770 2.336536e-02 -0.0452590659 0.0588493087 -0.0442220726
#> [54,] 0.0321476410 4.344545e-03 -0.1072919723 0.0446683501 -0.0438232449
#> [55,] 0.0404498762 7.025164e-02 0.0986229655 0.0243401407 -0.0493128627
#> [56,] -0.0188884364 3.485793e-02 -0.0633120663 0.0516649375 -0.0305288684
#> [57,] 0.1438942261 -2.922097e-01 -0.2487522278 0.1878343368 0.0880538835
#> [58,] 0.0091004708 2.723583e-02 -0.0569561706 0.0664620846 -0.0208367770
#> [59,] 0.0881333636 9.354251e-02 0.0058170309 0.0429191698 0.0710920385
#> [60,] 0.1126481808 1.879270e-02 0.0831718915 0.3164338063 -0.0537747504
#> [61,] 0.0387658504 -9.413197e-03 -0.0710769728 -0.0651831050 -0.0256980641
#> [62,] 0.0889589125 -5.139547e-02 -0.0703897367 -0.0604056034 -0.0596675990
#> [63,] 0.1678009465 -1.177955e-02 -0.0443380833 -0.2011215064 -0.0994761493
#> [64,] 0.1151844934 -3.067110e-02 0.2492941582 -0.1375567204 -0.0624134264
#> [65,] 0.0244754667 1.735743e-02 -0.0350143719 0.0796422503 -0.0330602176
#> [66,] 0.1056211361 4.834840e-02 0.0506117875 0.0814127630 0.0657425919
#> [67,] 0.0401060086 2.372619e-02 -0.0657800635 -0.0371657376 0.0461537892
#> [68,] 0.0303080907 1.027638e-01 0.0413353841 0.0100017463 -0.0151455793
#> [69,] 0.1147273315 2.143274e-01 -0.1383956209 0.0398317935 0.0580743442
#> [70,] 0.0224416041 3.709760e-02 0.0131977328 -0.1139083689 -0.0268965336
#> [71,] 0.0372539557 6.624342e-02 -0.0142656288 -0.1048656592 -0.0549175282
#> [72,] 0.2416317360 3.027947e-01 -0.0809408230 0.1453269456 0.0310899091
#> [73,] 0.0415510289 3.147302e-02 0.0160845051 -0.0984404385 -0.1100811349
#> [74,] 0.1867250409 -2.039713e-01 0.1228426297 0.0890582923 -0.2283637017
#> [75,] 0.0281007989 -1.294362e-02 -0.0665732715 0.0623690218 -0.0491946894
#> [76,] 0.2124574120 1.284490e-01 0.1546760809 0.2458556387 -0.1314332438
#> [77,] 0.1170150274 4.374834e-02 -0.0698689896 -0.0770929413 0.0197378109
#> [78,] 0.1496022375 1.203048e-01 0.0930899996 0.1297441365 -0.2636931190
#> [79,] 0.0383910007 1.395806e-01 0.1708355939 0.1079151463 0.0389262885
#> [80,] 0.0475253807 -2.097583e-02 0.1505831601 0.0316147923 0.0245105832
#> [81,] 0.1206395620 -1.463422e-01 -0.1088588030 -0.0470422781 0.0154885425
#> [82,] -0.0331112575 1.630393e-02 0.0034341302 0.0425238146 -0.0325672418
#> [83,] 0.1774825031 1.987985e-01 -0.0212239187 0.1796320754 0.0291993305
#> [84,] 0.0991928015 6.892928e-02 0.0062164983 0.0312663157 0.0215914195
#> [85,] 0.0447199794 7.454482e-02 -0.1255961428 0.0636132262 0.0691665125
#> [86,] 0.0328986092 -2.332612e-02 0.1180202524 0.1494714557 -0.0556888092
#> [87,] 0.1529226087 5.017280e-02 -0.1887237306 0.0548737319 -0.0512544022
#> [88,] 0.1729656718 1.180506e-01 0.1499932557 0.0411830573 -0.0681602828
#> [89,] 0.2621162536 -5.000792e-01 -0.1106533400 -0.2152392859 0.0013071658
#> [90,] 0.1535433038 1.572316e-01 -0.0942827752 0.0516211273 -0.0963399642
#> [91,] 0.1195379493 7.183507e-02 -0.0745787403 0.1934899467 -0.0692007911
#> [92,] 0.0140816701 -3.118873e-02 -0.0332242656 0.0166702593 -0.0371422785
#> [93,] 0.1111031164 3.210763e-02 -0.1593914109 0.0013792903 0.0275452972
#> [94,] 0.0634513237 8.022598e-02 -0.0993934927 0.0596453752 0.0882019987
#> [95,] 0.0575641990 5.256782e-02 -0.1216738017 -0.1150217454 0.1046296464
#> [96,] 0.1154559521 2.719382e-01 0.1449791775 -0.2147402383 -0.0467345114
#> [97,] 0.0888308958 -1.219501e-01 0.0233159236 0.0238150619 -0.0410215206
#> [98,] 0.0145436194 5.322205e-02 0.0425575130 0.0685137711 0.0298201816
#> [99,] 0.1295013915 -6.433418e-02 -0.1916729790 0.1399459422 0.1570209984
#> [100,] 0.0458464909 1.155659e-02 -0.0196672297 -0.0096317892 -0.1274901282
#> [,6] [,7] [,8] [,9] [,10]
#> [1,] 0.208176331 -0.4244605407 -0.1619187888 0.0493685146 1.579639e-01
#> [2,] 0.022628988 -0.0482215840 0.0392023634 0.0850043012 4.986152e-02
#> [3,] 0.013176348 0.0328130149 -0.1097599359 0.1661078603 5.372837e-02
#> [4,] -0.018397104 -0.3434033499 -0.2987413563 -0.1057473747 -1.285573e-02
#> [5,] 0.059126072 0.0910942357 -0.0435812612 0.0783955045 -9.985797e-03
#> [6,] 0.097532640 0.1930207003 -0.0379461082 0.1969013801 2.720401e-01
#> [7,] 0.061465838 0.0743903522 -0.0353751456 -0.0754749343 1.566203e-02
#> [8,] 0.080783375 0.0357875574 -0.3653210065 -0.0044682710 5.350230e-02
#> [9,] -0.030165658 0.0879251969 -0.0696500587 0.0140491269 -5.705238e-02
#> [10,] 0.110955589 0.0691715212 -0.0466325887 -0.0422039952 -9.962702e-02
#> [11,] 0.094739339 -0.0432288019 -0.1064740491 0.0370479313 -3.840687e-03
#> [12,] 0.161758560 0.1314541297 -0.0683662313 0.0273256807 -3.181473e-02
#> [13,] 0.051407686 0.1031704491 0.1543568298 0.1388908789 -7.203436e-04
#> [14,] -0.045625988 0.0996786724 -0.0066474949 0.1607933154 2.951463e-02
#> [15,] 0.110548586 0.0529091025 -0.1110524685 -0.1014682465 -7.313268e-02
#> [16,] -0.211690206 -0.1297453968 0.0091113458 0.1695411473 -1.184478e-02
#> [17,] 0.183715713 -0.1167408687 -0.1533352669 -0.1095224540 9.392686e-02
#> [18,] 0.022726447 0.1181676014 -0.1201887429 0.1025726325 5.755983e-02
#> [19,] 0.071116190 0.1884076650 0.0251634116 0.1130796369 2.098902e-01
#> [20,] 0.039465340 0.0126455634 -0.0563474716 0.0412239412 3.823639e-02
#> [21,] 0.147839784 0.1147806112 -0.1210666075 0.0313981738 2.147407e-02
#> [22,] -0.059758309 0.0696322904 -0.0017680137 -0.0498824430 2.820375e-01
#> [23,] 0.049099396 0.0047580132 -0.0926146699 0.0036741697 2.459304e-02
#> [24,] 0.186231397 -0.1626310548 -0.2937054169 0.0558043487 5.938487e-02
#> [25,] 0.105874280 0.1160972520 0.0277934980 0.1218299740 1.391006e-01
#> [26,] 0.085374701 0.0751815459 -0.0393976495 0.0578134752 6.629325e-04
#> [27,] 0.161411169 -0.0205577328 0.0435181184 0.0199850339 -8.226140e-02
#> [28,] -0.109466189 0.0239142936 -0.0970156849 -0.1662879840 2.921723e-02
#> [29,] -0.109619161 -0.0551339634 -0.1346530204 0.0231836143 5.614589e-03
#> [30,] 0.067285441 -0.0180484693 -0.1020909814 0.0014935568 6.372262e-02
#> [31,] 0.188788253 -0.0749112020 -0.1175861988 -0.4485472861 -4.690547e-02
#> [32,] 0.089148261 0.0637398502 -0.0039819427 0.0510979652 1.151467e-01
#> [33,] -0.077914665 0.0876528378 -0.0493457718 -0.0291340190 7.959136e-02
#> [34,] 0.085709356 0.0919373431 0.0242801466 0.0678066472 1.629041e-01
#> [35,] 0.040867260 0.0552961247 0.1069108843 0.1287971075 -7.154997e-02
#> [36,] -0.048461916 0.0806186548 -0.0899894313 0.0233954870 6.883652e-02
#> [37,] 0.028019191 -0.5250050832 -0.1732366467 0.0979486427 -7.213828e-02
#> [38,] 0.093240187 0.1526073426 -0.0464657864 -0.2177966328 3.411686e-02
#> [39,] 0.062957316 0.0959364963 0.0220498249 0.0510679190 1.201799e-01
#> [40,] 0.066308421 0.0276974485 -0.0887417370 0.0377701550 3.740978e-02
#> [41,] 0.017592053 -0.3760306726 0.0929267817 -0.0562214915 -2.732554e-02
#> [42,] -0.010704172 0.0492120645 -0.0941979540 -0.0486691439 -7.304380e-04
#> [43,] 0.142774625 0.1752291069 -0.2556335775 -0.2038525245 1.965002e-01
#> [44,] 0.032122681 0.1046700939 -0.1190775121 -0.1233854249 3.097333e-02
#> [45,] 0.032284105 0.0740106936 -0.0278407799 0.0126585070 3.943443e-02
#> [46,] 0.084926424 0.1558729044 -0.0829627838 0.0371797190 7.595525e-02
#> [47,] 0.101415488 0.1955115827 -0.1505089851 0.1156619380 -7.052376e-02
#> [48,] -0.055893946 0.0199581196 0.2102985337 0.2366544906 7.243967e-02
#> [49,] -0.118296471 -0.0732045005 -0.0710165121 0.1320079472 -8.825296e-02
#> [50,] -0.177257921 0.0425706282 -0.0703064361 -0.1301330588 1.087575e-01
#> [51,] -0.003517092 0.0092585609 0.0138239696 0.0430983276 -4.077510e-04
#> [52,] 0.033208017 0.0212280298 0.0310855027 0.1084585074 5.785034e-02
#> [53,] 0.024575246 -0.0368760802 -0.0402375122 0.0669341338 2.022732e-02
#> [54,] 0.136207992 0.0935482372 -0.0702770232 -0.0877565360 -1.016640e-02
#> [55,] 0.060770488 0.0476462320 0.0887371505 0.0707682592 -5.541385e-02
#> [56,] 0.074185002 0.0305312849 -0.0496613999 0.0240096636 -5.364790e-05
#> [57,] 0.138870806 0.0084158113 0.0297186760 0.2210063530 2.154089e-01
#> [58,] 0.048976837 0.0542520460 -0.0732097184 0.0201069106 7.232679e-02
#> [59,] 0.074801675 0.1354008581 -0.1559022550 0.0009583226 -2.825072e-02
#> [60,] 0.132141808 -0.1957834906 -0.0342797758 0.0968653761 3.177423e-02
#> [61,] -0.136166917 -0.0091143881 -0.0578239464 0.0629451631 -1.237207e-02
#> [62,] -0.069428151 0.0873483531 0.1131870717 0.0978109422 -5.863740e-02
#> [63,] -0.177585908 -0.0481565336 -0.2279521737 0.1734434752 4.402363e-02
#> [64,] -0.164015598 -0.0374658755 0.0332035893 0.1298113889 2.753941e-02
#> [65,] 0.058900742 -0.0028157879 -0.0745807220 0.0573134765 -5.584940e-02
#> [66,] 0.195280169 0.1767338816 0.0667018760 0.0596501579 3.519800e-02
#> [67,] -0.067374395 0.0817759191 -0.0531172069 0.1049032653 -4.063321e-02
#> [68,] 0.115095365 0.1060956908 -0.0085870128 0.0631923895 1.445944e-01
#> [69,] 0.083026341 -0.0069618021 0.0094252366 0.0934708541 1.107691e-01
#> [70,] 0.075473585 0.0003144115 0.0176762498 0.0583964234 -3.113601e-02
#> [71,] 0.036741062 -0.0768436815 -0.1270287995 0.0726802253 4.865493e-02
#> [72,] 0.107630362 -0.0455540484 0.0127211180 0.0646582175 -1.851277e-01
#> [73,] 0.090981520 -0.0445829503 -0.0100775780 -0.0281960496 9.884348e-02
#> [74,] -0.157084160 -0.3604296255 -0.1265321512 -0.1545001392 1.091936e-01
#> [75,] 0.048468607 0.0312810231 -0.0337865907 -0.0155399955 1.471336e-02
#> [76,] 0.184830734 0.1807541477 -0.2931806752 0.0801507714 4.271267e-02
#> [77,] 0.134697012 0.0908029684 -0.0412413069 0.0694495939 -6.035475e-02
#> [78,] -0.046389239 0.0299368232 0.1319093421 -0.0921905942 5.607704e-02
#> [79,] -0.011905204 0.1625075195 -0.0718334845 0.1029744775 1.162091e-01
#> [80,] 0.170106802 0.0524242477 0.0721909691 0.0620962307 3.998930e-02
#> [81,] 0.186706604 0.0327082410 -0.1539707091 0.1181008444 1.591435e-01
#> [82,] 0.090065184 0.0275427870 -0.0067536644 0.0356092341 4.033550e-02
#> [83,] -0.176001362 0.1795017330 0.1520085964 0.2404368460 1.844633e-01
#> [84,] -0.018750299 0.1141314299 -0.2474429173 0.0075622993 8.455163e-02
#> [85,] -0.077659827 -0.0128479198 -0.0367031117 -0.0815153687 6.587492e-02
#> [86,] 0.122186093 -0.0828534461 -0.0889982215 0.0363676598 4.995994e-03
#> [87,] -0.006931496 -0.0902201640 0.0226358992 -0.2278934623 3.200014e-02
#> [88,] 0.169108342 -0.1667514666 -0.1606780731 0.1260788089 1.709060e-02
#> [89,] 0.060337546 0.1228745640 0.1149764035 0.2887513923 1.423535e-01
#> [90,] 0.043276389 0.0947810704 0.2051392901 0.2120636180 6.012679e-02
#> [91,] -0.140439750 -0.0079812843 0.2450546830 0.1205803358 1.423457e-01
#> [92,] 0.014545082 0.0509030824 -0.0440998758 0.0613515966 5.196905e-03
#> [93,] 0.184132302 0.1110476920 0.0037206042 0.0973691931 -2.936501e-01
#> [94,] 0.151048221 0.0082513143 0.1084265301 -0.0163743472 -1.734043e-01
#> [95,] 0.062644157 0.1343520112 -0.0777500884 0.1144203604 -8.201987e-02
#> [96,] 0.233089665 -0.1185246196 -0.2227373588 0.1158359163 3.163908e-02
#> [97,] -0.235914116 -0.0872562720 0.1003636033 0.1336223890 -9.587557e-02
#> [98,] 0.167519158 0.1023879834 -0.0132203565 0.0309516046 -1.766326e-02
#> [99,] 0.021517058 -0.0151734950 -0.0924849608 0.0614425357 2.965564e-01
#> [100,] -0.051733088 -0.1831777984 -0.1022405821 0.0628142978 6.648866e-03
#> [ reached getOption("max.print") -- omitted 900 rows ]
#>
#> $v
#> [,1] [,2] [,3] [,4] [,5]
#> [1,] 0.2149260 -0.076190416 -0.0822264127 -0.2680622085 -0.088755782
#> [2,] 0.2974099 -0.232937362 -0.0901852067 -0.0483096043 -0.052899648
#> [3,] 0.3803527 0.489017270 0.5391790434 -0.0310228671 0.004716606
#> [4,] 0.2849800 0.400455240 0.1658514778 -0.0294860829 0.516942557
#> [5,] 0.4885513 0.677955132 0.6125165970 -1.1562645144 -0.155175908
#> [6,] 0.4722826 0.026631083 -0.2578209026 -0.4143583341 -0.466522173
#> [7,] 0.2505574 0.064087150 0.0925318622 0.2576670305 0.067305827
#> [8,] 0.2253700 -0.059300085 0.0150463412 -0.0335798270 0.177201830
#> [9,] 0.3175433 0.399589016 0.0537158239 -0.2465391217 0.069136695
#> [10,] 0.5258952 0.962911199 -0.8886627745 0.3508963388 0.078743600
#> [11,] 0.2303544 -0.065997430 0.1947274535 -0.2671454711 0.060202358
#> [12,] 0.3282651 -0.160493664 0.1385537558 -0.2120140026 0.509570564
#> [13,] 0.3228640 -0.306730119 0.1888055571 0.4611237092 0.365528059
#> [14,] 0.2342243 -0.185433574 -0.2301694613 -0.0001633888 0.227240111
#> [15,] 0.2014615 -0.180938351 0.0900545166 -0.0856514023 -0.108630519
#> [16,] 0.2572740 -0.067160079 0.0663337931 0.1048991230 -0.038117931
#> [17,] 0.3856626 -0.083690259 0.3803492442 0.6825630433 -0.119483416
#> [18,] 0.3021844 0.034948908 -0.2266511422 -0.0474070970 0.084740279
#> [19,] 0.2851221 -0.098416965 0.1545435813 -0.0923970828 0.254524320
#> [20,] 0.5328422 -1.949104402 -0.9771492112 -0.0619439768 0.369542679
#> [21,] 0.3222199 0.192036181 -0.1492608182 -0.1117176473 0.172440298
#> [22,] 0.1540626 -0.080286063 -0.0018381742 -0.0773085134 0.190073355
#> [23,] 0.2874883 -0.020398321 -0.0655621858 0.0510557254 -0.281658075
#> [24,] 0.2739418 0.035421543 0.0900500387 0.2209373124 0.001857070
#> [25,] 0.4008466 0.470269547 -0.0411605156 -1.0289755144 0.537932658
#> [26,] 0.3837962 -0.462649264 0.0422903344 0.0878602156 0.455321338
#> [27,] 0.4235156 -0.602883724 0.5915828841 0.3673574262 -0.556106952
#> [28,] 0.2568168 0.003390262 0.1157055351 0.0343647311 -0.052906957
#> [29,] 0.3607412 0.769415099 -0.3902604980 0.6635494374 -0.102940001
#> [30,] 0.2475474 -0.483461138 -0.1607114216 -0.2968399920 -0.354084025
#> [31,] 0.3075315 0.029731818 0.0576052587 0.3089014205 -0.429483726
#> [32,] 0.2188245 -0.209851108 0.2381509207 0.0389135957 0.020540955
#> [33,] 0.3533810 -0.056417296 0.2064956405 -0.4435013706 0.037010439
#> [34,] 0.1998890 -0.043688528 0.1664296844 0.0206046153 0.257260288
#> [35,] 0.2004848 -0.130836632 0.0711531102 -0.0780530306 -0.052179616
#> [36,] 0.2482743 -0.162133436 -0.0314601953 -0.0970668644 0.062552320
#> [37,] 0.3794031 -0.171107197 -0.1409103124 0.2956899868 1.107064676
#> [38,] 0.3068617 -0.008584008 -0.2006073302 0.0916998906 0.326934694
#> [39,] 0.3774921 0.297789679 -0.4535301836 0.0500693418 0.585463349
#> [40,] 0.2820791 -0.135704393 0.5190805015 0.0709882989 0.105457133
#> [41,] 0.2773272 -0.371826026 -0.0336403497 -0.1442956766 -0.048735280
#> [42,] 0.2947680 0.039454913 -0.1330100071 0.0981748056 -0.197113608
#> [43,] 0.2237271 -0.146204894 -0.0123229801 -0.0722456542 0.051779133
#> [44,] 0.4893929 0.041773064 -0.4866427394 0.3658038681 0.693396063
#> [45,] 0.3617168 -0.567190372 0.9087851351 -0.0545295751 0.211781264
#> [46,] 0.2503159 -0.186543493 0.0620409823 -0.0568035794 -0.092979131
#> [47,] 0.2323670 -0.098880133 -0.0350956006 -0.0345656037 0.132960984
#> [48,] 0.2834040 0.134723128 -0.2263320209 -0.4001687401 0.359097062
#> [49,] 0.2260726 0.123392451 -0.0943935892 -0.0521986671 -0.035178679
#> [50,] 0.2670545 -0.308342775 -0.0008648042 -0.1166895752 0.072783302
#> [51,] 0.2439850 -0.021424675 -0.1123268374 -0.1087985453 0.237570737
#> [52,] 0.2974427 -0.074175177 -0.1260432002 0.1693843777 0.244866871
#> [53,] 0.2007577 -0.118086624 0.3713457552 -0.0560206711 -0.197539814
#> [54,] 0.3850171 0.133816261 0.2917843135 0.0132542783 -0.023005848
#> [55,] 0.2666920 -0.115228898 -0.0359589458 0.3973397956 -0.175941945
#> [56,] 0.5422125 0.103843196 -0.5977387652 -0.0642093520 -1.947392760
#> [57,] 0.4237825 0.663122789 -0.4043623298 -0.5131672259 -0.013179763
#> [58,] 0.3148849 0.185369973 0.0122641053 0.3429166455 0.084489146
#> [59,] 0.2836933 0.035987587 -0.0782909331 0.0772479146 -0.154080314
#> [60,] 0.3713090 0.145371071 0.1528781717 -0.1680665462 0.024712810
#> [61,] 0.2825946 -0.110375344 -0.1301149892 -0.2375007506 0.009582661
#> [62,] 0.3029739 -0.315281663 0.3716817926 -0.0428740086 -0.005890357
#> [63,] 0.2113266 0.212236609 0.0770854037 0.3184069149 -0.142519887
#> [64,] 0.1941621 -0.148110851 0.0209731956 -0.1705038062 0.186098778
#> [65,] 0.2413661 -0.374593814 -0.0241057686 -0.0349413847 -0.011133682
#> [66,] 0.3652099 -0.233862183 -0.3143100465 -0.2447641855 0.009260882
#> [67,] 0.2280040 0.134207919 -0.1773819032 -0.1748581789 0.340656830
#> [68,] 0.3756730 0.261653746 0.0385259025 -0.1414398182 0.215887496
#> [69,] 0.2802568 -0.297344742 -0.0872736295 -0.2018010609 0.062441997
#> [70,] 0.4792029 0.528826766 0.0080841026 -0.7257836954 -0.025143970
#> [71,] 0.2677620 0.396922978 -0.0740489969 0.1718078696 0.295805542
#> [72,] 0.3049306 -0.029236244 -0.2893005208 -0.1787351194 -0.128458190
#> [73,] 0.4742746 -0.209666201 -0.0959889713 0.1239931896 -0.471218648
#> [74,] 0.5117096 0.423306042 1.6612708828 0.8917260353 0.448592988
#> [75,] 0.2981851 -0.132351141 -0.4053110756 0.2350450072 0.284131203
#> [76,] 0.3732082 -0.011811016 -0.0659619601 -0.6645208779 0.158884685
#> [77,] 0.3351017 0.068172313 0.2240278277 0.0473466200 -0.296940247
#> [78,] 0.2840126 -0.476961395 -0.2619424656 -0.0046694786 0.206718753
#> [79,] 0.4099535 -0.446932996 0.7843617008 1.1441859323 -0.144452308
#> [80,] 0.4414757 -0.182042783 0.5907650958 -0.5053725611 -0.315021549
#> [81,] 0.3489320 0.003042576 -0.1442799473 -0.0411653764 0.187536740
#> [82,] 0.4525799 0.246623366 -0.1751828678 0.7654963894 -0.865706229
#> [83,] 0.2656481 -0.400646732 -0.1181721326 -0.0762009136 0.034375926
#> [84,] 0.2830458 0.157693752 -0.1909134851 0.0782065966 0.153281751
#> [85,] 0.1930024 -0.182612498 0.1924962674 -0.1471814320 0.083300173
#> [86,] 0.2924796 -0.242241368 0.5356802558 -0.4947953212 0.446963949
#> [87,] 0.2570445 -0.187105456 0.0586041631 -0.0842090344 -0.131231020
#> [88,] 0.2804309 0.158442242 0.2147863508 -0.1679680407 0.265124886
#> [89,] 0.3057172 0.126014457 0.4608058605 -0.4738780869 -0.229220991
#> [90,] 0.2973654 -0.233879948 0.0429093674 -0.2162018343 -0.559440661
#> [91,] 0.4056711 0.216237703 0.6641400556 -0.2061471000 -0.011149722
#> [92,] 0.3025589 -0.128097667 0.4446207105 -0.2293002293 0.011846570
#> [93,] 0.5172331 0.572953234 -0.7714388093 0.9818085601 0.192697226
#> [94,] 0.2785030 -0.404921073 -0.4513148827 -0.1058019288 0.043498628
#> [95,] 0.1972177 -0.120601454 0.1169447614 0.0823344899 0.255632481
#> [96,] 0.3859856 -0.462559975 0.1431492172 -0.4855087484 -0.346314778
#> [97,] 0.4239506 -0.234623366 -0.0163490644 -0.6229517371 0.021633756
#> [98,] 0.2754687 0.298160466 0.0735441924 -0.2754830599 -0.011155356
#> [99,] 0.2839010 -0.383913076 0.0116048572 0.1047186716 0.046960720
#> [100,] 0.3702165 -0.273941741 0.4472915457 -0.2277925109 -0.233898219
#> [,6] [,7] [,8] [,9] [,10]
#> [1,] 0.2547440286 0.168580889 0.238120369 -0.178895832 -0.038999015
#> [2,] -0.0164694431 -0.079552584 -0.190327378 0.145703678 0.336695347
#> [3,] 0.3149674463 -0.159119998 0.311928987 0.150242263 -0.266908220
#> [4,] -0.2733627187 0.160516283 -0.267054246 -0.486495219 -0.167115839
#> [5,] 1.1960926641 -0.604500148 -0.942524266 0.078738985 0.277218956
#> [6,] 0.4849638784 0.586423244 0.024577359 -0.471299361 -0.228243604
#> [7,] 0.2602688165 -0.565803266 0.186190125 0.186915626 -0.049198341
#> [8,] -0.0314835697 -0.124186011 0.241863182 -0.280811722 -0.100573117
#> [9,] 0.4107510602 -0.379783993 -0.127956552 0.197622360 0.584510871
#> [10,] 0.0365730128 0.176126557 0.108033824 -0.362374070 0.097757108
#> [11,] -0.2974601837 -0.198393453 0.202273877 -0.007740457 -0.117793583
#> [12,] -0.2062595940 0.163015120 0.047647233 -0.181176253 0.244200423
#> [13,] 0.3177397804 -0.001065707 0.174220431 0.200680911 -0.266671030
#> [14,] -0.0307815163 0.218336491 -0.284536782 -0.158758656 -0.229542818
#> [15,] -0.1538896826 -0.179274634 0.103436551 0.273090505 -0.238357492
#> [16,] 0.0609852669 -0.281550116 0.298061968 -0.080507801 0.008839671
#> [17,] -0.2594970515 -0.275185404 0.574175116 -0.142554181 0.362411750
#> [18,] -0.4702621745 -0.076632518 0.365183457 0.103563270 -0.046762562
#> [19,] -0.1670904550 0.071321475 0.150762397 -0.206908581 0.231860148
#> [20,] 0.7595153547 0.287018207 -0.376119394 0.898159118 0.971702426
#> [21,] -0.0112900200 -0.500625069 0.448005550 -0.591598877 0.108089440
#> [22,] -0.2157457153 0.102106354 0.029750065 -0.176875075 -0.156582849
#> [23,] -0.0361053609 0.163938825 0.304640298 -0.072222067 0.108088967
#> [24,] -0.1711443390 -0.232283941 0.023201688 -0.212317961 0.437878253
#> [25,] 0.2394987388 -0.082703160 0.719735215 -0.335206216 0.603122823
#> [26,] -0.1676885446 -0.347194001 -0.278206789 -0.476664067 -0.021677305
#> [27,] -0.1447884096 0.244841883 -0.004201765 -1.267264098 0.298086623
#> [28,] 0.0689683297 -0.181751316 0.116788936 -0.105583195 0.053546737
#> [29,] 0.0336695323 -0.227913020 0.167505585 0.746992221 -0.038482186
#> [30,] -0.2595428668 0.374450384 0.011246945 0.129106641 -0.332461203
#> [31,] -0.1029592202 -0.307485506 -0.322067741 0.223347247 -0.411762696
#> [32,] -0.1860044600 -0.004671251 0.207380768 -0.311526677 -0.193561335
#> [33,] 0.3539427930 -0.412844422 -0.644798642 -0.288636145 0.014714085
#> [34,] -0.1368834375 -0.072060078 0.062925977 -0.124409054 -0.053840953
#> [35,] -0.3776745738 0.007010752 0.283885004 -0.059463311 -0.158191376
#> [36,] -0.0961602521 -0.118665908 0.167970399 -0.118466411 -0.218569160
#> [37,] 0.0374307193 -0.036637713 0.378914472 -0.520700884 -0.637846004
#> [38,] 0.4980481581 0.134463430 0.139236384 -0.480705127 -1.066487750
#> [39,] -0.6832876911 -0.289318188 0.173715387 -0.238475216 -0.290694719
#> [40,] 0.0289828001 0.040361844 0.179857022 0.076731262 -0.077058849
#> [41,] -0.5220171478 0.071087116 -0.097549725 0.023060073 -0.214726073
#> [42,] 0.3983257880 0.290219349 -0.262429721 -0.293749449 -0.351761845
#> [43,] 0.0024515738 -0.142922097 0.051711640 -0.009982002 -0.072819497
#> [44,] -0.5919924640 -0.231619309 -0.468293711 -0.572436096 0.992497128
#> [45,] 0.8383392720 0.049944049 0.668682258 -0.147813152 -0.087509893
#> [46,] 0.0496007839 -0.324037755 0.240942797 -0.090076234 -0.191135054
#> [47,] 0.0070510247 0.339936240 -0.215660176 -0.202270043 -0.150606032
#> [48,] -0.1939965206 0.334817856 -0.095900041 0.288311497 -0.258437472
#> [49,] 0.0172618638 -0.176193502 -0.289336137 0.086308545 -0.125265960
#> [50,] -0.2594380857 -0.061810135 0.188894283 0.412231821 -0.315673331
#> [51,] 0.2592854332 -0.054411829 0.075324259 0.021315835 0.211293029
#> [52,] 0.0449375164 -0.233382911 -0.031949867 -0.843878300 0.080491983
#> [53,] 0.0874539944 0.116820712 -0.058167031 0.115272545 -0.380116057
#> [54,] -0.2056009760 -0.432731728 0.038136142 0.232318697 -0.073828605
#> [55,] -0.0355122828 -0.179369598 -0.228635531 0.012359695 -0.059577505
#> [56,] 0.1532658687 0.306487067 0.770451072 -0.733076046 0.590510479
#> [57,] -0.3928205541 0.292812326 -0.021120457 -0.160528553 0.135058144
#> [58,] 0.4466652448 0.254694259 -0.250682967 -0.037018717 -0.611712455
#> [59,] -0.1594395643 0.044090082 -0.173840979 -0.219341187 0.190444076
#> [60,] 0.1817297170 0.236726281 0.337898608 0.331766578 0.076412823
#> [61,] 0.2727139418 -0.012716707 0.106912796 0.024947834 -0.032353838
#> [62,] -0.6897258823 -0.324831717 0.108435400 0.035004568 -0.291684712
#> [63,] -0.1770494839 0.055883012 -0.070065786 -0.087216390 -0.078540769
#> [64,] -0.0631795611 -0.096487292 0.058892649 0.086868857 0.228468395
#> [65,] -0.3705698490 -0.009465416 0.204554390 -0.097344226 -0.109352641
#> [66,] -0.7164747403 0.094197103 -0.100050951 0.036431516 -0.227853278
#> [67,] 0.0817635333 -0.018733244 0.428717114 0.244464277 -0.138709386
#> [68,] 0.2724349460 0.141130282 -0.075127829 0.039483747 0.149668453
#> [69,] -0.0659565565 -0.202052472 0.171284135 0.185663837 0.065886137
#> [70,] -0.4237009593 -0.411551633 -0.830976686 0.399866542 0.140388109
#> [71,] 0.3617910991 0.369864221 -0.200138439 -0.019184583 -0.363277377
#> [72,] 0.0009644288 0.317020790 0.020923623 -0.084949793 -0.090057873
#> [73,] -0.1565175605 -2.187194337 -0.323979144 0.051364162 -0.191720340
#> [74,] -0.2147188097 0.968993066 -0.394521616 0.447211449 0.517589343
#> [75,] 0.1176618066 -0.181010845 0.346515551 0.043034665 -0.049500501
#> [76,] 0.2957110099 0.207529294 0.013925144 0.427481222 -0.505042521
#> [77,] -0.4592205233 0.455723260 -0.308032753 -0.314869650 -0.452676346
#> [78,] 0.0602985740 -0.090332568 -0.460301454 0.100296124 -0.157818881
#> [79,] 0.4155500305 -0.492572007 -0.321668655 -0.182259051 -0.034729687
#> [80,] -0.9660313344 0.092002950 -0.501336661 0.305617663 -0.166487994
#> [81,] -0.1240598016 0.024665103 -0.885718585 -0.399502010 -0.168867232
#> [82,] -0.0335121833 0.519769689 -0.773172617 0.582306358 0.031319478
#> [83,] -0.0600854641 -0.268652726 -0.048556436 0.183298659 -0.257411071
#> [84,] 0.4118478557 -0.076813221 0.171620384 -0.522752493 -0.388878617
#> [85,] 0.1004374689 0.095970843 0.133387604 0.041416046 -0.292599274
#> [86,] -0.1331831871 -0.079646144 0.758173477 0.200967719 0.099787271
#> [87,] -0.0423677222 0.268038429 -0.071506221 0.266059779 -0.448278787
#> [88,] -0.5055875092 0.308970888 -0.382332009 -0.134452903 0.165640712
#> [89,] -0.3869736138 -0.165452687 0.207904279 0.295123996 -0.195371633
#> [90,] -0.2183914331 -0.298787367 0.448240324 -0.133528592 -0.357831679
#> [91,] 0.3258080498 0.434638587 0.541499734 0.133174438 0.726353936
#> [92,] -0.1009565750 -0.361110662 0.453477377 -0.302787257 -0.092732110
#> [93,] -0.2091130163 -0.062193785 0.854567797 1.252197921 0.237930779
#> [94,] 0.4531727762 0.234990628 -0.197789206 -0.149451490 -0.063951215
#> [95,] -0.0779829632 0.024509581 0.009517371 0.115410983 -0.026071212
#> [96,] -0.2259513981 0.443588806 0.336673079 0.634519716 -0.386170314
#> [97,] -0.9963902668 0.161206413 0.071057732 -0.105762544 0.273122691
#> [98,] 0.0353323918 -0.059380730 0.127655377 0.208348778 -0.056791752
#> [99,] -0.2968287439 -0.012681463 -0.206885586 -0.251720110 0.556237895
#> [100,] -0.5021861944 -0.247749187 0.480515033 0.131095213 0.067244273
#>
#> $penal_lik
#> [1] -100578.5
#>