if-the-statement for 3 class categorical covariate
Dear all
I am shoked because NONMEM sents a crash message when i factor in
covariate modeling using a 4 class categorical variable genotype:
IF (GENE.EQ.3) THEN ;GG
CL = THETA(1)* EXP(ETA(1))
ELSE IF (GENE.EQ.2)THEN ;GT
CL = THETA(6)* EXP(ETA(3))
IF (GENE.EQ.1) THEN ;TT
CL = THETA(7)* EXP(ETA(4))
IF (GENE.EQ.4) THEN ;MISSING
CL = THETA(7)* EXP(ETA(4))
END IF
However if i binarize the code, nonmem runs successfully
IF (GENE.EQ.3) THEN
CL = THETA(1)* EXP(ETA(1))
ELSE
CL = THETA(6)* EXP(ETA(3))
END IF
Thanks in advance
--
Bernard Ngara
Trainee Biostatistician
Africa Institute of Biomedical Science and Technology
P. O Box 2294
Harare
Zimbabwe
+263 772 160 896