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;