Hi all,
I'm running into an issue with my code. Briefly, I'm looking at how certain psychological constructs impact the effectiveness of a behavioral treatment. One of the treatments is time-out, where observers watch the kid for 30 seconds, and then record behavior for 10 seconds, and repeats this until the time-out ends. Naturally, it's an unbalanced dataset with some kids getting a few time-outs, and others getting many time-outs. I'm approaching it in a multi-level design:
Level 1 = the interval where data are collected (ranges from 1-197)
Level 2 = the unique ID # to identify the time-out (there are 545 unique time-outs)
Level 3 = the unique ID # to identify the child
So, interval data is nested in time-out, which is nested in child.
I'm running this in a stepwise approach. My first two models converged (see end of message). But model 3, where I want to include the full nesting term provides the following error:
1 OPTIONS NONOTES NOSTIMER NOSOURCE NOSYNTAXCHECK;
SYMBOLGEN: Macro variable _SASWSTEMP_ resolves to /home/u60849988/.sasstudio/.images/495a0ef2-2874-447e-a1c1-b64f1cce7501
SYMBOLGEN: Some characters in the above value which were subject to macro quoting have been unquoted for printing.
SYMBOLGEN: Macro variable GRAPHINIT resolves to GOPTIONS RESET=ALL GSFNAME=_GSFNAME;
68
69 PROC GLIMMIX DATA=impstdmrpsra2 METHOD=LAPLACE;
70 CLASS MRPSID UnqTONum IntervalNum Meds TOTxCond;
71 MODEL Effective = meds TOTxCond week ptadhd ptcp pticu24 pticu24*TOTxCond / SOLUTION CL;
72 RANDOM intercept / SUBJECT=MRPSID;
73 RANDOM intercept / SUBJECT= UnqTONum(MRPSID);
74 RANDOM intercept / SUBJECT=IntervalNum(UnqTONum(MRPSID)) TYPE=UN;
_
22
200
ERROR 22-322: Syntax error, expecting one of the following: a name, ), *.
ERROR 200-322: The symbol is not recognized and will be ignored.
Importantly, this occurs whether I use a singular random statement (i.e., RANDOM intercept / SUBJECT = IntervalNum(UnqTONum(MRPSID)) or multiple random statements as shown below. This is happening across both continuous and categorical variables. Any feedback or thoughts would be greatly appreciated. Thanks all.
Model 3 - not converging, receiving above error.
PROC GLIMMIX DATA=impstdmrpsra2 METHOD=LAPLACE;
CLASS MRPSID UnqTONum IntervalNum Meds TOTxCond;
MODEL Effective = meds TOTxCond week ptadhd ptcp pticu24 pticu24*TOTxCond / SOLUTION CL;
RANDOM intercept / SUBJECT=MRPSID;
RANDOM intercept / SUBJECT= UnqTONum(MRPSID);
RANDOM intercept / SUBJECT=IntervalNum(UnqTONum(MRPSID)) TYPE=UN;
options mprint mlogic symbolgen ;
RUN;
Model 1 - converges.
PROC GLIMMIX DATA=impstdmrpsra2 METHOD=LAPLACE;
CLASS MRPSID Meds TOTxCond;
MODEL Effective = meds TOTxCond week ptadhd ptcp pticu24 pticu24*TOTxCond / SOLUTION CL;
RANDOM intercept / SUBJECT=MRPSID TYPE=UN;
RUN;
Model 2 - converges.
PROC GLIMMIX DATA=impstdmrpsra2 METHOD=LAPLACE;
CLASS MRPSID UnqTONum Meds TOTxCond;
MODEL Effective = meds TOTxCond week ptadhd ptcp pticu24 pticu24*TOTxCond / SOLUTION CL;
RANDOM intercept / SUBJECT=UnqTONum(MRPSID) TYPE=UN;
RUN;
