"INPUT"  P, Q, P * [ Y, X1, X2, X3, X4, X5, N ];

"MODEL"  Y = C + N1055 * INDICATOR (1055, N, 1055) +
                 N1091 * INDICATOR (1.0910000000#+03, N, 1.0910000000#+03) +
                 N1105 * INDICATOR (1105, ROUND(N), 1105) +
                 N1136 * ( 1 - SIGN(ABS(ENTIER(N+.5)-1136)) ) +
      BED * X1 + TAS * X2 + PER * X3 + MED * X4 + FAC * X5 +
      BED2 * (X1*X1) + TAS2 * (X2*X2) + PER2 * (X3*X3) +
      MED2 * (X4*X4) + FAC2 * (X5*X5) +
      BEDTAS * (X1*X2) + BEDPER * (X1*X3) + BEDMED * (X1*X4) +
      BEDFAC * (X1*X5) + TASPER * (X2*X3) + TASMED * (X2*X4) +
      TASFAC * (X2*X5) + PERMED * (X3*X4) + PERFAC * (X3*X5) +
      MEDFAC * (X4*X5);

"OPTIONS"  PROCESS SUBMODELS(10,15),
           RESIDUAL ANALYSIS,
           TRANSFORMED DATA MATRIX,
           CORRELATION MATRIX;

"RUN"
"EXIT"