$MIX

Marks the beginning of abbreviated code for the MIX routine.

1
2
3
4
5
6
7
8


 $MIX
 abbreviated code

;# Example
 $MIX
 NSPOP=2
 P(1)=THETA(5)
 P(2)=1-THETA(5)

The $MIX record is used to describe the mixture probabilities of a mixture model. It is evaluated with individual records.

See abbreviated code for general rules. Specific rules follow.

NSPOP: an integer variable. The number of sub-populations. Must be given a value when ICALL=1.
P(i): for each i (i=1, … , NSPOP), P(i) is the modeled fraction of the population in the ith subpopulation. The sum of the P(i) should equal 1. In principle, the P(i) can change from individual to individual. If for a given individual, the second (for example) subpopulation doesn't apply, then set P(2)=0 for that individual. MIXP may be coded instead of P; they are synonymous in the context of the $MIX block.

THETA(n).
MIX-defined items that appeared earlier as left-hand quantities.
ICALL
- ICALL=0: Run initialization.
- ICALL=1: Problem initialization.
- ICALL=2: Normal call.
- ICALL=4: Simulation.
Data items listed in DATA option of $CONTR record. E.g., assume that the following records appear in the control stream prior to the $MIX block:
1 2

$INPUT ... STDY ... $CONTR DATA=(STDY)
Then STDY may be used on the right in $MIX. STDY and STDY(1) both refer to the value of STDY on the first observation record of the individual record. STDY(i) refers to the value of STDY on the i'th observation record of the individual record. TEMPLT may be used. (See MIX CONTR: TEMPLT).

Certain variables in FORTRAN modules can be used. (See Variables)

May not include:

EXIT, CALL , DO WHILE / ENDDO statements

COMRES, CALLFL pseudo-variables

Forbidden Variable Names: NEWIND,ETA(i),EPS(i),ERR(i) COM(i)

Variables defined in $MIX are not listed in module NMPRD4 and may not be displayed in $TABLE and $SCATTER.

Follows $SUBROUTINES.