R : Copyright 2005, The R Foundation for Statistical Computing Version 2.1.1 (2005-06-20), ISBN 3-900051-07-0 R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for a HTML browser interface to help. Type 'q()' to quit R. [Previously saved workspace restored] > # chex99.s J F Monahan, June 2001 > # IRWLS for Poisson regression using Frome smoking example > ni <- 9 # age groups > nj <- 7 # smoking rate group, nonsmoker = 1 > All <- scan("frome.dat") > All <- t(matrix(All,4,63)) > I <- matrix(All[,1],ni,nj) > J <- matrix(All[,2],ni,nj) > Nij <- matrix(All[,3],ni,nj) > y <- matrix(All[,4],ni,nj) > # beta = ( age effects (1:ni), smoke effects (2:nj) ) > # so nonsmoke effect is zero > # initial beta > bage <- rep(log(sum(y)/sum(Nij)),ni) > bsmok <- rep( 0, nj) > beta <- c( bage, bsmok[2:nj] ) > for( k in 1:12 ) { + print("betas") + print(beta) + bage <- beta[1:ni] + bsmok <- c( 0, beta[(ni-1)+(2:nj)] ) + Xb <- matrix( rep(bage,nj),ni,nj) + + + t(matrix( rep(bsmok,ni),nj,ni) ) + lam <- exp(Xb) + rll <- sum( y*log(Nij*lam) - Nij*lam ) + print(k); print(rll) + d <- y - Nij*lam + dell <- c( (d%*%rep(1,nj)), (rep(1,ni) %*% d)[2:nj] ) + print("gradient") + print(dell) + lam <- Nij*lam + H <- matrix(0,ni+nj-1,ni+nj-1) + H[1:ni,1:ni] <- diag( (lam%*%rep(1,nj))[,1],ni ) + H[(ni-1)+(2:nj),(ni-1)+(2:nj)] <- diag( (rep(1,ni)%*%lam[,2:nj])[1,], nj-1 ) + H[1:ni,(ni-1)+(2:nj)] <- lam[,2:nj] + H[(ni-1)+(2:nj),1:ni] <- t(lam[,2:nj]) + print(H) + step <- solve(H,dell) + print("step") + print(step) + beta <- beta + step + } # end of iteration loop on k [1] "betas" [1] -6.801519 -6.801519 -6.801519 -6.801519 -6.801519 -6.801519 -6.801519 [8] -6.801519 -6.801519 0.000000 0.000000 0.000000 0.000000 0.000000 [15] 0.000000 [1] 1 [1] -97.15273 [1] "gradient" [1] -33.2255309 -30.5849698 -22.0723640 0.7287428 1.6608402 24.6660997 [7] 21.4180393 19.7064750 17.7026677 -11.3062682 -4.8620753 -6.2430495 [13] 13.3752306 25.3555532 23.0221370 [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 34.225531 0.000000 0.000000 0.000000 0.000000 0.000000 0.0000000 [2,] 0.000000 33.584970 0.000000 0.000000 0.000000 0.000000 0.0000000 [3,] 0.000000 0.000000 28.072364 0.000000 0.000000 0.000000 0.0000000 [4,] 0.000000 0.000000 0.000000 24.271257 0.000000 0.000000 0.0000000 [5,] 0.000000 0.000000 0.000000 0.000000 19.339160 0.000000 0.0000000 [6,] 0.000000 0.000000 0.000000 0.000000 0.000000 13.333900 0.0000000 [7,] 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 8.5819607 [8,] 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.0000000 [9,] 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.0000000 [10,] 3.470818 3.266194 2.544451 2.240851 1.832716 1.456831 1.0309029 [11,] 3.977928 3.654312 2.831369 2.467717 2.030667 1.541350 1.0987401 [12,] 4.800871 4.686327 3.541991 2.846938 2.105177 1.483522 0.9441602 [13,] 6.319980 7.100663 6.097563 5.212343 4.054662 2.681237 1.7426373 [14,] 3.382963 4.503945 4.770845 4.746379 3.924548 2.695694 1.5669279 [15,] 0.745097 1.296691 1.648110 1.757094 1.485746 1.027567 0.6183193 [,8] [,9] [,10] [,11] [,12] [,13] [1,] 0.0000000 0.0000000 3.4708176 3.9779284 4.8008714 6.3199796 [2,] 0.0000000 0.0000000 3.2661939 3.6543116 4.6863266 7.1006633 [3,] 0.0000000 0.0000000 2.5444507 2.8313686 3.5419910 6.0975626 [4,] 0.0000000 0.0000000 2.2408515 2.4677168 2.8469378 5.2123428 [5,] 0.0000000 0.0000000 1.8327162 2.0306674 2.1051771 4.0546623 [6,] 0.0000000 0.0000000 1.4568315 1.5413499 1.4835215 2.6812372 [7,] 0.0000000 0.0000000 1.0309029 1.0987401 0.9441602 1.7426373 [8,] 5.2935250 0.0000000 0.7895804 0.7606662 0.5226800 0.9530569 [9,] 0.0000000 3.2973323 0.6739236 0.4993262 0.3113838 0.4626274 [10,] 0.7895804 0.6739236 17.3062682 0.0000000 0.0000000 0.0000000 [11,] 0.7606662 0.4993262 0.0000000 18.8620753 0.0000000 0.0000000 [12,] 0.5226800 0.3113838 0.0000000 0.0000000 21.2430495 0.0000000 [13,] 0.9530569 0.4626274 0.0000000 0.0000000 0.0000000 34.6247694 [14,] 0.7373124 0.3158322 0.0000000 0.0000000 0.0000000 0.0000000 [15,] 0.2835817 0.1156568 0.0000000 0.0000000 0.0000000 0.0000000 [,14] [,15] [1,] 3.3829629 0.7450970 [2,] 4.5039446 1.2966912 [3,] 4.7708451 1.6481101 [4,] 4.7463792 1.7570944 [5,] 3.9245483 1.4857457 [6,] 2.6956943 1.0275666 [7,] 1.5669279 0.6183193 [8,] 0.7373124 0.2835817 [9,] 0.3158322 0.1156568 [10,] 0.0000000 0.0000000 [11,] 0.0000000 0.0000000 [12,] 0.0000000 0.0000000 [13,] 0.0000000 0.0000000 [14,] 26.6444468 0.0000000 [15,] 0.0000000 8.9778630 [1] "step" [1] -1.57341224 -1.65688558 -1.65312221 -0.91554304 -0.87595847 0.89396456 [7] 1.58252046 2.97286329 4.80424147 -0.06292814 0.44498659 0.63602475 [13] 1.24989813 1.69686771 3.19472699 [1] "betas" [1] -8.37493080 -8.45840414 -8.45464077 -7.71706160 -7.67747703 -5.90755401 [7] -5.21899810 -3.82865527 -1.99727710 -0.06292814 0.44498659 0.63602475 [13] 1.24989813 1.69686771 3.19472699 [1] 2 [1] -1383.044 [1] "gradient" [1] -17.40380 -17.54200 -14.61355 -16.35215 -14.23094 -104.20655 [7] -144.08898 -335.31152 -1086.90860 -97.04768 -124.59157 -100.26508 [13] -293.06677 -347.32037 -632.20973 [,1] [,2] [,3] [,4] [,5] [,6] [1,] 18.4038013 0.0000000 0.0000000 0.0000000 0.0000000 0.000000 [2,] 0.0000000 20.5420040 0.0000000 0.0000000 0.0000000 0.000000 [3,] 0.0000000 0.0000000 20.6135544 0.0000000 0.0000000 0.000000 [4,] 0.0000000 0.0000000 0.0000000 41.3521540 0.0000000 0.000000 [5,] 0.0000000 0.0000000 0.0000000 0.0000000 35.2309433 0.000000 [6,] 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 142.206551 [7,] 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.000000 [8,] 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.000000 [9,] 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.000000 [10,] 0.6757378 0.5849737 0.4574281 0.8423032 0.7167075 3.344443 [11,] 1.2870278 1.0876384 0.8458820 1.5414706 1.3196844 5.880312 [12,] 1.8802598 1.6884116 1.2809366 2.1527061 1.6561007 6.851091 [13,] 4.5731483 4.7265746 4.0741618 7.2818469 5.8932466 22.877205 [14,] 3.8274829 4.6876703 4.9841803 10.3678194 8.9187919 35.962907 [15,] 3.7699958 6.0354910 7.7001055 17.1645381 15.0998577 61.306450 [,7] [,8] [,9] [,10] [,11] [,12] [1,] 0.000000 0.00000 0.00000 0.6757378 1.287028 1.880260 [2,] 0.000000 0.00000 0.00000 0.5849737 1.087638 1.688412 [3,] 0.000000 0.00000 0.00000 0.4574281 0.845882 1.280937 [4,] 0.000000 0.00000 0.00000 0.8423032 1.541471 2.152706 [5,] 0.000000 0.00000 0.00000 0.7167075 1.319684 1.656101 [6,] 0.000000 0.00000 0.00000 3.3444433 5.880312 6.851091 [7,] 174.088983 0.00000 0.00000 4.7115989 8.345072 8.680558 [8,] 0.000000 360.31152 0.00000 14.4932287 23.203178 19.299944 [9,] 0.000000 0.00000 1107.90860 77.2212551 95.081307 71.775074 [10,] 4.711599 14.49323 77.22126 103.0476763 0.000000 0.000000 [11,] 8.345072 23.20318 95.08131 0.0000000 138.591571 0.000000 [12,] 8.680558 19.29994 71.77507 0.0000000 0.000000 115.265082 [13,] 29.601300 65.01908 197.02021 0.0000000 0.000000 0.000000 [14,] 41.616849 78.64841 210.30626 0.0000000 0.000000 0.000000 [15,] 73.442088 135.27845 344.41275 0.0000000 0.000000 0.000000 [,13] [,14] [,15] [1,] 4.573148 3.827483 3.769996 [2,] 4.726575 4.687670 6.035491 [3,] 4.074162 4.984180 7.700105 [4,] 7.281847 10.367819 17.164538 [5,] 5.893247 8.918792 15.099858 [6,] 22.877205 35.962907 61.306450 [7,] 29.601300 41.616849 73.442088 [8,] 65.019081 78.648411 135.278452 [9,] 197.020207 210.306260 344.412750 [10,] 0.000000 0.000000 0.000000 [11,] 0.000000 0.000000 0.000000 [12,] 0.000000 0.000000 0.000000 [13,] 341.066771 0.000000 0.000000 [14,] 0.000000 399.320372 0.000000 [15,] 0.000000 0.000000 664.209729 [1] "step" [1] -0.99842883 -0.90162640 -0.75030264 -0.43351764 -0.44071204 -0.76947506 [7] -0.86524112 -0.97060882 -1.02379425 0.04808181 0.08058757 0.09332274 [13] 0.09559311 0.07067225 -0.01279945 [1] "betas" [1] -9.37335964 -9.36003054 -9.20494341 -8.15057924 -8.11818907 -6.67702907 [7] -6.08423922 -4.79926409 -3.02107135 -0.01484633 0.52557416 0.72934749 [13] 1.34549124 1.76753996 3.18192754 [1] 3 [1] -361.0926 [1] "gradient" [1] -6.155487 -5.754602 -4.156640 -2.876885 -2.549311 -30.415354 [7] -46.172342 -117.229349 -394.806647 -34.368371 -42.740379 -33.707600 [13] -97.799242 -117.194902 -227.269462 [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 7.1554869 0.0000000 0.0000000 0.000000 0.0000000 0.000000 0.000000 [2,] 0.0000000 8.7546018 0.0000000 0.000000 0.0000000 0.000000 0.000000 [3,] 0.0000000 0.0000000 10.1566400 0.000000 0.0000000 0.000000 0.000000 [4,] 0.0000000 0.0000000 0.0000000 27.876885 0.0000000 0.000000 0.000000 [5,] 0.0000000 0.0000000 0.0000000 0.000000 23.5493111 0.000000 0.000000 [6,] 0.0000000 0.0000000 0.0000000 0.000000 0.0000000 68.415354 0.000000 [7,] 0.0000000 0.0000000 0.0000000 0.000000 0.0000000 0.000000 76.172342 [8,] 0.0000000 0.0000000 0.0000000 0.000000 0.0000000 0.000000 0.000000 [9,] 0.0000000 0.0000000 0.0000000 0.000000 0.0000000 0.000000 0.000000 [10,] 0.2612448 0.2491418 0.2266482 0.572896 0.4839771 1.625649 2.081041 [11,] 0.5140135 0.4785330 0.4329684 1.083078 0.9205980 2.952707 3.807671 [12,] 0.7605631 0.7523787 0.6640563 1.531934 1.1700852 3.484260 4.011509 [13,] 1.8540383 2.1110120 2.1169058 5.193771 4.1732209 11.661106 13.710614 [14,] 1.5135393 2.0421059 2.5260039 7.212830 6.1602707 17.880044 18.801492 [15,] 1.3714185 2.4187028 3.5899272 10.984975 9.5943447 28.039397 30.522271 [,8] [,9] [,10] [,11] [,12] [,13] [1,] 0.000000 0.00000 0.2612448 0.5140135 0.7605631 1.854038 [2,] 0.000000 0.00000 0.2491418 0.4785330 0.7523787 2.111012 [3,] 0.000000 0.00000 0.2266482 0.4329684 0.6640563 2.116906 [4,] 0.000000 0.00000 0.5728960 1.0830778 1.5319337 5.193771 [5,] 0.000000 0.00000 0.4839771 0.9205980 1.1700852 4.173221 [6,] 0.000000 0.00000 1.6256488 2.9527067 3.4842599 11.661106 [7,] 0.000000 0.00000 2.0810405 3.8076707 4.0115089 13.710614 [8,] 142.229349 0.00000 5.7612505 9.5283159 8.0270420 27.103560 [9,] 0.000000 415.80665 29.1065231 37.0224950 28.3057721 77.875014 [10,] 5.761250 29.10652 40.3683708 0.0000000 0.0000000 0.000000 [11,] 9.528316 37.02249 0.0000000 56.7403790 0.0000000 0.000000 [12,] 8.027042 28.30577 0.0000000 0.0000000 48.7075999 0.000000 [13,] 27.103560 77.87501 0.0000000 0.0000000 0.0000000 145.799242 [14,] 31.978086 81.08053 0.0000000 0.0000000 0.0000000 0.000000 [15,] 50.598756 122.14967 0.0000000 0.0000000 0.0000000 0.000000 [,14] [,15] [1,] 1.513539 1.371419 [2,] 2.042106 2.418703 [3,] 2.526004 3.589927 [4,] 7.212830 10.984975 [5,] 6.160271 9.594345 [6,] 17.880044 28.039397 [7,] 18.801492 30.522271 [8,] 31.978086 50.598756 [9,] 81.080529 122.149670 [10,] 0.000000 0.000000 [11,] 0.000000 0.000000 [12,] 0.000000 0.000000 [13,] 0.000000 0.000000 [14,] 169.194902 0.000000 [15,] 0.000000 259.269462 [1] "step" [1] -0.98877590 -0.77497041 -0.51309923 -0.19977859 -0.20192808 -0.53810668 [7] -0.70170564 -0.92512131 -1.05640968 0.11951961 0.19297801 0.21616393 [13] 0.21861684 0.16500238 -0.02202024 [1] "betas" [1] -10.3621355 -10.1350010 -9.7180426 -8.3503578 -8.3201171 -7.2151358 [7] -6.7859449 -5.7243854 -4.0774810 0.1046733 0.7185522 0.9455114 [13] 1.5641081 1.9325423 3.1599073 [1] 4 [1] -29.62573 [1] "gradient" [1] -2.0418396 -1.5607090 -0.7825962 -0.2828016 -0.2500570 [6] -6.1026953 -11.7776265 -37.7171175 -140.7255368 -11.5779305 [11] -13.3769188 -10.1802501 -29.1769263 -35.9678623 -80.1476831 [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 3.0418396 0.0000000 0.0000000 0.0000000 0.0000000 0.000000 0.000000 [2,] 0.0000000 4.5607090 0.0000000 0.0000000 0.0000000 0.000000 0.000000 [3,] 0.0000000 0.0000000 6.7825962 0.0000000 0.0000000 0.000000 0.000000 [4,] 0.0000000 0.0000000 0.0000000 25.2828016 0.0000000 0.000000 0.000000 [5,] 0.0000000 0.0000000 0.0000000 0.0000000 21.2500570 0.000000 0.000000 [6,] 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 44.102695 0.000000 [7,] 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.000000 41.777627 [8,] 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.000000 0.000000 [9,] 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.000000 0.000000 [10,] 0.1095304 0.1293564 0.1529054 0.5287127 0.4456924 1.069638 1.162627 [11,] 0.2319337 0.2673966 0.3143615 1.0757373 0.9123953 2.090898 2.289401 [12,] 0.3512321 0.4302788 0.4934550 1.5572417 1.1868615 2.525181 2.468537 [13,] 0.8583074 1.2102342 1.5769193 5.2925398 4.2434512 8.472022 8.457736 [14,] 0.6640999 1.1096155 1.7834358 6.9663062 5.9369469 12.312081 10.992717 [15,] 0.4990991 1.0900687 2.1022590 8.7998053 7.6693012 16.014325 14.801533 [,8] [,9] [,10] [,11] [,12] [,13] [1,] 0.000000 0.00000 0.1095304 0.2319337 0.3512321 0.8583074 [2,] 0.000000 0.00000 0.1293564 0.2673966 0.4302788 1.2102342 [3,] 0.000000 0.00000 0.1529054 0.3143615 0.4934550 1.5769193 [4,] 0.000000 0.00000 0.5287127 1.0757373 1.5572417 5.2925398 [5,] 0.000000 0.00000 0.4456924 0.9123953 1.1868615 4.2434512 [6,] 0.000000 0.00000 1.0696376 2.0908979 2.5251814 8.4720220 [7,] 0.000000 0.00000 1.1626269 2.2894005 2.4685372 8.4577355 [8,] 62.717117 0.00000 2.5742361 4.5819501 3.9505654 13.3719683 [9,] 0.000000 161.72554 11.4052325 15.6128459 12.2168971 33.6937485 [10,] 2.574236 11.40523 17.5779305 0.0000000 0.0000000 0.0000000 [11,] 4.581950 15.61285 0.0000000 27.3769188 0.0000000 0.0000000 [12,] 3.950565 12.21690 0.0000000 0.0000000 25.1802501 0.0000000 [13,] 13.371968 33.69375 0.0000000 0.0000000 0.0000000 77.1769263 [14,] 14.953299 33.24936 0.0000000 0.0000000 0.0000000 0.0000000 [15,] 19.624631 41.54666 0.0000000 0.0000000 0.0000000 0.0000000 [,14] [,15] [1,] 0.6640999 0.4990991 [2,] 1.1096155 1.0900687 [3,] 1.7834358 2.1022590 [4,] 6.9663062 8.7998053 [5,] 5.9369469 7.6693012 [6,] 12.3120805 16.0143248 [7,] 10.9927174 14.8015327 [8,] 14.9532988 19.6246306 [9,] 33.2493612 41.5466617 [10,] 0.0000000 0.0000000 [11,] 0.0000000 0.0000000 [12,] 0.0000000 0.0000000 [13,] 0.0000000 0.0000000 [14,] 87.9678623 0.0000000 [15,] 0.0000000 112.1476831 [1] "step" [1] -0.9333161548 -0.5885469898 -0.3403509014 -0.2247660676 -0.2206723467 [6] -0.3473263005 -0.4943432510 -0.8221300065 -1.1000967466 0.2548152334 [11] 0.3779687639 0.3957738235 0.3910830439 0.3111303498 0.0005713531 [1] "betas" [1] -11.2954517 -10.7235479 -10.0583935 -8.5751239 -8.5407895 -7.5624621 [7] -7.2802881 -6.5465154 -5.1775778 0.3594885 1.0965209 1.3412852 [13] 1.9551911 2.2436727 3.1604787 [1] 5 [1] 70.45016 [1] "gradient" [1] -0.5746115 -0.2834900 -0.1316956 -0.3722847 -0.3140119 -0.9723892 [7] -1.9828243 -9.8761645 -47.7054099 -3.5422510 -3.7137607 -2.8504286 [13] -8.1835712 -10.1259765 -26.3066513 [,1] [,2] [,3] [,4] [,5] [,6] [1,] 1.57461152 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 [2,] 0.00000000 3.2834900 0.0000000 0.0000000 0.0000000 0.0000000 [3,] 0.00000000 0.0000000 6.1316956 0.0000000 0.0000000 0.0000000 [4,] 0.00000000 0.0000000 0.0000000 25.3722847 0.0000000 0.0000000 [5,] 0.00000000 0.0000000 0.0000000 0.0000000 21.3140119 0.0000000 [6,] 0.00000000 0.0000000 0.0000000 0.0000000 0.0000000 38.9723892 [7,] 0.00000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 [8,] 0.00000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 [9,] 0.00000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 [10,] 0.05557319 0.0926510 0.1403703 0.5448412 0.4611724 0.9751235 [11,] 0.13310047 0.2166221 0.3264123 1.2538376 1.0678150 2.1559599 [12,] 0.20518353 0.3548376 0.5215758 1.8476672 1.4139877 2.6505321 [13,] 0.49906146 0.9933720 1.6589838 6.2502110 5.0318494 8.8509602 [14,] 0.35646881 0.8407986 1.7320768 7.5946907 6.4990306 11.8744020 [15,] 0.19638165 0.6054793 1.4966554 7.0324481 6.1541368 11.3217842 [,7] [,8] [,9] [,10] [,11] [,12] [1,] 0.0000000 0.000000 0.000000 0.05557319 0.1331005 0.2051835 [2,] 0.0000000 0.000000 0.000000 0.09265100 0.2166221 0.3548376 [3,] 0.0000000 0.000000 0.000000 0.14037025 0.3264123 0.5215758 [4,] 0.0000000 0.000000 0.000000 0.54484120 1.2538376 1.8476672 [5,] 0.0000000 0.000000 0.000000 0.46117243 1.0678150 1.4139877 [6,] 0.0000000 0.000000 0.000000 0.97512353 2.1559599 2.6505321 [7,] 31.9828243 0.000000 0.000000 0.91498644 2.0378912 2.2368225 [8,] 0.0000000 34.876165 0.000000 1.45971042 2.9386886 2.5792600 [9,] 0.0000000 0.000000 68.705410 4.89782254 7.5834334 6.0405621 [10,] 0.9149864 1.459710 4.897823 9.54225100 0.0000000 0.0000000 [11,] 2.0378912 2.938689 7.583433 0.00000000 17.7137607 0.0000000 [12,] 2.2368225 2.579260 6.040562 0.00000000 0.0000000 17.8504286 [13,] 7.6279662 8.689485 16.581682 0.00000000 0.0000000 0.0000000 [14,] 9.1524362 8.970418 15.105655 0.00000000 0.0000000 0.0000000 [15,] 9.0336712 8.629847 13.836248 0.00000000 0.0000000 0.0000000 [,13] [,14] [,15] [1,] 0.4990615 0.3564688 0.1963816 [2,] 0.9933720 0.8407986 0.6054793 [3,] 1.6589838 1.7320768 1.4966554 [4,] 6.2502110 7.5946907 7.0324481 [5,] 5.0318494 6.4990306 6.1541368 [6,] 8.8509602 11.8744020 11.3217842 [7,] 7.6279662 9.1524362 9.0336712 [8,] 8.6894847 8.9704180 8.6298471 [9,] 16.5816824 15.1056547 13.8362475 [10,] 0.0000000 0.0000000 0.0000000 [11,] 0.0000000 0.0000000 0.0000000 [12,] 0.0000000 0.0000000 0.0000000 [13,] 56.1835712 0.0000000 0.0000000 [14,] 0.0000000 62.1259765 0.0000000 [15,] 0.0000000 0.0000000 58.3066513 [1] "step" [1] -0.7382572 -0.4494917 -0.3681063 -0.3524647 -0.3486711 -0.3595152 [7] -0.3998240 -0.6272443 -1.0451198 0.3873027 0.4953802 0.4803792 [13] 0.4665362 0.4094559 0.1173394 [1] "betas" [1] -12.0337089 -11.1730396 -10.4264998 -8.9275886 -8.8894606 -7.9219773 [7] -7.6801121 -7.1737597 -6.2226976 0.7467912 1.5919011 1.8216645 [13] 2.4217274 2.6531286 3.2778180 [1] 6 [1] 95.7756 [1] "gradient" [1] -0.10609207 -0.04709406 -0.07474224 -0.31141834 -0.26378365 [6] -0.48152237 -0.43247063 -1.60708798 -13.75437276 -0.90426859 [11] -1.04891082 -0.89747669 -2.58969754 -2.75407153 -6.49790232 [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 1.1060921 0.00000000 0.0000000 0.0000000 0.0000000 0.000000 0.000000 [2,] 0.0000000 3.04709406 0.0000000 0.0000000 0.0000000 0.000000 0.000000 [3,] 0.0000000 0.00000000 6.0747422 0.0000000 0.0000000 0.000000 0.000000 [4,] 0.0000000 0.00000000 0.0000000 25.3114183 0.0000000 0.000000 0.000000 [5,] 0.0000000 0.00000000 0.0000000 0.0000000 21.2637837 0.000000 0.000000 [6,] 0.0000000 0.00000000 0.0000000 0.0000000 0.0000000 38.481522 0.000000 [7,] 0.0000000 0.00000000 0.0000000 0.0000000 0.0000000 0.000000 30.432471 [8,] 0.0000000 0.00000000 0.0000000 0.0000000 0.0000000 0.000000 0.000000 [9,] 0.0000000 0.00000000 0.0000000 0.0000000 0.0000000 0.000000 0.000000 [10,] 0.0391244 0.08706464 0.1430909 0.5641569 0.4793369 1.002600 0.903601 [11,] 0.1043998 0.22679422 0.3707157 1.4464676 1.2365480 2.469711 2.242232 [12,] 0.1585432 0.36596866 0.5835485 2.0997923 1.6130424 2.991050 2.424466 [13,] 0.3803182 1.01044871 1.8305844 7.0054385 5.6612956 9.850744 8.154203 [14,] 0.2565813 0.80780155 1.8051990 8.0400929 6.9063265 12.482476 9.241017 [15,] 0.1055455 0.43435755 1.1647029 5.5589518 4.8831613 8.886674 6.810559 [,8] [,9] [,10] [,11] [,12] [,13] [1,] 0.000000 0.000000 0.03912440 0.1043998 0.1585432 0.3803182 [2,] 0.000000 0.000000 0.08706464 0.2267942 0.3659687 1.0104487 [3,] 0.000000 0.000000 0.14309090 0.3707157 0.5835485 1.8305844 [4,] 0.000000 0.000000 0.56415690 1.4464676 2.0997923 7.0054385 [5,] 0.000000 0.000000 0.47933689 1.2365480 1.6130424 5.6612956 [6,] 0.000000 0.000000 1.00259978 2.4697108 2.9910498 9.8507440 [7,] 0.000000 0.000000 0.90360104 2.2422317 2.4244665 8.1542027 [8,] 26.607088 0.000000 1.14831595 2.5756432 2.2269600 7.3994491 [9,] 0.000000 34.754373 2.53697809 4.3763998 3.4341055 9.2972164 [10,] 1.148316 2.536978 6.90426859 0.0000000 0.0000000 0.0000000 [11,] 2.575643 4.376400 0.00000000 15.0489108 0.0000000 0.0000000 [12,] 2.226960 3.434105 0.00000000 0.0000000 15.8974767 0.0000000 [13,] 7.399449 9.297216 0.00000000 0.0000000 0.0000000 50.5896975 [14,] 7.214868 7.999710 0.00000000 0.0000000 0.0000000 0.0000000 [15,] 5.182678 5.471273 0.00000000 0.0000000 0.0000000 0.0000000 [,14] [,15] [1,] 0.2565813 0.1055455 [2,] 0.8078015 0.4343575 [3,] 1.8051990 1.1647029 [4,] 8.0400929 5.5589518 [5,] 6.9063265 4.8831613 [6,] 12.4824755 8.8866739 [7,] 9.2410173 6.8105587 [8,] 7.2148678 5.1826776 [9,] 7.9997098 5.4712731 [10,] 0.0000000 0.0000000 [11,] 0.0000000 0.0000000 [12,] 0.0000000 0.0000000 [13,] 0.0000000 0.0000000 [14,] 54.7540715 0.0000000 [15,] 0.0000000 38.4979023 [1] "step" [1] -0.4077980 -0.3261214 -0.3183460 -0.3160203 -0.3149133 -0.3156789 [7] -0.3185722 -0.3659322 -0.7023772 0.3362196 0.3682993 0.3514631 [13] 0.3443307 0.3296149 0.2095304 [1] "betas" [1] -12.441507 -11.499161 -10.744846 -9.243609 -9.204374 -8.237656 [7] -7.998684 -7.539692 -6.925075 1.083011 1.960200 2.173128 [13] 2.766058 2.982744 3.487348 [1] 7 [1] 99.49502 [1] "gradient" [1] -0.008373134 -0.008826940 -0.016291261 -0.064354329 -0.054074893 [6] -0.095171049 -0.076672769 -0.118470179 -2.474077739 -0.152287090 [11] -0.230018988 -0.195342968 -0.546125461 -0.501576760 -0.783203307 [,1] [,2] [,3] [,4] [,5] [,6] [1,] 1.00837313 0.00000000 0.0000000 0.0000000 0.0000000 0.000000 [2,] 0.00000000 3.00882694 0.0000000 0.0000000 0.0000000 0.000000 [3,] 0.00000000 0.00000000 6.0162913 0.0000000 0.0000000 0.000000 [4,] 0.00000000 0.00000000 0.0000000 25.0643543 0.0000000 0.000000 [5,] 0.00000000 0.00000000 0.0000000 0.0000000 21.0540749 0.000000 [6,] 0.00000000 0.00000000 0.0000000 0.0000000 0.0000000 38.095171 [7,] 0.00000000 0.00000000 0.0000000 0.0000000 0.0000000 0.000000 [8,] 0.00000000 0.00000000 0.0000000 0.0000000 0.0000000 0.000000 [9,] 0.00000000 0.00000000 0.0000000 0.0000000 0.0000000 0.000000 [10,] 0.03642182 0.08794829 0.1456714 0.5756683 0.4896594 1.023407 [11,] 0.10035647 0.23656451 0.3897045 1.5240989 1.3043562 2.603148 [12,] 0.14985856 0.37536143 0.6031974 2.1755492 1.6730894 3.100020 [13,] 0.35693048 1.02901682 1.8787750 7.2065990 5.8303101 10.137068 [14,] 0.23728512 0.81062857 1.8256567 8.1501410 7.0086111 12.657650 [15,] 0.08656318 0.38655612 1.0446171 4.9974093 4.3947470 7.991707 [,7] [,8] [,9] [,10] [,11] [,12] [1,] 0.0000000 0.000000 0.000000 0.03642182 0.1003565 0.1498586 [2,] 0.0000000 0.000000 0.000000 0.08794829 0.2365645 0.3753614 [3,] 0.0000000 0.000000 0.000000 0.14567143 0.3897045 0.6031974 [4,] 0.0000000 0.000000 0.000000 0.57566830 1.5240989 2.1755492 [5,] 0.0000000 0.000000 0.000000 0.48965936 1.3043562 1.6730894 [6,] 0.0000000 0.000000 0.000000 1.02340680 2.6031478 3.1000201 [7,] 30.0766728 0.000000 0.000000 0.91968870 2.3565500 2.5055351 [8,] 0.0000000 25.118470 0.000000 1.11469835 2.5817472 2.1949698 [9,] 0.0000000 0.000000 23.474078 1.75912405 3.1334934 2.4177620 [10,] 0.9196887 1.114698 1.759124 6.15228709 0.0000000 0.0000000 [11,] 2.3565500 2.581747 3.133493 0.00000000 14.2300190 0.0000000 [12,] 2.5055351 2.194970 2.417762 0.00000000 0.0000000 15.1953430 [13,] 8.3669709 7.241324 6.499131 0.00000000 0.0000000 0.0000000 [14,] 9.3436285 6.957544 5.510432 0.00000000 0.0000000 0.0000000 [15,] 6.1069797 4.432306 3.342318 0.00000000 0.0000000 0.0000000 [,13] [,14] [,15] [1,] 0.3569305 0.2372851 0.08656318 [2,] 1.0290168 0.8106286 0.38655612 [3,] 1.8787750 1.8256567 1.04461714 [4,] 7.2065990 8.1501410 4.99740928 [5,] 5.8303101 7.0086111 4.39474699 [6,] 10.1370679 12.6576495 7.99170678 [7,] 8.3669709 9.3436285 6.10697972 [8,] 7.2413239 6.9575444 4.43230624 [9,] 6.4991314 5.5104318 3.34231786 [10,] 0.0000000 0.0000000 0.00000000 [11,] 0.0000000 0.0000000 0.00000000 [12,] 0.0000000 0.0000000 0.00000000 [13,] 48.5461255 0.0000000 0.00000000 [14,] 0.0000000 52.5015768 0.00000000 [15,] 0.0000000 0.0000000 32.78320331 [1] "step" [1] -0.1257074 -0.1211118 -0.1207870 -0.1206744 -0.1205843 -0.1207180 [7] -0.1208813 -0.1226766 -0.2228144 0.1255589 0.1274447 0.1244579 [13] 0.1234820 0.1221752 0.1075272 [1] "betas" [1] -12.567214 -11.620273 -10.865633 -9.364283 -9.324958 -8.358374 [7] -8.119566 -7.662369 -7.147889 1.208570 2.087645 2.297585 [13] 2.889540 3.104919 3.594876 [1] 8 [1] 99.66071 [1] "gradient" [1] -0.0003273174 -0.0006297096 -0.0010163365 -0.0035679725 -0.0029241751 [6] -0.0049134473 -0.0040062506 -0.0048527891 -0.1313213264 -0.0080983884 [11] -0.0140263734 -0.0113962349 -0.0311603738 -0.0270363436 -0.0240666175 [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 1.00032732 0.00000000 0.0000000 0.000000 0.0000000 0.000000 0.0000000 [2,] 0.00000000 3.00062971 0.0000000 0.000000 0.0000000 0.000000 0.0000000 [3,] 0.00000000 0.00000000 6.0010163 0.000000 0.0000000 0.000000 0.0000000 [4,] 0.00000000 0.00000000 0.0000000 25.003568 0.0000000 0.000000 0.0000000 [5,] 0.00000000 0.00000000 0.0000000 0.000000 21.0029242 0.000000 0.0000000 [6,] 0.00000000 0.00000000 0.0000000 0.000000 0.0000000 38.004913 0.0000000 [7,] 0.00000000 0.00000000 0.0000000 0.000000 0.0000000 0.000000 30.0040063 [8,] 0.00000000 0.00000000 0.0000000 0.000000 0.0000000 0.000000 0.0000000 [9,] 0.00000000 0.00000000 0.0000000 0.000000 0.0000000 0.000000 0.0000000 [10,] 0.03641641 0.08834027 0.1463682 0.578487 0.4921013 1.028373 0.9240006 [11,] 0.10053097 0.23806740 0.3923076 1.534453 1.3133354 2.620717 2.3720677 [12,] 0.14967143 0.37661953 0.6054158 2.183796 1.6795828 3.111636 2.5145122 [13,] 0.35613707 1.03145876 1.8838451 7.226861 5.8472294 10.165126 8.3887592 [14,] 0.23644847 0.81149111 1.8281929 8.162382 7.0197703 12.676109 9.3557261 [15,] 0.08500367 0.38134044 1.0308572 4.932138 4.3377375 7.886983 6.0259686 [,8] [,9] [,10] [,11] [,12] [,13] [1,] 0.000000 0.000000 0.03641641 0.1005310 0.1496714 0.3561371 [2,] 0.000000 0.000000 0.08834027 0.2380674 0.3766195 1.0314588 [3,] 0.000000 0.000000 0.14636822 0.3923076 0.6054158 1.8838451 [4,] 0.000000 0.000000 0.57848701 1.5344525 2.1837959 7.2268611 [5,] 0.000000 0.000000 0.49210127 1.3133354 1.6795828 5.8472294 [6,] 0.000000 0.000000 1.02837302 2.6207175 3.1116356 10.1651263 [7,] 0.000000 0.000000 0.92400065 2.3720677 2.5145122 8.3887592 [8,] 25.004853 0.000000 1.11791585 2.5940865 2.1988831 7.2471588 [9,] 0.000000 21.131321 1.59609569 2.8484607 2.1912799 5.8845847 [10,] 1.117916 1.596096 6.00809839 0.0000000 0.0000000 0.0000000 [11,] 2.594087 2.848461 0.00000000 14.0140264 0.0000000 0.0000000 [12,] 2.198883 2.191280 0.00000000 0.0000000 15.0113962 0.0000000 [13,] 7.247159 5.884585 0.00000000 0.0000000 0.0000000 48.0311604 [14,] 6.954057 4.982859 0.00000000 0.0000000 0.0000000 0.0000000 [15,] 4.365666 2.978374 0.00000000 0.0000000 0.0000000 0.0000000 [,14] [,15] [1,] 0.2364485 0.08500367 [2,] 0.8114911 0.38134044 [3,] 1.8281929 1.03085717 [4,] 8.1623824 4.93213753 [5,] 7.0197703 4.33773750 [6,] 12.6761088 7.88698252 [7,] 9.3557261 6.02596855 [8,] 6.9540572 4.36566566 [9,] 4.9828591 2.97837357 [10,] 0.0000000 0.00000000 [11,] 0.0000000 0.00000000 [12,] 0.0000000 0.00000000 [13,] 0.0000000 0.00000000 [14,] 52.0270363 0.00000000 [15,] 0.0000000 32.02406662 [1] "step" [1] -0.01115782 -0.01115464 -0.01115855 -0.01116341 -0.01116280 -0.01116902 [7] -0.01117180 -0.01116746 -0.01709886 0.01139497 0.01137179 0.01127329 [13] 0.01124462 0.01121522 0.01096707 [1] "betas" [1] -12.578372 -11.631427 -10.876791 -9.375447 -9.336121 -8.369543 [7] -8.130738 -7.673536 -7.164988 1.219965 2.099017 2.308859 [13] 2.900785 3.116134 3.605843 [1] 9 [1] 99.66131 [1] "gradient" [1] -2.247634e-06 -4.567245e-06 -7.118566e-06 -2.412055e-05 -1.960370e-05 [6] -3.236381e-05 -2.653810e-05 -3.295757e-05 -4.466812e-04 -2.602936e-05 [11] -4.685907e-05 -3.718429e-05 -1.007712e-04 -8.613333e-05 -5.645754e-05 [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 1.00000225 0.0000000 0.0000000 0.000000 0.0000000 0.000000 0.0000000 [2,] 0.00000000 3.0000046 0.0000000 0.000000 0.0000000 0.000000 0.0000000 [3,] 0.00000000 0.0000000 6.0000071 0.000000 0.0000000 0.000000 0.0000000 [4,] 0.00000000 0.0000000 0.0000000 25.000024 0.0000000 0.000000 0.0000000 [5,] 0.00000000 0.0000000 0.0000000 0.000000 21.0000196 0.000000 0.0000000 [6,] 0.00000000 0.0000000 0.0000000 0.000000 0.0000000 38.000032 0.0000000 [7,] 0.00000000 0.0000000 0.0000000 0.000000 0.0000000 0.000000 30.0000265 [8,] 0.00000000 0.0000000 0.0000000 0.000000 0.0000000 0.000000 0.0000000 [9,] 0.00000000 0.0000000 0.0000000 0.000000 0.0000000 0.000000 0.0000000 [10,] 0.03642504 0.0883615 0.1464028 0.578621 0.4922155 1.028605 0.9242069 [11,] 0.10055248 0.2381191 0.3923913 1.534772 1.3136099 2.621249 2.3725421 [12,] 0.14968871 0.3766642 0.6054852 2.184036 1.6797684 3.111960 2.5147674 [13,] 0.35616798 1.0315516 1.8840073 7.227448 5.8477078 10.165895 8.3893700 [14,] 0.23646204 0.8115403 1.8282965 8.162805 7.0201383 12.676694 9.3561323 [15,] 0.08498745 0.3812689 1.0306598 4.931169 4.3368886 7.885390 6.0247350 [,8] [,9] [,10] [,11] [,12] [,13] [,14] [1,] 0.000000 0.000000 0.03642504 0.1005525 0.1496887 0.356168 0.2364620 [2,] 0.000000 0.000000 0.08836150 0.2381191 0.3766642 1.031552 0.8115403 [3,] 0.000000 0.000000 0.14640283 0.3923913 0.6054852 1.884007 1.8282965 [4,] 0.000000 0.000000 0.57862098 1.5347723 2.1840359 7.227448 8.1628053 [5,] 0.000000 0.000000 0.49221554 1.3136099 1.6797684 5.847708 7.0201383 [6,] 0.000000 0.000000 1.02860540 2.6212489 3.1119600 10.165895 12.6766944 [7,] 0.000000 0.000000 0.92420688 2.3725421 2.5147674 8.389370 9.3561323 [8,] 25.000033 0.000000 1.11817021 2.5946166 2.1991158 7.247718 6.9543893 [9,] 0.000000 21.000447 1.58701765 2.8321940 2.1785515 5.850236 4.9536278 [10,] 1.118170 1.587018 6.00002603 0.0000000 0.0000000 0.000000 0.0000000 [11,] 2.594617 2.832194 0.00000000 14.0000469 0.0000000 0.000000 0.0000000 [12,] 2.199116 2.178552 0.00000000 0.0000000 15.0000372 0.000000 0.0000000 [13,] 7.247718 5.850236 0.00000000 0.0000000 0.0000000 48.000101 0.0000000 [14,] 6.954389 4.953628 0.00000000 0.0000000 0.0000000 0.000000 52.0000861 [15,] 4.364791 2.960167 0.00000000 0.0000000 0.0000000 0.000000 0.0000000 [,15] [1,] 0.08498745 [2,] 0.38126892 [3,] 1.03065980 [4,] 4.93116925 [5,] 4.33688857 [6,] 7.88538990 [7,] 6.02473498 [8,] 4.36479090 [9,] 2.96016668 [10,] 0.00000000 [11,] 0.00000000 [12,] 0.00000000 [13,] 0.00000000 [14,] 0.00000000 [15,] 32.00005646 [1] "step" [1] -7.659084e-05 -7.673296e-05 -7.679047e-05 -7.683239e-05 -7.683700e-05 [6] -7.686055e-05 -7.686166e-05 -7.678560e-05 -9.603630e-05 7.757069e-05 [11] 7.736765e-05 7.714036e-05 7.707256e-05 7.700775e-05 7.684807e-05 [1] "betas" [1] -12.578449 -11.631504 -10.876868 -9.375524 -9.336198 -8.369620 [7] -8.130814 -7.673613 -7.165084 1.220042 2.099094 2.308936 [13] 2.900862 3.116211 3.605920 [1] 10 [1] 99.66131 [1] "gradient" [1] -1.048975e-10 -2.136322e-10 -3.327318e-10 -1.125863e-09 -9.148153e-10 [6] -1.508132e-09 -1.236406e-09 -1.537950e-09 -6.591769e-09 -2.717868e-10 [11] -4.951666e-10 -3.895466e-10 -1.053174e-09 -8.975061e-10 -5.449610e-10 [,1] [,2] [,3] [,4] [,5] [,6] [1,] 1.00000000 0.00000000 0.0000000 0.0000000 0.0000000 0.000000 [2,] 0.00000000 3.00000000 0.0000000 0.0000000 0.0000000 0.000000 [3,] 0.00000000 0.00000000 6.0000000 0.0000000 0.0000000 0.000000 [4,] 0.00000000 0.00000000 0.0000000 25.0000000 0.0000000 0.000000 [5,] 0.00000000 0.00000000 0.0000000 0.0000000 21.0000000 0.000000 [6,] 0.00000000 0.00000000 0.0000000 0.0000000 0.0000000 38.000000 [7,] 0.00000000 0.00000000 0.0000000 0.0000000 0.0000000 0.000000 [8,] 0.00000000 0.00000000 0.0000000 0.0000000 0.0000000 0.000000 [9,] 0.00000000 0.00000000 0.0000000 0.0000000 0.0000000 0.000000 [10,] 0.03642508 0.08836158 0.1464029 0.5786214 0.4922159 1.028606 [11,] 0.10055256 0.23811926 0.3923915 1.5347731 1.3136106 2.621250 [12,] 0.14968879 0.37666437 0.6054854 2.1840365 1.6797689 3.111961 [13,] 0.35616815 1.03155193 1.8840078 7.2274497 5.8477092 10.165897 [14,] 0.23646214 0.81154050 1.8282969 8.1628067 7.0201395 12.676696 [15,] 0.08498748 0.38126896 1.0306599 4.9311693 4.3368886 7.885390 [,7] [,8] [,9] [,10] [,11] [,12] [1,] 0.0000000 0.000000 0.000000 0.03642508 0.1005526 0.1496888 [2,] 0.0000000 0.000000 0.000000 0.08836158 0.2381193 0.3766644 [3,] 0.0000000 0.000000 0.000000 0.14640294 0.3923915 0.6054854 [4,] 0.0000000 0.000000 0.000000 0.57862141 1.5347731 2.1840365 [5,] 0.0000000 0.000000 0.000000 0.49221590 1.3136106 1.6797689 [6,] 0.0000000 0.000000 0.000000 1.02860614 2.6212503 3.1119609 [7,] 30.0000000 0.000000 0.000000 0.92420753 2.3725433 2.5147681 [8,] 0.0000000 25.000000 0.000000 1.11817109 2.5946182 2.1991166 [9,] 0.0000000 0.000000 21.000000 1.58698834 2.8321412 2.1785104 [10,] 0.9242075 1.118171 1.586988 6.00000000 0.0000000 0.0000000 [11,] 2.3725433 2.594618 2.832141 0.00000000 14.0000000 0.0000000 [12,] 2.5147681 2.199117 2.178510 0.00000000 0.0000000 15.0000000 [13,] 8.3893718 7.247720 5.850125 0.00000000 0.0000000 0.0000000 [14,] 9.3561337 6.954391 4.953534 0.00000000 0.0000000 0.0000000 [15,] 6.0247349 4.364791 2.960110 0.00000000 0.0000000 0.0000000 [,13] [,14] [,15] [1,] 0.3561682 0.2364621 0.08498748 [2,] 1.0315519 0.8115405 0.38126896 [3,] 1.8840078 1.8282969 1.03065986 [4,] 7.2274497 8.1628067 4.93116933 [5,] 5.8477092 7.0201395 4.33688862 [6,] 10.1658969 12.6766963 7.88538980 [7,] 8.3893718 9.3561337 6.02473490 [8,] 7.2477200 6.9543908 4.36479117 [9,] 5.8501246 4.9535335 2.96010988 [10,] 0.0000000 0.0000000 0.00000000 [11,] 0.0000000 0.0000000 0.00000000 [12,] 0.0000000 0.0000000 0.00000000 [13,] 48.0000000 0.0000000 0.00000000 [14,] 0.0000000 52.0000000 0.00000000 [15,] 0.0000000 0.0000000 32.00000000 [1] "step" [1] -3.251675e-09 -3.255407e-09 -3.256951e-09 -3.258000e-09 -3.258136e-09 [6] -3.258625e-09 -3.258579e-09 -3.256645e-09 -3.478315e-09 3.270828e-09 [11] 3.267057e-09 3.263900e-09 3.262812e-09 3.261715e-09 3.261365e-09 [1] "betas" [1] -12.578449 -11.631504 -10.876868 -9.375524 -9.336198 -8.369620 [7] -8.130814 -7.673613 -7.165084 1.220042 2.099094 2.308936 [13] 2.900862 3.116211 3.605920 [1] 11 [1] 99.66131 [1] "gradient" [1] -3.469447e-16 -4.773959e-15 -4.884981e-15 -7.549517e-15 -2.398082e-14 [6] 8.881784e-15 8.437695e-15 -9.769963e-15 -8.881784e-15 3.108624e-15 [11] 8.881784e-15 -3.996803e-15 -1.731948e-14 -1.953993e-14 -1.465494e-14 [,1] [,2] [,3] [,4] [,5] [,6] [1,] 1.00000000 0.00000000 0.0000000 0.0000000 0.0000000 0.000000 [2,] 0.00000000 3.00000000 0.0000000 0.0000000 0.0000000 0.000000 [3,] 0.00000000 0.00000000 6.0000000 0.0000000 0.0000000 0.000000 [4,] 0.00000000 0.00000000 0.0000000 25.0000000 0.0000000 0.000000 [5,] 0.00000000 0.00000000 0.0000000 0.0000000 21.0000000 0.000000 [6,] 0.00000000 0.00000000 0.0000000 0.0000000 0.0000000 38.000000 [7,] 0.00000000 0.00000000 0.0000000 0.0000000 0.0000000 0.000000 [8,] 0.00000000 0.00000000 0.0000000 0.0000000 0.0000000 0.000000 [9,] 0.00000000 0.00000000 0.0000000 0.0000000 0.0000000 0.000000 [10,] 0.03642508 0.08836158 0.1464029 0.5786214 0.4922159 1.028606 [11,] 0.10055256 0.23811926 0.3923915 1.5347731 1.3136106 2.621250 [12,] 0.14968879 0.37666437 0.6054854 2.1840365 1.6797689 3.111961 [13,] 0.35616815 1.03155193 1.8840078 7.2274497 5.8477092 10.165897 [14,] 0.23646214 0.81154050 1.8282969 8.1628067 7.0201395 12.676696 [15,] 0.08498748 0.38126896 1.0306599 4.9311693 4.3368886 7.885390 [,7] [,8] [,9] [,10] [,11] [,12] [1,] 0.0000000 0.000000 0.000000 0.03642508 0.1005526 0.1496888 [2,] 0.0000000 0.000000 0.000000 0.08836158 0.2381193 0.3766644 [3,] 0.0000000 0.000000 0.000000 0.14640294 0.3923915 0.6054854 [4,] 0.0000000 0.000000 0.000000 0.57862141 1.5347731 2.1840365 [5,] 0.0000000 0.000000 0.000000 0.49221590 1.3136106 1.6797689 [6,] 0.0000000 0.000000 0.000000 1.02860614 2.6212503 3.1119609 [7,] 30.0000000 0.000000 0.000000 0.92420753 2.3725433 2.5147681 [8,] 0.0000000 25.000000 0.000000 1.11817109 2.5946182 2.1991166 [9,] 0.0000000 0.000000 21.000000 1.58698834 2.8321412 2.1785104 [10,] 0.9242075 1.118171 1.586988 6.00000000 0.0000000 0.0000000 [11,] 2.3725433 2.594618 2.832141 0.00000000 14.0000000 0.0000000 [12,] 2.5147681 2.199117 2.178510 0.00000000 0.0000000 15.0000000 [13,] 8.3893718 7.247720 5.850125 0.00000000 0.0000000 0.0000000 [14,] 9.3561337 6.954391 4.953534 0.00000000 0.0000000 0.0000000 [15,] 6.0247349 4.364791 2.960110 0.00000000 0.0000000 0.0000000 [,13] [,14] [,15] [1,] 0.3561682 0.2364621 0.08498748 [2,] 1.0315519 0.8115405 0.38126896 [3,] 1.8840078 1.8282969 1.03065986 [4,] 7.2274497 8.1628067 4.93116933 [5,] 5.8477092 7.0201395 4.33688862 [6,] 10.1658969 12.6766963 7.88538980 [7,] 8.3893718 9.3561337 6.02473490 [8,] 7.2477200 6.9543908 4.36479117 [9,] 5.8501246 4.9535335 2.96010988 [10,] 0.0000000 0.0000000 0.00000000 [11,] 0.0000000 0.0000000 0.00000000 [12,] 0.0000000 0.0000000 0.00000000 [13,] 48.0000000 0.0000000 0.00000000 [14,] 0.0000000 52.0000000 0.00000000 [15,] 0.0000000 0.0000000 32.00000000 [1] "step" [1] 1.387533e-16 -1.060901e-15 -2.514931e-16 2.739929e-16 -5.646918e-16 [6] 8.020028e-16 8.336061e-16 1.142546e-16 8.806408e-18 2.694071e-16 [11] 3.670714e-16 -5.318749e-16 -6.355196e-16 -6.793956e-16 -7.742584e-16 [1] "betas" [1] -12.578449 -11.631504 -10.876868 -9.375524 -9.336198 -8.369620 [7] -8.130814 -7.673613 -7.165084 1.220042 2.099094 2.308936 [13] 2.900862 3.116211 3.605920 [1] 12 [1] 99.66131 [1] "gradient" [1] 2.636780e-16 4.496403e-15 -4.440892e-16 1.221245e-14 -8.881784e-15 [6] 4.085621e-14 2.886580e-14 5.773160e-15 2.220446e-15 3.330669e-15 [11] -3.996803e-15 1.110223e-14 2.930989e-14 2.930989e-14 1.598721e-14 [,1] [,2] [,3] [,4] [,5] [,6] [1,] 1.00000000 0.00000000 0.0000000 0.0000000 0.0000000 0.000000 [2,] 0.00000000 3.00000000 0.0000000 0.0000000 0.0000000 0.000000 [3,] 0.00000000 0.00000000 6.0000000 0.0000000 0.0000000 0.000000 [4,] 0.00000000 0.00000000 0.0000000 25.0000000 0.0000000 0.000000 [5,] 0.00000000 0.00000000 0.0000000 0.0000000 21.0000000 0.000000 [6,] 0.00000000 0.00000000 0.0000000 0.0000000 0.0000000 38.000000 [7,] 0.00000000 0.00000000 0.0000000 0.0000000 0.0000000 0.000000 [8,] 0.00000000 0.00000000 0.0000000 0.0000000 0.0000000 0.000000 [9,] 0.00000000 0.00000000 0.0000000 0.0000000 0.0000000 0.000000 [10,] 0.03642508 0.08836158 0.1464029 0.5786214 0.4922159 1.028606 [11,] 0.10055256 0.23811926 0.3923915 1.5347731 1.3136106 2.621250 [12,] 0.14968879 0.37666437 0.6054854 2.1840365 1.6797689 3.111961 [13,] 0.35616815 1.03155193 1.8840078 7.2274497 5.8477092 10.165897 [14,] 0.23646214 0.81154050 1.8282969 8.1628067 7.0201395 12.676696 [15,] 0.08498748 0.38126896 1.0306599 4.9311693 4.3368886 7.885390 [,7] [,8] [,9] [,10] [,11] [,12] [1,] 0.0000000 0.000000 0.000000 0.03642508 0.1005526 0.1496888 [2,] 0.0000000 0.000000 0.000000 0.08836158 0.2381193 0.3766644 [3,] 0.0000000 0.000000 0.000000 0.14640294 0.3923915 0.6054854 [4,] 0.0000000 0.000000 0.000000 0.57862141 1.5347731 2.1840365 [5,] 0.0000000 0.000000 0.000000 0.49221590 1.3136106 1.6797689 [6,] 0.0000000 0.000000 0.000000 1.02860614 2.6212503 3.1119609 [7,] 30.0000000 0.000000 0.000000 0.92420753 2.3725433 2.5147681 [8,] 0.0000000 25.000000 0.000000 1.11817109 2.5946182 2.1991166 [9,] 0.0000000 0.000000 21.000000 1.58698834 2.8321412 2.1785104 [10,] 0.9242075 1.118171 1.586988 6.00000000 0.0000000 0.0000000 [11,] 2.3725433 2.594618 2.832141 0.00000000 14.0000000 0.0000000 [12,] 2.5147681 2.199117 2.178510 0.00000000 0.0000000 15.0000000 [13,] 8.3893718 7.247720 5.850125 0.00000000 0.0000000 0.0000000 [14,] 9.3561337 6.954391 4.953534 0.00000000 0.0000000 0.0000000 [15,] 6.0247349 4.364791 2.960110 0.00000000 0.0000000 0.0000000 [,13] [,14] [,15] [1,] 0.3561682 0.2364621 0.08498748 [2,] 1.0315519 0.8115405 0.38126896 [3,] 1.8840078 1.8282969 1.03065986 [4,] 7.2274497 8.1628067 4.93116933 [5,] 5.8477092 7.0201395 4.33688862 [6,] 10.1658969 12.6766963 7.88538980 [7,] 8.3893718 9.3561337 6.02473490 [8,] 7.2477200 6.9543908 4.36479117 [9,] 5.8501246 4.9535335 2.96010988 [10,] 0.0000000 0.0000000 0.00000000 [11,] 0.0000000 0.0000000 0.00000000 [12,] 0.0000000 0.0000000 0.00000000 [13,] 48.0000000 0.0000000 0.00000000 [14,] 0.0000000 52.0000000 0.00000000 [15,] 0.0000000 0.0000000 32.00000000 [1] "step" [1] -7.189268e-17 1.150613e-15 -4.253449e-16 1.387076e-16 -7.692204e-16 [6] 7.337557e-16 6.272614e-16 -8.559533e-17 -1.886284e-16 4.421415e-16 [11] -4.253097e-16 5.776423e-16 4.468305e-16 3.807283e-16 3.128744e-16 > rm(list=ls()) >