Arguments for VPC runs
NlmeVpcParams-class.RdClass initializer for arguments of visual predictive check (VPC) runs
Arguments
- numReplicates
Integer; Number of replicates to simulate the model
- seed
Integer; Random number generator seed
- predCorrection
Character; Type of correction to use when calculating a prediction-corrected observation. Options are
"none", "proportional", "additive". This option is ignored for discontinuous observed variables (categorical, count, and time-to-event).- predVarCorr
Logical; Set to
TRUEto use Prediction Variance Correction. Only applicable to the case wherepredCorrectionis set to either"proportional"or"additive".- outputPRED
Logical; Set to
TRUEto include population prediction (PRED) results for continuous observed variables in output.- stratifyColumns
Character or character vector; Names of categorical covariates (up to 3) used to stratify modeling simulation results.
- observationVars
NlmeObservationVar class instance or list of these instances
- simulationTables
Optional list of simulation tables.
NlmeSimTableDefclass instance or a list of such instances. Could be generated bytableParamswrapper function or byNlmeSimTableDefclass instance initializing directly.
Examples
# \donttest{
model <- pkmodel(
parameterization = "Clearance",
numCompartments = 2,
data = pkData,
ID = "Subject",
Time = "Act_Time",
A1 = "Amount",
CObs = "Conc",
workingDir = tempdir()
)
# Define the host
host <- hostParams(parallelMethod = "NONE",
hostName = "local",
numCores = 1,
sharedDirectory = tempdir())
job <- fitmodel(model,
hostPlatform = host)
#>
#> NLME Job
#>
#> Compiling 1 of 1 NLME models
#> TDL5 version: 25.7.1.1
#>
#> Status: OK
#> License expires: 2026-08-06
#> Refresh until: 2025-09-05 07:01:25
#> Current Date: 2025-08-06
#> The model compiled
#>
#>
#> Iteration -2LL tvV tvCl tvV2 tvCl2 nSubj nObs
#> 1 1734.38 1.07796 1.05372 1.10597 1.07380 16 112
#> 2 1640.76 1.19901 1.16447 1.21074 1.19472 16 112
#> 3 1587.09 1.35133 1.31271 1.35521 1.34724 16 112
#> 4 1559.94 1.51691 1.47592 1.51570 1.51267 16 112
#> 5 1545.53 1.69253 1.65038 1.68761 1.68812 16 112
#> 6 1537.46 1.87322 1.83076 1.86587 1.86828 16 112
#> 7 1532.81 2.04832 2.00610 2.03978 2.04280 16 112
#> 8 1530.43 2.19847 2.15746 2.19006 2.19264 16 112
#> 9 1529.44 2.28870 2.24956 2.28142 2.28275 16 112
#> 10 1528.44 2.36768 2.33174 2.36244 2.36178 16 112
#> 11 1525.84 2.52046 2.49445 2.52134 2.51486 16 112
#> 12 1519.67 2.80483 2.80466 2.82032 2.80038 16 112
#> 13 1515.10 3.21194 3.25592 3.24936 3.20963 16 112
#> 14 1510.03 3.05238 3.08750 3.07917 3.04942 16 112
#> 15 1504.03 3.18762 3.25244 3.21772 3.18563 16 112
#> 16 1503.29 3.75212 3.93610 3.79491 3.75353 16 112
#> 17 1495.99 3.57164 3.73550 3.60451 3.57215 16 112
#> 18 1493.74 3.72773 3.93150 3.76212 3.72926 16 112
#> 19 1486.86 4.46486 4.87109 4.50720 4.47162 16 112
#> 20 1482.33 4.54962 4.99801 4.59027 4.55754 16 112
#> 21 1465.23 5.63303 6.49837 5.67268 5.65262 16 112
#> 22 1449.27 6.59887 7.94730 6.63178 6.63436 16 112
#> 23 1448.20 6.63453 8.02035 6.66403 6.67397 16 112
#> 24 1447.20 6.48840 7.81175 6.50995 6.53610 16 112
#> 25 1445.80 6.66308 8.03269 6.67198 6.72787 16 112
#> 26 1443.68 7.05368 8.50609 7.03508 7.15239 16 112
#> 27 1442.28 7.73989 9.30359 7.67427 7.88704 16 112
#> 28 1439.90 8.00197 9.57586 7.91965 8.15561 16 112
#> 29 1438.04 8.71977 10.37810 8.58723 8.90172 16 112
#> 30 1433.56 10.51470 12.38290 10.24620 10.76480 16 112
#> 31 1430.80 13.20840 15.39390 12.73360 13.55560 16 112
#> 32 1426.24 13.39910 15.55390 12.89450 13.73680 16 112
#> 33 1423.81 13.37820 15.51600 12.86180 13.71630 16 112
#> 34 1420.55 13.70890 15.83190 13.11970 14.05080 16 112
#> 35 1419.51 13.04010 15.05450 12.47940 13.35950 16 112
#> 36 1418.94 13.20170 15.18740 12.60690 13.52300 16 112
#> 37 1418.52 13.34050 15.25530 12.69030 13.65650 16 112
#> 38 1418.36 13.37370 15.22970 12.68750 13.68070 16 112
#> 39 1418.06 13.40110 15.13910 12.64770 13.68500 16 112
#> 40 1417.51 13.43880 14.93790 12.55190 13.66830 16 112
#> 41 1416.28 13.54310 14.48840 12.34650 13.63650 16 112
#> 42 1412.63 13.81930 13.15790 11.73050 13.50180 16 112
#> 43 1407.84 14.02130 11.77100 11.06010 13.29600 16 112
#> 44 1396.03 14.26380 9.51721 9.94067 12.90450 16 112
#> 45 1395.16 14.29720 9.34148 9.85834 12.88540 16 112
#> 46 1393.35 14.74660 9.26154 9.97137 13.20540 16 112
#> 47 1392.09 15.37510 9.54412 10.34430 13.74600 16 112
#> 48 1391.28 15.54000 9.56985 10.42690 13.86320 16 112
#> 49 1389.36 15.99990 9.88680 10.78810 14.24390 16 112
#> 50 1383.23 16.82190 10.23320 11.36040 14.80120 16 112
#> 51 1369.52 18.30030 10.88910 12.42670 15.78270 16 112
#> 52 1369.33 17.92850 10.52810 12.09090 15.44410 16 112
#> 53 1368.42 17.88140 10.44130 12.04870 15.36520 16 112
#> 54 1366.48 17.61510 10.38310 11.94490 15.13790 16 112
#> 55 1363.06 16.78580 10.04080 11.47630 14.45060 16 112
#> 56 1362.31 13.82850 8.69583 9.76333 11.93010 16 112
#> 57 1357.67 15.93530 9.66199 10.99520 13.73650 16 112
#> 58 1355.79 15.74690 9.56148 10.89120 13.57230 16 112
#> 59 1355.65 15.58150 9.41209 10.79040 13.40820 16 112
#> 60 1355.00 15.68560 9.44333 10.83690 13.49060 16 112
#> 61 1354.90 15.68560 9.41837 10.83000 13.49110 16 112
#> 62 1354.80 15.68740 9.37022 10.82650 13.48340 16 112
#> 63 1354.61 15.70510 9.29862 10.83830 13.48970 16 112
#> 64 1353.66 15.72660 8.80199 10.87760 13.46370 16 112
#> 65 1352.13 15.77960 8.26341 10.97590 13.48090 16 112
#> 66 1349.69 15.87630 7.49375 11.15190 13.53640 16 112
#> 67 1348.18 16.29090 6.02434 11.68150 13.86540 16 112
#> 68 1346.34 16.71400 6.71746 11.93300 14.29090 16 112
#> 69 1345.37 16.43280 6.82777 11.72850 14.04620 16 112
#> 70 1345.16 16.42540 6.64103 11.76870 14.03860 16 112
#> 71 1345.13 16.45730 6.57512 11.81440 14.06960 16 112
#> 72 1345.11 16.48070 6.56530 11.85440 14.09680 16 112
#> 73 1345.06 16.47850 6.55807 11.89300 14.10530 16 112
#> 74 1344.80 16.45650 6.57700 12.10850 14.14470 16 112
#> 75 1344.41 16.39200 6.57185 12.45800 14.18110 16 112
#> 76 1343.45 16.19350 6.59046 13.30500 14.23230 16 112
#> 77 1341.39 15.78050 6.60390 15.30820 14.37450 16 112
#> 78 1337.19 15.05110 6.65227 19.61790 14.76860 16 112
#> 79 1327.88 13.93180 6.70938 28.02470 15.71720 16 112
#> 80 1302.82 12.80220 6.76827 38.86270 17.14170 16 112
#> 81 1284.74 12.33980 6.79452 43.29190 17.72490 16 112
#> 82 1282.34 12.76410 6.77007 41.87900 17.70010 16 112
#> 83 1279.32 13.15660 6.72522 39.98650 17.56910 16 112
#> 84 1277.51 13.16020 6.62288 41.20580 17.80480 16 112
#> 85 1276.86 13.41570 6.58475 40.84230 17.89380 16 112
#> 86 1276.27 13.89050 6.52270 40.70480 18.16860 16 112
#> 87 1276.13 14.08000 6.52653 40.79190 18.30560 16 112
#> 88 1276.05 14.20080 6.54900 40.96770 18.41330 16 112
#> 89 1276.01 14.20620 6.56174 41.01130 18.41290 16 112
#> 90 1275.87 14.31720 6.60695 41.11430 18.43160 16 112
#> 91 1275.68 14.34790 6.62197 41.15470 18.37020 16 112
#> 92 1273.84 14.55290 6.75451 41.40800 17.47950 16 112
#> 93 1271.14 14.86770 6.88971 41.58780 15.89640 16 112
#> 94 1269.45 15.41650 6.99641 41.70710 12.99820 16 112
#> 95 1267.17 15.67980 6.82707 41.49850 14.43090 16 112
#> 96 1266.02 15.77030 6.69588 41.32660 14.35680 16 112
#> 97 1265.63 15.43280 6.57327 41.14070 14.04780 16 112
#> 98 1265.59 15.43120 6.59821 41.17400 14.03770 16 112
#> 99 1265.59 15.40390 6.59947 41.18690 14.02430 16 112
#> 100 1265.59 15.39370 6.59984 41.19000 14.02580 16 112
#>
#> Trying to generate job results...
#>
#> Generating Overall.csv
#> Generating EtaEta.csv
#> Generating Eta.csv
#> Generating EtaStacked.csv
#> Generating bluptable.dat
#> Generating ConvergenceData.csv
#> Generating initest.csv
#> Generating doses.csv
#> Generating omega.csv
#> Generating omega_stderr.csv
#> Generating theta.csv
#> Generating thetaCorrelation.csv
#> Generating thetaCovariance.csv
#> Generating Covariance.csv
#> Generating Residuals.csv
#> Generating posthoc.csv
#>
#> Finished summarizing results. Transferring data and loading the results...
#> Done generating job results.
# View estimation results
print(job)
#> $Overall
#> Scenario RetCode LogLik -2LL AIC BIC nParm nObs nSub
#> <char> <int> <num> <num> <num> <num> <int> <int> <int>
#> 1: WorkFlow 1 -632.7953 1265.591 1283.591 1308.057 9 112 16
#> EpsShrinkage Condition
#> <num> <num>
#> 1: 0.17312 4.86919
#>
#> $theta
#> Scenario Parameter Estimate Units Stderr CV% 2.5%CI
#> <char> <char> <num> <lgcl> <num> <num> <num>
#> 1: WorkFlow tvV 15.393731 NA 1.13948176 7.402246 13.1338373
#> 2: WorkFlow tvCl 6.599838 NA 0.72057416 10.918058 5.1707492
#> 3: WorkFlow tvV2 41.190009 NA 1.06233995 2.579120 39.0831087
#> 4: WorkFlow tvCl2 14.025770 NA 1.03509504 7.379952 11.9729030
#> 5: WorkFlow CEps 0.161206 NA 0.01964795 12.188101 0.1222389
#> 97.5%CI Var.Inf.factor
#> <num> <lgcl>
#> 1: 17.6536240 NA
#> 2: 8.0289269 NA
#> 3: 43.2969102 NA
#> 4: 16.0786368 NA
#> 5: 0.2001731 NA
#>
#> $residuals
#> Scenario Subject IVAR TAD PRED IPRED DV IRES PREDSE
#> <char> <int> <num> <num> <num> <num> <num> <num> <lgcl>
#> 1: WorkFlow 1 0.00 0.00 1624.0400 1868.5100 2010 141.49100 NA
#> 2: WorkFlow 1 0.26 0.26 1159.4700 1260.1900 1330 69.80530 NA
#> 3: WorkFlow 1 1.10 1.10 479.8060 446.3790 565 118.62100 NA
#> 4: WorkFlow 1 2.10 2.10 270.6230 233.9580 216 -17.95780 NA
#> 5: WorkFlow 1 4.13 4.13 186.4900 159.1890 180 20.81080 NA
#> ---
#> 108: WorkFlow 16 0.85 0.85 600.5140 663.4330 535 -128.43300 NA
#> 109: WorkFlow 16 1.95 1.95 286.6770 366.3280 333 -33.32780 NA
#> 110: WorkFlow 16 3.91 3.91 191.0610 300.9670 341 40.03280 NA
#> 111: WorkFlow 16 7.98 7.98 129.5720 239.5410 273 33.45880 NA
#> 112: WorkFlow 16 24.19 24.19 29.2072 97.1843 101 3.81572 NA
#> Weight IWRES WRES CWRES PPRED PCWRES CdfPCWRES NPDE
#> <num> <num> <num> <num> <lgcl> <lgcl> <lgcl> <lgcl>
#> 1: 2.86424e-07 0.469733 0.770502 0.700128 NA NA NA NA
#> 2: 6.29687e-07 0.343613 0.164720 0.185784 NA NA NA NA
#> 3: 5.01872e-06 1.648460 0.727621 0.669790 NA NA NA NA
#> 4: 1.82694e-05 -0.476139 -1.233770 -1.425860 NA NA NA NA
#> 5: 3.94614e-05 0.810950 0.224090 0.349749 NA NA NA NA
#> ---
#> 108: 2.27199e-06 -1.200880 -1.007460 -1.059690 NA NA NA NA
#> 109: 7.45179e-06 -0.564359 1.094610 0.712562 NA NA NA NA
#> 110: 1.10398e-05 0.825116 2.424550 1.902710 NA NA NA NA
#> 111: 1.74277e-05 0.866460 1.738750 1.239370 NA NA NA NA
#> 112: 1.05879e-04 0.243556 2.906660 0.195257 NA NA NA NA
#> CdfDV NPD WhichDose ObsName WhichReset
#> <lgcl> <lgcl> <int> <char> <int>
#> 1: NA NA 1 CObs 0
#> 2: NA NA 1 CObs 0
#> 3: NA NA 1 CObs 0
#> 4: NA NA 1 CObs 0
#> 5: NA NA 1 CObs 0
#> ---
#> 108: NA NA 1 CObs 0
#> 109: NA NA 1 CObs 0
#> 110: NA NA 1 CObs 0
#> 111: NA NA 1 CObs 0
#> 112: NA NA 1 CObs 0
#>
#> $thetaCorrelation
#> Scenario tvV tvCl tvV2 tvCl2 CEps
#> <char> <num> <num> <num> <num> <int>
#> 1: WorkFlow 1.0000000 NA NA NA NA
#> 2: WorkFlow 0.8593004 1.0000000 NA NA NA
#> 3: WorkFlow 0.2894919 0.2079076 1.0000000 NA NA
#> 4: WorkFlow -0.1389020 -0.3488953 0.1590860 1.00000000 NA
#> 5: WorkFlow -0.2822157 -0.2928308 0.1160611 -0.07855954 1
#>
#> $thetaCovariance
#> Scenario tvV tvCl tvV2 tvCl2 CEps
#> <char> <num> <num> <num> <num> <num>
#> 1: WorkFlow 1.298418676 NA NA NA NA
#> 2: WorkFlow 0.705555361 0.519227114 NA NA NA
#> 3: WorkFlow 0.350434817 0.159152149 1.128566179 NA NA
#> 4: WorkFlow -0.163831025 -0.260228006 0.174934558 1.071421748 NA
#> 5: WorkFlow -0.006318379 -0.004145842 0.002422521 -0.001597704 0.0003860419
#>
#> $Covariance
#> Scenario Var Name tvV tvCl tvV2 tvCl2
#> <char> <char> <num> <num> <num> <num>
#> 1: WorkFlow tvV 1.298419e+00 NA NA NA
#> 2: WorkFlow tvCl 7.055554e-01 5.192271e-01 NA NA
#> 3: WorkFlow tvV2 3.504348e-01 1.591521e-01 1.128566e+00 NA
#> 4: WorkFlow tvCl2 -1.638310e-01 -2.602280e-01 1.749346e-01 1.071422e+00
#> 5: WorkFlow stdev0 -6.318379e-03 -4.145841e-03 2.422520e-03 -1.597704e-03
#> 6: WorkFlow nV-nV 3.295851e-03 -3.974253e-03 3.460331e-03 -7.978365e-04
#> 7: WorkFlow nCl-nCl -2.824936e-02 -2.128393e-02 1.239261e-02 3.031259e-03
#> 8: WorkFlow nV2-nV2 -2.012206e-06 -1.664998e-06 -1.420576e-06 1.872002e-06
#> 9: WorkFlow nCl2-nCl2 1.078400e-02 5.464257e-03 -3.682945e-03 -1.141343e-02
#> stdev0 nV-nV nCl-nCl nV2-nV2 nCl2-nCl2
#> <num> <num> <num> <num> <num>
#> 1: NA NA NA NA NA
#> 2: NA NA NA NA NA
#> 3: NA NA NA NA NA
#> 4: NA NA NA NA NA
#> 5: 3.860419e-04 NA NA NA NA
#> 6: 1.938038e-04 5.942643e-04 NA NA NA
#> 7: 5.552742e-04 8.076182e-04 2.810335e-03 NA NA
#> 8: -2.092889e-08 1.640274e-08 8.133134e-08 1.812036e-11 NA
#> 9: -2.736106e-04 -2.149110e-04 -6.455762e-04 8.716363e-09 0.001201008
#>
#> $ConvergenceData
#> Scenario Iter Parameter Value
#> <char> <int> <char> <num>
#> 1: WorkFlow 1 -2LL 1734.37818
#> 2: WorkFlow 1 tvV 1.07796
#> 3: WorkFlow 1 tvCl 1.05372
#> 4: WorkFlow 1 tvV2 1.10597
#> 5: WorkFlow 1 tvCl2 1.07380
#> ---
#> 496: WorkFlow 100 -2LL 1265.59061
#> 497: WorkFlow 100 tvV 15.39370
#> 498: WorkFlow 100 tvCl 6.59984
#> 499: WorkFlow 100 tvV2 41.19000
#> 500: WorkFlow 100 tvCl2 14.02580
#>
#> $posthocStacked
#> Scenario Subject time TableSource Params Value
#> <char> <int> <num> <char> <char> <num>
#> 1: WorkFlow 1 0.00 Main V 13.379650
#> 2: WorkFlow 1 0.00 Main Cl 7.672770
#> 3: WorkFlow 1 0.00 Main V2 41.189985
#> 4: WorkFlow 1 0.00 Main Cl2 13.208453
#> 5: WorkFlow 1 0.26 Main V 13.379650
#> ---
#> 444: WorkFlow 16 7.98 Main Cl2 16.627960
#> 445: WorkFlow 16 24.19 Main V 11.524632
#> 446: WorkFlow 16 24.19 Main Cl 3.300209
#> 447: WorkFlow 16 24.19 Main V2 41.190005
#> 448: WorkFlow 16 24.19 Main Cl2 16.627960
#>
#> $Eta
#> Scenario Subject nV nCl nV2 nCl2
#> <char> <int> <num> <num> <num> <num>
#> 1: WorkFlow 1 -0.140225458 0.15063262 -5.910225e-07 -0.06003932
#> 2: WorkFlow 2 0.479233952 0.50712172 2.710213e-07 -0.13042943
#> 3: WorkFlow 3 0.143254391 0.09407487 1.126426e-07 0.09854059
#> 4: WorkFlow 4 0.210215296 0.40147284 1.100524e-07 -0.04484961
#> 5: WorkFlow 5 -0.007981039 0.15095118 1.645053e-08 0.07616206
#> 6: WorkFlow 6 0.019609501 0.56872533 4.714352e-07 -0.16320367
#> 7: WorkFlow 7 0.109707888 0.13319750 3.831503e-07 0.14145686
#> 8: WorkFlow 8 -0.255086214 -0.36509625 -3.621856e-08 -0.09254025
#> 9: WorkFlow 9 0.203815159 0.35342619 -2.997510e-07 0.13252957
#> 10: WorkFlow 10 -0.307170037 -0.62454337 -4.546653e-07 0.12483107
#> 11: WorkFlow 11 0.057787931 0.17780777 -7.101335e-08 -0.37751129
#> 12: WorkFlow 12 -0.438271913 -0.82453717 1.020711e-08 -0.03548031
#> 13: WorkFlow 13 0.066650844 0.29173113 -1.066623e-07 -0.10274682
#> 14: WorkFlow 14 0.189123315 -0.01438328 2.045132e-08 0.18575368
#> 15: WorkFlow 15 -0.103658405 -0.25226211 1.289549e-07 0.02649993
#> 16: WorkFlow 16 -0.289473668 -0.69305946 -1.137622e-07 0.17018928
#> nV_shrinkage nCl_shrinkage nV2_shrinkage nCl2_shrinkage
#> <num> <num> <num> <num>
#> 1: 0.199262 0.017544 0.999997 0.496458
#> 2: 0.188192 0.015425 0.999997 0.587153
#> 3: 0.220227 0.019892 0.999997 0.567375
#> 4: 0.223739 0.015565 0.999997 0.522552
#> 5: 0.220407 0.017856 0.999997 0.519462
#> 6: 0.215600 0.013255 0.999996 0.396486
#> 7: 0.254336 0.018879 0.999997 0.525748
#> 8: 0.230788 0.030799 0.999998 0.456568
#> 9: 0.229899 0.015715 0.999997 0.541134
#> 10: 0.230452 0.044376 0.999998 0.494883
#> 11: 0.202502 0.018919 0.999998 0.481259
#> 12: 0.258034 0.055488 0.999998 0.383910
#> 13: 0.238111 0.016567 0.999997 0.464855
#> 14: 0.240825 0.022845 0.999998 0.572972
#> 15: 0.251156 0.027781 0.999998 0.481207
#> 16: 0.253619 0.049125 0.999998 0.454066
#>
#> $EtaStacked
#> Scenario Subject Eta Value
#> <char> <int> <char> <num>
#> 1: WorkFlow 1 nV -1.402255e-01
#> 2: WorkFlow 1 nCl 1.506326e-01
#> 3: WorkFlow 1 nV2 -5.910225e-07
#> 4: WorkFlow 1 nCl2 -6.003932e-02
#> 5: WorkFlow 2 nV 4.792340e-01
#> 6: WorkFlow 2 nCl 5.071217e-01
#> 7: WorkFlow 2 nV2 2.710213e-07
#> 8: WorkFlow 2 nCl2 -1.304294e-01
#> 9: WorkFlow 3 nV 1.432544e-01
#> 10: WorkFlow 3 nCl 9.407487e-02
#> 11: WorkFlow 3 nV2 1.126426e-07
#> 12: WorkFlow 3 nCl2 9.854059e-02
#> 13: WorkFlow 4 nV 2.102153e-01
#> 14: WorkFlow 4 nCl 4.014728e-01
#> 15: WorkFlow 4 nV2 1.100524e-07
#> 16: WorkFlow 4 nCl2 -4.484961e-02
#> 17: WorkFlow 5 nV -7.981039e-03
#> 18: WorkFlow 5 nCl 1.509512e-01
#> 19: WorkFlow 5 nV2 1.645053e-08
#> 20: WorkFlow 5 nCl2 7.616206e-02
#> 21: WorkFlow 6 nV 1.960950e-02
#> 22: WorkFlow 6 nCl 5.687253e-01
#> 23: WorkFlow 6 nV2 4.714352e-07
#> 24: WorkFlow 6 nCl2 -1.632037e-01
#> 25: WorkFlow 7 nV 1.097079e-01
#> 26: WorkFlow 7 nCl 1.331975e-01
#> 27: WorkFlow 7 nV2 3.831503e-07
#> 28: WorkFlow 7 nCl2 1.414569e-01
#> 29: WorkFlow 8 nV -2.550862e-01
#> 30: WorkFlow 8 nCl -3.650963e-01
#> 31: WorkFlow 8 nV2 -3.621856e-08
#> 32: WorkFlow 8 nCl2 -9.254025e-02
#> 33: WorkFlow 9 nV 2.038152e-01
#> 34: WorkFlow 9 nCl 3.534262e-01
#> 35: WorkFlow 9 nV2 -2.997510e-07
#> 36: WorkFlow 9 nCl2 1.325296e-01
#> 37: WorkFlow 10 nV -3.071700e-01
#> 38: WorkFlow 10 nCl -6.245434e-01
#> 39: WorkFlow 10 nV2 -4.546653e-07
#> 40: WorkFlow 10 nCl2 1.248311e-01
#> 41: WorkFlow 11 nV 5.778793e-02
#> 42: WorkFlow 11 nCl 1.778078e-01
#> 43: WorkFlow 11 nV2 -7.101335e-08
#> 44: WorkFlow 11 nCl2 -3.775113e-01
#> 45: WorkFlow 12 nV -4.382719e-01
#> 46: WorkFlow 12 nCl -8.245372e-01
#> 47: WorkFlow 12 nV2 1.020711e-08
#> 48: WorkFlow 12 nCl2 -3.548031e-02
#> 49: WorkFlow 13 nV 6.665084e-02
#> 50: WorkFlow 13 nCl 2.917311e-01
#> 51: WorkFlow 13 nV2 -1.066623e-07
#> 52: WorkFlow 13 nCl2 -1.027468e-01
#> 53: WorkFlow 14 nV 1.891233e-01
#> 54: WorkFlow 14 nCl -1.438328e-02
#> 55: WorkFlow 14 nV2 2.045132e-08
#> 56: WorkFlow 14 nCl2 1.857537e-01
#> 57: WorkFlow 15 nV -1.036584e-01
#> 58: WorkFlow 15 nCl -2.522621e-01
#> 59: WorkFlow 15 nV2 1.289549e-07
#> 60: WorkFlow 15 nCl2 2.649993e-02
#> 61: WorkFlow 16 nV -2.894737e-01
#> 62: WorkFlow 16 nCl -6.930595e-01
#> 63: WorkFlow 16 nV2 -1.137622e-07
#> 64: WorkFlow 16 nCl2 1.701893e-01
#> Scenario Subject Eta Value
#>
#> $EtaEta
#> Scenario Subject Eta1Name Eta2Name Eta1 Eta2
#> <char> <int> <char> <char> <num> <num>
#> 1: WorkFlow 1 nV nCl -1.402255e-01 1.506326e-01
#> 2: WorkFlow 1 nV nV2 -1.402255e-01 -5.910225e-07
#> 3: WorkFlow 1 nV nCl2 -1.402255e-01 -6.003932e-02
#> 4: WorkFlow 1 nCl nV 1.506326e-01 -1.402255e-01
#> 5: WorkFlow 1 nCl nV2 1.506326e-01 -5.910225e-07
#> ---
#> 188: WorkFlow 16 nV2 nCl -1.137622e-07 -6.930595e-01
#> 189: WorkFlow 16 nV2 nCl2 -1.137622e-07 1.701893e-01
#> 190: WorkFlow 16 nCl2 nV 1.701893e-01 -2.894737e-01
#> 191: WorkFlow 16 nCl2 nCl 1.701893e-01 -6.930595e-01
#> 192: WorkFlow 16 nCl2 nV2 1.701893e-01 -1.137622e-07
#>
#> $omega
#> Label nV nCl nV2 nCl2
#> <char> <num> <num> <num> <num>
#> 1: nV 0.06914959 0.0000000 0.000000e+00 0.00000000
#> 2: nCl 0.00000000 0.1806635 0.000000e+00 0.00000000
#> 3: nV2 0.00000000 0.0000000 4.649898e-08 0.00000000
#> 4: nCl2 0.00000000 0.0000000 0.000000e+00 0.04239193
#>
#> $omega_Correlation
#> Label nV nCl nV2 nCl2
#> <char> <num> <num> <num> <num>
#> 1: nV 1 0 0 0
#> 2: nCl 0 1 0 0
#> 3: nV2 0 0 1 0
#> 4: nCl2 0 0 0 1
#>
#> $Eta_Shrinkage
#> Label nV nCl nV2 nCl2
#> <char> <num> <num> <num> <num>
#> 1: Shrinkage 0.1196793 0.01532516 0.9987514 0.2869761
#>
#> $omega_stderr
#> Label nV nCl nV2 nCl2
#> <char> <num> <num> <num> <num>
#> 1: nV 0.02437754 0.0000000 0.000000e+00 0.00000000
#> 2: nCl 0.00000000 0.0530126 0.000000e+00 0.00000000
#> 3: nV2 0.00000000 0.0000000 4.256802e-06 0.00000000
#> 4: nCl2 0.00000000 0.0000000 0.000000e+00 0.03465556
#>
#> $bluptable
#> Scenario Subject Eta ETA_VAL ETA_SE OMEGA_SD ETA_SE/OMEGA_SD
#> <char> <int> <char> <num> <num> <num> <num>
#> 1: WorkFlow 1 nV -0.140225 0.117383 0.262963 0.446387
#> 2: WorkFlow 2 nV 0.479234 0.114076 0.262963 0.433811
#> 3: WorkFlow 3 nV 0.143254 0.123404 0.262963 0.469284
#> 4: WorkFlow 4 nV 0.210215 0.124384 0.262963 0.473010
#> 5: WorkFlow 5 nV -0.007981 0.123455 0.262963 0.469475
#> 6: WorkFlow 6 nV 0.019610 0.122101 0.262963 0.464327
#> 7: WorkFlow 7 nV 0.109708 0.132617 0.262963 0.504317
#> 8: WorkFlow 8 nV -0.255086 0.126329 0.262963 0.480404
#> 9: WorkFlow 9 nV 0.203815 0.126085 0.262963 0.479478
#> 10: WorkFlow 10 nV -0.307170 0.126236 0.262963 0.480054
#> 11: WorkFlow 11 nV 0.057788 0.118334 0.262963 0.450002
#> 12: WorkFlow 12 nV -0.438272 0.133577 0.262963 0.507970
#> 13: WorkFlow 13 nV 0.066651 0.128317 0.262963 0.487966
#> 14: WorkFlow 14 nV 0.189123 0.129046 0.262963 0.490739
#> 15: WorkFlow 15 nV -0.103658 0.131785 0.262963 0.501154
#> 16: WorkFlow 16 nV -0.289474 0.132430 0.262963 0.503606
#> 17: WorkFlow 1 nCl 0.150633 0.056299 0.425045 0.132454
#> 18: WorkFlow 2 nCl 0.507122 0.052790 0.425045 0.124198
#> 19: WorkFlow 3 nCl 0.094075 0.059948 0.425045 0.141039
#> 20: WorkFlow 4 nCl 0.401473 0.053028 0.425045 0.124758
#> 21: WorkFlow 5 nCl 0.150951 0.056798 0.425045 0.133628
#> 22: WorkFlow 6 nCl 0.568725 0.048936 0.425045 0.115130
#> 23: WorkFlow 7 nCl 0.133198 0.058401 0.425045 0.137400
#> 24: WorkFlow 8 nCl -0.365096 0.074594 0.425045 0.175497
#> 25: WorkFlow 9 nCl 0.353426 0.053283 0.425045 0.125358
#> 26: WorkFlow 10 nCl -0.624543 0.089538 0.425045 0.210655
#> 27: WorkFlow 11 nCl 0.177808 0.058463 0.425045 0.137545
#> 28: WorkFlow 12 nCl -0.824537 0.100123 0.425045 0.235558
#> 29: WorkFlow 13 nCl 0.291731 0.054709 0.425045 0.128712
#> 30: WorkFlow 14 nCl -0.014383 0.064244 0.425045 0.151146
#> 31: WorkFlow 15 nCl -0.252262 0.070845 0.425045 0.166675
#> 32: WorkFlow 16 nCl -0.693059 0.094208 0.425045 0.221643
#> 33: WorkFlow 1 nV2 -0.000001 0.000216 0.000216 0.999998
#> 34: WorkFlow 2 nV2 0.000000 0.000216 0.000216 0.999999
#> 35: WorkFlow 3 nV2 0.000000 0.000216 0.000216 0.999999
#> 36: WorkFlow 4 nV2 0.000000 0.000216 0.000216 0.999998
#> 37: WorkFlow 5 nV2 0.000000 0.000216 0.000216 0.999999
#> 38: WorkFlow 6 nV2 0.000000 0.000216 0.000216 0.999998
#> 39: WorkFlow 7 nV2 0.000000 0.000216 0.000216 0.999999
#> 40: WorkFlow 8 nV2 0.000000 0.000216 0.000216 0.999999
#> 41: WorkFlow 9 nV2 0.000000 0.000216 0.000216 0.999998
#> 42: WorkFlow 10 nV2 0.000000 0.000216 0.000216 0.999999
#> 43: WorkFlow 11 nV2 0.000000 0.000216 0.000216 0.999999
#> 44: WorkFlow 12 nV2 0.000000 0.000216 0.000216 0.999999
#> 45: WorkFlow 13 nV2 0.000000 0.000216 0.000216 0.999998
#> 46: WorkFlow 14 nV2 0.000000 0.000216 0.000216 0.999999
#> 47: WorkFlow 15 nV2 0.000000 0.000216 0.000216 0.999999
#> 48: WorkFlow 16 nV2 0.000000 0.000216 0.000216 0.999999
#> 49: WorkFlow 1 nCl2 -0.060039 0.145072 0.205893 0.704598
#> 50: WorkFlow 2 nCl2 -0.130429 0.157767 0.205893 0.766259
#> 51: WorkFlow 3 nCl2 0.098541 0.155087 0.205893 0.753243
#> 52: WorkFlow 4 nCl2 -0.044850 0.148835 0.205893 0.722877
#> 53: WorkFlow 5 nCl2 0.076162 0.148395 0.205893 0.720737
#> 54: WorkFlow 6 nCl2 -0.163204 0.129645 0.205893 0.629671
#> 55: WorkFlow 7 nCl2 0.141457 0.149290 0.205893 0.725085
#> 56: WorkFlow 8 nCl2 -0.092540 0.139122 0.205893 0.675698
#> 57: WorkFlow 9 nCl2 0.132530 0.151459 0.205893 0.735618
#> 58: WorkFlow 10 nCl2 0.124831 0.144841 0.205893 0.703479
#> 59: WorkFlow 11 nCl2 -0.377511 0.142834 0.205893 0.693728
#> 60: WorkFlow 12 nCl2 -0.035480 0.127572 0.205893 0.619605
#> 61: WorkFlow 13 nCl2 -0.102747 0.140378 0.205893 0.681803
#> 62: WorkFlow 14 nCl2 0.185754 0.155850 0.205893 0.756949
#> 63: WorkFlow 15 nCl2 0.026500 0.142826 0.205893 0.693691
#> 64: WorkFlow 16 nCl2 0.170189 0.138740 0.205893 0.673844
#> Scenario Subject Eta ETA_VAL ETA_SE OMEGA_SD ETA_SE/OMEGA_SD
#> SHRINK_SUB_VAR SHRINK_SUB_SD
#> <num> <num>
#> 1: 0.199262 0.105160
#> 2: 0.188192 0.098996
#> 3: 0.220227 0.116953
#> 4: 0.223739 0.118943
#> 5: 0.220407 0.117054
#> 6: 0.215600 0.114336
#> 7: 0.254336 0.136482
#> 8: 0.230788 0.122953
#> 9: 0.229899 0.122446
#> 10: 0.230452 0.122761
#> 11: 0.202502 0.106973
#> 12: 0.258034 0.138625
#> 13: 0.238111 0.127137
#> 14: 0.240825 0.128693
#> 15: 0.251156 0.134642
#> 16: 0.253619 0.136067
#> 17: 0.017544 0.008811
#> 18: 0.015425 0.007743
#> 19: 0.019892 0.009996
#> 20: 0.015565 0.007813
#> 21: 0.017856 0.008968
#> 22: 0.013255 0.006650
#> 23: 0.018879 0.009484
#> 24: 0.030799 0.015520
#> 25: 0.015715 0.007888
#> 26: 0.044376 0.022440
#> 27: 0.018919 0.009505
#> 28: 0.055488 0.028140
#> 29: 0.016567 0.008318
#> 30: 0.022845 0.011489
#> 31: 0.027781 0.013988
#> 32: 0.049125 0.024872
#> 33: 0.999997 0.998213
#> 34: 0.999997 0.998411
#> 35: 0.999997 0.998396
#> 36: 0.999997 0.998232
#> 37: 0.999997 0.998271
#> 38: 0.999996 0.998074
#> 39: 0.999997 0.998296
#> 40: 0.999998 0.998536
#> 41: 0.999997 0.998250
#> 42: 0.999998 0.998530
#> 43: 0.999998 0.998485
#> 44: 0.999998 0.998644
#> 45: 0.999997 0.998238
#> 46: 0.999998 0.998491
#> 47: 0.999998 0.998465
#> 48: 0.999998 0.998566
#> 49: 0.496458 0.290393
#> 50: 0.587153 0.357468
#> 51: 0.567375 0.342258
#> 52: 0.522552 0.309024
#> 53: 0.519462 0.306792
#> 54: 0.396486 0.223138
#> 55: 0.525748 0.311340
#> 56: 0.456568 0.262822
#> 57: 0.541134 0.322603
#> 58: 0.494883 0.289284
#> 59: 0.481259 0.279763
#> 60: 0.383910 0.215086
#> 61: 0.464855 0.268464
#> 62: 0.572972 0.346526
#> 63: 0.481207 0.279727
#> 64: 0.454066 0.261126
#> SHRINK_SUB_VAR SHRINK_SUB_SD
#>
#> $posthoc
#> Scenario Subject time V Cl V2 Cl2 TableSource
#> <char> <int> <num> <num> <num> <num> <num> <char>
#> 1: WorkFlow 1 0.00 13.37965 7.672770 41.18999 13.20845 Main
#> 2: WorkFlow 1 0.26 13.37965 7.672770 41.18999 13.20845 Main
#> 3: WorkFlow 1 1.10 13.37965 7.672770 41.18999 13.20845 Main
#> 4: WorkFlow 1 2.10 13.37965 7.672770 41.18999 13.20845 Main
#> 5: WorkFlow 1 4.13 13.37965 7.672770 41.18999 13.20845 Main
#> ---
#> 108: WorkFlow 16 0.85 11.52463 3.300209 41.19000 16.62796 Main
#> 109: WorkFlow 16 1.95 11.52463 3.300209 41.19000 16.62796 Main
#> 110: WorkFlow 16 3.91 11.52463 3.300209 41.19000 16.62796 Main
#> 111: WorkFlow 16 7.98 11.52463 3.300209 41.19000 16.62796 Main
#> 112: WorkFlow 16 24.19 11.52463 3.300209 41.19000 16.62796 Main
#>
#> $nlme7engine.log
#> [1] " ELS FOCE/Laplace engine log file "
#> [2] " Closed-form model can give incorrect results"
#> [3] " if time-varying model parameters are present."
#> [4] ""
#> [5] " MAX ITERS = 1000"
#> [6] ""
#> [7] ""
#> [8] " Model evaluation ODE level used:"
#> [9] " 0=none"
#> [10] ""
#> [11] " Other flags and run conditions/tolerances:"
#> [12] ""
#> [13] " 16 NSUB"
#> [14] " 9 NUMFREEPARAM"
#> [15] " 4 NUMFREEFIXEF"
#> [16] " 1 NUMFREEEPS "
#> [17] " 4 NUMFREEOMEGAPARAM"
#> [18] " 4 NUMRANEF"
#> [19] " 1 Allow Gaussian Fit"
#> [20] " 5 IDENGINE"
#> [21] " 0 IFLAGRESTART"
#> [22] " 1 IFLAGSTDERR "
#> [23] " 0 IFLAGNP "
#> [24] " 1 IFLAGFOCEHESS"
#> [25] " 1 NORDERAGQ"
#> [26] " 1 NUMPROCS "
#> [27] " 13 NDIGITblup "
#> [28] " 7 NDIGITlagl "
#> [29] " 0.200E-02 tolmodlinz "
#> [30] " 0.100E-01 tolstderr "
#> [31] " 0.200E-03 GradTolOuter "
#> [32] " 0.100E-03 StepTolOuter "
#> [33] " 0.171E-04 GradTolInner "
#> [34] " 0.707E-07 StepTolInner "
#> [35] " 0.100E-02 RefDeltaLAGL"
#> [36] " 0 NREP_PCWRES"
#> [37] " 1 NPRESAMPLE "
#> [38] " 0 NGETMAPNP "
#> [39] " 0 iflaganagrad"
#> [40] " 0 iodelevelused"
#> [41] " 1 iflagstderr "
#> [42] " 1 iflagwhichstderr"
#> [43] ""
#> [44] ""
#> [45] " Initial parameter values:"
#> [46] ""
#> [47] " 0.10000000E+01 1 1 THETA"
#> [48] " 0.10000000E+01 2 1 THETA"
#> [49] " 0.10000000E+01 3 1 THETA"
#> [50] " 0.10000000E+01 4 1 THETA"
#> [51] " 0.10000000E+00 1 1 EPS_STD_DEV"
#> [52] " 0.10000000E-01 1 1 EPS_VARIANCE"
#> [53] " 0.10000000E+01 1 1 OMEGA"
#> [54] " 0.00000000E+00 2 1 OMEGA"
#> [55] " 0.10000000E+01 2 2 OMEGA"
#> [56] " 0.00000000E+00 3 1 OMEGA"
#> [57] " 0.00000000E+00 3 2 OMEGA"
#> [58] " 0.10000000E+01 3 3 OMEGA"
#> [59] " 0.00000000E+00 4 1 OMEGA"
#> [60] " 0.00000000E+00 4 2 OMEGA"
#> [61] " 0.00000000E+00 4 3 OMEGA"
#> [62] " 0.10000000E+01 4 4 OMEGA"
#> [63] " 0.96154314E+03 1 1 -LOGLIKE"
#> [64] " 0.17172440E+04 1 1 ELSOBJ "
#> [65] ""
#> [66] " Total # Subjects = 16"
#> [67] " Total # Observations = 112"
#> [68] ""
#> [69] ""
#> [70] " FOCE Hessian approximation used"
#> [71] "iteration = 1 LL value = 867.18909 "
#> [72] "iteration = 2 LL value = 820.37827 "
#> [73] "iteration = 3 LL value = 793.54515 "
#> [74] "iteration = 4 LL value = 779.96933 "
#> [75] "iteration = 5 LL value = 772.76678 "
#> [76] "iteration = 6 LL value = 768.72801 "
#> [77] "iteration = 7 LL value = 766.40380 "
#> [78] "iteration = 8 LL value = 765.21363 "
#> [79] "iteration = 9 LL value = 764.72055 "
#> [80] "iteration = 10 LL value = 764.21784 "
#> [81] "iteration = 11 LL value = 762.92153 "
#> [82] "iteration = 12 LL value = 759.83450 "
#> [83] "iteration = 13 LL value = 757.54976 "
#> [84] "iteration = 14 LL value = 755.01567 "
#> [85] "iteration = 15 LL value = 752.01467 "
#> [86] "iteration = 16 LL value = 751.64331 "
#> [87] "iteration = 17 LL value = 747.99571 "
#> [88] "iteration = 18 LL value = 746.87217 "
#> [89] "iteration = 19 LL value = 743.42994 "
#> [90] "iteration = 20 LL value = 741.16621 "
#> [91] "iteration = 21 LL value = 732.61665 "
#> [92] "iteration = 22 LL value = 724.63673 "
#> [93] "iteration = 23 LL value = 724.09957 "
#> [94] "iteration = 24 LL value = 723.59835 "
#> [95] "iteration = 25 LL value = 722.89780 "
#> [96] "iteration = 26 LL value = 721.83993 "
#> [97] "iteration = 27 LL value = 721.14183 "
#> [98] "iteration = 28 LL value = 719.95062 "
#> [99] "iteration = 29 LL value = 719.02061 "
#> [100] "iteration = 30 LL value = 716.77849 "
#> [101] "iteration = 31 LL value = 715.40198 "
#> [102] "iteration = 32 LL value = 713.12037 "
#> [103] "iteration = 33 LL value = 711.90441 "
#> [104] "iteration = 34 LL value = 710.27424 "
#> [105] "iteration = 35 LL value = 709.75611 "
#> [106] "iteration = 36 LL value = 709.47099 "
#> [107] "iteration = 37 LL value = 709.26146 "
#> [108] "iteration = 38 LL value = 709.17819 "
#> [109] "iteration = 39 LL value = 709.03148 "
#> [110] "iteration = 40 LL value = 708.75429 "
#> [111] "iteration = 41 LL value = 708.14185 "
#> [112] "iteration = 42 LL value = 706.31437 "
#> [113] "iteration = 43 LL value = 703.92112 "
#> [114] "iteration = 44 LL value = 698.01360 "
#> [115] "iteration = 45 LL value = 697.57812 "
#> [116] "iteration = 46 LL value = 696.67465 "
#> [117] "iteration = 47 LL value = 696.04595 "
#> [118] "iteration = 48 LL value = 695.63985 "
#> [119] "iteration = 49 LL value = 694.67991 "
#> [120] "iteration = 50 LL value = 691.61401 "
#> [121] "iteration = 51 LL value = 684.75854 "
#> [122] "iteration = 52 LL value = 684.66293 "
#> [123] "iteration = 53 LL value = 684.21021 "
#> [124] "iteration = 54 LL value = 683.24017 "
#> [125] "iteration = 55 LL value = 681.53076 "
#> [126] "iteration = 56 LL value = 681.15502 "
#> [127] "iteration = 57 LL value = 678.83650 "
#> [128] "iteration = 58 LL value = 677.89262 "
#> [129] "iteration = 59 LL value = 677.82653 "
#> [130] "iteration = 60 LL value = 677.50112 "
#> [131] "iteration = 61 LL value = 677.45029 "
#> [132] "iteration = 62 LL value = 677.39873 "
#> [133] "iteration = 63 LL value = 677.30299 "
#> [134] "iteration = 64 LL value = 676.82981 "
#> [135] "iteration = 65 LL value = 676.06654 "
#> [136] "iteration = 66 LL value = 674.84407 "
#> [137] "iteration = 67 LL value = 674.08786 "
#> [138] "iteration = 68 LL value = 673.17017 "
#> [139] "iteration = 69 LL value = 672.68646 "
#> [140] "iteration = 70 LL value = 672.58149 "
#> [141] "iteration = 71 LL value = 672.56561 "
#> [142] "iteration = 72 LL value = 672.55297 "
#> [143] "iteration = 73 LL value = 672.52843 "
#> [144] "iteration = 74 LL value = 672.40163 "
#> [145] "iteration = 75 LL value = 672.20549 "
#> [146] "iteration = 76 LL value = 671.72710 "
#> [147] "iteration = 77 LL value = 670.69505 "
#> [148] "iteration = 78 LL value = 668.59442 "
#> [149] "iteration = 79 LL value = 663.94055 "
#> [150] "iteration = 80 LL value = 651.40834 "
#> [151] "iteration = 81 LL value = 642.36917 "
#> [152] "iteration = 82 LL value = 641.17190 "
#> [153] "iteration = 83 LL value = 639.66183 "
#> [154] "iteration = 84 LL value = 638.75342 "
#> [155] "iteration = 85 LL value = 638.42781 "
#> [156] "iteration = 86 LL value = 638.13694 "
#> [157] "iteration = 87 LL value = 638.06704 "
#> [158] "iteration = 88 LL value = 638.02520 "
#> [159] "iteration = 89 LL value = 638.00360 "
#> [160] "iteration = 90 LL value = 637.93683 "
#> [161] "iteration = 91 LL value = 637.84198 "
#> [162] "iteration = 92 LL value = 636.92085 "
#> [163] "iteration = 93 LL value = 635.57006 "
#> [164] "iteration = 94 LL value = 634.72493 "
#> [165] "iteration = 95 LL value = 633.58526 "
#> [166] "iteration = 96 LL value = 633.00975 "
#> [167] "iteration = 97 LL value = 632.81534 "
#> [168] "iteration = 98 LL value = 632.79614 "
#> [169] "iteration = 99 LL value = 632.79535 "
#> [170] "iteration = 100 LL value = 632.79530 "
#> [171] " Main umdriver optimization for engine 5 terminated"
#> [172] " ITRMCD exit code from umdriver = 1"
#> [173] " 1,2,3 = probable success, 4 = maxiterations reached"
#> [174] " 7 = Linesearch failed to find any good gradient to move"
#> [175] ""
#> [176] " Starting std error computation"
#> [177] ""
#> [178] " Hessian based std error successful"
#> [179] " Standard errors of estimated parameters"
#> [180] ""
#> [181] "internal_coords Param Param_val Std_err rel_std_err"
#> [182] ""
#> [183] " 1 tvV 0.153937E+02 0.114953E+01 0.746749E-01"
#> [184] " 2 tvCl 0.659984E+01 0.709606E+00 0.107519E+00"
#> [185] " 3 tvV2 0.411900E+02 0.138707E+01 0.336749E-01"
#> [186] " 4 tvCl2 0.140258E+02 0.105920E+01 0.755181E-01"
#> [187] " 5 CEps 0.161206E+00 0.146372E-01 0.907979E-01"
#> [188] " 6 nV 0.262963E+00 0.591546E-01 0.224954E+00"
#> [189] " 7 nCl 0.425045E+00 0.779119E-01 0.183303E+00"
#> [190] " 8 nV2 0.215636E-03 0.199261E-01 0.924063E+02"
#> [191] " 9 nCl2 0.205893E+00 0.761848E-01 0.370021E+00"
#> [192] ""
#> [193] ""
#> [194] "external_coords Param Param_val Std_err rel_std_err"
#> [195] ""
#> [196] " 1 tvV 0.153937E+02 0.114953E+01 0.746749E-01"
#> [197] " 2 tvCl 0.659984E+01 0.709606E+00 0.107519E+00"
#> [198] " 3 tvV2 0.411900E+02 0.138707E+01 0.336749E-01"
#> [199] " 4 tvCl2 0.140258E+02 0.105920E+01 0.755181E-01"
#> [200] " 5 CEps 0.161206E+00 0.146372E-01 0.907979E-01"
#> [201] " 6 nV 0.691496E-01 0.311110E-01 0.449908E+00"
#> [202] " 7 nCl 0.180664E+00 0.662322E-01 0.366605E+00"
#> [203] " 8 nV2 0.464990E-07 0.859359E-05 0.184813E+03"
#> [204] " 9 nCl2 0.423919E-01 0.313718E-01 0.740042E+00"
#> [205] ""
#> [206] ""
#> [207] " Score based std error successful"
#> [208] " Standard errors of estimated parameters"
#> [209] ""
#> [210] "internal_coords Param Param_val Std_err rel_std_err"
#> [211] ""
#> [212] " 1 tvV 0.153937E+02 0.710233E+01 0.461378E+00"
#> [213] " 2 tvCl 0.659984E+01 0.373966E+01 0.566629E+00"
#> [214] " 3 tvV2 0.411900E+02 0.261175E+01 0.634073E-01"
#> [215] " 4 tvCl2 0.140258E+02 0.288048E+01 0.205370E+00"
#> [216] " 5 CEps 0.161206E+00 0.252833E-01 0.156838E+00"
#> [217] " 6 nV 0.262963E+00 0.269659E+00 0.102546E+01"
#> [218] " 7 nCl 0.425045E+00 0.264192E+00 0.621562E+00"
#> [219] " 8 nV2 0.215636E-03 0.912650E-01 0.423236E+03"
#> [220] " 9 nCl2 0.205893E+00 0.172252E+00 0.836607E+00"
#> [221] ""
#> [222] ""
#> [223] "external_coords Param Param_val Std_err rel_std_err"
#> [224] ""
#> [225] " 1 tvV 0.153937E+02 0.710233E+01 0.461378E+00"
#> [226] " 2 tvCl 0.659984E+01 0.373966E+01 0.566629E+00"
#> [227] " 3 tvV2 0.411900E+02 0.261175E+01 0.634073E-01"
#> [228] " 4 tvCl2 0.140258E+02 0.288048E+01 0.205370E+00"
#> [229] " 5 CEps 0.161206E+00 0.252833E-01 0.156838E+00"
#> [230] " 6 nV 0.691496E-01 0.141821E+00 0.205093E+01"
#> [231] " 7 nCl 0.180664E+00 0.224587E+00 0.124312E+01"
#> [232] " 8 nV2 0.464990E-07 0.393601E-04 0.846472E+03"
#> [233] " 9 nCl2 0.423919E-01 0.709308E-01 0.167321E+01"
#> [234] ""
#> [235] ""
#> [236] " Sandwich based std error successful"
#> [237] " Standard errors of estimated parameters"
#> [238] ""
#> [239] "internal_coords Param Param_val Std_err rel_std_err"
#> [240] ""
#> [241] " 1 tvV 0.153937E+02 0.113948E+01 0.740225E-01"
#> [242] " 2 tvCl 0.659984E+01 0.720574E+00 0.109181E+00"
#> [243] " 3 tvV2 0.411900E+02 0.106234E+01 0.257912E-01"
#> [244] " 4 tvCl2 0.140258E+02 0.103510E+01 0.737995E-01"
#> [245] " 5 CEps 0.161206E+00 0.196479E-01 0.121881E+00"
#> [246] " 6 nV 0.262963E+00 0.463516E-01 0.176267E+00"
#> [247] " 7 nCl 0.425045E+00 0.623611E-01 0.146716E+00"
#> [248] " 8 nV2 0.215636E-03 0.987033E-02 0.457731E+02"
#> [249] " 9 nCl2 0.205893E+00 0.841591E-01 0.408752E+00"
#> [250] ""
#> [251] ""
#> [252] "external_coords Param Param_val Std_err rel_std_err"
#> [253] ""
#> [254] " 1 tvV 0.153937E+02 0.113948E+01 0.740225E-01"
#> [255] " 2 tvCl 0.659984E+01 0.720574E+00 0.109181E+00"
#> [256] " 3 tvV2 0.411900E+02 0.106234E+01 0.257912E-01"
#> [257] " 4 tvCl2 0.140258E+02 0.103510E+01 0.737995E-01"
#> [258] " 5 CEps 0.161206E+00 0.196479E-01 0.121881E+00"
#> [259] " 6 nV 0.691496E-01 0.243775E-01 0.352533E+00"
#> [260] " 7 nCl 0.180664E+00 0.530126E-01 0.293433E+00"
#> [261] " 8 nV2 0.464990E-07 0.425680E-05 0.915461E+02"
#> [262] " 9 nCl2 0.423919E-01 0.346556E-01 0.817504E+00"
#> [263] ""
#> [264] ""
#> [265] " Engine Convergence Summary:"
#> [266] " Convergence achieved"
#> [267] ""
#> [268] " ITRMCD value on last iteration = 1"
#> [269] " NLME7 return code= 1"
#> [270] ""
#> [271] " Optimal parameters:"
#> [272] ""
#> [273] " 0.15393731E+02 1 1 THETA BOUNDS: -0.100+101 0.100+101"
#> [274] " 0.65998381E+01 2 1 THETA BOUNDS: -0.100+101 0.100+101"
#> [275] " 0.41190009E+02 3 1 THETA BOUNDS: -0.100+101 0.100+101"
#> [276] " 0.14025770E+02 4 1 THETA BOUNDS: -0.100+101 0.100+101"
#> [277] " 0.16120598E+00 1 1 EPS_STD_DEV BOUNDS: 0.100E-02 0.100+101"
#> [278] " 0.25987369E-01 1 1 EPS_VARIANCE"
#> [279] " 0.69149588E-01 1 1 OMEGA"
#> [280] " 0.00000000E+00 2 1 OMEGA"
#> [281] " 0.18066353E+00 2 2 OMEGA"
#> [282] " 0.00000000E+00 3 1 OMEGA"
#> [283] " 0.00000000E+00 3 2 OMEGA"
#> [284] " 0.46498975E-07 3 3 OMEGA"
#> [285] " 0.00000000E+00 4 1 OMEGA"
#> [286] " 0.00000000E+00 4 2 OMEGA"
#> [287] " 0.00000000E+00 4 3 OMEGA"
#> [288] " 0.42391935E-01 4 4 OMEGA"
#> [289] " 0.63279530E+03 1 1 -LOGLIKE"
#> [290] " 0.10597484E+04 1 1 ELSOBJ "
#> [291] " "
#> [292] " EPSILON SHRINKAGES"
#> [293] " Epsilon shrinkage for CEps = 0.173121"
#> [294] ""
#> [295] " ETA SHRINKAGES"
#> [296] " eta-shrinkage( 1) = 0.11968"
#> [297] " eta-shrinkage( 2) = 0.01533"
#> [298] " eta-shrinkage( 3) = 0.99875"
#> [299] " eta-shrinkage( 4) = 0.28698"
#> [300] ""
#> [301] " RUNTIMES"
#> [302] ""
#> [303] " engine runtime (secs) = 0.141"
#> [304] " stderr runtime (secs) = 0.047"
#> [305] ""
#> [306] " EXITCODE 1"
#>
#> $dmp.txt
#> $coefficients
#> $coefficients$fixed
#> tvV tvCl tvV2 tvCl2 CEps
#> 15.393731 6.599838 41.190009 14.025770 0.161206
#>
#> $coefficients$random
#> $coefficients$random$Subject
#> nV nCl nV2 nCl2
#> 1 -0.140225458 0.15063262 -5.910225e-07 -0.06003932
#> 2 0.479233952 0.50712172 2.710213e-07 -0.13042943
#> 3 0.143254391 0.09407487 1.126426e-07 0.09854059
#> 4 0.210215296 0.40147284 1.100524e-07 -0.04484961
#> 5 -0.007981039 0.15095118 1.645053e-08 0.07616206
#> 6 0.019609501 0.56872533 4.714352e-07 -0.16320367
#> 7 0.109707888 0.13319750 3.831503e-07 0.14145686
#> 8 -0.255086214 -0.36509625 -3.621856e-08 -0.09254025
#> 9 0.203815159 0.35342619 -2.997510e-07 0.13252957
#> 10 -0.307170037 -0.62454337 -4.546653e-07 0.12483107
#> 11 0.057787931 0.17780777 -7.101335e-08 -0.37751129
#> 12 -0.438271913 -0.82453717 1.020711e-08 -0.03548031
#> 13 0.066650844 0.29173113 -1.066623e-07 -0.10274682
#> 14 0.189123315 -0.01438328 2.045132e-08 0.18575368
#> 15 -0.103658405 -0.25226211 1.289549e-07 0.02649993
#> 16 -0.289473668 -0.69305946 -1.137622e-07 0.17018928
#>
#>
#>
#> $logLik
#> [1] -632.795
#>
#> $nObs
#> [1] 112
#>
#> $nSubj
#> [1] 16
#>
#> $nParm
#> [1] 9
#>
#> $returnCode
#> [1] 1
#>
#> $omega
#> nV nCl nV2 nCl2
#> nV 0.06914959 0.0000000 0.000000e+00 0.00000000
#> nCl 0.00000000 0.1806635 0.000000e+00 0.00000000
#> nV2 0.00000000 0.0000000 4.649898e-08 0.00000000
#> nCl2 0.00000000 0.0000000 0.000000e+00 0.04239193
#>
#> $eta
#> mean var sd shrinkage
#> nV -3.904279e-03 5.358847e-02 2.314918e-01 0.11967934
#> nCl 3.453719e-03 1.751686e-01 4.185314e-01 0.01532516
#> nV2 -9.295596e-09 7.249374e-14 2.692466e-07 0.99875138
#> nCl2 -3.177354e-03 2.155219e-02 1.468067e-01 0.28697606
#>
#> $eta_shrinkage
#> $eta_shrinkage$Subject
#> nV_shrinkage nCl_shrinkage nV2_shrinkage nCl2_shrinkage
#> 1 0.1992617 0.01754401 0.9999968 0.4964577
#> 2 0.1881922 0.01542515 0.9999975 0.5871525
#> 3 0.2202274 0.01989199 0.9999974 0.5673752
#> 4 0.2237389 0.01556465 0.9999969 0.5225516
#> 5 0.2204068 0.01785632 0.9999970 0.5194621
#> 6 0.2155996 0.01325500 0.9999963 0.3964855
#> 7 0.2543358 0.01887873 0.9999971 0.5257477
#> 8 0.2307882 0.03079931 0.9999979 0.4565681
#> 9 0.2298989 0.01571451 0.9999969 0.5411340
#> 10 0.2304517 0.04437570 0.9999978 0.4948831
#> 11 0.2025019 0.01891869 0.9999977 0.4812590
#> 12 0.2580338 0.05548771 0.9999982 0.3839103
#> 13 0.2381110 0.01656682 0.9999969 0.4648553
#> 14 0.2408250 0.02284508 0.9999977 0.5729716
#> 15 0.2511555 0.02778062 0.9999976 0.4812069
#> 16 0.2536193 0.04912547 0.9999979 0.4540657
#>
#>
#> $sigma
#> CEps
#> CEps 0.02598737
#>
#> $method
#> [1] "FOCE-ELS"
#>
#> $varFix
#> tvV tvCl tvV2 tvCl2 CEps
#> tvV 1.298418676 0.705555361 0.350434817 -0.163831025 -0.0063183795
#> tvCl 0.705555361 0.519227114 0.159152149 -0.260228006 -0.0041458415
#> tvV2 0.350434817 0.159152149 1.128566179 0.174934558 0.0024225205
#> tvCl2 -0.163831025 -0.260228006 0.174934558 1.071421748 -0.0015977041
#> CEps -0.006318379 -0.004145842 0.002422521 -0.001597704 0.0003860419
#>
#> $omegaSE
#> nV nCl nV2 nCl2
#> nV 0.02437754 0.0000000 0.000000e+00 0.00000000
#> nCl 0.00000000 0.0530126 0.000000e+00 0.00000000
#> nV2 0.00000000 0.0000000 4.256802e-06 0.00000000
#> nCl2 0.00000000 0.0000000 0.000000e+00 0.03465556
#>
#> $Covariance
#> tvV tvCl tvV2 tvCl2 stdev0
#> tvV 1.298419e+00 7.055554e-01 3.504348e-01 -1.638310e-01 -6.318379e-03
#> tvCl 7.055554e-01 5.192271e-01 1.591521e-01 -2.602280e-01 -4.145841e-03
#> tvV2 3.504348e-01 1.591521e-01 1.128566e+00 1.749346e-01 2.422520e-03
#> tvCl2 -1.638310e-01 -2.602280e-01 1.749346e-01 1.071422e+00 -1.597704e-03
#> stdev0 -6.318379e-03 -4.145841e-03 2.422520e-03 -1.597704e-03 3.860419e-04
#> nV-nV 3.295851e-03 -3.974253e-03 3.460331e-03 -7.978365e-04 1.938038e-04
#> nCl-nV 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00
#> nCl-nCl -2.824936e-02 -2.128393e-02 1.239261e-02 3.031259e-03 5.552742e-04
#> nV2-nV 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00
#> nV2-nCl 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00
#> nV2-nV2 -2.012206e-06 -1.664998e-06 -1.420576e-06 1.872002e-06 -2.092889e-08
#> nCl2-nV 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00
#> nCl2-nCl 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00
#> nCl2-nV2 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00
#> nCl2-nCl2 1.078400e-02 5.464257e-03 -3.682945e-03 -1.141343e-02 -2.736106e-04
#> nV-nV nCl-nV nCl-nCl nV2-nV nV2-nCl nV2-nV2
#> tvV 3.295851e-03 0 -2.824936e-02 0 0 -2.012206e-06
#> tvCl -3.974253e-03 0 -2.128393e-02 0 0 -1.664998e-06
#> tvV2 3.460331e-03 0 1.239261e-02 0 0 -1.420576e-06
#> tvCl2 -7.978365e-04 0 3.031259e-03 0 0 1.872002e-06
#> stdev0 1.938038e-04 0 5.552742e-04 0 0 -2.092889e-08
#> nV-nV 5.942643e-04 0 8.076182e-04 0 0 1.640274e-08
#> nCl-nV 0.000000e+00 0 0.000000e+00 0 0 0.000000e+00
#> nCl-nCl 8.076182e-04 0 2.810335e-03 0 0 8.133134e-08
#> nV2-nV 0.000000e+00 0 0.000000e+00 0 0 0.000000e+00
#> nV2-nCl 0.000000e+00 0 0.000000e+00 0 0 0.000000e+00
#> nV2-nV2 1.640274e-08 0 8.133134e-08 0 0 1.812036e-11
#> nCl2-nV 0.000000e+00 0 0.000000e+00 0 0 0.000000e+00
#> nCl2-nCl 0.000000e+00 0 0.000000e+00 0 0 0.000000e+00
#> nCl2-nV2 0.000000e+00 0 0.000000e+00 0 0 0.000000e+00
#> nCl2-nCl2 -2.149110e-04 0 -6.455762e-04 0 0 8.716363e-09
#> nCl2-nV nCl2-nCl nCl2-nV2 nCl2-nCl2
#> tvV 0 0 0 1.078400e-02
#> tvCl 0 0 0 5.464257e-03
#> tvV2 0 0 0 -3.682945e-03
#> tvCl2 0 0 0 -1.141343e-02
#> stdev0 0 0 0 -2.736106e-04
#> nV-nV 0 0 0 -2.149110e-04
#> nCl-nV 0 0 0 0.000000e+00
#> nCl-nCl 0 0 0 -6.455762e-04
#> nV2-nV 0 0 0 0.000000e+00
#> nV2-nCl 0 0 0 0.000000e+00
#> nV2-nV2 0 0 0 8.716363e-09
#> nCl2-nV 0 0 0 0.000000e+00
#> nCl2-nCl 0 0 0 0.000000e+00
#> nCl2-nV2 0 0 0 0.000000e+00
#> nCl2-nCl2 0 0 0 1.201008e-03
#>
#> $Correlation
#> tvV tvCl tvV2 tvCl2 stdev0 nV-nV
#> tvV 1.1394818 0.8593004 0.2894919 -0.13890201 -0.28221567 0.11865073
#> tvCl 0.8593004 0.7205742 0.2079076 -0.34889530 -0.29283084 -0.22624916
#> tvV2 0.2894919 0.2079076 1.0623400 0.15908597 0.11606111 0.13361779
#> tvCl2 -0.1389020 -0.3488953 0.1590860 1.03509500 -0.07855953 -0.03161869
#> stdev0 -0.2822157 -0.2928308 0.1160611 -0.07855953 0.01964795 0.40462729
#> nV-nV 0.1186507 -0.2262492 0.1336178 -0.03161869 0.40462729 0.02437754
#> nCl-nV 0.0000000 0.0000000 0.0000000 0.00000000 0.00000000 0.00000000
#> nCl-nCl -0.4676514 -0.5571781 0.2200494 0.05524128 0.53310308 0.62493838
#> nV2-nV 0.0000000 0.0000000 0.0000000 0.00000000 0.00000000 0.00000000
#> nV2-nCl 0.0000000 0.0000000 0.0000000 0.00000000 0.00000000 0.00000000
#> nV2-nV2 -0.4148410 -0.5428148 -0.3141359 0.42485688 -0.25023354 0.15806771
#> nCl2-nV 0.0000000 0.0000000 0.0000000 0.00000000 0.00000000 0.00000000
#> nCl2-nCl 0.0000000 0.0000000 0.0000000 0.00000000 0.00000000 0.00000000
#> nCl2-nV2 0.0000000 0.0000000 0.0000000 0.00000000 0.00000000 0.00000000
#> nCl2-nCl2 0.2730860 0.2188162 -0.1000366 -0.31817271 -0.40183034 -0.25438763
#> nCl-nV nCl-nCl nV2-nV nV2-nCl nV2-nV2 nCl2-nV nCl2-nCl
#> tvV 0 -0.46765136 0 0 -4.148410e-01 0 0
#> tvCl 0 -0.55717807 0 0 -5.428148e-01 0 0
#> tvV2 0 0.22004937 0 0 -3.141359e-01 0 0
#> tvCl2 0 0.05524128 0 0 4.248569e-01 0 0
#> stdev0 0 0.53310308 0 0 -2.502335e-01 0 0
#> nV-nV 0 0.62493838 0 0 1.580677e-01 0 0
#> nCl-nV 0 0.00000000 0 0 0.000000e+00 0 0
#> nCl-nCl 0 0.05301260 0 0 3.604089e-01 0 0
#> nV2-nV 0 0.00000000 0 0 0.000000e+00 0 0
#> nV2-nCl 0 0.00000000 0 0 0.000000e+00 0 0
#> nV2-nV2 0 0.36040886 0 0 4.256802e-06 0 0
#> nCl2-nV 0 0.00000000 0 0 0.000000e+00 0 0
#> nCl2-nCl 0 0.00000000 0 0 0.000000e+00 0 0
#> nCl2-nV2 0 0.00000000 0 0 0.000000e+00 0 0
#> nCl2-nCl2 0 -0.35139494 0 0 5.908524e-02 0 0
#> nCl2-nV2 nCl2-nCl2
#> tvV 0 0.27308600
#> tvCl 0 0.21881623
#> tvV2 0 -0.10003658
#> tvCl2 0 -0.31817271
#> stdev0 0 -0.40183034
#> nV-nV 0 -0.25438763
#> nCl-nV 0 0.00000000
#> nCl-nCl 0 -0.35139494
#> nV2-nV 0 0.00000000
#> nV2-nCl 0 0.00000000
#> nV2-nV2 0 0.05908524
#> nCl2-nV 0 0.00000000
#> nCl2-nCl 0 0.00000000
#> nCl2-nV2 0 0.00000000
#> nCl2-nCl2 0 0.03465556
#>
#> $Inverse
#> tvV tvCl tvV2 tvCl2 stdev0
#> tvV 30.669119 -3.966410e+01 -1.472446e+00 -9.664792e+00 1.125657e+02
#> tvCl -39.664100 5.609381e+01 7.581991e-01 1.369820e+01 -1.164470e+02
#> tvV2 -1.472446 7.581991e-01 2.246036e+00 -9.930480e-01 3.508040e+01
#> tvCl2 -9.664792 1.369820e+01 -9.930480e-01 5.562615e+00 -6.494345e+01
#> stdev0 112.565750 -1.164470e+02 3.508040e+01 -6.494345e+01 7.303467e+03
#> nV-nV -788.509640 9.964689e+02 4.209836e+01 2.484849e+02 -3.346073e+03
#> nCl-nV 0.000000 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00
#> nCl-nCl 148.440890 -1.515144e+02 -5.078641e+01 -6.339149e+00 -1.363358e+03
#> nV2-nV 0.000000 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00
#> nV2-nCl 0.000000 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00
#> nV2-nV2 934123.710000 -1.112535e+06 4.188411e+05 -7.897050e+05 2.926216e+07
#> nCl2-nV 0.000000 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00
#> nCl2-nCl 0.000000 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00
#> nCl2-nV2 0.000000 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00
#> nCl2-nCl2 -233.724270 3.118519e+02 -7.591698e+00 1.062688e+02 -8.706470e+02
#> nV-nV nCl-nV nCl-nCl nV2-nV nV2-nCl nV2-nV2
#> tvV -7.885096e+02 0 1.484409e+02 0 0 9.341237e+05
#> tvCl 9.964689e+02 0 -1.515144e+02 0 0 -1.112535e+06
#> tvV2 4.209836e+01 0 -5.078641e+01 0 0 4.188411e+05
#> tvCl2 2.484849e+02 0 -6.339149e+00 0 0 -7.897050e+05
#> stdev0 -3.346073e+03 0 -1.363358e+03 0 0 2.926216e+07
#> nV-nV 2.320232e+04 0 -4.723265e+03 0 0 -2.495784e+07
#> nCl-nV 0.000000e+00 0 0.000000e+00 0 0 0.000000e+00
#> nCl-nCl -4.723265e+03 0 2.744943e+03 0 0 -1.016001e+07
#> nV2-nV 0.000000e+00 0 0.000000e+00 0 0 0.000000e+00
#> nV2-nCl 0.000000e+00 0 0.000000e+00 0 0 0.000000e+00
#> nV2-nV2 -2.495784e+07 0 -1.016001e+07 0 0 2.802425e+11
#> nCl2-nV 0.000000e+00 0 0.000000e+00 0 0 0.000000e+00
#> nCl2-nCl 0.000000e+00 0 0.000000e+00 0 0 0.000000e+00
#> nCl2-nV2 0.000000e+00 0 0.000000e+00 0 0 0.000000e+00
#> nCl2-nCl2 6.068787e+03 0 -4.660676e+02 0 0 -1.484096e+07
#> nCl2-nV nCl2-nCl nCl2-nV2 nCl2-nCl2
#> tvV 0 0 0 -2.337243e+02
#> tvCl 0 0 0 3.118519e+02
#> tvV2 0 0 0 -7.591698e+00
#> tvCl2 0 0 0 1.062688e+02
#> stdev0 0 0 0 -8.706470e+02
#> nV-nV 0 0 0 6.068787e+03
#> nCl-nV 0 0 0 0.000000e+00
#> nCl-nCl 0 0 0 -4.660676e+02
#> nV2-nV 0 0 0 0.000000e+00
#> nV2-nCl 0 0 0 0.000000e+00
#> nV2-nV2 0 0 0 -1.484096e+07
#> nCl2-nV 0 0 0 0.000000e+00
#> nCl2-nCl 0 0 0 0.000000e+00
#> nCl2-nV2 0 0 0 0.000000e+00
#> nCl2-nCl2 0 0 0 3.243844e+03
#>
#> $LLTable
#> id LL
#> 1 1 -39.17393
#> 2 2 -38.95108
#> 3 3 -36.39699
#> 4 4 -35.17740
#> 5 5 -37.60798
#> 6 6 -36.09242
#> 7 7 -38.38033
#> 8 8 -41.61708
#> 9 9 -39.13651
#> 10 10 -42.74386
#> 11 11 -40.97003
#> 12 12 -51.55430
#> 13 13 -36.23927
#> 14 14 -37.40598
#> 15 15 -37.94734
#> 16 16 -43.40081
#>
#> $`.\\test.mdl`
#> [1] "test(){"
#> [2] " cfMicro(A1,Cl/V, Cl2/V, Cl2/V2)"
#> [3] " dosepoint(A1)"
#> [4] " C = A1 / V"
#> [5] " error(CEps=0.1)"
#> [6] " observe(CObs=C * ( 1 + CEps))"
#> [7] " stparm(V = tvV * exp(nV))"
#> [8] " stparm(Cl = tvCl * exp(nCl))"
#> [9] " stparm(V2 = tvV2 * exp(nV2))"
#> [10] " stparm(Cl2 = tvCl2 * exp(nCl2))"
#> [11] " fixef( tvV = c(,1,))"
#> [12] " fixef( tvCl = c(,1,))"
#> [13] " fixef( tvV2 = c(,1,))"
#> [14] " fixef( tvCl2 = c(,1,))"
#> [15] " ranef(diag(nV,nCl,nV2,nCl2) = c(1,1,1,1))"
#> [16] "}"
#>
#> $cols1.txt
#> [1] "id(\"Subject\")"
#> [2] "time(\"Act_Time\")"
#> [3] "dose(A1<-\"Amount\")"
#> [4] "obs(CObs<-\"Conc\")"
#> [5] "table(file=\"posthoc.csv\", V, Cl, V2, Cl2, mode=keep)"
#>
#> attr(,"class")
#> [1] "nlme"
#>
finalModelVPC <- copyModel(model,
acceptAllEffects = TRUE,
modelName = "model_VPC",
workingDir = tempdir())
# View the model
print(finalModelVPC)
#>
#> Model Overview
#> -------------------------------------------
#> Model Name : model_VPC
#> Working Directory : C:\Users\jcraig\AppData\Local\Temp\RtmpI3HPMy
#> Is population : TRUE
#> Model Type : PK
#>
#> PK
#> -------------------------------------------
#> Parameterization : Clearance
#> Absorption : Intravenous
#> Num Compartments : 2
#> Dose Tlag? : FALSE
#> Elimination Comp ?: FALSE
#> Infusion Allowed ?: FALSE
#> Sequential : FALSE
#> Freeze PK : FALSE
#>
#> PML
#> -------------------------------------------
#> test(){
#> cfMicro(A1,Cl/V, Cl2/V, Cl2/V2)
#> dosepoint(A1)
#> C = A1 / V
#> error(CEps=0.161205984597858)
#> observe(CObs=C * ( 1 + CEps))
#> stparm(V = tvV * exp(nV))
#> stparm(Cl = tvCl * exp(nCl))
#> stparm(V2 = tvV2 * exp(nV2))
#> stparm(Cl2 = tvCl2 * exp(nCl2))
#> fixef( tvV = c(,15.393730662165,))
#> fixef( tvCl = c(,6.59983805751661,))
#> fixef( tvV2 = c(,41.1900094327716,))
#> fixef( tvCl2 = c(,14.0257699237193,))
#> ranef(diag(nV,nCl,nV2,nCl2) = c(0.0691495878371252,0.180663528781695,4.64989751058228e-08,0.0423919347837654))
#>
#> }
#>
#> Structural Parameters
#> -------------------------------------------
#> V Cl V2 Cl2
#> -------------------------------------------
#> Observations:
#> Observation Name : CObs
#> Effect Name : C
#> Epsilon Name : CEps
#> Epsilon Type : Multiplicative
#> Epsilon frozen : FALSE
#> is BQL : FALSE
#> -------------------------------------------
#> Column Mappings
#> -------------------------------------------
#> Model Variable Name : Data Column name
#> id : Subject
#> time : Act_Time
#> A1 : Amount
#> CObs : Conc
#>
# Set up VPC arguments to have PRED outputted to simulation output dataset "predout.csv"
vpcSetup <- NlmeVpcParams(outputPRED = TRUE)
# Run VPC using the default host, default values for the relevant NLME engine arguments
finalVPCJob <-
vpcmodel(model = finalModelVPC, vpcParams = vpcSetup, )
#> Using localhost without parallelization.
#>
#> NLME Job
#> sharedDirectory is not given in the host class instance.
#> Valid NLME_ROOT_DIRECTORY is not given, using current working directory:
#> C:/Repos/R-RsNLME/docs/reference
#>
#> Compiling 1 of 1 NLME models
#> TDL5 version: 25.7.1.1
#>
#> Status: OK
#> License expires: 2026-08-06
#> Refresh until: 2025-09-05 07:01:25
#> Current Date: 2025-08-06
#> The model compiled
#>
#> Trying to generate job results...
#>
#> Copying tables predcheck0.csv, predcheck0_cat.csv, predcheck1.csv, predcheck1_cat.csv, predcheck1_npd.csv, predcheck2.csv, predcheck2_cat.csv, predcheck2_npd.csv, predcheck_bql.csv, predout.csv
#>
#> Finished summarizing results. Transferring data and loading the results...
#> Done generating job results.
#>
#> VPC/Simulation results are ready in C:/Users/jcraig/AppData/Local/Temp/RtmpI3HPMy
#> Loading the results
#> Loading predcheck_bql.csv
#> Loading predcheck0.csv
#> Loading predcheck0_cat.csv
#> Loading predcheck1.csv
#> Loading predcheck1_cat.csv
#> Loading predcheck2.csv
#> Loading predcheck2_cat.csv
#> Loading predout.csv
# }