RE: [NMusers] $PRIOR with normal for OMEGA?

2016-11-11 Thread Mats Karlsson
Hi Mark, As I’m sure you know NONMEM has the TNPRI functionality for that. If you want to use it without the TNPRI functionality, you can estimate OMEGA as THETA and use a multivariate normal for that. If you have off-diagonal elements, you may want to do a Cholesky transformation (you can get

Re: [NMusers] $PRIOR with normal for OMEGA?

2016-11-11 Thread Jakob Ribbing
Hi Mark, Indeed there is: As an alternative to NWPRI, there is the TNPRI subroutine that you can use with $PRIOR (frequentist prior). This functionality is tripple normal, with regards to thetas, omegas and sigma(s). I will describe this more in detail than Mark would need (hopefully for the be

RE: [NMusers] $PRIOR with normal for OMEGA?

2016-11-11 Thread Åstrand , Magnus
Hi To add on Martin’s suggestion, one item to think of when using theta for estimating the variances of eta is to use log scaled STD as your model parameter, so instead of THETA*ETA (with OMEGA fixed) use exp(THETA)*ETA. This way you would assume a log-normal prior for the standard deviation of E

Re: [NMusers] $PRIOR with normal for OMEGA?

2016-11-11 Thread Mark Sale
thanks Mats, I've never had much luck getting TNPRI to work, pretty limited examples in the docs. I actually not not want to run NONMEM with normal OMEGA PRIOR, just want to validate some other work where I calculate the PRIOR penalty from a normal, and thought this might be a practical way t