wnl-package.RdThis is a minimization tool to solve the examples in the book Gabrielsson J, Weiner D. 'Pharmacokinetic and Pharmacodynamic Data Analysis - Concepts and Applications' 5th ed. 2016. (ISBN:9198299107).
This is a set of minimization tools to solve all the examples in the book 'Pharmacokinetic and Pharmacodynamic Data Analysis - Concepts and Applications' 5th ed. 2016.
Kyun-Seop Bae <k@acr.kr>
Gabrielsson J, Weiner D. Pharmacokinetic and Pharmacodynamic Data Analysis - Concepts and Applications. 5th ed. 2016.
tData = Theoph colnames(tData) = c("ID", "BWT", "DOSE", "TIME", "DV") fPK = function(THETA) # Prediction function { DOSE = 320000 # in microgram TIME = e$DATA[,"TIME"] # use data in e$DATA K = THETA[1] Ka = THETA[2] V = THETA[3] Cp = DOSE/V*Ka/(Ka - K)*(exp(-K*TIME) - exp(-Ka*TIME)) return(Cp) } IDs = unique(tData[,"ID"]) nID = length(IDs) for (i in 1:nID) { Data = tData[tData$ID == IDs[i],] Res = nlr(fPK, Data, pNames=c("k", "ka", "V"), IE=c(0.1, 3, 500), SecNames=c("CL", "Thalf", "MRT"), SecForms=c(~V*k, ~log(2)/k, ~1/k)) print(paste("## ID =", i, "##")) print(Res) }#> [1] "## ID = 1 ##" #> $Est #> k ka V AddErrVar AddErrSD CL #> PE 0.053954312 1.7774192 29394.222169 0.3896372 0.6242092 1585.94503 #> SE 0.007784005 0.2299941 1449.530064 0.1661418 0.1330819 172.85019 #> RSE 14.427031240 12.9397799 4.931343 42.6401419 21.3200709 10.89888 #> Thalf MRT #> PE 12.84693 18.534200 #> SE 1.85343 2.673935 #> RSE 14.42703 14.427031 #> #> $Cov #> k ka V AddErrVar #> k 6.059074e-05 -9.280489e-04 -9.013867e+00 -1.664726e-08 #> ka -9.280489e-04 5.289730e-02 2.119645e+02 3.840233e-07 #> V -9.013867e+00 2.119645e+02 2.101137e+06 3.776607e-03 #> AddErrVar -1.664726e-08 3.840233e-07 3.776607e-03 2.760311e-02 #> #> $run #> $run$m #> [1] 5 #> #> $run$n #> [1] 6 #> #> $run$run #> [1] 5 #> #> $run$p.value #> [1] 0.2619048 #> #> #> $`Objective Function Value` #> [1] 0.632068 #> #> $`-2LL` #> [1] 20.84872 #> #> $AIC #> [1] 28.84872 #> #> $AICc #> [1] 35.51538 #> #> $BIC #> [1] 30.4403 #> #> $Convergence #> NULL #> #> $Message #> [1] "CONVERGENCE: REL_REDUCTION_OF_F <= FACTR*EPSMCH" #> #> $Prediction #> [1] 0.000000 3.877507 6.810852 9.035310 9.758253 9.123553 8.525224 7.683263 #> [9] 6.889936 5.838199 3.014645 #> #> $Residual #> [1] 0.7400000000 -1.0375074250 -0.2408522030 1.4646901948 -0.0982525982 #> [6] -0.5435530326 -0.1652236347 -0.2132625771 0.0000635282 0.1018013770 #> [11] 0.2653553178 #> #> $`Elapsed Time` #> Time difference of 0.01168203 secs #> #> [1] "## ID = 2 ##" #> $Est #> k ka V AddErrVar AddErrSD CL #> PE 0.1016610 1.9426639 32024.771537 0.8134795 0.9019310 3255.67032 #> SE 0.0207964 0.4307417 3064.675035 0.3468676 0.1922917 451.30587 #> RSE 20.4566171 22.1727339 9.569701 42.6399995 21.3199997 13.86215 #> Thalf MRT #> PE 6.818221 9.836614 #> SE 1.394777 2.012238 #> RSE 20.456617 20.456617 #> #> $Cov #> k ka V AddErrVar #> k 4.324903e-04 -5.234795e-03 -5.174777e+01 -2.618365e-08 #> ka -5.234795e-03 1.855384e-01 9.264106e+02 1.125743e-07 #> V -5.174777e+01 9.264106e+02 9.392233e+06 4.884969e-03 #> AddErrVar -2.618365e-08 1.125743e-07 4.884969e-03 1.203172e-01 #> #> $run #> $run$m #> [1] 5 #> #> $run$n #> [1] 5 #> #> $run$run #> [1] 6 #> #> $run$p.value #> [1] 0.6428571 #> #> #> $`Objective Function Value` #> [1] 8.729257 #> #> $`-2LL` #> [1] 28.9459 #> #> $AIC #> [1] 36.9459 #> #> $AICc #> [1] 43.61257 #> #> $BIC #> [1] 38.53749 #> #> $Convergence #> NULL #> #> $Message #> [1] "CONVERGENCE: REL_REDUCTION_OF_F <= FACTR*EPSMCH" #> #> $Prediction #> [1] 0.0000000 4.0181763 6.1614927 8.0136206 8.4213596 7.3754388 6.3288963 #> [8] 5.1597171 4.2232793 3.1131310 0.8915432 #> #> $Residual #> [1] 0.000000000 -2.298176259 1.748507295 0.296379432 -0.091359642 #> [6] -0.525438800 -0.248896251 0.240282904 0.326720701 -0.103131003 #> [11] 0.008456754 #> #> $`Elapsed Time` #> Time difference of 0.005707026 secs #> #> [1] "## ID = 3 ##" #> $Est #> k ka V AddErrVar AddErrSD CL #> PE 0.081424641 2.4535753 34319.365330 0.03966119 0.19915117 2794.44201 #> SE 0.003809873 0.1437188 684.951604 0.01691155 0.04245909 92.85260 #> RSE 4.679017016 5.8575265 1.995817 42.64005756 21.32002878 3.32276 #> Thalf MRT #> PE 8.5127447 12.2812945 #> SE 0.3983128 0.5746439 #> RSE 4.6790170 4.6790170 #> #> $Cov #> k ka V AddErrVar #> k 1.451513e-05 -3.094035e-04 -2.072882e+00 -2.225382e-09 #> ka -3.094035e-04 2.065510e-02 6.848231e+01 6.715511e-08 #> V -2.072882e+00 6.848231e+01 4.691587e+05 4.890908e-04 #> AddErrVar -2.225382e-09 6.715510e-08 4.890907e-04 2.860006e-04 #> #> $run #> $run$m #> [1] 5 #> #> $run$n #> [1] 5 #> #> $run$run #> [1] 7 #> #> $run$p.value #> [1] 0.3571429 #> #> #> $`Objective Function Value` #> [1] -24.50118 #> #> $`-2LL` #> [1] -4.284535 #> #> $AIC #> [1] 3.715465 #> #> $AICc #> [1] 10.38213 #> #> $BIC #> [1] 5.307046 #> #> $Convergence #> NULL #> #> $Message #> [1] "CONVERGENCE: REL_REDUCTION_OF_F <= FACTR*EPSMCH" #> #> $Prediction #> [1] 0.000000 4.462162 6.875370 8.086058 8.113669 7.180843 6.377126 5.423212 #> [9] 4.634548 3.586049 1.347602 #> #> $Residual #> [1] 0.00000000 -0.06216214 0.02463040 0.11394202 -0.31366860 0.31915699 #> [7] -0.17712563 -0.12321179 0.26545182 0.11395080 -0.29760226 #> #> $`Elapsed Time` #> Time difference of 0.006668091 secs #> #> [1] "## ID = 4 ##" #> $Est #> k ka V AddErrVar AddErrSD CL #> PE 0.08746650 1.1714812 31097.469810 0.5210868 0.7218634 2719.98679 #> SE 0.01605189 0.2035763 2607.290790 0.2221923 0.1539019 328.29334 #> RSE 18.35203984 17.3776875 8.384254 42.6401720 21.3200860 12.06967 #> Thalf MRT #> PE 7.924716 11.432949 #> SE 1.454347 2.098179 #> RSE 18.352040 18.352040 #> #> $Cov #> k ka V AddErrVar #> k 2.576631e-04 -2.094472e-03 -3.555243e+01 -3.661761e-08 #> ka -2.094472e-03 4.144333e-02 3.931963e+02 4.001435e-07 #> V -3.555243e+01 3.931963e+02 6.797965e+06 6.916608e-03 #> AddErrVar -3.661761e-08 4.001436e-07 6.916609e-03 4.936942e-02 #> #> $run #> $run$m #> [1] 4 #> #> $run$n #> [1] 6 #> #> $run$run #> [1] 5 #> #> $run$p.value #> [1] 0.4047619 #> #> #> $`Objective Function Value` #> [1] 3.829767 #> #> $`-2LL` #> [1] 24.04641 #> #> $AIC #> [1] 32.04641 #> #> $AICc #> [1] 38.71308 #> #> $BIC #> [1] 33.638 #> #> $Convergence #> NULL #> #> $Message #> [1] "CONVERGENCE: REL_REDUCTION_OF_F <= FACTR*EPSMCH" #> #> $Prediction #> [1] 0.000000 3.405243 5.045608 6.951979 8.313101 8.003652 7.137545 6.015151 #> [9] 5.052017 3.899862 1.287552 #> #> $Residual #> [1] 0.00000000 -1.51524257 -0.44560789 1.64802135 0.06689894 -0.46365195 #> [7] -0.25754508 -0.23515116 0.27798341 0.29013842 -0.13755210 #> #> $`Elapsed Time` #> Time difference of 0.005770922 secs #> #> [1] "## ID = 5 ##" #> $Est #> k ka V AddErrVar AddErrSD CL #> PE 0.08843508 1.4715011 26925.05354 1.2239519 1.106324 2381.11936 #> SE 0.02009436 0.3192602 2739.87065 0.5218949 0.235869 366.26415 #> RSE 22.72216233 21.6962239 10.17592 42.6401459 21.320073 15.38202 #> Thalf MRT #> PE 7.837921 11.307729 #> SE 1.780945 2.569361 #> RSE 22.722162 22.722162 #> #> $Cov #> k ka V AddErrVar #> k 4.037834e-04 -3.825702e-03 -4.562695e+01 -7.384261e-08 #> ka -3.825702e-03 1.019271e-01 6.136298e+02 1.049592e-06 #> V -4.562695e+01 6.136298e+02 7.506891e+06 1.180135e-02 #> AddErrVar -7.384261e-08 1.049592e-06 1.180135e-02 2.723742e-01 #> #> $run #> $run$m #> [1] 4 #> #> $run$n #> [1] 6 #> #> $run$run #> [1] 6 #> #> $run$p.value #> [1] 0.5952381 #> #> #> $`Objective Function Value` #> [1] 13.22293 #> #> $`-2LL` #> [1] 33.43958 #> #> $AIC #> [1] 41.43958 #> #> $AICc #> [1] 48.10625 #> #> $BIC #> [1] 43.03116 #> #> $Convergence #> NULL #> #> $Message #> [1] "CONVERGENCE: REL_REDUCTION_OF_F <= FACTR*EPSMCH" #> #> $Prediction #> [1] 0.000000 4.181820 6.193476 8.671558 9.929046 9.205385 8.103808 6.796237 #> [9] 5.654668 4.375509 1.467925 #> #> $Residual #> [1] 0.0000000 -2.1618201 -0.5634761 2.7284423 -0.5990462 -0.4653845 #> [7] -0.5438080 0.2937634 0.2453322 -0.0055095 0.1020751 #> #> $`Elapsed Time` #> Time difference of 0.006057024 secs #> #> [1] "## ID = 6 ##" #> $Est #> k ka V AddErrVar AddErrSD CL #> PE 0.09952552 1.1637350 41104.663654 0.2222050 0.4713862 4090.9629 #> SE 0.01617973 0.1850828 3235.606751 0.0947491 0.1005005 420.7268 #> RSE 16.25686601 15.9042025 7.871629 42.6403998 21.3201999 10.2843 #> Thalf MRT #> PE 6.964517 10.047674 #> SE 1.132212 1.633437 #> RSE 16.256866 16.256866 #> #> $Cov #> k ka V AddErrVar #> k 2.617837e-04 -1.988203e-03 -4.509907e+01 -3.223842e-08 #> ka -1.988203e-03 3.425563e-02 4.527839e+02 3.985651e-07 #> V -4.509907e+01 4.527839e+02 1.046915e+07 6.770589e-03 #> AddErrVar -3.223842e-08 3.985651e-07 6.770589e-03 8.977391e-03 #> #> $run #> $run$m #> [1] 5 #> #> $run$n #> [1] 5 #> #> $run$run #> [1] 4 #> #> $run$p.value #> [1] 0.1666667 #> #> #> $`Objective Function Value` #> [1] -5.54577 #> #> $`-2LL` #> [1] 14.67088 #> #> $AIC #> [1] 22.67088 #> #> $AICc #> [1] 29.33754 #> #> $BIC #> [1] 24.26246 #> #> $Convergence #> NULL #> #> $Message #> [1] "CONVERGENCE: REL_REDUCTION_OF_F <= FACTR*EPSMCH" #> #> $Prediction #> [1] 0.0000000 2.0696889 3.7009358 5.3594859 6.1538318 5.8337105 5.1504012 #> [8] 4.2390588 3.4004975 2.5531822 0.7928813 #> #> $Residual #> [1] 0.00000000 -0.77968892 -0.62093582 1.08051409 0.16616820 -0.30371051 #> [7] -0.21040119 -0.21905884 0.05950247 0.22681781 0.12711872 #> #> $`Elapsed Time` #> Time difference of 0.008221149 secs #> #> [1] "## ID = 7 ##" #> $Est #> k ka V AddErrVar AddErrSD CL #> PE 0.10224535 0.67974242 32621.575146 0.09059623 0.30099207 3335.404346 #> SE 0.01161466 0.07256341 1986.491058 0.03863041 0.06417181 207.010940 #> RSE 11.35959330 10.67513321 6.089501 42.64020001 21.32010001 6.206472 #> Thalf MRT #> PE 6.7792539 9.780396 #> SE 0.7700957 1.111013 #> RSE 11.3595933 11.359593 #> #> $Cov #> k ka V AddErrVar #> k 1.349002e-04 -6.670885e-04 -2.128023e+01 -1.436814e-08 #> ka -6.670885e-04 5.265448e-03 1.234222e+02 8.046788e-08 #> V -2.128023e+01 1.234222e+02 3.946147e+06 2.638606e-03 #> AddErrVar -1.436814e-08 8.046787e-08 2.638605e-03 1.492309e-03 #> #> $run #> $run$m #> [1] 4 #> #> $run$n #> [1] 7 #> #> $run$run #> [1] 5 #> #> $run$p.value #> [1] 0.3333333 #> #> #> $`Objective Function Value` #> [1] -15.41478 #> #> $`-2LL` #> [1] 4.801863 #> #> $AIC #> [1] 12.80186 #> #> $AICc #> [1] 19.46853 #> #> $BIC #> [1] 14.39344 #> #> $Convergence #> NULL #> #> $Message #> [1] "CONVERGENCE: REL_REDUCTION_OF_F <= FACTR*EPSMCH" #> #> $Prediction #> [1] 0.0000000 1.5130557 2.7514642 4.6307109 6.4667083 7.0051008 6.5391191 #> [8] 5.5553852 4.5749861 3.3647396 0.9704208 #> #> $Residual #> [1] 0.15000000 -0.66305575 -0.40146424 0.38928907 0.11329169 0.08489924 #> [7] 0.12088088 -0.30538517 -0.18498611 0.16526043 0.17957924 #> #> $`Elapsed Time` #> Time difference of 0.006792068 secs #> #> [1] "## ID = 8 ##" #> $Est #> k ka V AddErrVar AddErrSD CL #> PE 0.09195638 1.3755270 35692.013652 0.3348497 0.578662 3282.10853 #> SE 0.01528989 0.2316089 2770.809865 0.1427802 0.123371 355.81477 #> RSE 16.62732923 16.8378316 7.763109 42.6400962 21.320048 10.84104 #> Thalf MRT #> PE 7.537782 10.874721 #> SE 1.253332 1.808176 #> RSE 16.627329 16.627329 #> #> $Cov #> k ka V AddErrVar #> k 2.337808e-04 -2.274929e-03 -3.597285e+01 -2.495078e-08 #> ka -2.274929e-03 5.364269e-02 4.735930e+02 3.317225e-07 #> V -3.597285e+01 4.735930e+02 7.677387e+06 5.009122e-03 #> AddErrVar -2.495078e-08 3.317225e-07 5.009122e-03 2.038620e-02 #> #> $run #> $run$m #> [1] 5 #> #> $run$n #> [1] 5 #> #> $run$run #> [1] 7 #> #> $run$p.value #> [1] 0.3571429 #> #> #> $`Objective Function Value` #> [1] -1.034796 #> #> $`-2LL` #> [1] 19.18185 #> #> $AIC #> [1] 27.18185 #> #> $AICc #> [1] 33.84852 #> #> $BIC #> [1] 28.77343 #> #> $Convergence #> NULL #> #> $Message #> [1] "CONVERGENCE: REL_REDUCTION_OF_F <= FACTR*EPSMCH" #> #> $Prediction #> [1] 0.000000 2.577418 4.460464 6.284227 7.382246 6.869909 6.029557 4.977828 #> [9] 4.172570 3.157918 1.045607 #> #> $Residual #> [1] 0.0000000 0.4725823 -1.4104642 1.0257728 0.1777542 -0.2799094 #> [7] -0.1495575 -0.2478285 0.3974305 -0.1579175 0.2043934 #> #> $`Elapsed Time` #> Time difference of 0.005850077 secs #> #> [1] "## ID = 9 ##" #> $Est #> k ka V AddErrVar AddErrSD CL #> PE 0.086631747 8.865626 38948.216859 0.22625893 0.4756668 3374.152060 #> SE 0.009328053 2.846690 1504.712606 0.09647691 0.1014123 286.814500 #> RSE 10.767476217 32.109297 3.863367 42.64004647 21.3200232 8.500343 #> Thalf MRT #> PE 8.0010759 11.543113 #> SE 0.8615139 1.242902 #> RSE 10.7674762 10.767476 #> #> $Cov #> k ka V AddErrVar #> k 8.701257e-05 -9.311313e-03 -9.887669e+00 -7.334334e-09 #> ka -9.311313e-03 8.103646e+00 1.962685e+03 7.025506e-07 #> V -9.887669e+00 1.962685e+03 2.264160e+06 1.090413e-03 #> AddErrVar -7.334334e-09 7.025506e-07 1.090413e-03 9.307795e-03 #> #> $run #> $run$m #> [1] 5 #> #> $run$n #> [1] 5 #> #> $run$run #> [1] 6 #> #> $run$p.value #> [1] 0.6428571 #> #> #> $`Objective Function Value` #> [1] -5.346802 #> #> $`-2LL` #> [1] 14.86985 #> #> $AIC #> [1] 22.86985 #> #> $AICc #> [1] 29.53651 #> #> $BIC #> [1] 24.46143 #> #> $Convergence #> NULL #> #> $Message #> [1] "CONVERGENCE: REL_REDUCTION_OF_F <= FACTR*EPSMCH" #> #> $Prediction #> [1] 0.0000000 7.5037042 7.8252750 7.5749377 6.9650972 6.1110441 5.3710120 #> [8] 4.4582579 3.8711373 3.0373318 0.9994943 #> #> $Residual #> [1] 0.0000000 -0.1337042 1.2047250 -0.4349377 -0.6350972 -0.4510441 #> [7] 0.2989880 -0.2182579 0.2388627 0.1226682 0.1205057 #> #> $`Elapsed Time` #> Time difference of 0.007771969 secs #> #> [1] "## ID = 10 ##" #> $Est #> k ka V AddErrVar AddErrSD CL #> PE 0.073965665 0.69550537 25519.76585 0.12285496 0.35050671 1887.586459 #> SE 0.006923395 0.05965454 1127.00909 0.05238562 0.07472841 108.080728 #> RSE 9.360281733 8.57715071 4.41622 42.64021517 21.32010759 5.725869 #> Thalf MRT #> PE 9.3712019 13.519786 #> SE 0.8771709 1.265490 #> RSE 9.3602817 9.360282 #> #> $Cov #> k ka V AddErrVar #> k 4.793339e-05 -3.108263e-04 -7.015432e+00 -1.223972e-08 #> ka -3.108263e-04 3.558665e-03 5.633714e+01 9.243685e-08 #> V -7.015432e+00 5.633714e+01 1.270149e+06 2.133122e-03 #> AddErrVar -1.223972e-08 9.243685e-08 2.133122e-03 2.744253e-03 #> #> $run #> $run$m #> [1] 5 #> #> $run$n #> [1] 6 #> #> $run$run #> [1] 6 #> #> $run$p.value #> [1] 0.7380952 #> #> #> $`Objective Function Value` #> [1] -12.06428 #> #> $`-2LL` #> [1] 8.15237 #> #> $AIC #> [1] 16.15237 #> #> $AICc #> [1] 22.81904 #> #> $BIC #> [1] 17.74395 #> #> $Convergence #> NULL #> #> $Message #> [1] "CONVERGENCE: REL_REDUCTION_OF_F <= FACTR*EPSMCH" #> #> $Prediction #> [1] 0.000000 2.804848 5.041307 6.109314 8.685332 9.603184 9.239402 8.209423 #> [9] 6.990600 5.730370 2.431042 #> #> $Residual #> [1] 0.24000000 0.08515194 0.17869281 0.30068631 -0.85533217 0.60681616 #> [7] -0.05940207 -0.18942298 0.14939969 -0.05037035 -0.01104203 #> #> $`Elapsed Time` #> Time difference of 0.008107901 secs #> #> [1] "## ID = 11 ##" #> $Est #> k ka V AddErrVar AddErrSD CL #> PE 0.098123028 3.8490542 37945.347628 0.03874697 0.19684251 3723.312416 #> SE 0.004603223 0.2483307 768.521663 0.01652178 0.04196701 126.928376 #> RSE 4.691277227 6.4517342 2.025338 42.64019331 21.32009666 3.409018 #> Thalf MRT #> PE 7.0640623 10.1912876 #> SE 0.3313947 0.4781016 #> RSE 4.6912772 4.6912772 #> #> $Cov #> k ka V AddErrVar #> k 2.118966e-05 -5.781500e-04 -2.697293e+00 -1.807690e-09 #> ka -5.781500e-04 6.166816e-02 1.233536e+02 7.806698e-08 #> V -2.697293e+00 1.233536e+02 5.906255e+05 3.550745e-04 #> AddErrVar -1.807690e-09 7.806697e-08 3.550745e-04 2.729694e-04 #> #> $run #> $run$m #> [1] 4 #> #> $run$n #> [1] 6 #> #> $run$run #> [1] 4 #> #> $run$p.value #> [1] 0.1904762 #> #> #> $`Objective Function Value` #> [1] -24.75774 #> #> $`-2LL` #> [1] -4.541094 #> #> $AIC #> [1] 3.458906 #> #> $AICc #> [1] 10.12557 #> #> $BIC #> [1] 5.050487 #> #> $Convergence #> NULL #> #> $Message #> [1] "CONVERGENCE: REL_REDUCTION_OF_F <= FACTR*EPSMCH" #> #> $Prediction #> [1] 0.0000000 5.1381059 6.9764937 7.6613204 7.1215015 6.0784708 5.2878927 #> [8] 4.3413812 3.5677904 2.6346394 0.8148061 #> #> $Residual #> [1] 0.00000000 -0.27810591 0.26350627 0.33867961 -0.31150152 -0.20847084 #> [7] -0.06789273 0.10861884 0.05220955 0.05536060 0.04519389 #> #> $`Elapsed Time` #> Time difference of 0.007643938 secs #> #> [1] "## ID = 12 ##" #> $Est #> k ka V AddErrVar AddErrSD CL #> PE 0.10557578 0.83291039 24017.514932 0.2553549 0.5053266 2535.667942 #> SE 0.01289139 0.09971057 1557.293843 0.1088724 0.1077248 176.024973 #> RSE 12.21055245 11.97134456 6.483992 42.6356972 21.3178486 6.941957 #> Thalf MRT #> PE 6.5653994 9.471869 #> SE 0.8016715 1.156568 #> RSE 12.2105525 12.210552 #> #> $Cov #> k ka V AddErrVar #> k 1.661878e-04 -9.662099e-04 -1.812356e+01 4.339652e-09 #> ka -9.662099e-04 9.942198e-03 1.276843e+02 4.984468e-07 #> V -1.812356e+01 1.276843e+02 2.425164e+06 9.254574e-04 #> AddErrVar 4.339652e-09 4.984468e-07 9.254574e-04 1.185319e-02 #> #> $run #> $run$m #> [1] 5 #> #> $run$n #> [1] 5 #> #> $run$run #> [1] 5 #> #> $run$p.value #> [1] 0.3571429 #> #> #> $`Objective Function Value` #> [1] -4.014962 #> #> $`-2LL` #> [1] 16.20169 #> #> $AIC #> [1] 24.20169 #> #> $AICc #> [1] 30.86835 #> #> $BIC #> [1] 25.79327 #> #> $Convergence #> NULL #> #> $Message #> [1] "CONVERGENCE: REL_REDUCTION_OF_F <= FACTR*EPSMCH" #> #> $Prediction #> [1] 0.000000 2.470648 4.412512 7.095149 9.469103 9.708602 8.709870 7.190744 #> [9] 5.872744 4.274780 1.191749 #> #> $Residual #> [1] 0.00000000 -1.22064823 -0.45251234 0.72485114 0.25089707 0.04139793 #> [7] -0.13987006 -0.60074407 0.23725635 0.29521982 -0.02174890 #> #> $`Elapsed Time` #> Time difference of 0.008006096 secs #>