Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

alag problems #143

Closed
mattfidler opened this issue Jun 14, 2023 · 4 comments
Closed

alag problems #143

mattfidler opened this issue Jun 14, 2023 · 4 comments

Comments

@mattfidler
Copy link
Member

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 filedonereading in phi filedonereading in lst fileabbreviated list parsingdonedonesplitting control stream by recordsdoneProcessing record $INPUTProcessing record $MODELProcessing record $THETAProcessing record $OMEGAProcessing record $SIGMAProcessing record $PROBLEMProcessing record $DATAProcessing record $SUBROUTINESProcessing record $PKProcessing record $DESProcessing record $ERRORProcessing record $ESTIMATIONIgnore record $ESTIMATIONProcessing record $COVARIANCEIgnore record $COVARIANCEProcessing record $TABLEchange 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 directlychange 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 directlychange 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.csvignoring lines that begin with 'C'applying names specified by $INPUTsubsetting accept/ignore filters code: .data[-which((.data$FLAG == 1 | .data$FED == 0)),]
ℹ doneread in nonmem IPRED data (for model validation): E:\MicuRX\poppk1\run401f.nm7\run401f.fitdonechanging most variables to lower casedonereplace theta namesdonereplace eta namesdonerenaming compartmentsdonesolving ipred problem
Error : 'alag(.)'/'rate(.)'/'dur(.)' cannot depend on the state valuesdonesolving pred problem
Error : 'alag(.)'/'rate(.)'/'dur(.)' cannot depend on the state valuesdone

Originally posted by @sliao999 in #125 (comment)

@mattfidler
Copy link
Member Author

From @sliao999

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

@mattfidler
Copy link
Member Author

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
Copy link
Member Author

Note that this may be related:

nlmixr2/rxode2#526

@mattfidler mattfidler mentioned this issue Jun 14, 2023
@mattfidler
Copy link
Member Author

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

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant