Thanks! ... View more

It worked! Thank you, could you tell me why I cannot have class statement? ... View more

Here is my code: data wafdata;   set wafdata;   PARM_fit = PARM;     if Tukey1 = 'fail' then PARM = .; run;   proc glm data = wafdata noprint;   by F_ParmName;   class x y;   model PARM = x                            y                            x * x                            x * y                            y * y                            x * x * x                            x * x * y                            x * y * y                            y * y * y                            x * x * x * x                            x * x * x * y                            x * x * y * y                            x * y * y * y                            y * y * y * y /p;   output out = poly p = prediction r = residual;   ods graphics off; run;   data poly; set poly; if missing(residual) then residual = PARM_fit - prediction; run;   Part of table wafdata is as this, also attached a xlsx file for whole data PARMNAME                PARM                x                y                    Field                  F_ParmName                Filter MR5                            5.702313        -11.475        9.28                 10A                   10A_MR5                       pass MR5                            5.701248         -9.945         9.28                 10A                   10A_MR5                      pass MR5                            5.687179         -8.415         9.28                 11A                    11A_MR5                      pass MR5                            -7777.7            -6.885         9.28                 11A                    11A_MR5                      fail MR5                            5.632577         -5.355         9.28                 11A                    11A_MR5                      pass MR5                            5.627558         -3.825         9.28                 11A                    11A_MR5                      pass MR5                            5.653947         -2.295         9.28                 12A                    12A_MR5                      pass MR5                            5.648395         -0.765         9.28                 12A                    12A_MR5                      pass MR5                            5.740413          -22.185      8.7                   14A                     14A_MR5                    pass MR5                            5.70026            -20.655      8.7                   14A                    14A_MR5                     pass MR5                            5.715317         -19.125       8.7                   14A                     14A_MR5                    pass MRSUB                       7295.616           18.36       39.73               11A                     11A_MRSUB                pass MRSUB                       7404.475           19.89       39.73               11A                      11A_MRSUB               pass MRSUB                       7431.816           21.42        39.73               11A                     11A_MRSUB               pass MRSUB                      1500000000       0.765        51.91               13A                     13A_MRSUB               fail MRSUB                       7061.483           2.295        51.91               13A                     13A_MRSUB               pass MRSUB                       7110.674           3.825        51.91               13A                     13A_MRSUB                pass MRSUB                      7063.209            5.355        51.91               14A                     14A_MRSUB                pass MRSUB                      7122.559            6.885        51.91               14A                     14A_MRSUB                pass MRSUB                      7108.894            8.415        51.91               14A                     14A_MRSUB                pass   Thank you! ... View more

I also thought By statement should not affect the solution, however, case 1. use your solution with by Field              fitting works, have prediction and residual but fitting result is wired.  case 2. use your solution with by F_ParmName or ParmName              fitting works, but missing part has no prediction, thus cannot calculate residual. ... View more

My actual dataset is much more conplicated, there are many other columns beside the ones I listed, for more detail in the "by" group, let me describe as below: PARM       x         y          Field        ParmName      F_ParmName               Filter 3               x0      y0           1A              AA                    1A_AA                      pass 3.1            x1      y1           1A              BB                    1A_BB                      pass 2.9            x2      y2           1B              AA                    1B_AA                      pass 5               x3      y3           1B              BB                    1B_BB                      fail   Field and ParmName is independent, I contact them as a new column F_ParmName, PARM is the actual value of each ParmName. The dataset is much larger so there are much more than one row of data for each category as I showed. I used by Field fitting, PaigeMiller's way is working but fitting result is wired since by group is not right, but when I try by F_ParmName fitting, the PARM '.' row will have prediction '.' result.  ... View more

I there any requirement for "by" in fitting in your way? I changed "by" group and it works... both "by" groups are character bu the working one is simpler (not working one is included in working one) ... View more

I tried your way but it seems not work, predictions for PARM = '.' points are also '.' ... View more

thx, let me try ... View more

if I use WHERE and do the fit, will SAS calculate residual for me on those excluded points? ... View more