B Theoph의 nlr의 전체 실행 결과
5장 구획분석의 자료해석의 5.7절에서 언급된 wnl::nlr() 실행결과를 별첨으로 수록하였다.
library(wnl)
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],]
Theoph_nlr = 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(Theoph_nlr)
}## [1] "## ID = 1 ##"
## $Est
## k ka V AddErrVar AddErrSD
## PE 0.053954 1.777 29394.222 0.389637 0.6242
## SE 0.007784 0.230 1449.530 0.166142 0.1331
## RSE 14.427031 12.940 4.931 42.640142 21.3201
## LL 0.035548 1.234 25966.628 -0.003226 0.3095
## UL 0.072361 2.321 32821.816 0.782500 0.9389
## CL Thalf MRT
## PE 1585.9 12.847 18.534
## SE 172.9 1.853 2.674
## RSE 10.9 14.427 14.427
## LL 1177.2 8.464 12.211
## UL 1994.7 17.230 24.857
##
## $LI
## k ka V AddErrVar
## LL 0.04471 1.443 27826 0.1784
## UL 0.06448 2.200 31150 1.1189
## attr(,"k")
## [1] 9.18
##
## $Skewness
## [1] 0.2770 0.5231 0.0884
##
## $Cov
## k ka V AddErrVar
## k 6.059e-05 -9.28e-04 -9.014e+00 -1.665e-08
## ka -9.280e-04 5.29e-02 2.120e+02 3.840e-07
## V -9.014e+00 2.12e+02 2.101e+06 3.777e-03
## AddErrVar -1.665e-08 3.84e-07 3.777e-03 2.760e-02
##
## $run
## $run$m
## [1] 4
##
## $run$n
## [1] 7
##
## $run$run
## [1] 5
##
## $run$p.value
## [1] 0.3333
##
##
## $`Objective Function Value`
## [1] 0.6321
##
## $`-2LL`
## [1] 20.85
##
## $AIC
## [1] 28.85
##
## $AICc
## [1] 35.52
##
## $BIC
## [1] 30.44
##
## $Convergence
## NULL
##
## $Message
## [1] "CONVERGENCE: REL_REDUCTION_OF_F <= FACTR*EPSMCH"
##
## $Prediction
## [1] 0.000 3.878 6.811 9.035 9.758 9.124 8.525 7.683
## [9] 6.890 5.838 3.015
##
## $Residual
## [1] 0.74000 -1.08659 -0.32707 1.35032 -0.22177
## [6] -0.65904 -0.27314 -0.31052 -0.08715 0.02790
## [11] 0.22720
##
## $`Elapsed Time`
## Time difference of 0.03354 secs
##
## [1] "## ID = 2 ##"
## $Est
## k ka V AddErrVar AddErrSD
## PE 0.10166 1.9427 32024.77 0.813479 0.9019
## SE 0.02080 0.4307 3064.68 0.346868 0.1923
## RSE 20.45662 22.1727 9.57 42.639999 21.3200
## LL 0.05249 0.9241 24777.97 -0.006732 0.4472
## UL 0.15084 2.9612 39271.58 1.633691 1.3566
## CL Thalf MRT
## PE 3255.67 6.818 9.837
## SE 451.31 1.395 2.012
## RSE 13.86 20.457 20.457
## LL 2188.50 3.520 5.078
## UL 4322.84 10.116 14.595
##
## $LI
## k ka V AddErrVar
## LL 0.07872 1.394 29029 0.3724
## UL 0.13044 2.724 35710 2.3360
## attr(,"k")
## [1] 9.18
##
## $Skewness
## [1] 0.5245 0.8626 0.1441
##
## $Cov
## k ka V AddErrVar
## k 4.325e-04 -5.235e-03 -5.175e+01 -2.618e-08
## ka -5.235e-03 1.855e-01 9.264e+02 1.126e-07
## V -5.175e+01 9.264e+02 9.392e+06 4.885e-03
## AddErrVar -2.618e-08 1.126e-07 4.885e-03 1.203e-01
##
## $run
## $run$m
## [1] 4
##
## $run$n
## [1] 6
##
## $run$run
## [1] 5
##
## $run$p.value
## [1] 0.4048
##
##
## $`Objective Function Value`
## [1] 8.729
##
## $`-2LL`
## [1] 28.95
##
## $AIC
## [1] 36.95
##
## $AICc
## [1] 43.61
##
## $BIC
## [1] 38.54
##
## $Convergence
## NULL
##
## $Message
## [1] "CONVERGENCE: REL_REDUCTION_OF_F <= FACTR*EPSMCH"
##
## $Prediction
## [1] 0.0000 4.0182 6.1615 8.0136 8.4214 7.3754 6.3289
## [8] 5.1597 4.2233 3.1131 0.8915
##
## $Residual
## [1] 0.000000 -2.349039 1.670514 0.194941 -0.197959
## [6] -0.618799 -0.329009 0.174970 0.273261 -0.142538
## [11] -0.002829
##
## $`Elapsed Time`
## Time difference of 0.006591 secs
##
## [1] "## ID = 3 ##"
## $Est
## k ka V AddErrVar AddErrSD
## PE 0.08142 2.4536 34319.365 0.0396612 0.19915
## SE 0.00381 0.1437 684.952 0.0169116 0.04246
## RSE 4.67902 5.8575 1.996 42.6400576 21.32003
## LL 0.07242 2.1137 32699.712 -0.0003283 0.09875
## UL 0.09043 2.7934 35939.019 0.0796507 0.29955
## CL Thalf MRT
## PE 2794.442 8.5127 12.2813
## SE 92.853 0.3983 0.5746
## RSE 3.323 4.6790 4.6790
## LL 2574.881 7.5709 10.9225
## UL 3014.004 9.4546 13.6401
##
## $LI
## k ka V AddErrVar
## LL 0.07666 2.248 33573 0.01815
## UL 0.08643 2.684 35099 0.11389
## attr(,"k")
## [1] 9.18
##
## $Skewness
## [1] 0.10515 0.22728 0.03572
##
## $Cov
## k ka V AddErrVar
## k 1.452e-05 -3.094e-04 -2.073e+00 -2.225e-09
## ka -3.094e-04 2.066e-02 6.848e+01 6.716e-08
## V -2.073e+00 6.848e+01 4.692e+05 4.891e-04
## AddErrVar -2.225e-09 6.716e-08 4.891e-04 2.860e-04
##
## $run
## $run$m
## [1] 4
##
## $run$n
## [1] 6
##
## $run$run
## [1] 7
##
## $run$p.value
## [1] 0.3095
##
##
## $`Objective Function Value`
## [1] -24.5
##
## $`-2LL`
## [1] -4.285
##
## $AIC
## [1] 3.715
##
## $AICc
## [1] 10.38
##
## $BIC
## [1] 5.307
##
## $Convergence
## NULL
##
## $Message
## [1] "CONVERGENCE: REL_REDUCTION_OF_F <= FACTR*EPSMCH"
##
## $Prediction
## [1] 0.000 4.462 6.875 8.086 8.114 7.181 6.377 5.423
## [9] 4.635 3.586 1.348
##
## $Residual
## [1] 0.00000 -0.11865 -0.06240 0.01159 -0.41637
## [6] 0.22826 -0.25785 -0.19186 0.20679 0.06856
## [11] -0.31466
##
## $`Elapsed Time`
## Time difference of 0.006918 secs
##
## [1] "## ID = 4 ##"
## $Est
## k ka V AddErrVar AddErrSD
## PE 0.08747 1.1715 31097.470 0.521087 0.7219
## SE 0.01605 0.2036 2607.291 0.222192 0.1539
## RSE 18.35204 17.3777 8.384 42.640172 21.3201
## LL 0.04951 0.6901 24932.207 -0.004315 0.3579
## UL 0.12542 1.6529 37262.733 1.046488 1.0858
## CL Thalf MRT
## PE 2719.99 7.925 11.433
## SE 328.29 1.454 2.098
## RSE 12.07 18.352 18.352
## LL 1943.70 4.486 6.472
## UL 3496.28 11.364 16.394
##
## $LI
## k ka V AddErrVar
## LL 0.07108 0.9143 28755 0.2385
## UL 0.10698 1.4968 33856 1.4963
## attr(,"k")
## [1] 9.18
##
## $Skewness
## [1] 0.53716 0.56163 0.03693
##
## $Cov
## k ka V AddErrVar
## k 2.577e-04 -2.094e-03 -3.555e+01 -3.662e-08
## ka -2.094e-03 4.144e-02 3.932e+02 4.001e-07
## V -3.555e+01 3.932e+02 6.798e+06 6.917e-03
## AddErrVar -3.662e-08 4.001e-07 6.917e-03 4.937e-02
##
## $run
## $run$m
## [1] 3
##
## $run$n
## [1] 7
##
## $run$run
## [1] 5
##
## $run$p.value
## [1] 0.7167
##
##
## $`Objective Function Value`
## [1] 3.83
##
## $`-2LL`
## [1] 24.05
##
## $AIC
## [1] 32.05
##
## $AICc
## [1] 38.71
##
## $BIC
## [1] 33.64
##
## $Convergence
## NULL
##
## $Message
## [1] "CONVERGENCE: REL_REDUCTION_OF_F <= FACTR*EPSMCH"
##
## $Prediction
## [1] 0.000 3.405 5.046 6.952 8.313 8.004 7.138 6.015
## [9] 5.052 3.900 1.288
##
## $Residual
## [1] 0.00000 -1.55835 -0.50948 1.56002 -0.03833
## [6] -0.56496 -0.34789 -0.31129 0.21403 0.24077
## [11] -0.15385
##
## $`Elapsed Time`
## Time difference of 0.009309 secs
##
## [1] "## ID = 5 ##"
## $Est
## k ka V AddErrVar AddErrSD
## PE 0.08844 1.4715 26925.05 1.22395 1.1063
## SE 0.02009 0.3193 2739.87 0.52189 0.2359
## RSE 22.72216 21.6962 10.18 42.64015 21.3201
## LL 0.04092 0.7166 20446.29 -0.01013 0.5486
## UL 0.13595 2.2264 33403.82 2.45804 1.6641
## CL Thalf MRT
## PE 2381.12 7.838 11.308
## SE 366.26 1.781 2.569
## RSE 15.38 22.722 22.722
## LL 1515.04 3.627 5.232
## UL 3247.20 12.049 17.383
##
## $LI
## k ka V AddErrVar
## LL 0.06723 1.057 24336 0.5603
## UL 0.11513 2.036 30131 3.5147
## attr(,"k")
## [1] 9.18
##
## $Skewness
## [1] 0.5985 0.7966 0.1248
##
## $Cov
## k ka V AddErrVar
## k 4.038e-04 -3.826e-03 -4.563e+01 -7.384e-08
## ka -3.826e-03 1.019e-01 6.136e+02 1.050e-06
## V -4.563e+01 6.136e+02 7.507e+06 1.180e-02
## AddErrVar -7.384e-08 1.050e-06 1.180e-02 2.724e-01
##
## $run
## $run$m
## [1] 4
##
## $run$n
## [1] 6
##
## $run$run
## [1] 6
##
## $run$p.value
## [1] 0.5952
##
##
## $`Objective Function Value`
## [1] 13.22
##
## $`-2LL`
## [1] 33.44
##
## $AIC
## [1] 41.44
##
## $AICc
## [1] 48.11
##
## $BIC
## [1] 43.03
##
## $Convergence
## NULL
##
## $Message
## [1] "CONVERGENCE: REL_REDUCTION_OF_F <= FACTR*EPSMCH"
##
## $Prediction
## [1] 0.000 4.182 6.193 8.672 9.929 9.205 8.104 6.796
## [9] 5.655 4.376 1.468
##
## $Residual
## [1] 0.00000 -2.21475 -0.64187 2.61868 -0.72473
## [6] -0.58191 -0.64639 0.20774 0.17375 -0.06090
## [11] 0.08349
##
## $`Elapsed Time`
## Time difference of 0.006563 secs
##
## [1] "## ID = 6 ##"
## $Est
## k ka V AddErrVar AddErrSD
## PE 0.09953 1.1637 41104.664 0.222205 0.4714
## SE 0.01618 0.1851 3235.607 0.094749 0.1005
## RSE 16.25687 15.9042 7.872 42.640400 21.3202
## LL 0.06127 0.7261 33453.669 -0.001841 0.2337
## UL 0.13778 1.6014 48755.658 0.446251 0.7090
## CL Thalf MRT
## PE 4090.96 6.965 10.048
## SE 420.73 1.132 1.633
## RSE 10.28 16.257 16.257
## LL 3096.10 4.287 6.185
## UL 5085.82 9.642 13.910
##
## $LI
## k ka V AddErrVar
## LL 0.08349 0.9343 38286 0.1017
## UL 0.11828 1.4495 44371 0.6381
## attr(,"k")
## [1] 9.18
##
## $Skewness
## [1] 0.50175 0.52432 0.01456
##
## $Cov
## k ka V AddErrVar
## k 2.618e-04 -1.988e-03 -4.510e+01 -3.224e-08
## ka -1.988e-03 3.426e-02 4.528e+02 3.986e-07
## V -4.510e+01 4.528e+02 1.047e+07 6.771e-03
## AddErrVar -3.224e-08 3.986e-07 6.771e-03 8.977e-03
##
## $run
## $run$m
## [1] 5
##
## $run$n
## [1] 5
##
## $run$run
## [1] 4
##
## $run$p.value
## [1] 0.1667
##
##
## $`Objective Function Value`
## [1] -5.546
##
## $`-2LL`
## [1] 14.67
##
## $AIC
## [1] 22.67
##
## $AICc
## [1] 29.34
##
## $BIC
## [1] 24.26
##
## $Convergence
## NULL
##
## $Message
## [1] "CONVERGENCE: REL_REDUCTION_OF_F <= FACTR*EPSMCH"
##
## $Prediction
## [1] 0.0000 2.0697 3.7009 5.3595 6.1538 5.8337 5.1504
## [8] 4.2391 3.4005 2.5532 0.7929
##
## $Residual
## [1] 0.00000 -0.80589 -0.66778 1.01267 0.08827
## [6] -0.37755 -0.27560 -0.27272 0.01646 0.19450
## [11] 0.11708
##
## $`Elapsed Time`
## Time difference of 0.006665 secs
##
## [1] "## ID = 7 ##"
## $Est
## k ka V AddErrVar AddErrSD
## PE 0.10225 0.67974 32621.58 0.0905962 0.30099
## SE 0.01161 0.07256 1986.49 0.0386304 0.06417
## RSE 11.35959 10.67513 6.09 42.6402000 21.32010
## LL 0.07478 0.50816 27924.27 -0.0007502 0.14925
## UL 0.12971 0.85133 37318.88 0.1819426 0.45273
## CL Thalf MRT
## PE 3335.404 6.7793 9.780
## SE 207.011 0.7701 1.111
## RSE 6.206 11.3596 11.360
## LL 2845.901 4.9583 7.153
## UL 3824.907 8.6002 12.408
##
## $LI
## k ka V AddErrVar
## LL 0.09314 0.6051 31310 0.04149
## UL 0.11212 0.7633 34048 0.26016
## attr(,"k")
## [1] 9.18
##
## $Skewness
## [1] 0.4748 0.1809 -0.1357
##
## $Cov
## k ka V AddErrVar
## k 1.349e-04 -6.671e-04 -2.128e+01 -1.437e-08
## ka -6.671e-04 5.265e-03 1.234e+02 8.047e-08
## V -2.128e+01 1.234e+02 3.946e+06 2.639e-03
## AddErrVar -1.437e-08 8.047e-08 2.639e-03 1.492e-03
##
## $run
## $run$m
## [1] 5
##
## $run$n
## [1] 6
##
## $run$run
## [1] 7
##
## $run$p.value
## [1] 0.4784
##
##
## $`Objective Function Value`
## [1] -15.41
##
## $`-2LL`
## [1] 4.802
##
## $AIC
## [1] 12.8
##
## $AICc
## [1] 19.47
##
## $BIC
## [1] 14.39
##
## $Convergence
## NULL
##
## $Message
## [1] "CONVERGENCE: REL_REDUCTION_OF_F <= FACTR*EPSMCH"
##
## $Prediction
## [1] 0.0000 1.5131 2.7515 4.6307 6.4667 7.0051 6.5391
## [8] 5.5554 4.5750 3.3647 0.9704
##
## $Residual
## [1] 0.150000 -0.682208 -0.436293 0.330672 0.031435
## [6] -0.003773 0.038107 -0.375707 -0.242897 0.122669
## [11] 0.167295
##
## $`Elapsed Time`
## Time difference of 0.006954 secs
##
## [1] "## ID = 8 ##"
## $Est
## k ka V AddErrVar AddErrSD
## PE 0.09196 1.3755 35692.014 0.334850 0.5787
## SE 0.01529 0.2316 2770.810 0.142780 0.1234
## RSE 16.62733 16.8378 7.763 42.640096 21.3200
## LL 0.05580 0.8279 29140.089 -0.002772 0.2869
## UL 0.12811 1.9232 42243.938 0.672471 0.8704
## CL Thalf MRT
## PE 3282.11 7.538 10.875
## SE 355.81 1.253 1.808
## RSE 10.84 16.627 16.627
## LL 2440.74 4.574 6.599
## UL 4123.48 10.501 15.150
##
## $LI
## k ka V AddErrVar
## LL 0.07624 1.082 33172 0.1533
## UL 0.11043 1.749 38626 0.9616
## attr(,"k")
## [1] 9.18
##
## $Skewness
## [1] 0.44072 0.55360 0.07714
##
## $Cov
## k ka V AddErrVar
## k 2.338e-04 -2.275e-03 -3.597e+01 -2.495e-08
## ka -2.275e-03 5.364e-02 4.736e+02 3.317e-07
## V -3.597e+01 4.736e+02 7.677e+06 5.009e-03
## AddErrVar -2.495e-08 3.317e-07 5.009e-03 2.039e-02
##
## $run
## $run$m
## [1] 5
##
## $run$n
## [1] 5
##
## $run$run
## [1] 7
##
## $run$p.value
## [1] 0.3571
##
##
## $`Objective Function Value`
## [1] -1.035
##
## $`-2LL`
## [1] 19.18
##
## $AIC
## [1] 27.18
##
## $AICc
## [1] 33.85
##
## $BIC
## [1] 28.77
##
## $Convergence
## NULL
##
## $Message
## [1] "CONVERGENCE: REL_REDUCTION_OF_F <= FACTR*EPSMCH"
##
## $Prediction
## [1] 0.000 2.577 4.460 6.284 7.382 6.870 6.030 4.978
## [9] 4.173 3.158 1.046
##
## $Residual
## [1] 0.00000 0.43996 -1.46693 0.94623 0.08431
## [6] -0.36687 -0.22588 -0.31084 0.34461 -0.19789
## [11] 0.19116
##
## $`Elapsed Time`
## Time difference of 0.009353 secs
##
## [1] "## ID = 9 ##"
## $Est
## k ka V AddErrVar AddErrSD
## PE 0.086632 8.866 38948.217 0.226259 0.4757
## SE 0.009328 2.847 1504.713 0.096477 0.1014
## RSE 10.767476 32.109 3.863 42.640046 21.3200
## LL 0.064574 2.134 35390.137 -0.001873 0.2359
## UL 0.108689 15.597 42506.297 0.454391 0.7155
## CL Thalf MRT
## PE 3374.2 8.0011 11.543
## SE 286.8 0.8615 1.243
## RSE 8.5 10.7675 10.767
## LL 2695.9 5.9639 8.604
## UL 4052.4 10.0382 14.482
##
## $LI
## k ka V AddErrVar
## LL 0.07359 5.721 36927 0.1036
## UL 0.10157 Inf 41204 0.6497
## attr(,"k")
## [1] 9.18
##
## $Skewness
## [1] 0.2066 3.0078 0.1045
##
## $Cov
## k ka V AddErrVar
## k 8.701e-05 -9.311e-03 -9.888e+00 -7.334e-09
## ka -9.311e-03 8.104e+00 1.963e+03 7.026e-07
## V -9.888e+00 1.963e+03 2.264e+06 1.090e-03
## AddErrVar -7.334e-09 7.026e-07 1.090e-03 9.308e-03
##
## $run
## $run$m
## [1] 5
##
## $run$n
## [1] 5
##
## $run$run
## [1] 6
##
## $run$p.value
## [1] 0.6429
##
##
## $`Objective Function Value`
## [1] -5.347
##
## $`-2LL`
## [1] 14.87
##
## $AIC
## [1] 22.87
##
## $AICc
## [1] 29.54
##
## $BIC
## [1] 24.46
##
## $Convergence
## NULL
##
## $Message
## [1] "CONVERGENCE: REL_REDUCTION_OF_F <= FACTR*EPSMCH"
##
## $Prediction
## [1] 0.0000 7.5037 7.8253 7.5749 6.9651 6.1110 5.3710
## [8] 4.4583 3.8711 3.0373 0.9995
##
## $Residual
## [1] 0.00000 -0.22869 1.10567 -0.53082 -0.72326
## [6] -0.52840 0.23100 -0.27469 0.18986 0.08422
## [11] 0.10785
##
## $`Elapsed Time`
## Time difference of 0.007591 secs
##
## [1] "## ID = 10 ##"
## $Est
## k ka V AddErrVar AddErrSD
## PE 0.073966 0.69551 25519.766 0.122855 0.35051
## SE 0.006923 0.05965 1127.009 0.052386 0.07473
## RSE 9.360282 8.57715 4.416 42.640215 21.32011
## LL 0.057594 0.55444 22854.813 -0.001017 0.17380
## UL 0.090337 0.83657 28184.719 0.246727 0.52721
## CL Thalf MRT
## PE 1887.586 9.3712 13.520
## SE 108.081 0.8772 1.265
## RSE 5.726 9.3603 9.360
## LL 1632.016 7.2970 10.527
## UL 2143.157 11.4454 16.512
##
## $LI
## k ka V AddErrVar
## LL 0.06778 0.6302 24687 0.05626
## UL 0.08059 0.7676 26411 0.35280
## attr(,"k")
## [1] 9.18
##
## $Skewness
## [1] 0.28945 0.18789 -0.03627
##
## $Cov
## k ka V AddErrVar
## k 4.793e-05 -3.108e-04 -7.015e+00 -1.224e-08
## ka -3.108e-04 3.559e-03 5.634e+01 9.244e-08
## V -7.015e+00 5.634e+01 1.270e+06 2.133e-03
## AddErrVar -1.224e-08 9.244e-08 2.133e-03 2.744e-03
##
## $run
## $run$m
## [1] 5
##
## $run$n
## [1] 6
##
## $run$run
## [1] 6
##
## $run$p.value
## [1] 0.7381
##
##
## $`Objective Function Value`
## [1] -12.06
##
## $`-2LL`
## [1] 8.152
##
## $AIC
## [1] 16.15
##
## $AICc
## [1] 22.82
##
## $BIC
## [1] 17.74
##
## $Convergence
## NULL
##
## $Message
## [1] "CONVERGENCE: REL_REDUCTION_OF_F <= FACTR*EPSMCH"
##
## $Prediction
## [1] 0.000 2.805 5.041 6.109 8.685 9.603 9.239 8.209
## [9] 6.991 5.730 2.431
##
## $Residual
## [1] 0.24000 0.04965 0.11488 0.22335 -0.96527
## [6] 0.48526 -0.17636 -0.29334 0.06091 -0.12291
## [11] -0.04181
##
## $`Elapsed Time`
## Time difference of 0.007993 secs
##
## [1] "## ID = 11 ##"
## $Est
## k ka V AddErrVar AddErrSD
## PE 0.098123 3.8491 37945.348 0.0387470 0.19684
## SE 0.004603 0.2483 768.522 0.0165218 0.04197
## RSE 4.691277 6.4517 2.025 42.6401933 21.32010
## LL 0.087238 3.2618 36128.083 -0.0003208 0.09761
## UL 0.109008 4.4363 39762.613 0.0778148 0.29608
## CL Thalf MRT
## PE 3723.312 7.0641 10.1913
## SE 126.928 0.3314 0.4781
## RSE 3.409 4.6913 4.6913
## LL 3423.174 6.2804 9.0608
## UL 4023.450 7.8477 11.3218
##
## $LI
## k ka V AddErrVar
## LL 0.09201 3.477 37042 0.01774
## UL 0.10459 4.279 38894 0.11127
## attr(,"k")
## [1] 9.18
##
## $Skewness
## [1] 0.10154 0.29127 0.04615
##
## $Cov
## k ka V AddErrVar
## k 2.119e-05 -5.782e-04 -2.697e+00 -1.808e-09
## ka -5.782e-04 6.167e-02 1.234e+02 7.807e-08
## V -2.697e+00 1.234e+02 5.906e+05 3.551e-04
## AddErrVar -1.808e-09 7.807e-08 3.551e-04 2.730e-04
##
## $run
## $run$m
## [1] 4
##
## $run$n
## [1] 6
##
## $run$run
## [1] 4
##
## $run$p.value
## [1] 0.1905
##
##
## $`Objective Function Value`
## [1] -24.76
##
## $`-2LL`
## [1] -4.541
##
## $AIC
## [1] 3.459
##
## $AICc
## [1] 10.13
##
## $BIC
## [1] 5.05
##
## $Convergence
## NULL
##
## $Message
## [1] "CONVERGENCE: REL_REDUCTION_OF_F <= FACTR*EPSMCH"
##
## $Prediction
## [1] 0.0000 5.1381 6.9765 7.6613 7.1215 6.0785 5.2879
## [8] 4.3414 3.5678 2.6346 0.8148
##
## $Residual
## [1] 0.000000 -0.343145 0.175196 0.241701 -0.401647
## [6] -0.285414 -0.134828 0.053665 0.007048 0.022011
## [11] 0.034880
##
## $`Elapsed Time`
## Time difference of 0.01042 secs
##
## [1] "## ID = 12 ##"
## $Est
## k ka V AddErrVar AddErrSD
## PE 0.10558 0.83291 24017.515 0.255355 0.5053
## SE 0.01289 0.09971 1557.294 0.108872 0.1077
## RSE 12.21055 11.97134 6.484 42.635697 21.3178
## LL 0.07509 0.59713 20335.100 -0.002087 0.2506
## UL 0.13606 1.06869 27699.930 0.512797 0.7601
## CL Thalf MRT
## PE 2535.668 6.5654 9.472
## SE 176.025 0.8017 1.157
## RSE 6.942 12.2106 12.211
## LL 2119.435 4.6697 6.737
## UL 2951.901 8.4611 12.207
##
## $LI
## k ka V AddErrVar
## LL 0.09449 0.7212 22860 0.1169
## UL 0.11788 0.9608 25299 0.7334
## attr(,"k")
## [1] 9.18
##
## $Skewness
## [1] 0.46551 0.26478 -0.07838
##
## $Cov
## k ka V AddErrVar
## k 1.662e-04 -9.662e-04 -1.812e+01 4.340e-09
## ka -9.662e-04 9.942e-03 1.277e+02 4.984e-07
## V -1.812e+01 1.277e+02 2.425e+06 9.255e-04
## AddErrVar 4.340e-09 4.984e-07 9.255e-04 1.185e-02
##
## $run
## $run$m
## [1] 4
##
## $run$n
## [1] 6
##
## $run$run
## [1] 5
##
## $run$p.value
## [1] 0.4048
##
##
## $`Objective Function Value`
## [1] -4.015
##
## $`-2LL`
## [1] 16.2
##
## $AIC
## [1] 24.2
##
## $AICc
## [1] 30.87
##
## $BIC
## [1] 25.79
##
## $Convergence
## NULL
##
## $Message
## [1] "CONVERGENCE: REL_REDUCTION_OF_F <= FACTR*EPSMCH"
##
## $Prediction
## [1] 0.000 2.471 4.413 7.095 9.469 9.709 8.710 7.191
## [9] 5.873 4.275 1.192
##
## $Residual
## [1] 0.00000 -1.25192 -0.50837 0.63504 0.13104
## [6] -0.08150 -0.25012 -0.69177 0.16292 0.24111
## [11] -0.03683
##
## $`Elapsed Time`
## Time difference of 0.006656 secs