alag problems

[mattfidler]

It seems to have problem with between occasion variability (BOV) of the model. When I tried the same model without BOV, I got different errors.

> setwd("E:\\MicuRX\\poppk1\\run401f.nm7")
> mod <- nonmem2rx("run401f.mod")
ℹ getting information from  'run401f.mod'
ℹ reading in xml file
ℹ done
ℹ reading in phi file
ℹ done
ℹ reading in lst file
ℹ abbreviated list parsing
ℹ done
ℹ done
ℹ splitting control stream by records
ℹ done
ℹ Processing record $INPUT
ℹ Processing record $MODEL
ℹ Processing record $THETA
ℹ Processing record $OMEGA
ℹ Processing record $SIGMA
ℹ Processing record $PROBLEM
ℹ Processing record $DATA
ℹ Processing record $SUBROUTINES
ℹ Processing record $PK
ℹ Processing record $DES
ℹ Processing record $ERROR
ℹ Processing record $ESTIMATION
ℹ Ignore record $ESTIMATION
ℹ Processing record $COVARIANCE
ℹ Ignore record $COVARIANCE
ℹ Processing record $TABLE
ℹ change initial estimate of `theta1` to `16.9184465874669`
ℹ change initial estimate of `theta2` to `29.0692797324232`
ℹ change initial estimate of `theta3` to `2.43235177034855`
ℹ change initial estimate of `theta4` to `1.34653279917032`
ℹ change initial estimate of `theta5` to `0.304490361221458`
ℹ change initial estimate of `theta6` to `2.5`
ℹ some correlations may have been dropped for the variables: `eta1`, `eta2`
ℹ the piping should specify the needed covariances directly
ℹ change initial estimate of `eta1` to `0.0722553408614083`
ℹ add covariance between `eta2` and `eta1` with initial estimate `0.0650039024597345`
ℹ change initial estimate of `eta2` to `0.0700334194752202`
ℹ some correlations may have been dropped for the variables: `eta3`, `eta4`
ℹ the piping should specify the needed covariances directly
ℹ change initial estimate of `eta3` to `0.16268513883389`
ℹ add covariance between `eta4` and `eta3` with initial estimate `0.0352032351476214`
ℹ change initial estimate of `eta4` to `0.376883716036017`
ℹ read in nonmem input data (for model validation): E:\MicuRX\poppk1\nm1.csv
ℹ ignoring lines that begin with 'C'
ℹ applying names specified by $INPUT
ℹ subsetting accept/ignore filters code: .data[-which((.data$FLAG == 1 | .data$FED == 0)),]
ℹ done
ℹ read in nonmem IPRED data (for model validation): E:\MicuRX\poppk1\run401f.nm7\run401f.fit
ℹ done
ℹ changing most variables to lower case
ℹ done
ℹ replace theta names
ℹ done
ℹ replace eta names
ℹ done
ℹ renaming compartments
ℹ done
ℹ solving ipred problem
Error : 'alag(.)'/'rate(.)'/'dur(.)' cannot depend on the state values
ℹ done
ℹ solving pred problem
Error : 'alag(.)'/'rate(.)'/'dur(.)' cannot depend on the state values
ℹ done

Originally posted by @sliao999 in https://github.com/nlmixr2/nonmem2rx/issues/125#issuecomment-1591788488

[mattfidler]

From @sliao999

'alag(.)'/'rate(.)'/'dur(.)' none of these parameters were defined in the mod file. Any ideal to resolve this problem?

[mattfidler]

Unfortunately, without the data and model, I cannot figure out what is happening here.

I know you cannot share proprietary data and models (and I cannot take them either), but if you have permission to share an open source version (even with things redacted simulated etc), then I can help with this issue.

[mattfidler]

Note that this may be related:

https://github.com/nlmixr2/rxode2/issues/526

[mattfidler]

For now I will close this issue since I cannot figure out how to reproduce it.