Resources

Variable Star Data

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

                    SAS Sample Library

        Name: ucmex02.sas
 Description: Example program from SAS/ETS User's Guide,
              The UCM Procedure
       Title: Variable Star Data
     Product: SAS/ETS Software
        Keys: equally spaced univariate time series data
        PROC: UCM
       Notes:

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

ods graphics on;

data star;
   input magnitude @@;
   day = _n_;
datalines;
  25  28  31  32  33  33  32  31  28   25   22   18
  14  10   7   4   2   0   0   0   2    4    8   11
  15  19  23  26  29  32  33  34  33   32   30   27
  24  20  17  13  10   7   5   3   3    3    4    5
   7  10  13  16  19  22  24  26  27   28   29   28
  27  25  24  21  19  17  15  13  12   11   11   10
  10  11  12  12  13  14  15  16  17   18   19   19
  19  19  20  20  20  20  20  20  20   20   21   20
  20  20  20  19  18  17  16  15  13   12   11   10
   9   9  10  10  11  12  14  16  19   21   24   25
  27  28  29  29  28  27  25  23  20   17   14   11
   8   5   4   2   2   2   4   6   9   12   16   19
  23  27  30  32  33  34  33  32  30   27   24   20
  16  12   9   5   3   1   0   0   1    3    6    9
  13  17  21  24  27  30  32  33  33   32   31   28
  25  22  19  15  12   9   7   5   4    4    5    5
   7   9  12  14  17  20  22  24  25   26   27   27
  26  25  24  22  20  18  17  15  14   13   13   12
  12  12  13  13  13  14  14  15  15   16   17   17
  17  17  18  18  19  19  20  20  21   21   22   22
  22  22  22  21  20  19  17  16  14   12   11    9
   8   7   8   8   9  10  12  14  17   20   23   25
  27  29  30  30  30  29  27  25  22   19   16   12
   9   6   4   2   1   1   2   4   7   10   14   17
  21  25  29  31  33  34  34  33  31   29   26   22
  19  14  11   7   4   2   1   0   1    2    5    7
  11  15  19  22  25  28  30  32  32   32   31   29
  26  23  21  17  14  11   9   7   6    5    6    6
   7   9  11  13  15  18  20  22  23   24   25   25
  25  24  24  22  21  19  18  17  16   15   15   14
  14  14  14  14  14  14  14  14  14   14   15   15
  15  15  16  16  17  18  19  20  21   22   23   23
  24  24  24  23  22  21  19  17  15   13   11    9
   7   6   6   6   7   8  10  12  15   18   22   24
  27  29  31  31  31  31  29  27  24   21   18   14
  10   7   5   2   1   0   1   2   5    8   12   15
  19  23  27  30  32  34  34  34  32   30   28   24
  20  16  13   9   6   3   2   1   1    2    4    6
   9  13  17  20  23  26  28  30  31   31   31   29
  27  24  22  19  16  13  11   9   8    7    7    7
   8   9  11  12  14  16  18  20  21   22   23   23
  23  23  23  22  21  20  19  18  18   17   17   16
  16  16  16  15  15  15  14  14  13   13   13   13
  13  13  14  14  15  16  18  19  21   22   24   24
  25  26  26  25  24  23  21  19  16   14   12    9
   7   5   5   4   5   6   8  10  13   16   20   23
  26  29  31  32  32  32  31  29  26   23   20   16
  12   8   6   3   1   0   0   1   3    6   10   13
  17  21  25  28  31  33  34  34  33   31   29   26
  22  18  15  11   8   5   3   2   2    2    4    5
;

proc timeseries data=star plot=series;
   var magnitude;
run;

ods output CycleDescription = cycledesc;
proc ucm data=star;
   model magnitude;
   irregular;
   level var=0 noest;
   cycle;
   cycle;
   estimate;
run;

proc print data=cycledesc noobs;
   var name type period rho  ErrorVar;
run;