Multivariate Modeling: Analysis of Long-Term Temperature Trends

/*--------------------------------------------------------------
                    SAS Sample Library

       Name: ssmex12.sas
       Description: Example program from SAS/ETS User's Guide,
              The SSM Procedure
       Title: Multivariate Modeling: Analysis of Long-Term Temperature Trends
       Product: SAS/ETS Software
        Keys: Multivariate Modeling,  Long-Term Temperature Trends
        PROC: SSM
        Notes:

--------------------------------------------------------------*/

 title;

data Temp;
input UAH CRU GISS @@;
date = intnx('month', '01jan1850'd, _n_-1);
format date date.;
datalines;
.    8.243    .
.    9.733    .
.    9.591    .
.    9.221    .
.    9.448    .
.    10.13    .
.    9.348    .
.    9.709    .
.    9.771    .
.    9.781    .
.    9.621    .
.    10.304    .
.    10.57    .
.    10.269    .
.    9.289    .
.    10.402    .
.    9.643    .
.    9.852    .
.    9.991    .
.    10.424    .
.    10.031    .
.    10.655    .
.    9.47    .
.    9.267    .
.    9.774    .
.    10.166    .
.    9.985    .
.    9.021    .
.    9.902    .
.    9.86    .
.    10.029    .
.    9.601    .
.    9.822    .
.    9.205    .
.    9.499    .
.    10.634    .
.    10.055    .
.    9.078    .
.    8.567    .
.    9.35    .
.    9.675    .
.    9.637    .
.    10.32    .
.    10.167    .
.    9.599    .
.    9.948    .
.    9.65    .
.    9.176    .
.    9.746    .
.    10.34    .
.    9.979    .
.    9.956    .
.    10.489    .
.    10.034    .
.    10.577    .
.    10.423    .
.    10.193    .
.    10.425    .
.    9.282    .
.    10.131    .
.    10.328    .
.    8.82    .
.    9.524    .
.    10.53    .
.    10.337    .
.    10.219    .
.    10.309    .
.    10.266    .
.    9.809    .
.    10.226    .
.    9.464    .
.    8.84    .
.    9.585    .
.    9.42    .
.    9.279    .
.    10.325    .
.    9.576    .
.    10.276    .
.    9.644    .
.    9.98    .
.    9.671    .
.    9.499    .
.    8.636    .
.    9.678    .
.    8.946    .
.    10.45    .
.    9.475    .
.    9.039    .
.    9.211    .
.    9.949    .
.    10.129    .
.    10.214    .
.    9.861    .
.    9.651    .
.    9.416    .
.    10.695    .
.    10.32    .
.    9.001    .
.    9.55    .
.    9.739    .
.    9.694    .
.    10.141    .
.    9.536    .
.    9.58    .
.    9.514    .
.    9.933    .
.    8.504    .
.    9.672    .
.    9.733    .
.    9.83    .
.    9.788    .
.    9.888    .
.    9.966    .
.    9.87    .
.    9.879    .
.    9.815    .
.    9.275    .
.    10.002    .
.    9.587    .
.    9.184    .
.    10.34    .
.    9.263    .
.    9.044    .
.    9.645    .
.    9.709    .
.    10.106    .
.    9.876    .
.    10.033    .
.    9.75    .
.    9.873    .
.    9.115    .
.    8.992    .
.    8.753    .
.    10.203    .
.    10.18    .
.    9.636    .
.    9.18    .
.    10.245    .
.    9.695    .
.    9.792    .
.    9.75    .
.    9.797    .
.    9.425    .
.    9.739    .
.    9.106    .
.    8.834    .
.    9.77    .
.    9.701    .
.    9.792    .
.    9.714    .
.    9.467    .
.    9.138    .
.    9.61    .
.    9.39    .
.    9.197    .
.    9.204    .
.    10.827    .
.    10.195    .
.    9.519    .
.    9.64    .
.    9.654    .
.    9.573    .
.    9.541    .
.    9.449    .
.    9.237    .
.    9.311    .
.    9.407    .
.    9.634    .
.    8.722    .
.    9.272    .
.    9.461    .
.    9.39    .
.    9.534    .
.    9.969    .
.    9.883    .
.    9.352    .
.    9.815    .
.    8.854    .
.    9.078    .
.    9.094    .
.    9.665    .
.    8.924    .
.    8.998    .
.    10.056    .
.    9.788    .
.    10.019    .
.    10.043    .
.    9.681    .
.    10.218    .
.    9.494    .
.    10.084    .
.    9.49    .
.    10.279    .
.    9.615    .
.    9.199    .
.    9.856    .
.    9.529    .
.    9.865    .
.    9.767    .
.    9.544    .
.    9.926    .
.    9.351    .
.    9.173    .
.    9.591    .
.    9.458    .
.    10.44    .
.    8.572    .
.    9.646    .
.    9.245    .
.    9.841    .
.    9.631    .
.    9.827    .
.    9.947    .
.    9.961    .
.    9.51    .
.    9.312    .
.    8.749    .
.    9.148    .
.    9.869    .
.    9.411    .
.    10.211    .
.    10.066    .
.    10.238    .
.    9.929    .
.    9.657    .
.    9.916    .
.    9.244    .
.    9.764    .
.    10.062    .
.    10.464    .
.    9.091    .
.    9.702    .
.    9.613    .
.    9.748    .
.    9.665    .
.    10.071    .
.    9.956    .
.    9.044    .
.    9.485    .
.    9.828    .
.    10.295    .
.    9.568    .
.    8.996    .
.    9.687    .
.    10.011    .
.    10.394    .
.    10.093    .
.    9.302    .
.    9.83    .
.    9.812    .
.    9.882    .
.    8.628    .
.    9.358    .
.    9.159    .
.    9.881    .
.    9.764    .
.    9.696    .
.    9.82    .
.    10.102    .
.    10.211    .
.    9.548    .
.    9.304    .
.    8.912    .
.    9.185    .
.    10.134    .
.    9.752    .
.    9.637    .
.    9.888    .
.    9.942    .
.    10.103    .
.    9.973    .
.    9.462    .
.    9.991    .
.    9.918    .
.    9.709    .
.    9.516    .
.    9.912    .
.    9.365    .
.    9.534    .
.    9.164    .
.    9.735    .
.    10.514    .
.    9.78    .
.    10.179    .
.    9.87    .
.    9.762    .
.    9.268    .
.    10.127    .
.    10.188    .
.    9.383    .
.    9.057    .
.    9.608    .
.    9.18    .
.    9.747    .
.    9.793    .
.    9.835    .
.    9.853    .
.    9.718    .
.    9.155    .
.    9.787    .
.    9.106    .
.    8.834    .
.    8.761    .
.    9.364    .
.    9.823    .
.    9.821    .
.    9.784    .
.    9.735    .
.    9.532    .
.    9.375    .
.    9.1    .
.    9.494    .
.    9.931    .
.    10.004    .
.    9.778    .
.    9.84    .
.    9.654    .
.    9.91    .
.    10.241    .
.    9.805    .
.    10.013    .
.    9.734    .
.    8.962    .
.    9.003    .
.    9.997    .
.    10.198    .
.    9.63    .
.    9.539    .
.    9.205    .
.    10.018    .
.    10.409    .
.    10.355    .
.    9.765    .
.    9.726    .
.    10.115    .
.    10.216    .
.    9.909    .
.    10.293    .
.    10.669    .
.    10.51    .
.    9.593    .
.    9.87    .
.    10.112    .
.    10.085    .
.    10.404    .
.    10.265    .
.    10.175    .
.    9.321    .
.    9.661    .
.    9.619    .
.    9.657    .
.    9.562    .
.    9.507    .
.    9.493    .
.    9.788    .
.    9.75    .
.    9.835    .
.    9.839    .
.    9.375    .
.    9.088    .
.    10.258    13.45
.    9.904    13.76
.    9.723    13.69
.    9.838    13.53
.    9.945    13.81
.    10.061    13.67
.    9.911    13.8
.    10.346    14.28
.    9.936    13.76
.    9.168    13.6
.    9.203    13.81
.    9.734    13.83
.    8.925    13.52
.    9.553    13.75
.    9.435    13.93
.    9.662    14.01
.    10.074    14.16
.    9.463    13.14
.    9.7    13.72
.    9.929    13.89
.    10.069    13.84
.    9.546    13.74
.    9.533    13.68
.    10.378    14.12
.    10.722    14.32
.    10.182    14.09
.    10.314    14.06
.    9.533    13.58
.    9.429    13.8
.    9.47    13.19
.    9.605    13.53
.    10.106    14.01
.    10.154    14.09
.    9.913    13.88
.    9.57    13.81
.    9.39    13.55
.    9.286    13.52
.    9.312    13.25
.    9.253    13.77
.    9.483    13.94
.    9.553    13.85
.    10.287    14.57
.    9.836    14.2
.    9.676    14.08
.    9.498    13.79
.    9.541    13.58
.    9.528    13.56
.    9.743    13.85
.    9.461    13.63
.    9.659    13.9
.    9.407    13.79
.    9.124    13.27
.    9.058    13.03
.    9.371    13.44
.    9.4    13.33
.    10.08    14.31
.    9.738    13.8
.    9.602    13.36
.    9.208    13.39
.    9.379    13.22
.    8.961    13.29
.    9.324    13.42
.    9.117    13.29
.    9.186    13.36
.    9.282    13.68
.    9.391    13.24
.    9.65    13.56
.    9.594    14
.    9.847    14
.    9.73    14.13
.    9.74    13.86
.    10.131    14.26
.    9.459    13.51
.    9.206    13.45
.    9.36    13.55
.    9.863    13.61
.    9.644    13.67
.    9.448    13.47
.    9.721    13.82
.    9.68    13.49
.    9.864    13.81
.    9.343    13.35
.    9.657    13.66
.    9.707    13.46
.    9.572    13.28
.    9.256    12.94
.    9.564    13.57
.    9.453    13.36
.    9.407    13.28
.    9.606    13.45
.    9.694    13.66
.    9.653    13.78
.    9.574    13.5
.    9.225    13.67
.    9.32    13.26
.    9.472    13.47
.    9.062    13.41
.    9.111    13.34
.    8.672    13.25
.    9.769    13.84
.    9.267    13.4
.    9.709    13.6
.    9.663    13.8
.    9.639    13.8
.    9.89    13.73
.    9.705    13.9
.    9.753    14.07
.    9.945    14.02
.    9.682    13.83
.    9.725    14.3
.    9.877    14.2
.    10.023    14.25
.    10.123    14.23
.    9.808    13.82
.    9.67    13.68
.    9.496    13.5
.    9.42    13.72
.    9.724    13.89
.    9.408    13.68
.    10.181    13.71
.    9.902    13.45
.    9.533    13.47
.    9.51    13.51
.    9.762    13.6
.    9.351    13.21
.    9.769    13.54
.    9.396    13.51
.    9.426    13.49
.    9.674    13.74
.    9.395    13.82
.    9.25    13.36
.    9.367    13.56
.    9.01    13.11
.    9.031    13.02
.    9.275    13.57
.    9.28    13.45
.    9.596    13.64
.    9.565    13.53
.    9.41    13.38
.    9.38    13.52
.    9.681    13.57
.    9.468    13.52
.    9.263    13.55
.    10.024    14.06
.    9.342    13.42
.    9.974    13.71
.    9.41    13.48
.    9.169    13.34
.    9.162    13.37
.    9.517    13.47
.    9.494    13.68
.    9.672    13.79
.    9.553    13.54
.    9.453    13.69
.    9.327    13.45
.    8.778    13.41
.    8.209    12.63
.    8.638    13.06
.    9.373    13.46
.    9.294    13.4
.    9.253    13.33
.    9.354    13.36
.    9.794    13.89
.    9.573    13.7
.    9.442    13.77
.    9.681    13.94
.    9.474    13.87
.    9.634    13.59
.    9.494    13.26
.    9.434    13.38
.    9.667    13.59
.    9.468    13.36
.    9.26    13.37
.    9.305    13.57
.    9.42    13.72
.    9.488    13.71
.    9.11    13.57
.    9.56    13.95
.    9.658    13.92
.    9.569    13.88
.    8.964    13.38
.    8.619    13.25
.    9.255    13.56
.    9.501    13.58
.    9.47    13.68
.    9.669    13.75
.    9.214    13.51
.    9.686    13.86
.    9.796    13.89
.    9.732    13.98
.    9.565    13.93
.    9.837    13.98
.    9.751    13.84
.    9.658    13.76
.    9.289    13.52
.    9.278    13.58
.    9.539    13.62
.    9.611    13.6
.    9.689    13.86
.    9.61    13.66
.    9.719    13.94
.    10.016    14.03
.    9.254    13.61
.    10.024    13.86
.    9.414    13.85
.    9.648    13.63
.    9.308    13.78
.    10.035    13.96
.    9.825    13.82
.    9.93    13.88
.    9.838    14.01
.    9.637    13.77
.    9.845    13.88
.    9.667    13.87
.    9.659    14.1
.    9.608    13.85
.    10.215    13.94
.    9.573    13.61
.    8.914    13.33
.    9.329    13.62
.    9.224    13.44
.    9.832    13.8
.    9.628    13.75
.    9.776    13.81
.    9.845    14.06
.    9.43    13.69
.    9.592    13.62
.    9.897    13.88
.    9.689    13.7
.    9.393    13.52
.    9.381    13.73
.    9.808    13.84
.    9.481    13.65
.    9.42    13.53
.    9.695    13.84
.    9.732    13.83
.    10.014    14.1
.    9.669    13.83
.    10.168    14.11
.    9.602    13.93
.    9.643    13.62
.    9.767    14.09
.    9.579    13.95
.    9.63    13.94
.    9.732    13.91
.    10.077    13.95
.    9.746    13.74
.    9.807    13.91
.    9.743    14
.    10.228    14.27
.    9.895    14.02
.    10.118    14.06
.    9.754    13.85
.    9.546    14.04
.    9.734    14.01
.    9.834    13.98
.    9.862    13.89
.    10.066    14.11
.    9.892    13.8
.    9.876    13.85
.    9.948    14.08
.    9.799    13.93
.    9.81    13.98
.    9.699    13.93
.    9.932    13.82
.    9.915    14.1
.    9.747    13.94
.    9.656    13.94
.    9.701    13.82
.    9.783    13.75
.    9.62    13.71
.    9.422    13.6
.    9.583    13.72
.    9.504    13.65
.    9.568    13.62
.    9.309    13.56
.    9.858    13.79
.    10.046    13.99
.    9.88    13.95
.    9.459    13.59
.    9.45    13.68
.    9.463    13.48
.    9.524    13.8
.    9.526    13.69
.    9.628    13.66
.    9.444    13.62
.    9.382    13.67
.    9.197    13.49
.    9.137    13.31
.    9.2    13.38
.    9.077    13.59
.    9.561    13.68
.    9.453    13.44
.    9.54    13.65
.    9.523    13.53
.    9.504    13.62
.    9.46    13.66
.    9.517    13.58
.    9.684    13.88
.    9.496    13.75
.    9.442    13.72
.    8.797    13.16
.    9.496    13.87
.    9.387    13.79
.    9.592    13.78
.    9.847    13.77
.    9.604    13.69
.    9.623    13.87
.    9.688    13.87
.    9.446    13.75
.    9.877    14.04
.    9.919    13.86
.    10.108    13.96
.    9.687    13.79
.    9.499    13.81
.    10.135    14.18
.    10.002    13.9
.    9.856    13.9
.    9.815    13.9
.    9.759    13.92
.    9.629    13.81
.    9.728    13.92
.    9.331    13.6
.    9.74    13.86
.    9.364    13.66
.    9.301    13.46
.    9.454    13.66
.    9.131    13.52
.    9.078    13.46
.    9.334    13.47
.    9.504    13.72
.    9.391    13.57
.    9.573    13.72
.    9.804    13.78
.    9.253    13.5
.    9.429    13.47
.    9.795    13.77
.    9.594    13.8
.    9.254    13.62
.    9.457    13.68
.    9.555    13.8
.    9.601    13.72
.    9.604    13.76
.    9.41    13.73
.    9.61    13.94
.    9.32    13.73
.    9.333    13.53
.    9.608    13.72
.    9.377    13.39
.    9.297    13.54
.    9.263    13.66
.    9.23    13.53
.    9.198    13.56
.    9.678    13.72
.    9.605    13.81
.    10.009    14.02
.    9.826    13.91
.    9.841    13.94
.    9.827    13.96
.    9.28    13.59
.    9.832    13.71
.    9.604    13.79
.    9.605    13.72
.    9.745    13.85
.    9.696    13.9
.    9.8    13.82
.    9.798    13.95
.    9.866    14
.    9.871    14.03
.    9.528    13.81
.    9.213    13.64
.    9.338    13.5
.    9.397    13.52
.    9.102    13.54
.    9.248    13.56
.    9.277    13.64
.    9.532    13.78
.    9.445    13.76
.    9.701    13.97
.    9.689    13.88
.    9.633    13.87
.    9.614    13.95
.    9.855    13.97
.    9.98    13.87
.    9.688    13.95
.    10.015    14.01
.    9.609    13.71
.    9.643    13.77
.    9.69    13.85
.    9.919    13.89
.    9.414    13.66
.    9.388    13.51
.    9.408    13.74
.    9.335    13.61
.    9.463    13.66
.    9.773    13.8
.    9.742    13.69
.    9.424    13.6
.    9.358    13.53
.    9.782    13.83
.    9.158    13.51
.    9.369    13.58
.    9.779    13.92
.    9.803    13.93
.    9.661    13.92
.    9.648    13.81
.    10.121    13.96
.    10.191    14.13
.    10.436    14.25
.    10.14    14.06
.    9.899    14.07
.    9.778    13.94
.    9.918    14.11
.    10.009    14.07
.    9.864    13.93
.    9.879    14
.    9.82    13.99
.    9.985    14.19
.    9.727    13.76
.    9.749    13.88
.    9.893    13.92
.    10.394    14.15
.    9.523    14
.    10.106    14.09
.    9.701    13.92
.    9.765    13.88
.    10.105    14.16
.    9.861    13.94
.    10.042    14.09
.    9.832    14
.    9.996    14.12
.    9.87    13.96
.    10.028    13.99
.    9.946    13.98
.    9.438    13.66
.    9.691    13.81
.    9.69    13.73
.    9.482    13.62
.    9.704    13.77
.    9.717    14
.    9.775    13.96
.    9.667    13.94
.    9.449    13.71
.    9.121    13.3
.    9.427    13.6
.    8.753    13.37
.    8.882    13.32
.    9.227    13.57
.    8.797    13.29
.    9.58    13.66
.    9.896    14.12
.    9.538    13.73
.    9.754    13.77
.    9.378    13.6
.    9.567    13.74
.    8.89    13.2
.    9.114    13.46
.    9.327    13.45
.    9.383    13.58
.    9.373    13.55
.    9.219    13.34
.    9.649    13.68
.    9.479    13.63
.    9.594    13.65
.    9.642    13.91
.    10.024    14.12
.    9.692    13.86
.    9.773    13.78
.    9.966    13.87
.    9.798    13.88
.    9.578    13.95
.    10.009    13.97
.    9.709    13.93
.    9.951    13.91
.    9.788    14.04
.    9.736    14.03
.    9.817    13.9
.    9.682    14.04
.    9.264    13.67
.    9.643    13.79
.    9.829    13.98
.    9.554    13.81
.    9.799    13.86
.    9.65    13.91
.    9.742    13.73
.    9.737    13.74
.    9.822    13.92
.    9.724    13.92
.    9.88    13.92
.    9.595    13.66
.    9.501    13.76
.    9.534    13.64
.    10.15    14.07
.    9.805    13.91
.    9.869    13.97
.    9.866    13.93
.    10.09    13.88
.    10.058    14
.    9.938    14.11
.    9.598    13.64
.    9.797    13.94
.    9.703    14
.    9.422    13.82
.    9.8    13.92
.    9.447    13.74
.    9.588    13.72
.    9.671    13.81
.    9.978    14.04
.    9.781    13.76
.    9.759    13.81
.    9.733    14.07
.    9.692    13.86
.    9.833    13.97
.    9.706    13.85
.    9.874    13.98
.    9.81    13.97
.    10.012    13.87
.    9.521    13.87
.    9.723    13.88
.    9.389    13.63
.    9.625    13.79
.    9.748    13.81
.    9.428    13.6
.    9.445    13.7
.    9.712    13.76
.    9.76    13.97
.    10.024    14.04
.    10.122    14.01
.    9.566    13.82
.    9.494    13.72
.    9.528    13.98
.    9.319    13.61
.    9.601    13.82
.    9.826    13.93
.    9.785    13.95
.    9.736    13.96
.    9.85    13.99
.    9.77    13.92
.    9.703    14
.    9.325    13.72
.    9.479    13.64
.    9.824    13.73
.    9.705    13.84
.    9.824    13.92
.    9.556    13.77
.    9.675    13.71
.    9.528    13.75
.    9.868    13.98
.    9.597    13.77
.    9.436    13.78
.    9.893    14.1
.    10.19    14.23
.    10.328    14.33
.    10.552    14.43
.    10.102    14.29
.    9.766    13.97
.    9.439    13.77
.    9.882    13.87
.    9.789    13.95
.    9.975    14.11
.    9.987    14.12
.    9.956    13.98
.    10.052    14.14
.    9.553    13.85
.    9.795    13.95
.    9.829    13.96
.    9.676    13.77
.    9.593    13.75
.    9.495    13.81
.    9.754    13.83
.    9.884    14.06
.    9.781    13.83
.    9.957    14.02
.    10.155    14.23
.    10.003    14.21
.    9.256    13.86
.    9.959    14.13
.    9.824    14.07
.    9.601    13.98
.    9.868    13.99
.    9.614    13.7
.    9.458    13.54
.    9.888    14.02
.    9.897    14.01
.    9.996    14.03
.    9.959    14
.    10.123    14.16
.    9.937    14.05
.    9.483    13.83
.    8.954    13.65
.    9.5    13.73
.    9.382    13.67
.    9.496    13.73
.    9.607    13.69
.    9.557    13.76
.    9.751    13.74
.    9.717    13.89
.    9.97    14.12
.    9.971    14.07
.    9.346    13.52
.    9.608    13.87
.    9.915    13.95
.    9.917    14.04
.    9.835    13.94
.    9.713    13.98
.    10.04    13.91
.    9.959    13.94
.    9.876    13.94
.    9.79    13.94
.    9.707    13.84
.    10.062    14.22
.    9.913    14.06
.    9.999    13.98
.    9.697    13.84
.    9.772    14.04
.    9.608    13.84
.    9.645    13.92
.    9.985    13.98
.    10.034    14.19
.    10.022    14.19
.    9.845    13.91
.    10.067    14.12
.    9.687    13.94
.    10.072    14.19
.    10.761    14.41
.    9.653    13.81
.    9.516    13.88
.    10.041    14.18
.    9.769    13.94
.    9.808    13.67
.    9.979    13.9
.    9.854    13.9
.    10.029    14.1
.    10.028    14.17
.    9.832    13.99
.    9.892    13.9
.    9.604    13.76
.    9.522    13.84
.    9.629    13.85
.    9.748    13.9
.    9.803    13.95
.    9.814    13.8
.    9.847    13.98
.    9.882    13.88
.    9.903    13.96
.    9.992    13.99
.    9.607    13.72
.    9.183    13.55
.    9.895    13.93
.    9.918    14
.    9.661    13.74
.    9.621    13.77
.    10.161    14.11
.    10.02    13.95
.    10.056    14.11
.    10.026    14.12
.    9.877    14.01
.    9.896    14.06
.    10.021    14.18
.    10.03    14.2
.    9.622    13.79
.    10.52    14.32
.    9.979    14.12
.    9.436    13.64
.    9.557    13.77
.    9.818    13.89
.    9.934    13.85
.    9.884    13.85
.    9.77    13.66
.    10.1    14.04
.    9.582    13.79
.    9.743    13.91
.    9.787    13.8
.    9.292    13.64
.    9.697    13.88
.    9.67    14.02
.    9.832    13.79
.    9.918    13.92
.    10.292    14.27
.    10.068    14.06
.    9.905    14.01
.    10.009    14.05
.    10.032    14.08
.    10.128    14.14
.    9.704    13.93
.    10.044    14.09
.    9.461    13.86
.    9.641    13.88
.    9.883    14.05
.    10.06    14.04
.    9.931    13.99
.    10.174    14.11
.    10.168    14.27
.    10.22    14.27
.    10.118    14.37
.    9.781    14.02
.    10.016    14.08
.    10.084    14.08
.    10.305    14.26
.    10.241    14.32
.    10.113    14.19
.    9.97    13.96
.    10.063    14.09
.    9.988    13.97
.    10.186    14.13
.    10.408    14.37
.    10.209    14.17
.    9.736    13.98
.    10.16    14.03
.    9.981    13.92
.    9.519    13.71
.    9.764    13.9
.    10.103    14.06
.    10.149    13.97
.    9.938    13.94
.    10.007    13.8
.    9.762    13.92
.    9.668    13.94
.    9.8    13.9
.    10.424    14.39
.    9.376    13.87
.    9.869    14.09
.    9.92    14.13
.    9.889    14.03
.    9.802    14.08
.    10.188    14.18
.    10.158    14.25
.    10.019    14.1
.    10.148    14.19
.    9.981    14.09
.    9.882    14.07
.    10.366    14.25
.    9.839    14.09
.    10.095    14.25
.    9.759    13.95
.    10.046    14.04
.    9.967    14.13
.    10.144    14.14
.    10.204    14.28
.    10.051    14.09
.    9.761    13.83
.    10.153    14.12
.    9.851    14.01
.    9.956    14.07
.    10.115    14.28
.    9.709    14.03
.    9.908    14.02
.    9.956    14.03
.    9.894    14.06
.    10.078    14.09
.    9.864    14.06
.    10.057    14.09
.    10.115    14.07
.    10.007    14.03
.    10.045    14.17
.    9.896    14.09
.    9.648    13.84
.    10.109    14.08
.    9.72    13.88
.    9.995    14.04
.    9.925    13.96
.    9.701    13.92
.    9.995    14.13
.    9.674    13.84
.    9.959    14.07
.    10.348    14.31
.    9.866    14.1
.    10.194    14.22
.    10.505    14.35
.    10.147    14.16
.    9.979    14.14
.    9.83    13.94
.    9.837    13.95
.    9.994    13.9
.    10.12    14.05
.    10.184    14.07
.    10.416    14.2
.    10.307    14.22
.    9.972    14.06
.    9.688    13.89
.    9.795    14.09
.    9.621    13.93
.    9.857    14.03
.    10.151    14.14
.    9.638    13.86
.    9.859    13.75
.    9.793    13.83
.    10.452    14.23
.    10.135    14.09
.    10.069    14.09
.    9.881    14
.    9.425    13.71
.    10.029    14.06
.    10.143    14.02
.    9.711    13.9
.    10.076    14.1
.    9.945    14.07
.    9.715    13.9
.    10.013    14.08
.    9.963    13.95
.    10.082    14.04
.    9.803    13.85
.    9.989    13.97
.    9.543    13.75
.    9.845    13.96
.    9.822    14.1
.    9.994    14.23
.    10.047    14.17
.    9.984    14.04
.    10.147    14
.    9.946    14.06
.    10.043    14.07
.    9.962    14.03
.    10.186    14.28
.    10.045    14.25
.    9.723    13.92
.    10.256    14.23
.    9.873    13.89
.    9.531    13.77
.    9.875    13.93
.    10.056    14.06
.    10.127    14.12
.    9.892    13.87
.    9.889    13.85
.    9.998    13.99
.    10.066    14.02
.    9.885    13.87
.    9.838    13.89
.    10.212    14.08
.    9.584    13.71
.    9.617    13.96
.    9.763    13.79
.    9.826    13.87
.    9.774    13.81
.    9.888    13.94
.    9.94    13.88
.    9.907    13.97
.    10.027    14.04
.    9.984    13.95
.    9.706    13.9
.    9.392    13.75
.    9.624    13.7
.    9.743    13.93
.    9.703    13.79
.    9.813    13.8
.    9.902    13.9
.    9.8    13.91
.    9.674    13.68
.    9.814    13.84
.    9.691    13.82
.    9.33    13.66
.    9.675    13.86
.    9.463    13.67
.    9.266    13.57
.    9.672    13.88
.    9.818    13.98
.    9.932    14.1
.    9.925    13.91
.    9.94    13.93
.    10.039    14.13
.    10.11    14.03
.    10.047    14.08
.    9.982    14.05
.    10.248    14.18
.    10.299    14.13
.    10.083    14.15
.    9.754    13.87
.    9.937    13.96
.    9.955    13.91
.    9.918    13.93
.    10.078    14.08
.    10.068    14.13
.    10.088    14.15
.    9.905    14
.    9.565    13.82
.    9.884    13.98
.    10.077    14.14
.    10.109    14.09
.    10.183    14.18
.    10.123    14.12
.    9.984    14.01
.    10.098    14.06
.    9.85    14
.    10.075    14.13
.    10.066    14.1
.    10.088    14.07
.    9.884    14.05
.    10.163    14.19
.    9.543    13.74
.    9.798    13.94
.    9.728    13.9
.    9.717    13.88
.    9.678    13.82
.    9.809    13.81
.    9.762    13.76
.    9.897    13.84
.    9.939    13.83
.    9.991    14.03
.    10.122    14.21
.    9.767    13.92
.    10.243    14.23
.    9.912    13.93
.    9.352    13.66
.    9.622    13.85
.    9.738    13.87
.    9.859    13.86
.    9.701    13.88
.    9.976    14.15
.    9.957    13.89
.    9.938    14.04
.    9.666    13.77
.    9.568    13.74
.    9.568    13.78
.    9.331    13.71
.    9.459    13.68
.    9.512    13.75
.    9.517    13.74
.    9.692    13.91
.    9.741    13.9
.    9.605    13.68
.    9.583    13.9
.    9.718    13.78
.    9.507    13.8
.    9.568    13.84
.    9.715    13.88
.    9.766    13.97
.    9.68    13.97
.    9.87    14.1
.    9.927    14.07
.    10.064    14.19
.    9.761    13.96
.    10.008    14.16
.    10.032    14.12
.    9.983    13.99
.    10.215    14.18
.    10.368    14.19
.    10.424    14.4
.    10.256    14.22
.    10.08    14.19
.    10.126    14.11
.    10.085    14.08
.    9.963    13.85
.    10.085    14.1
.    9.899    13.92
.    9.845    13.94
.    10.035    14.1
.    10.142    14.1
.    10.001    13.95
.    10.071    14.06
.    10.045    14.06
.    10.154    14.23
.    10.027    14.15
.    9.814    13.99
.    10.083    14.09
.    10.062    14.07
.    9.976    14
.    10.001    13.86
.    9.845    13.93
.    9.821    13.91
.    9.968    14.1
.    10.057    14.07
.    10.317    14.28
.    9.39    13.57
.    9.676    13.79
.    9.551    13.87
.    9.974    13.95
.    9.828    13.87
.    9.949    14.09
.    10.01    14.05
.    10.058    14.1
.    9.802    13.92
.    10.186    14.19
.    10.068    14.08
.    10.207    14.22
.    10.057    14.04
.    10.09    14.1
.    10.097    14.23
.    10.067    14.1
.    9.961    13.99
.    10.001    14.04
.    10.095    14.14
.    10.039    14.14
.    9.99    14.09
.    9.842    13.86
.    10.077    14.08
.    10.163    14.2
.    9.964    14.1
.    9.928    14.11
.    9.837    13.87
.    9.932    14.11
.    9.877    13.87
.    9.929    13.92
.    10.028    14.08
.    10.125    14.04
.    10.1    14.11
.    10.046    14.01
.    9.848    14
.    10.314    14.25
.    9.769    13.8
.    9.786    13.8
.    9.848    13.87
.    9.897    14
.    10.09    14.1
.    10.045    14.16
.    10.144    14.24
.    10.304    14.06
.    10.105    14.03
.    9.915    13.98
.    10.001    13.95
.    9.707    13.84
.    9.618    13.71
.    9.721    13.62
.    9.79    13.62
.    9.703    13.85
.    9.782    13.92
.    9.729    13.7
.    9.733    13.62
.    9.574    13.66
.    9.72    13.79
.    9.471    13.68
.    9.887    13.95
.    9.689    13.81
.    9.707    13.9
.    9.537    13.7
.    9.798    13.88
.    9.967    13.86
.    9.69    13.73
.    9.818    13.91
.    9.847    13.77
.    9.945    13.89
.    9.81    13.87
.    9.893    13.79
.    9.852    13.75
.    9.942    13.93
.    9.969    14.13
.    9.803    13.88
.    9.811    13.92
.    10.098    13.94
.    9.998    14.11
.    9.904    13.89
.    9.926    13.96
.    9.821    13.75
.    9.855    13.95
.    9.561    13.93
.    9.806    13.93
.    9.597    13.73
.    9.78    14.02
.    9.926    13.92
.    10.102    14.12
.    9.722    13.86
.    9.828    14.07
.    9.91    14
.    9.866    14.03
.    10.276    14.16
.    9.921    14.04
.    9.837    13.99
.    9.669    13.81
.    9.689    13.83
.    10.195    14.33
.    9.717    13.92
.    9.598    13.86
.    9.71    13.93
.    9.849    13.88
.    9.817    13.94
.    9.776    13.77
.    9.887    14.08
.    9.876    13.94
.    9.685    13.8
.    9.546    13.85
.    9.536    13.74
.    9.858    13.88
.    10.068    14.13
.    10.035    14.12
.    9.958    14.09
.    10.056    13.94
.    10.001    13.92
.    9.913    14
.    9.964    14.08
.    10.155    14.09
.    10.201    14.34
.    9.981    14.06
.    10.336    14.3
.    9.819    14.07
.    10.086    14.06
.    9.918    13.94
.    10.144    14.09
.    10.001    14.07
.    9.954    13.89
.    10.016    14.14
.    9.908    14.08
.    9.883    14.03
.    9.703    13.88
.    9.896    13.97
.    9.596    13.85
.    9.611    13.78
.    9.647    13.95
.    9.658    13.91
.    9.655    13.8
.    9.87    13.93
.    9.914    14.04
.    10.017    14.02
.    9.934    13.98
.    9.889    13.87
.    9.861    13.94
.    9.329    13.68
.    9.383    13.67
.    9.688    13.89
.    9.803    13.91
.    9.867    13.96
.    10.045    14.05
.    9.957    14.03
.    9.959    14.17
.    9.897    14.03
.    9.886    13.93
.    9.892    13.91
.    10.208    14.12
.    10.3    14.22
.    10.456    14.33
.    10.408    14.36
.    10.249    14.34
.    10.19    14.31
.    10.306    14.26
.    10.106    14.2
.    10.025    14.05
.    9.965    14.1
.    10.05    14.17
.    9.86    14.06
.    9.913    13.96
.    9.525    13.93
.    9.396    13.77
.    9.743    13.96
.    9.769    13.9
.    9.748    14
.    9.787    14
.    9.931    13.98
.    9.823    14.11
.    9.669    13.86
.    9.652    13.91
.    9.695    13.89
.    9.625    13.94
.    9.997    14.07
.    9.902    14.07
.    9.939    14.15
.    9.953    14.02
.    10.011    14.27
.    9.987    14.02
.    9.966    14.05
.    9.788    13.74
.    10.036    13.96
.    9.729    13.94
.    9.681    13.9
.    9.636    13.76
.    9.773    13.94
.    9.606    13.9
.    9.262    13.65
.    9.701    13.75
.    9.648    13.72
.    9.693    13.79
.    9.718    13.81
.    9.684    13.76
.    9.719    13.85
.    9.37    13.63
.    9.58    13.79
.    9.805    13.95
.    9.627    14.11
.    10.161    14.21
.    10.282    14.27
.    10.202    14.3
.    10.037    14.3
.    10.173    14.3
.    10.067    14.26
.    10.025    14.27
.    9.991    13.89
.    10.069    13.96
.    10.311    14.21
.    10.021    14.11
.    10.025    14.11
.    10.068    14.17
.    10.139    14.24
.    10.009    14.24
.    9.951    14.2
.    9.823    14.03
.    9.974    14.2
.    9.716    13.85
.    9.904    14.15
.    9.862    14.02
.    10.034    14.1
9.72    9.892    14.04
9.75    9.943    14.13
9.75    9.807    13.88
9.75    10.089    14.22
9.73    9.83    14.21
9.76    9.888    14.01
9.8    10.08    14.18
9.88    10.006    14.01
9.77    10.118    14.11
9.9    10.066    14.24
10.01    10.192    14.23
9.89    10.14    14.26
10.03    10.599    14.49
9.93    10.097    14.31
10.02    10.195    14.38
9.9    10.021    14.38
10.04    10.22    14.43
10.09    10.315    14.45
10.05    10.118    14.19
10    10.044    14.32
10.05    10.088    14.29
10.07    10.108    14.3
9.97    10.08    14.17
10    10.297    14.26
9.85    10.131    14.15
10.02    10.55    14.59
10.07    10.493    14.53
10.03    10.449    14.63
9.89    10.388    14.45
9.97    10.147    14.28
9.92    10.093    14.34
9.97    10.117    14.48
9.98    10.23    14.49
9.9    10.101    14.22
9.87    9.94    14.05
9.87    10.033    14.2
10.04    10.444    14.44
9.8    9.949    14.08
9.8    10.05    14.22
9.65    9.828    13.94
9.76    10    13.98
9.77    10.07    14.21
9.84    9.885    13.98
9.67    9.996    14.14
9.74    10.063    14.04
9.73    10.069    14.04
9.65    9.879    13.98
9.76    10.002    14.09
9.91    10.41    14.36
10.02    10.64    14.47
9.88    10.543    14.36
10.14    10.36    14.41
10.06    10.208    14.36
10.11    10.242    14.45
9.87    10.123    14.22
10.05    10.204    14.16
9.99    10.296    14.44
9.96    10.369    14.5
9.79    10.204    14.27
9.85    10.473    14.29
9.61    10.131    14.13
9.58    10.015    14.2
9.71    9.925    14.1
9.77    9.991    14.34
9.68    9.907    14.12
9.88    10.19    14.43
9.72    9.969    13.99
9.67    9.949    14.14
9.72    10.013    14.21
9.39    9.874    14.22
9.69    10.076    14.14
9.48    9.841    14.05
9.53    9.558    13.98
9.74    10.07    14.26
9.68    9.773    14
9.77    9.985    14.21
9.69    10.02    14.12
9.67    10.073    14.15
9.71    9.946    14.22
9.55    9.915    14
9.73    10.156    14.31
9.68    9.984    14.18
9.59    10.01    14.1
9.74    9.937    14.02
9.78    10.03    14.1
9.87    10.389    14.33
9.74    10.183    14.46
9.74    10.212    14.34
9.87    10.296    14.36
9.88    10.21    14.27
9.8    10.162    14.11
9.74    9.961    14.1
9.7    10.029    14.16
9.62    9.992    14
9.62    9.966    14.07
9.78    9.911    14.07
9.79    10.096    14.15
10.05    10.259    14.38
10.09    10.642    14.5
9.82    9.955    14.2
10.01    10.301    14.35
9.87    10.146    14.24
10.1    10.262    14.49
10.05    10.387    14.54
9.96    10.161    14.2
9.93    10.267    14.28
10.09    10.25    14.3
9.98    10.249    14.28
10.28    10.566    14.52
10.17    10.595    14.62
9.93    10.254    14.43
10.12    10.438    14.54
9.95    10.34    14.42
10.02    10.317    14.56
10.04    10.39    14.54
10.12    10.291    14.36
10.09    10.335    14.43
10.16    10.391    14.46
10.01    10.383    14.45
9.81    10.102    14.08
9.79    10.339    14.28
9.57    10.142    14.2
9.72    10.322    14.38
9.7    10.256    14.4
9.79    10.21    14.32
9.7    10.117    14.17
9.73    10.119    14.14
9.83    10.158    14.41
9.87    10.193    14.4
9.94    10.21    14.47
9.93    10.304    14.36
9.82    10.087    14.15
9.99    10.27    14.31
9.9    10.429    14.45
9.76    10.414    14.37
10    10.91    14.8
9.9    10.452    14.55
10.02    10.339    14.47
10.04    10.347    14.4
9.97    10.27    14.58
9.95    10.312    14.38
9.88    10.158    14.21
10.02    10.478    14.49
10.21    10.668    14.49
10.14    10.399    14.41
10.03    10.371    14.44
10.07    10.523    14.54
10.18    10.33    14.42
10.03    10.508    14.64
10.09    10.316    14.36
10.28    10.554    14.62
10.12    10.394    14.56
10.14    10.32    14.51
9.95    10.355    14.53
9.84    10.267    14.27
9.79    10.106    14.27
9.8    10.073    14.28
9.87    10.559    14.46
9.78    10.416    14.35
9.89    10.447    14.49
9.72    10.15    14.2
9.74    10.167    14.22
9.75    10.165    14.25
9.6    9.788    14
9.54    9.961    14.07
9.53    9.814    13.89
9.75    9.831    14
9.74    9.752    13.97
9.7    10.056    14.18
9.66    10.362    14.33
9.71    10.219    14.34
9.55    10.295    14.31
9.68    10.266    14.25
9.74    10.286    14.24
9.89    10.198    14.18
9.89    10.203    14.28
9.75    10.157    14.13
9.59    9.904    14.04
9.82    10.128    14.16
9.83    9.811    14.04
10    10.271    14.16
9.91    10.371    14.36
9.77    9.834    14.04
9.77    10.259    14.24
9.85    10.386    14.48
9.86    10.436    14.24
10    10.419    14.52
9.97    10.267    14.33
9.87    10.262    14.19
9.94    10.36    14.4
9.76    10.369    14.32
10.03    10.521    14.41
10.04    10.517    14.41
10.03    10.682    14.53
9.99    10.81    14.85
9.88    10.36    14.47
10.11    10.344    14.55
9.99    10.293    14.33
10.07    10.475    14.45
10.03    10.435    14.59
10.21    10.564    14.53
10.11    10.446    14.37
9.99    10.536    14.63
10.04    10.437    14.49
9.8    10.231    14.36
9.82    10.145    14.38
10.01    10.319    14.56
9.98    10.09    14.41
9.85    9.997    14.39
9.81    10.265    14.35
9.82    10.272    14.29
9.91    10.261    14.51
10.05    10.301    14.72
10.03    10.157    14.35
9.97    10.187    14.28
10.01    10.112    14.45
9.93    10.342    14.48
9.83    10.346    14.43
9.91    10.526    14.5
9.87    10.3    14.6
9.69    10.201    14.34
9.87    10.277    14.38
9.95    10.53    14.61
10.01    10.453    14.27
10.02    10.457    14.37
9.97    10.609    14.41
10.02    10.528    14.5
10.09    10.56    14.63
10.23    10.773    14.63
10.48    10.726    14.67
10.66    11.369    15.11
10.43    10.779    14.79
10.66    10.956    14.76
10.57    10.873    14.88
10.52    10.85    14.95
10.44    10.928    14.83
10.45    10.896    14.8
10.34    10.661    14.52
10.3    10.706    14.62
10.08    10.37    14.61
10.2    10.729    14.74
10.03    10.697    14.62
10.14    10.946    14.78
9.9    10.353    14.41
9.97    10.354    14.37
9.94    10.394    14.43
9.83    10.44    14.57
9.93    10.45    14.54
9.88    10.429    14.49
10.01    10.56    14.51
9.95    10.438    14.52
9.92    10.268    14.47
9.92    10.54    14.46
9.67    10.218    14.25
9.91    10.659    14.72
9.94    10.529    14.62
10.03    10.705    14.68
10.05    10.327    14.49
10    10.345    14.61
9.91    10.288    14.56
9.88    10.427    14.53
10    10.389    14.48
9.98    10.225    14.27
10    10.048    14.42
9.96    10.177    14.42
9.96    10.611    14.52
10.06    10.499    14.56
10.04    10.66    14.6
10.18    10.523    14.51
10.19    10.524    14.68
9.99    10.469    14.55
10.06    10.523    14.61
10.25    10.738    14.61
10.09    10.515    14.63
10.19    10.53    14.55
10.17    10.672    14.79
10.17    10.365    14.67
10.21    10.934    14.87
10.19    11    14.83
10.25    10.898    15.07
10.24    10.686    14.75
10.31    10.595    14.75
10.32    10.513    14.55
10.22    10.685    14.77
10.17    10.57    14.59
10.24    10.656    14.68
10.13    10.461    14.57
10.25    10.616    14.61
10.15    10.349    14.47
10.29    10.811    14.86
10.17    10.573    14.68
10.08    10.448    14.66
10.15    10.532    14.62
10.22    10.749    14.8
10.03    10.774    14.57
10.11    10.527    14.57
10.11    10.588    14.8
10.19    10.644    14.74
10.3    10.663    14.79
10.22    10.606    14.6
10.37    10.843    14.8
10.19    10.73    14.66
10.21    10.906    14.86
10.32    10.822    14.8
10.15    10.714    14.68
10.05    10.228    14.45
10.02    10.539    14.58
9.81    10.329    14.21
9.94    10.419    14.57
10.12    10.476    14.56
10.24    10.8    14.76
10.15    10.839    14.81
10.09    10.53    14.61
10.3    10.767    14.91
10.18    10.563    14.81
10.21    10.783    14.87
10.36    10.882    14.79
10.21    10.673    14.69
10.22    10.825    14.73
10.27    10.801    14.7
10.13    10.643    14.68
10.32    10.761    14.86
10.35    10.854    14.91
10.31    10.865    14.82
10.18    10.542    14.84
10.19    10.657    14.69
10.22    10.74    14.86
10.19    10.638    14.75
10.13    10.551    14.63
9.97    10.269    14.47
10.12    10.655    14.69
10.15    10.797    14.59
10.19    10.638    14.74
10.24    10.629    14.65
10.31    10.891    14.76
10.19    10.712    14.81
10.21    10.851    14.87
10.41    11.104    15.16
10.22    10.867    14.86
10.24    10.749    14.85
10.18    10.851    14.87
10.19    10.623    14.86
10.16    10.412    14.58
10.18    10.403    14.71
10.17    10.555    14.8
10.14    10.78    14.75
10.15    10.792    14.75
10.06    10.505    14.62
9.94    10.49    14.63
9.7    10.234    14.41
9.75    10.337    14.37
9.93    10.907    14.82
9.98    10.377    14.58
9.81    10.297    14.52
9.88    10.428    14.42
9.99    10.69    14.66
9.94    10.45    14.36
10.14    10.509    14.72
10.13    10.8    14.67
10.17    10.773    14.76
10.1    10.532    14.68
10.15    10.576    14.7
10.16    10.731    14.7
10.09    10.615    14.61
10.05    10.694    14.68
10.06    10.536    14.74
9.99    10.566    14.75
10.36    10.514    14.78
10.17    10.896    14.8
10.38    10.801    14.86
10.25    10.608    14.68
10.39    10.73    14.78
10.21    10.439    14.72
10.55    10.692    14.91
10.52    10.731    15
10.55    10.946    15.09
10.4    10.806    14.94
10.46    10.746    14.93
10.39    10.84    14.82
10.42    10.747    14.67
10.44    10.668    14.76
10.48    10.68    14.79
10.31    10.613    14.89
10.27    10.821    14.96
10.18    10.263    14.63
10    10.24    14.66
9.99    10.254    14.61
9.9    10.315    14.7
10.12    10.692    14.76
10.14    10.383    14.59
10.32    10.669    14.67
10.38    10.68    14.84
10.33    10.692    14.89
10.29    10.736    14.75
10.12    10.684    14.76
10.12    10.577    14.73
10.13    10.529    14.74
9.91    10.428    14.57
;
run;

 proc sgplot data=Temp;
       title "Scatter Plots of the Temperature Series";
       scatter x=date y=cru;
       scatter x=date y=uah ;
       scatter x=date y=giss;
  run;

 title;

data append(keep=date cru giss uah);
    do i=1 to 1200;
        cru = .; giss=.; uah=.;
        date = intnx('month','01jan2012'D, i);
    format date monyy7.;
    output;
    end;
run;

 proc append base=Temp data=Append; run;

data temp;
   set temp;
   tindex = _n_;
run;

 ods output RegressionEstimates=regEst;
 proc ssm data=Temp;
     id date interval=month;
     parms loga1 logr1 logr3  logsigma;
     parms b=0 c=0;
     parms rhoParm;
     rho = (exp(rhoParm)-1)/(exp(rhoParm)+1);
     sigmaSq = exp(2*logsigma);
     initSigmaSq = sigmaSq/(1-rho*rho);
     a1 = exp(loga1);
     a1Sq = a1*a1;
     r1sq  = exp(2*logr1);
     r3sq  = exp(2*logr3);
     giss_var = a1Sq*r1sq*sigmaSq;
     nobs=1945;
     if tindex <= nobs then
         cru_var = a1Sq*exp(2*b + 2*c*tindex/nobs)*sigmaSq;
     else cru_var = a1Sq*exp(2*b + 2*c)*sigmaSq;
     uah_var = a1Sq*r3sq*sigmaSq;
     UAH_Intercept=1.0;
     CRU_Intercept=1.0;
     trend level(ll) variance=0 slopevar=sigmaSq;
     state auto(1) T(g)=(rho) cov(g)=(sigmaSq)  cov1(g)=(initSigmaSq);
     comp auto_common = auto*(a1);
     state wn(3) type=wn cov(d)=(giss_var cru_var uah_var);
     comp wn_giss = wn[1];
     comp wn_cru = wn[2];
     comp wn_uah = wn[3];
     model GISS = level auto_common  wn_giss;
     model CRU = CRU_Intercept level auto_common  wn_cru;
     model UAH = UAH_Intercept level auto_common  wn_uah;
     comp slope = level_state_*(0 1);
     output out=For pdv press;
 run;

data _NULL_;
   set regEst;
   if _n_ = 1 then call symput('intercept1',trim(left(estimate)));
   else call symput('intercept2',trim(left(estimate)));
run;
data for;
    set For;
    cru_adj = cru - &intercept1;
    uah_adj = uah - &intercept2;
run;

 proc sgplot data=For;
    where date < '01feb2021'd;
    title "Fitted Trend for the Temp Series (Up to Year 2020) ";
    title2 "(CRU and UAH Adjusted by Their Estimated Intercepts) ";
    scatter x=date y=cru_adj / LEGENDLABEL="Adjusted CRU"
        MARKERATTRS=GraphData1(symbol=star size=3);
    scatter x=date y=uah_adj / LEGENDLABEL="Adjusted UAH"
        MARKERATTRS=GraphData2(symbol=plus size=3);
    scatter x=date y=giss /
        MARKERATTRS=GraphData3(symbol=triangle size=3);
    series x=date y=smoothed_level / MARKERATTRS=GraphData4;
 run;

 proc sgplot data=For;
    title "Temperature Projections for the Next 100 Years";
    band x=date lower=smoothed_lower_level upper=smoothed_upper_level;
    series x=date y=smoothed_level;
    refline '01feb2012'd / axis=x lineattrs=(pattern=shortdash)
       LEGENDLABEL= "Start of Multistep Forecasts"
       name="Forecast Reference Line";
 run;

 proc sgplot data=For;
    where date <=  '01feb2031'd;
    title "The Monthly Rate of Temperature Change  (Up to Year 2030)";
    band x=date lower=smoothed_lower_slope upper=smoothed_upper_slope;
    series x=date y=smoothed_slope;
    refline '01feb2012'd / axis=x lineattrs=(pattern=shortdash)
       LEGENDLABEL= "Start of Multistep Forecasts"
       name="Forecast Reference Line";
    refline 0 / axis=y lineattrs=(pattern=dash);
 run;

 title;