data one; input trial $ year mgha mse; cards; e502 1971 2.37 0.041 e502 1972 1.708 0.0308 e502 1973 . . e502 1974 1.9 0.0335 e502 1975 2.92 0.0387 e502 1976 2.51 0.0396 e502 1977 1.85 0.034 e502 1978 2.409 0.125 e502 1979 2.64 0.318 e502 1980 3.01 0.081 e502 1981 2.304 0.118 e502 1982 1.96 0.0502 e502 1983 3.03 0.084 e502 1984 2.758 0.0825 e502 1985 2.115 0.0325 e502 1986 2.91 0.023 e502 1987 2.589 0.078 e502 1988 3.68 0.172 e502 1989 2.46 0.065 e502 1990 3.018 0.0422 e502 1991 1.889 0.03206 e502 1992 2.392 0.0539 e502 1993 2.22 0.07037 e502 1994 1.968 0.078 e502 1995 2.65 0.0513 e502 1996 2.011 0.198 e502 1997 2.525 0.264 e502 1998 3.09 0.0456 e502 1999 2.67 0.092 e502 2000 2.537 0.147 e502 2001 1.766 0.0926 e502 2002 2.89 0.157 e502 2003 4.98 0.228 e502 2004 3.346 0.098 e502 2005 2.418 0.155 e502 2006 2.64 0.359 e502 2007 2.982 0.2714 e502 2008 4.96 0.253 e502 2009 3.099 0.0789 e502 2010 1.544 0.0255 e502 2011 2.582 0.299 e502 2012 3.43 0.0733 e502 2013 2.426 0.1127 e502 2014 1.938 0.0313 e222 1969 2.51 0.119 e222 1970 1.62 0.071 e222 1971 2.04 0.04 e222 1972 1.52 0.072 e222 1973 3.46 0.089 e222 1974 . . e222 1975 . . e222 1976 1.29 0.112 e222 1977 1.3 0.209 e222 1978 0.738 0.052 e222 1979 3.059 0.075 e222 1980 2.135 0.043 e222 1981 2.26 0.017 e222 1982 2.22 0.0528 e222 1983 1.14 0.0507 e222 1984 2.89 0.128 e222 1985 1.48 0.174 e222 1986 0.889 0.034 e222 1987 0.845 0.0184 e222 1988 1.54 0.025 e222 1989 1.48 0.039 e222 1990 1.579 0.0605 e222 1991 1.86 0.0312 e222 1992 1.86 0.0344 e222 1993 1.64 0.0606 e222 1994 1.51 0.0217 e222 1995 0.27 0.0103 e222 1996 1.15 0.017 e222 1997 1.376 0.032 e222 1998 1.419 0.0285 e222 1999 1.53 0.072 e222 2000 2.15 0.0607 e222 2001 1.688 0.0807 e222 2002 2.201 0.0289 e222 2003 2.549 0.095 e222 2004 2.55 0.229 e222 2005 2.166 0.038 e222 2006 0.827 0.0136 e222 2007 0.447 0.015 e222 2008 1.77 0.241 e222 2009 . . e222 2010 2.12 0.015 e222 2011 1.24 0.029 e222 2012 2.14 0.0501 e222 2013 0.706 0.0325 e222 2014 2.39 0.046 proc print; data one; set one; if trial = 'e502' then delete; proc rsreg data = one out = two; by trial; model mse = mgha year /predict; proc g3grid data = two out = three; grid mgha*year=mse/spline; proc g3d data = three gout=new; plot mgha*year=mse; run;