Increasing the umber of sub-population in a mixture model
From: "HUTMACHER, MATTHEW [Non-Pharmacia/1825]" <matthew.hutmacher@pharmacia.com>
Subject: Increasing the umber of sub-population in a mixture model
Date: Fri, 25 May 2001 13:07:28 -0500
Hello,
I have a very complicated data pattern. Many of you have given great advice on how to use a mixture model at the IOV stage and I appreciate it. I have found, that I need more sub-populations to characterize the lag-time distribution in my model (the distribution of me predicted etas are still very skewed). I have tried to increase the number of sub-populations and I get this error message:
AN ERROR WAS FOUND IN THE CONTROL STATEMENTS.
AN ERROR WAS FOUND ON LINE 43 AT THE APPROXIMATE POSITION NOTED:
P(5)=THETA(.)
X
THE CHARACTERS IN ERROR ARE: 5
332 $MIX: NSPOP AND INDICES I OF P(I) MAY ONLY BE 1, 2, 3, OR 4.
STOP 4 statement executed
Is there a way to increase the number of allowable sub-populations.
I have extended the code Dr Sheiner supplied to me for two sub-populations (at IOV level) to three:
$MIX
P(1)=THETA(.)
P(2)=1.-THETA(.)
P(3)=THETA(.)
P(4)=1.-THETA(.)
P(5)=THETA(.)
P(6)=1.-THETA(.)
NSPOP=6
Any help in this would be appreciated.
Thanks again.
Matt