R Under development (unstable) (2024-06-07 r86704 ucrt) -- "Unsuffered Consequences" Copyright (C) 2024 The R Foundation for Statistical Computing Platform: x86_64-w64-mingw32/x64 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 an HTML browser interface to help. Type 'q()' to quit R. > # This file is part of the standard setup for testthat. > # It is recommended that you do not modify it. > # > # Where should you do additional test configuration? > # Learn more about the roles of various files in: > # * https://r-pkgs.org/tests.html > # * https://testthat.r-lib.org/reference/test_package.html#special-files > > library(testthat) > library(ETASbootstrap) > > test_check("ETASbootstrap") declustering: | | | 0% | |====== | 9% | |============= | 18% | |=================== | 27% | |========================= | 36% | |================================ | 45% | |====================================== | 55% | |============================================= | 64% | |=================================================== | 73% | |========================================================= | 82% | |================================================================ | 91% | |======================================================================| 100% iteration: 1 ====================================================== background seismicity rate: Min. 1st Qu. Median Mean 3rd Qu. Max. 0.0002489 0.0084420 0.0217383 0.0338881 0.0386664 0.1432186 probability of being a background event: Min. 1st Qu. Median Mean 3rd Qu. Max. 0.0000005 0.0002743 0.0390920 0.4376626 0.9756477 1.0000000 integral of background seismicity rate: 93.28025 ====================================================== estimating: start Davidon-Fletcher-Powell procedure ... Function Value = 206.3994 Gradient[1] = 47.06 theta[1] = 0.806226 Gradient[2] = 39.25 theta[2] = 0.489898 Gradient[3] = 82.49 theta[3] = 0.082462 Gradient[4] = 8.94 theta[4] = 0.984886 Gradient[5] = -17.71 theta[5] = 1.104536 Gradient[6] = -141.76 theta[6] = 0.057446 Gradient[7] = 16.14 theta[7] = 1.574802 Gradient[8] = -5.58 theta[8] = 0.412311 line search along the specified direction ... zeta = 0.000210 Function Value = 203.8967 Gradient[1] = 44.50 theta[1] = 0.796345 Gradient[2] = 33.79 theta[2] = 0.481655 Gradient[3] = 63.10 theta[3] = 0.065141 Gradient[4] = 5.20 theta[4] = 0.983008 Gradient[5] = 0.03 theta[5] = 1.108255 Gradient[6] = 71.41 theta[6] = 0.087214 Gradient[7] = -6.25 theta[7] = 1.571412 Gradient[8] = 3.09 theta[8] = 0.413483 line search along the specified direction ... zeta = 0.000413 Function Value = 201.6793 Gradient[1] = 36.59 theta[1] = 0.777667 Gradient[2] = 29.58 theta[2] = 0.467041 Gradient[3] = -33.84 theta[3] = 0.036728 Gradient[4] = 3.62 theta[4] = 0.980405 Gradient[5] = 37.63 theta[5] = 1.110409 Gradient[6] = 62.64 theta[6] = 0.083797 Gradient[7] = -4.99 theta[7] = 1.571255 Gradient[8] = 2.31 theta[8] = 0.413269 line search along the specified direction ... zeta = 0.001588 Function Value = 199.0465 Gradient[1] = 16.91 theta[1] = 0.714314 Gradient[2] = -11.73 theta[2] = 0.415286 Gradient[3] = 5.05 theta[3] = 0.031841 Gradient[4] = -15.32 theta[4] = 0.972035 Gradient[5] = -10.26 theta[5] = 1.078368 Gradient[6] = 64.05 theta[6] = 0.083525 Gradient[7] = -5.37 theta[7] = 1.568931 Gradient[8] = 2.48 theta[8] = 0.413424 line search along the specified direction ... zeta = 0.002626 Function Value = 198.2806 Gradient[1] = -1.29 theta[1] = 0.659930 Gradient[2] = 0.34 theta[2] = 0.422837 Gradient[3] = 20.77 theta[3] = 0.039724 Gradient[4] = -10.04 theta[4] = 1.000583 Gradient[5] = 5.96 theta[5] = 1.089333 Gradient[6] = 59.30 theta[6] = 0.082945 Gradient[7] = -4.72 theta[7] = 1.565257 Gradient[8] = 2.34 theta[8] = 0.413866 line search along the specified direction ... zeta = 0.005364 Function Value = 197.8575 Gradient[1] = -3.98 theta[1] = 0.643441 Gradient[2] = 3.80 theta[2] = 0.428627 Gradient[3] = 2.59 theta[3] = 0.028420 Gradient[4] = -6.88 theta[4] = 1.055347 Gradient[5] = -5.93 theta[5] = 1.060374 Gradient[6] = 65.77 theta[6] = 0.084269 Gradient[7] = -5.71 theta[7] = 1.556844 Gradient[8] = 2.85 theta[8] = 0.414845 line search along the specified direction ... zeta = 0.027070 Function Value = 196.7963 Gradient[1] = -3.67 theta[1] = 0.648943 Gradient[2] = 18.83 theta[2] = 0.327405 Gradient[3] = 4.69 theta[3] = 0.034073 Gradient[4] = 9.14 theta[4] = 1.249255 Gradient[5] = 8.68 theta[5] = 1.073558 Gradient[6] = 62.43 theta[6] = 0.080456 Gradient[7] = -5.94 theta[7] = 1.520099 Gradient[8] = 2.58 theta[8] = 0.414775 line search along the specified direction ... zeta = 0.160322 Function Value = 196.5812 Gradient[1] = -3.50 theta[1] = 0.652233 Gradient[2] = 9.44 theta[2] = 0.314742 Gradient[3] = 0.93 theta[3] = 0.034589 Gradient[4] = 3.80 theta[4] = 1.252200 Gradient[5] = 6.54 theta[5] = 1.078586 Gradient[6] = 77.89 theta[6] = 0.064290 Gradient[7] = -10.31 theta[7] = 1.398323 Gradient[8] = 2.49 theta[8] = 0.447850 line search along the specified direction ... zeta = 0.434424 Function Value = 196.1720 Gradient[1] = -3.09 theta[1] = 0.661169 Gradient[2] = -11.87 theta[2] = 0.274447 Gradient[3] = -5.50 theta[3] = 0.035522 Gradient[4] = -5.53 theta[4] = 1.281958 Gradient[5] = 3.50 theta[5] = 1.088432 Gradient[6] = 62.71 theta[6] = 0.051806 Gradient[7] = -7.18 theta[7] = 1.394747 Gradient[8] = 2.59 theta[8] = 0.614777 line search along the specified direction ... zeta = 0.572785 Function Value = 195.9907 Gradient[1] = -1.55 theta[1] = 0.665491 Gradient[2] = -9.07 theta[2] = 0.280076 Gradient[3] = -5.98 theta[3] = 0.034461 Gradient[4] = -4.20 theta[4] = 1.274245 Gradient[5] = 1.74 theta[5] = 1.085830 Gradient[6] = 3.11 theta[6] = 0.047396 Gradient[7] = -1.48 theta[7] = 1.386456 Gradient[8] = -0.62 theta[8] = 0.575558 line search along the specified direction ... zeta = 0.529547 Function Value = 195.9286 Gradient[1] = -0.48 theta[1] = 0.667890 Gradient[2] = -4.49 theta[2] = 0.290229 Gradient[3] = -5.41 theta[3] = 0.033986 Gradient[4] = -2.58 theta[4] = 1.262960 Gradient[5] = 1.54 theta[5] = 1.083946 Gradient[6] = -11.58 theta[6] = 0.040280 Gradient[7] = -0.24 theta[7] = 1.370372 Gradient[8] = -0.62 theta[8] = 0.685481 line search along the specified direction ... zeta = 0.609987 Function Value = 195.8928 Gradient[1] = 0.48 theta[1] = 0.668948 Gradient[2] = -0.95 theta[2] = 0.296707 Gradient[3] = -1.46 theta[3] = 0.033718 Gradient[4] = -1.04 theta[4] = 1.259098 Gradient[5] = -0.32 theta[5] = 1.080705 Gradient[6] = 2.81 theta[6] = 0.044730 Gradient[7] = -0.35 theta[7] = 1.403063 Gradient[8] = 0.06 theta[8] = 0.673120 line search along the specified direction ... zeta = 0.908554 Function Value = 195.8889 Gradient[1] = 0.14 theta[1] = 0.667991 Gradient[2] = -0.16 theta[2] = 0.292519 Gradient[3] = 0.70 theta[3] = 0.034126 Gradient[4] = -0.08 theta[4] = 1.268673 Gradient[5] = -0.81 theta[5] = 1.080707 Gradient[6] = 0.71 theta[6] = 0.044290 Gradient[7] = -0.11 theta[7] = 1.403259 Gradient[8] = 0.03 theta[8] = 0.681273 line search along the specified direction ... zeta = 0.934919 Function Value = 195.8886 Gradient[1] = -0.03 theta[1] = 0.667679 Gradient[2] = 0.06 theta[2] = 0.292349 Gradient[3] = -0.33 theta[3] = 0.034144 Gradient[4] = 0.04 theta[4] = 1.268757 Gradient[5] = 0.03 theta[5] = 1.081247 Gradient[6] = 0.15 theta[6] = 0.044337 Gradient[7] = -0.02 theta[7] = 1.404523 Gradient[8] = 0.01 theta[8] = 0.681748 line search along the specified direction ... zeta = 0.867637 Function Value = 195.8886 Gradient[1] = -0.01 theta[1] = 0.667775 Gradient[2] = 0.01 theta[2] = 0.292270 Gradient[3] = 0.02 theta[3] = 0.034239 Gradient[4] = -0.00 theta[4] = 1.268684 Gradient[5] = -0.00 theta[5] = 1.081391 Gradient[6] = -0.04 theta[6] = 0.044357 Gradient[7] = 0.00 theta[7] = 1.404796 Gradient[8] = -0.00 theta[8] = 0.681470 line search along the specified direction ... zeta = 1.013941 Function Value = 195.8886 Gradient[1] = 0.00 theta[1] = 0.667813 Gradient[2] = 0.00 theta[2] = 0.292228 Gradient[3] = 0.00 theta[3] = 0.034238 Gradient[4] = 0.00 theta[4] = 1.268745 Gradient[5] = 0.00 theta[5] = 1.081400 Gradient[6] = -0.00 theta[6] = 0.044363 Gradient[7] = -0.00 theta[7] = 1.404854 Gradient[8] = 0.00 theta[8] = 0.681510 line search along the specified direction ... zeta = 0.976690 Function Value = 195.8886 Gradient[1] = 0.00 theta[1] = 0.667811 Gradient[2] = 0.00 theta[2] = 0.292225 Gradient[3] = -0.00 theta[3] = 0.034238 Gradient[4] = 0.00 theta[4] = 1.268748 Gradient[5] = 0.00 theta[5] = 1.081399 Gradient[6] = -0.00 theta[6] = 0.044366 Gradient[7] = 0.00 theta[7] = 1.404868 Gradient[8] = -0.00 theta[8] = 0.681475 line search along the specified direction ... zeta = 0.877815 Function Value = 195.8886 Gradient[1] = -0.00 theta[1] = 0.667811 Gradient[2] = -0.00 theta[2] = 0.292225 Gradient[3] = -0.00 theta[3] = 0.034238 Gradient[4] = -0.00 theta[4] = 1.268748 Gradient[5] = -0.00 theta[5] = 1.081399 Gradient[6] = -0.00 theta[6] = 0.044366 Gradient[7] = 0.00 theta[7] = 1.404869 Gradient[8] = -0.00 theta[8] = 0.681478 line search along the specified direction ... zeta = 0.877815 Function Value = 195.8886 Gradient[1] = -0.00 theta[1] = 0.667811 Gradient[2] = -0.00 theta[2] = 0.292225 Gradient[3] = -0.00 theta[3] = 0.034238 Gradient[4] = -0.00 theta[4] = 1.268748 Gradient[5] = -0.00 theta[5] = 1.081399 Gradient[6] = -0.00 theta[6] = 0.044366 Gradient[7] = 0.00 theta[7] = 1.404869 Gradient[8] = -0.00 theta[8] = 0.681478 line search along the specified direction ... zeta = 0.877815 Function Value = 195.8886 Gradient[1] = -0.00 theta[1] = 0.667811 Gradient[2] = -0.00 theta[2] = 0.292225 Gradient[3] = -0.00 theta[3] = 0.034238 Gradient[4] = -0.00 theta[4] = 1.268748 Gradient[5] = -0.00 theta[5] = 1.081399 Gradient[6] = -0.00 theta[6] = 0.044366 Gradient[7] = 0.00 theta[7] = 1.404869 Gradient[8] = -0.00 theta[8] = 0.681478 line search along the specified direction ... zeta = 0.087782 Function Value = 195.8886 Gradient[1] = -0.00 theta[1] = 0.667811 Gradient[2] = -0.00 theta[2] = 0.292225 Gradient[3] = -0.00 theta[3] = 0.034238 Gradient[4] = -0.00 theta[4] = 1.268748 Gradient[5] = -0.00 theta[5] = 1.081399 Gradient[6] = -0.00 theta[6] = 0.044366 Gradient[7] = 0.00 theta[7] = 1.404869 Gradient[8] = -0.00 theta[8] = 0.681478 loglikelihood = -195.88857 AIC = 407.77714 theta[1] = 0.44597125 gradient[1] = -0.0000 theta[2] = 0.08539571 gradient[2] = -0.0000 theta[3] = 0.00117225 gradient[3] = -0.0000 theta[4] = 1.60972087 gradient[4] = -0.0000 theta[5] = 1.16942287 gradient[5] = -0.0000 theta[6] = 0.00196834 gradient[6] = -0.0000 theta[7] = 1.97365744 gradient[7] = 0.0000 theta[8] = 0.46441259 gradient[8] = -0.0000 ====================================================== MLE: mu A c alpha p D 0.445971251 0.085395710 0.001172247 1.609720868 1.169422871 0.001968337 q gamma 1.973657444 0.464412589 ====================================================== declustering: | | | 0% | |======= | 10% | |============== | 20% | |===================== | 30% | |============================ | 40% | |=================================== | 50% | |========================================== | 60% | |================================================= | 70% | |======================================================== | 80% | |=============================================================== | 90% | |======================================================================| 100% iteration: 2 ====================================================== background seismicity rate: Min. 1st Qu. Median Mean 3rd Qu. Max. 0.0002452 0.0084954 0.0221936 0.0342122 0.0396519 0.1449802 probability of being a background event: Min. 1st Qu. Median Mean 3rd Qu. Max. 0.0000002 0.0004388 0.0393883 0.4377489 0.9720255 1.0000000 integral of background seismicity rate: 93.36126 ====================================================== estimating: start Davidon-Fletcher-Powell procedure ... Function Value = 195.9294 Gradient[1] = 0.17 theta[1] = 0.667811 Gradient[2] = -0.13 theta[2] = 0.292225 Gradient[3] = -0.20 theta[3] = 0.034238 Gradient[4] = -0.06 theta[4] = 1.268748 Gradient[5] = 0.20 theta[5] = 1.081399 Gradient[6] = -1.09 theta[6] = 0.044366 Gradient[7] = 0.16 theta[7] = 1.404869 Gradient[8] = -0.04 theta[8] = 0.681478 line search along the specified direction ... zeta = 1.007774 Function Value = 195.9291 Gradient[1] = -0.00 theta[1] = 0.667138 Gradient[2] = 0.00 theta[2] = 0.292874 Gradient[3] = 0.00 theta[3] = 0.034156 Gradient[4] = -0.00 theta[4] = 1.268317 Gradient[5] = -0.01 theta[5] = 1.080975 Gradient[6] = -0.02 theta[6] = 0.044326 Gradient[7] = 0.00 theta[7] = 1.402885 Gradient[8] = -0.00 theta[8] = 0.680411 line search along the specified direction ... zeta = 1.036246 Function Value = 195.9291 Gradient[1] = 0.00 theta[1] = 0.667142 Gradient[2] = 0.00 theta[2] = 0.292860 Gradient[3] = -0.00 theta[3] = 0.034159 Gradient[4] = 0.00 theta[4] = 1.268333 Gradient[5] = 0.00 theta[5] = 1.080987 Gradient[6] = -0.00 theta[6] = 0.044327 Gradient[7] = -0.00 theta[7] = 1.402864 Gradient[8] = -0.00 theta[8] = 0.680376 line search along the specified direction ... zeta = 1.036246 Function Value = 195.9291 Gradient[1] = -0.00 theta[1] = 0.667142 Gradient[2] = -0.00 theta[2] = 0.292860 Gradient[3] = 0.00 theta[3] = 0.034159 Gradient[4] = -0.00 theta[4] = 1.268333 Gradient[5] = -0.00 theta[5] = 1.080987 Gradient[6] = 0.00 theta[6] = 0.044328 Gradient[7] = 0.00 theta[7] = 1.402865 Gradient[8] = 0.00 theta[8] = 0.680376 line search along the specified direction ... zeta = 2.072491 Function Value = 195.9291 Gradient[1] = 0.00 theta[1] = 0.667142 Gradient[2] = 0.00 theta[2] = 0.292860 Gradient[3] = -0.00 theta[3] = 0.034159 Gradient[4] = 0.00 theta[4] = 1.268333 Gradient[5] = 0.00 theta[5] = 1.080987 Gradient[6] = -0.00 theta[6] = 0.044328 Gradient[7] = -0.00 theta[7] = 1.402865 Gradient[8] = -0.00 theta[8] = 0.680376 line search along the specified direction ... zeta = 0.000000 Function Value = 195.9291 Gradient[1] = 0.00 theta[1] = 0.667142 Gradient[2] = 0.00 theta[2] = 0.292860 Gradient[3] = -0.00 theta[3] = 0.034159 Gradient[4] = 0.00 theta[4] = 1.268333 Gradient[5] = 0.00 theta[5] = 1.080987 Gradient[6] = -0.00 theta[6] = 0.044328 Gradient[7] = -0.00 theta[7] = 1.402865 Gradient[8] = -0.00 theta[8] = 0.680376 loglikelihood = -195.92914 AIC = 407.85828 theta[1] = 0.44507839 gradient[1] = 0.0000 theta[2] = 0.08576675 gradient[2] = 0.0000 theta[3] = 0.00116685 gradient[3] = -0.0000 theta[4] = 1.60866927 gradient[4] = 0.0000 theta[5] = 1.16853203 gradient[5] = 0.0000 theta[6] = 0.00196493 gradient[6] = -0.0000 theta[7] = 1.96802946 gradient[7] = -0.0000 theta[8] = 0.46291099 gradient[8] = -0.0000 ====================================================== MLE: mu A c alpha p D 0.445078385 0.085766755 0.001166852 1.608669266 1.168532034 0.001964928 q gamma 1.968029455 0.462910995 ====================================================== iteration 2 0.0043450080 0.0002071276 0.1246958106 declustering: | | | 0% | |======== | 11% | |================ | 22% | |======================= | 33% | |=============================== | 44% | |======================================= | 56% | |=============================================== | 67% | |====================================================== | 78% | |============================================================== | 89% | |======================================================================| 100% iteration: 3 ====================================================== background seismicity rate: Min. 1st Qu. Median Mean 3rd Qu. Max. 0.0002451 0.0084789 0.0221711 0.0341867 0.0395858 0.1449093 probability of being a background event: Min. 1st Qu. Median Mean 3rd Qu. Max. 0.0000002 0.0004366 0.0390968 0.4374084 0.9717464 1.0000000 integral of background seismicity rate: 93.28812 ====================================================== estimating: start Davidon-Fletcher-Powell procedure ... Function Value = 195.9270 Gradient[1] = -0.09 theta[1] = 0.667142 Gradient[2] = -0.02 theta[2] = 0.292860 Gradient[3] = -0.03 theta[3] = 0.034159 Gradient[4] = -0.01 theta[4] = 1.268333 Gradient[5] = 0.03 theta[5] = 1.080987 Gradient[6] = -0.06 theta[6] = 0.044328 Gradient[7] = 0.01 theta[7] = 1.402865 Gradient[8] = -0.00 theta[8] = 0.680376 line search along the specified direction ... zeta = 1.002053 Function Value = 195.9270 Gradient[1] = -0.00 theta[1] = 0.667384 Gradient[2] = 0.00 theta[2] = 0.292896 Gradient[3] = -0.00 theta[3] = 0.034156 Gradient[4] = -0.00 theta[4] = 1.268301 Gradient[5] = 0.00 theta[5] = 1.080968 Gradient[6] = 0.00 theta[6] = 0.044329 Gradient[7] = -0.00 theta[7] = 1.402805 Gradient[8] = 0.00 theta[8] = 0.680302 line search along the specified direction ... zeta = 0.991473 Function Value = 195.9270 Gradient[1] = 0.00 theta[1] = 0.667384 Gradient[2] = 0.00 theta[2] = 0.292895 Gradient[3] = -0.00 theta[3] = 0.034156 Gradient[4] = 0.00 theta[4] = 1.268303 Gradient[5] = 0.00 theta[5] = 1.080968 Gradient[6] = -0.00 theta[6] = 0.044329 Gradient[7] = 0.00 theta[7] = 1.402805 Gradient[8] = -0.00 theta[8] = 0.680303 line search along the specified direction ... zeta = 0.000010 Function Value = 195.9270 Gradient[1] = 0.00 theta[1] = 0.667384 Gradient[2] = 0.00 theta[2] = 0.292895 Gradient[3] = -0.00 theta[3] = 0.034156 Gradient[4] = 0.00 theta[4] = 1.268303 Gradient[5] = 0.00 theta[5] = 1.080968 Gradient[6] = -0.00 theta[6] = 0.044329 Gradient[7] = 0.00 theta[7] = 1.402805 Gradient[8] = -0.00 theta[8] = 0.680303 loglikelihood = -195.92701 AIC = 407.85401 theta[1] = 0.44540183 gradient[1] = 0.0000 theta[2] = 0.08578747 gradient[2] = 0.0000 theta[3] = 0.00116663 gradient[3] = -0.0000 theta[4] = 1.60859313 gradient[4] = 0.0000 theta[5] = 1.16849145 gradient[5] = 0.0000 theta[6] = 0.00196503 gradient[6] = -0.0000 theta[7] = 1.96786150 gradient[7] = 0.0000 theta[8] = 0.46281198 gradient[8] = -0.0000 ====================================================== MLE: mu A c alpha p D 0.445401828 0.085787471 0.001166625 1.608593134 1.168491449 0.001965034 q gamma 1.967861496 0.462811976 ====================================================== iteration 3 7.267099e-04 1.090124e-05 3.143269e-03 declustering: | | | 0% | |========= | 12% | |================== | 25% | |========================== | 38% | |=================================== | 50% | |============================================ | 62% | |==================================================== | 75% | |============================================================= | 88% | |======================================================================| 100% iteration: 4 ====================================================== background seismicity rate: Min. 1st Qu. Median Mean 3rd Qu. Max. 0.0002451 0.0084794 0.0221718 0.0341877 0.0395883 0.1449116 probability of being a background event: Min. 1st Qu. Median Mean 3rd Qu. Max. 0.0000002 0.0004368 0.0391177 0.4374183 0.9717563 1.0000000 integral of background seismicity rate: 93.29044 ====================================================== estimating: start Davidon-Fletcher-Powell procedure ... Function Value = 195.9271 Gradient[1] = 0.00 theta[1] = 0.667384 Gradient[2] = 0.00 theta[2] = 0.292895 Gradient[3] = 0.00 theta[3] = 0.034156 Gradient[4] = 0.00 theta[4] = 1.268303 Gradient[5] = -0.00 theta[5] = 1.080968 Gradient[6] = 0.00 theta[6] = 0.044329 Gradient[7] = -0.00 theta[7] = 1.402805 Gradient[8] = 0.00 theta[8] = 0.680303 line search along the specified direction ... zeta = 1.000124 Function Value = 195.9271 Gradient[1] = -0.00 theta[1] = 0.667376 Gradient[2] = -0.00 theta[2] = 0.292894 Gradient[3] = 0.00 theta[3] = 0.034156 Gradient[4] = -0.00 theta[4] = 1.268304 Gradient[5] = -0.00 theta[5] = 1.080968 Gradient[6] = 0.00 theta[6] = 0.044329 Gradient[7] = -0.00 theta[7] = 1.402805 Gradient[8] = 0.00 theta[8] = 0.680304 line search along the specified direction ... zeta = 1.166667 Function Value = 195.9271 Gradient[1] = 0.00 theta[1] = 0.667376 Gradient[2] = 0.00 theta[2] = 0.292894 Gradient[3] = -0.00 theta[3] = 0.034156 Gradient[4] = 0.00 theta[4] = 1.268304 Gradient[5] = 0.00 theta[5] = 1.080968 Gradient[6] = -0.00 theta[6] = 0.044329 Gradient[7] = 0.00 theta[7] = 1.402805 Gradient[8] = -0.00 theta[8] = 0.680304 loglikelihood = -195.92712 AIC = 407.85425 theta[1] = 0.44539126 gradient[1] = 0.0000 theta[2] = 0.08578713 gradient[2] = 0.0000 theta[3] = 0.00116663 gradient[3] = -0.0000 theta[4] = 1.60859412 gradient[4] = 0.0000 theta[5] = 1.16849225 gradient[5] = 0.0000 theta[6] = 0.00196503 gradient[6] = -0.0000 theta[7] = 1.96786272 gradient[7] = 0.0000 theta[8] = 0.46281317 gradient[8] = -0.0000 ====================================================== MLE: mu A c alpha p D 0.445391258 0.085787129 0.001166630 1.608594121 1.168492255 0.001965032 q gamma 1.967862719 0.462813174 ====================================================== iteration 4 3.862685e-06 6.035089e-07 7.875552e-05 Execution time: user system elapsed 148.47 0.00 148.46 declustering: | | | 0% | |====== | 9% | |============= | 18% | |=================== | 27% | |========================= | 36% | |================================ | 45% | |====================================== | 55% | |============================================= | 64% | |=================================================== | 73% | |========================================================= | 82% | |================================================================ | 91% | |======================================================================| 100% iteration: 1 ====================================================== background seismicity rate: Min. 1st Qu. Median Mean 3rd Qu. Max. 0.0001471 0.0070248 0.0151810 0.0299157 0.0353152 0.1334580 probability of being a background event: Min. 1st Qu. Median Mean 3rd Qu. Max. 0.0000016 0.9264104 0.9855878 0.8097728 0.9945622 1.0000000 integral of background seismicity rate: 91.61371 ====================================================== estimating: start Davidon-Fletcher-Powell procedure ... Function Value = 252.6672 Gradient[1] = -3.58 theta[1] = 0.667376 Gradient[2] = 73.24 theta[2] = 0.292894 Gradient[3] = -12.76 theta[3] = 0.034156 Gradient[4] = 35.66 theta[4] = 1.268304 Gradient[5] = 46.54 theta[5] = 1.080968 Gradient[6] = -78.53 theta[6] = 0.044329 Gradient[7] = 8.24 theta[7] = 1.402805 Gradient[8] = -5.32 theta[8] = 0.680304 line search along the specified direction ... zeta = 0.000405 Function Value = 250.0463 Gradient[1] = -2.93 theta[1] = 0.668825 Gradient[2] = 29.33 theta[2] = 0.263231 Gradient[3] = 22.46 theta[3] = 0.039324 Gradient[4] = 8.59 theta[4] = 1.253861 Gradient[5] = 18.77 theta[5] = 1.062119 Gradient[6] = 37.68 theta[6] = 0.076134 Gradient[7] = -5.79 theta[7] = 1.399467 Gradient[8] = 1.99 theta[8] = 0.682458 line search along the specified direction ... zeta = 0.000693 Function Value = 249.0368 Gradient[1] = -3.55 theta[1] = 0.670933 Gradient[2] = -1.99 theta[2] = 0.237568 Gradient[3] = 15.36 theta[3] = 0.028021 Gradient[4] = -6.70 theta[4] = 1.244623 Gradient[5] = -14.49 theta[5] = 1.045734 Gradient[6] = 15.60 theta[6] = 0.064094 Gradient[7] = -2.96 theta[7] = 1.401782 Gradient[8] = -0.19 theta[8] = 0.681966 line search along the specified direction ... zeta = 0.000359 Function Value = 248.9928 Gradient[1] = -3.45 theta[1] = 0.672231 Gradient[2] = 0.06 theta[2] = 0.234814 Gradient[3] = -3.87 theta[3] = 0.024425 Gradient[4] = -5.82 theta[4] = 1.244999 Gradient[5] = -7.35 theta[5] = 1.048382 Gradient[6] = 17.27 theta[6] = 0.064861 Gradient[7] = -3.17 theta[7] = 1.402079 Gradient[8] = -0.04 theta[8] = 0.682412 line search along the specified direction ... zeta = 0.005515 Function Value = 248.7233 Gradient[1] = 1.25 theta[1] = 0.691891 Gradient[2] = -5.16 theta[2] = 0.200422 Gradient[3] = 7.04 theta[3] = 0.031736 Gradient[4] = -7.30 theta[4] = 1.255864 Gradient[5] = -5.80 theta[5] = 1.077610 Gradient[6] = -0.05 theta[6] = 0.059452 Gradient[7] = -0.92 theta[7] = 1.408728 Gradient[8] = -1.33 theta[8] = 0.688263 line search along the specified direction ... zeta = 0.019695 Function Value = 248.5048 Gradient[1] = -1.35 theta[1] = 0.686453 Gradient[2] = -18.56 theta[2] = 0.162510 Gradient[3] = 8.67 theta[3] = 0.031149 Gradient[4] = -9.28 theta[4] = 1.330327 Gradient[5] = -11.36 theta[5] = 1.076932 Gradient[6] = 7.42 theta[6] = 0.063242 Gradient[7] = -1.90 theta[7] = 1.427562 Gradient[8] = -0.75 theta[8] = 0.712064 line search along the specified direction ... zeta = 0.005626 Function Value = 248.4802 Gradient[1] = 2.22 theta[1] = 0.695922 Gradient[2] = -17.50 theta[2] = 0.161451 Gradient[3] = 8.30 theta[3] = 0.030104 Gradient[4] = -8.65 theta[4] = 1.339845 Gradient[5] = -12.95 theta[5] = 1.072836 Gradient[6] = 10.30 theta[6] = 0.064628 Gradient[7] = -2.23 theta[7] = 1.433323 Gradient[8] = -0.51 theta[8] = 0.718039 line search along the specified direction ... zeta = 0.169405 Function Value = 248.2413 Gradient[1] = 0.80 theta[1] = 0.691343 Gradient[2] = -15.88 theta[2] = 0.162440 Gradient[3] = 7.45 theta[3] = 0.029658 Gradient[4] = -8.12 theta[4] = 1.342274 Gradient[5] = -12.83 theta[5] = 1.070860 Gradient[6] = 25.45 theta[6] = 0.075518 Gradient[7] = -2.82 theta[7] = 1.613655 Gradient[8] = 1.45 theta[8] = 0.877214 line search along the specified direction ... zeta = 1.355240 Function Value = 247.7069 Gradient[1] = -1.12 theta[1] = 0.682864 Gradient[2] = -6.02 theta[2] = 0.166069 Gradient[3] = 4.34 theta[3] = 0.029678 Gradient[4] = -3.73 theta[4] = 1.363308 Gradient[5] = -6.78 theta[5] = 1.068722 Gradient[6] = 14.18 theta[6] = 0.036916 Gradient[7] = -1.69 theta[7] = 1.509048 Gradient[8] = 0.33 theta[8] = 1.140858 line search along the specified direction ... zeta = 0.023211 Function Value = 247.6871 Gradient[1] = -1.22 theta[1] = 0.682412 Gradient[2] = -5.40 theta[2] = 0.166274 Gradient[3] = 4.15 theta[3] = 0.029690 Gradient[4] = -3.44 theta[4] = 1.364474 Gradient[5] = -6.39 theta[5] = 1.068657 Gradient[6] = 2.04 theta[6] = 0.033782 Gradient[7] = -1.02 theta[7] = 1.500427 Gradient[8] = -0.41 theta[8] = 1.162970 line search along the specified direction ... zeta = 0.497357 Function Value = 247.5881 Gradient[1] = -1.41 theta[1] = 0.681303 Gradient[2] = 1.57 theta[2] = 0.162072 Gradient[3] = 2.26 theta[3] = 0.030120 Gradient[4] = 0.09 theta[4] = 1.387017 Gradient[5] = -2.21 theta[5] = 1.070669 Gradient[6] = 13.39 theta[6] = 0.031034 Gradient[7] = -1.06 theta[7] = 1.570132 Gradient[8] = 1.09 theta[8] = 1.268808 line search along the specified direction ... zeta = 1.498758 Function Value = 247.5451 Gradient[1] = 0.85 theta[1] = 0.690008 Gradient[2] = -1.75 theta[2] = 0.149930 Gradient[3] = 0.09 theta[3] = 0.031442 Gradient[4] = -0.50 theta[4] = 1.401040 Gradient[5] = -0.85 theta[5] = 1.081620 Gradient[6] = -0.31 theta[6] = 0.033413 Gradient[7] = -0.28 theta[7] = 1.610429 Gradient[8] = 0.04 theta[8] = 1.249117 line search along the specified direction ... zeta = 1.546368 Function Value = 247.5374 Gradient[1] = 0.44 theta[1] = 0.689313 Gradient[2] = -1.44 theta[2] = 0.150482 Gradient[3] = -0.69 theta[3] = 0.031909 Gradient[4] = -0.57 theta[4] = 1.396980 Gradient[5] = 0.11 theta[5] = 1.084827 Gradient[6] = 7.77 theta[6] = 0.035922 Gradient[7] = -0.50 theta[7] = 1.656420 Gradient[8] = 0.70 theta[8] = 1.252232 line search along the specified direction ... zeta = 3.738758 Function Value = 247.5298 Gradient[1] = -0.06 theta[1] = 0.687554 Gradient[2] = 0.78 theta[2] = 0.152233 Gradient[3] = -0.48 theta[3] = 0.032016 Gradient[4] = 0.24 theta[4] = 1.396480 Gradient[5] = 0.80 theta[5] = 1.084317 Gradient[6] = -1.31 theta[6] = 0.037936 Gradient[7] = 0.01 theta[7] = 1.693162 Gradient[8] = -0.10 theta[8] = 1.238832 line search along the specified direction ... zeta = 1.810327 Function Value = 247.5288 Gradient[1] = -0.16 theta[1] = 0.687151 Gradient[2] = 0.34 theta[2] = 0.151118 Gradient[3] = -0.26 theta[3] = 0.031787 Gradient[4] = 0.23 theta[4] = 1.399852 Gradient[5] = 0.36 theta[5] = 1.083031 Gradient[6] = 0.37 theta[6] = 0.039308 Gradient[7] = -0.03 theta[7] = 1.710076 Gradient[8] = -0.00 theta[8] = 1.233453 line search along the specified direction ... zeta = 1.859335 Function Value = 247.5284 Gradient[1] = -0.04 theta[1] = 0.687366 Gradient[2] = 0.14 theta[2] = 0.152284 Gradient[3] = -0.11 theta[3] = 0.031677 Gradient[4] = 0.03 theta[4] = 1.396917 Gradient[5] = 0.10 theta[5] = 1.082323 Gradient[6] = 0.79 theta[6] = 0.039428 Gradient[7] = -0.04 theta[7] = 1.716854 Gradient[8] = 0.06 theta[8] = 1.236206 line search along the specified direction ... zeta = 1.304278 Function Value = 247.5284 Gradient[1] = 0.00 theta[1] = 0.687481 Gradient[2] = -0.02 theta[2] = 0.152100 Gradient[3] = -0.00 theta[3] = 0.031670 Gradient[4] = -0.01 theta[4] = 1.397269 Gradient[5] = 0.00 theta[5] = 1.082194 Gradient[6] = -0.01 theta[6] = 0.039387 Gradient[7] = 0.00 theta[7] = 1.717730 Gradient[8] = -0.01 theta[8] = 1.236350 line search along the specified direction ... zeta = 0.792885 Function Value = 247.5284 Gradient[1] = 0.00 theta[1] = 0.687464 Gradient[2] = 0.00 theta[2] = 0.152089 Gradient[3] = 0.00 theta[3] = 0.031664 Gradient[4] = 0.00 theta[4] = 1.397398 Gradient[5] = -0.00 theta[5] = 1.082142 Gradient[6] = 0.01 theta[6] = 0.039319 Gradient[7] = -0.00 theta[7] = 1.717406 Gradient[8] = 0.00 theta[8] = 1.237007 line search along the specified direction ... zeta = 1.106619 Function Value = 247.5284 Gradient[1] = -0.00 theta[1] = 0.687461 Gradient[2] = -0.00 theta[2] = 0.152094 Gradient[3] = -0.00 theta[3] = 0.031664 Gradient[4] = -0.00 theta[4] = 1.397381 Gradient[5] = 0.00 theta[5] = 1.082144 Gradient[6] = -0.00 theta[6] = 0.039311 Gradient[7] = 0.00 theta[7] = 1.717276 Gradient[8] = -0.00 theta[8] = 1.237024 line search along the specified direction ... zeta = 1.127930 Function Value = 247.5284 Gradient[1] = 0.00 theta[1] = 0.687461 Gradient[2] = 0.00 theta[2] = 0.152095 Gradient[3] = -0.00 theta[3] = 0.031663 Gradient[4] = 0.00 theta[4] = 1.397381 Gradient[5] = 0.00 theta[5] = 1.082143 Gradient[6] = -0.00 theta[6] = 0.039311 Gradient[7] = 0.00 theta[7] = 1.717267 Gradient[8] = -0.00 theta[8] = 1.237020 line search along the specified direction ... zeta = 0.982143 Function Value = 247.5284 Gradient[1] = -0.00 theta[1] = 0.687461 Gradient[2] = -0.00 theta[2] = 0.152095 Gradient[3] = -0.00 theta[3] = 0.031663 Gradient[4] = -0.00 theta[4] = 1.397381 Gradient[5] = 0.00 theta[5] = 1.082142 Gradient[6] = 0.00 theta[6] = 0.039311 Gradient[7] = -0.00 theta[7] = 1.717266 Gradient[8] = 0.00 theta[8] = 1.237020 line search along the specified direction ... zeta = 0.098214 Function Value = 247.5284 Gradient[1] = -0.00 theta[1] = 0.687461 Gradient[2] = -0.00 theta[2] = 0.152095 Gradient[3] = -0.00 theta[3] = 0.031663 Gradient[4] = -0.00 theta[4] = 1.397381 Gradient[5] = 0.00 theta[5] = 1.082142 Gradient[6] = 0.00 theta[6] = 0.039311 Gradient[7] = -0.00 theta[7] = 1.717266 Gradient[8] = 0.00 theta[8] = 1.237020 line search along the specified direction ... zeta = 0.196429 Function Value = 247.5284 Gradient[1] = -0.00 theta[1] = 0.687461 Gradient[2] = -0.00 theta[2] = 0.152095 Gradient[3] = -0.00 theta[3] = 0.031663 Gradient[4] = -0.00 theta[4] = 1.397381 Gradient[5] = 0.00 theta[5] = 1.082142 Gradient[6] = 0.00 theta[6] = 0.039311 Gradient[7] = -0.00 theta[7] = 1.717266 Gradient[8] = 0.00 theta[8] = 1.237020 line search along the specified direction ... zeta = 0.000001 Function Value = 247.5284 Gradient[1] = -0.00 theta[1] = 0.687461 Gradient[2] = -0.00 theta[2] = 0.152095 Gradient[3] = -0.00 theta[3] = 0.031663 Gradient[4] = -0.00 theta[4] = 1.397381 Gradient[5] = 0.00 theta[5] = 1.082142 Gradient[6] = 0.00 theta[6] = 0.039311 Gradient[7] = -0.00 theta[7] = 1.717266 Gradient[8] = 0.00 theta[8] = 1.237020 loglikelihood = -247.52839 AIC = 511.05678 theta[1] = 0.47260255 gradient[1] = -0.0000 theta[2] = 0.02313275 gradient[2] = -0.0000 theta[3] = 0.00100257 gradient[3] = -0.0000 theta[4] = 1.95267423 gradient[4] = -0.0000 theta[5] = 1.17103233 gradient[5] = 0.0000 theta[6] = 0.00154532 gradient[6] = 0.0000 theta[7] = 2.94900406 gradient[7] = -0.0000 theta[8] = 1.53021757 gradient[8] = 0.0000 ====================================================== MLE: mu A c alpha p D 0.472602548 0.023132755 0.001002569 1.952674229 1.171032328 0.001545319 q gamma 2.949004060 1.530217574 ====================================================== declustering: | | | 0% | |======= | 10% | |============== | 20% | |===================== | 30% | |============================ | 40% | |=================================== | 50% | |========================================== | 60% | |================================================= | 70% | |======================================================== | 80% | |=============================================================== | 90% | |======================================================================| 100% iteration: 2 ====================================================== background seismicity rate: Min. 1st Qu. Median Mean 3rd Qu. Max. 0.0001472 0.0074063 0.0164543 0.0309562 0.0360039 0.1377179 probability of being a background event: Min. 1st Qu. Median Mean 3rd Qu. Max. 0.0000022 0.9723211 0.9953023 0.8382518 0.9995039 1.0000000 integral of background seismicity rate: 95.2966 ====================================================== estimating: start Davidon-Fletcher-Powell procedure ... Function Value = 247.5072 Gradient[1] = 4.86 theta[1] = 0.687461 Gradient[2] = 0.94 theta[2] = 0.152095 Gradient[3] = 0.77 theta[3] = 0.031663 Gradient[4] = 0.28 theta[4] = 1.397381 Gradient[5] = -0.79 theta[5] = 1.082142 Gradient[6] = -0.92 theta[6] = 0.039311 Gradient[7] = 0.04 theta[7] = 1.717266 Gradient[8] = -0.06 theta[8] = 1.237020 line search along the specified direction ... zeta = 0.973040 Function Value = 247.4747 Gradient[1] = 0.00 theta[1] = 0.674359 Gradient[2] = 0.00 theta[2] = 0.151266 Gradient[3] = 0.02 theta[3] = 0.031805 Gradient[4] = 0.00 theta[4] = 1.398372 Gradient[5] = -0.02 theta[5] = 1.083134 Gradient[6] = 0.04 theta[6] = 0.039502 Gradient[7] = -0.00 theta[7] = 1.718871 Gradient[8] = 0.00 theta[8] = 1.236528 line search along the specified direction ... zeta = 1.013055 Function Value = 247.4747 Gradient[1] = -0.00 theta[1] = 0.674371 Gradient[2] = -0.00 theta[2] = 0.151255 Gradient[3] = 0.00 theta[3] = 0.031816 Gradient[4] = -0.00 theta[4] = 1.398317 Gradient[5] = -0.00 theta[5] = 1.083215 Gradient[6] = -0.00 theta[6] = 0.039536 Gradient[7] = 0.00 theta[7] = 1.719501 Gradient[8] = -0.00 theta[8] = 1.236436 line search along the specified direction ... zeta = 1.015072 Function Value = 247.4747 Gradient[1] = 0.00 theta[1] = 0.674371 Gradient[2] = 0.00 theta[2] = 0.151255 Gradient[3] = 0.00 theta[3] = 0.031816 Gradient[4] = -0.00 theta[4] = 1.398318 Gradient[5] = -0.00 theta[5] = 1.083215 Gradient[6] = 0.00 theta[6] = 0.039537 Gradient[7] = -0.00 theta[7] = 1.719502 Gradient[8] = 0.00 theta[8] = 1.236435 line search along the specified direction ... zeta = 0.000000 Function Value = 247.4747 Gradient[1] = 0.00 theta[1] = 0.674371 Gradient[2] = 0.00 theta[2] = 0.151255 Gradient[3] = 0.00 theta[3] = 0.031816 Gradient[4] = -0.00 theta[4] = 1.398318 Gradient[5] = -0.00 theta[5] = 1.083215 Gradient[6] = 0.00 theta[6] = 0.039537 Gradient[7] = -0.00 theta[7] = 1.719502 Gradient[8] = 0.00 theta[8] = 1.236435 loglikelihood = -247.47465 AIC = 510.94930 theta[1] = 0.45477616 gradient[1] = 0.0000 theta[2] = 0.02287811 gradient[2] = 0.0000 theta[3] = 0.00101223 gradient[3] = 0.0000 theta[4] = 1.95529427 gradient[4] = -0.0000 theta[5] = 1.17335524 gradient[5] = -0.0000 theta[6] = 0.00156315 gradient[6] = 0.0000 theta[7] = 2.95668877 gradient[7] = -0.0000 theta[8] = 1.52877174 gradient[8] = 0.0000 ====================================================== MLE: mu A c alpha p D 0.454776163 0.022878114 0.001012233 1.955294271 1.173355240 0.001563153 q gamma 2.956688773 1.528771745 ====================================================== iteration 2 0.0115407439 0.0002171021 0.1118169531 declustering: | | | 0% | |======== | 11% | |================ | 22% | |======================= | 33% | |=============================== | 44% | |======================================= | 56% | |=============================================== | 67% | |====================================================== | 78% | |============================================================== | 89% | |======================================================================| 100% iteration: 3 ====================================================== background seismicity rate: Min. 1st Qu. Median Mean 3rd Qu. Max. 0.0001472 0.0074029 0.0164381 0.0309389 0.0359985 0.1376676 probability of being a background event: Min. 1st Qu. Median Mean 3rd Qu. Max. 0.0000021 0.9716104 0.9952096 0.8378006 0.9994939 1.0000000 integral of background seismicity rate: 95.23894 ====================================================== estimating: start Davidon-Fletcher-Powell procedure ... Function Value = 247.4784 Gradient[1] = -0.07 theta[1] = 0.674371 Gradient[2] = -0.02 theta[2] = 0.151255 Gradient[3] = -0.02 theta[3] = 0.031816 Gradient[4] = -0.01 theta[4] = 1.398318 Gradient[5] = 0.02 theta[5] = 1.083215 Gradient[6] = 0.00 theta[6] = 0.039537 Gradient[7] = -0.00 theta[7] = 1.719502 Gradient[8] = -0.00 theta[8] = 1.236435 line search along the specified direction ... zeta = 0.992200 Function Value = 247.4784 Gradient[1] = 0.00 theta[1] = 0.674559 Gradient[2] = -0.00 theta[2] = 0.151284 Gradient[3] = -0.00 theta[3] = 0.031809 Gradient[4] = -0.00 theta[4] = 1.398305 Gradient[5] = 0.00 theta[5] = 1.083166 Gradient[6] = -0.00 theta[6] = 0.039532 Gradient[7] = 0.00 theta[7] = 1.719546 Gradient[8] = -0.00 theta[8] = 1.236516 line search along the specified direction ... zeta = 1.008345 Function Value = 247.4784 Gradient[1] = 0.00 theta[1] = 0.674559 Gradient[2] = 0.00 theta[2] = 0.151284 Gradient[3] = 0.00 theta[3] = 0.031809 Gradient[4] = 0.00 theta[4] = 1.398305 Gradient[5] = 0.00 theta[5] = 1.083166 Gradient[6] = 0.00 theta[6] = 0.039532 Gradient[7] = -0.00 theta[7] = 1.719548 Gradient[8] = 0.00 theta[8] = 1.236517 line search along the specified direction ... zeta = 0.001008 Function Value = 247.4784 Gradient[1] = 0.00 theta[1] = 0.674559 Gradient[2] = 0.00 theta[2] = 0.151284 Gradient[3] = 0.00 theta[3] = 0.031809 Gradient[4] = 0.00 theta[4] = 1.398305 Gradient[5] = 0.00 theta[5] = 1.083166 Gradient[6] = 0.00 theta[6] = 0.039532 Gradient[7] = -0.00 theta[7] = 1.719548 Gradient[8] = 0.00 theta[8] = 1.236517 loglikelihood = -247.47842 AIC = 510.95683 theta[1] = 0.45502959 gradient[1] = 0.0000 theta[2] = 0.02288686 gradient[2] = 0.0000 theta[3] = 0.00101182 gradient[3] = 0.0000 theta[4] = 1.95525678 gradient[4] = 0.0000 theta[5] = 1.17324777 gradient[5] = 0.0000 theta[6] = 0.00156276 gradient[6] = 0.0000 theta[7] = 2.95684689 gradient[7] = -0.0000 theta[8] = 1.52897498 gradient[8] = 0.0000 ====================================================== MLE: mu A c alpha p D 0.455029592 0.022886864 0.001011817 1.955256782 1.173247772 0.001562763 q gamma 2.956846893 1.528974976 ====================================================== iteration 3 5.572630e-04 1.521958e-05 2.745565e-03 declustering: | | | 0% | |========= | 12% | |================== | 25% | |========================== | 38% | |=================================== | 50% | |============================================ | 62% | |==================================================== | 75% | |============================================================= | 88% | |======================================================================| 100% iteration: 4 ====================================================== background seismicity rate: Min. 1st Qu. Median Mean 3rd Qu. Max. 0.0001472 0.0074028 0.0164379 0.0309387 0.0359984 0.1376668 probability of being a background event: Min. 1st Qu. Median Mean 3rd Qu. Max. 0.0000021 0.9716071 0.9952074 0.8377959 0.9994939 1.0000000 integral of background seismicity rate: 95.23828 ====================================================== estimating: start Davidon-Fletcher-Powell procedure ... Function Value = 247.4784 Gradient[1] = -0.00 theta[1] = 0.674559 Gradient[2] = -0.00 theta[2] = 0.151284 Gradient[3] = -0.00 theta[3] = 0.031809 Gradient[4] = -0.00 theta[4] = 1.398305 Gradient[5] = 0.00 theta[5] = 1.083166 Gradient[6] = 0.00 theta[6] = 0.039532 Gradient[7] = -0.00 theta[7] = 1.719548 Gradient[8] = 0.00 theta[8] = 1.236517 line search along the specified direction ... zeta = 1.000304 Function Value = 247.4784 Gradient[1] = 0.00 theta[1] = 0.674561 Gradient[2] = -0.00 theta[2] = 0.151284 Gradient[3] = 0.00 theta[3] = 0.031809 Gradient[4] = -0.00 theta[4] = 1.398305 Gradient[5] = -0.00 theta[5] = 1.083165 Gradient[6] = 0.00 theta[6] = 0.039532 Gradient[7] = -0.00 theta[7] = 1.719548 Gradient[8] = 0.00 theta[8] = 1.236518 line search along the specified direction ... zeta = 0.000100 Function Value = 247.4784 Gradient[1] = 0.00 theta[1] = 0.674561 Gradient[2] = -0.00 theta[2] = 0.151284 Gradient[3] = 0.00 theta[3] = 0.031809 Gradient[4] = -0.00 theta[4] = 1.398305 Gradient[5] = -0.00 theta[5] = 1.083165 Gradient[6] = 0.00 theta[6] = 0.039532 Gradient[7] = -0.00 theta[7] = 1.719548 Gradient[8] = 0.00 theta[8] = 1.236518 loglikelihood = -247.47843 AIC = 510.95686 theta[1] = 0.45503257 gradient[1] = 0.0000 theta[2] = 0.02288696 gradient[2] = -0.0000 theta[3] = 0.00101181 gradient[3] = 0.0000 theta[4] = 1.95525572 gradient[4] = -0.0000 theta[5] = 1.17324697 gradient[5] = -0.0000 theta[6] = 0.00156276 gradient[6] = 0.0000 theta[7] = 2.95684679 gradient[7] = -0.0000 theta[8] = 1.52897652 gradient[8] = 0.0000 ====================================================== MLE: mu A c alpha p D 0.455032566 0.022886958 0.001011814 1.955255718 1.173246970 0.001562758 q gamma 2.956846792 1.528976521 ====================================================== iteration 4 6.534604e-06 5.474307e-08 2.317690e-05 Execution time: user system elapsed 85.10 0.00 85.11 [ FAIL 0 | WARN 0 | SKIP 0 | PASS 12 ] > > proc.time() user system elapsed 236.37 0.15 236.53