如何使用proc模板在SAS中生成时间序列上的直方图图表

时间:2014-10-07 03:34:56

标签: plot sas proc

我有一个包含两个时间序列数据的数据集。一个是y,另一个是y。我想在一个图中为这两个系列绘制一个过度的直方图图形。

data a;
input y haty;
datalines;
2    2.1
2.12 2.24
2.3  2.5
3.1  3.1
1.23 0.98
;
run;

我有一个示例代码,用于定义覆盖两个直方图的模板。

proc template;
define statgraph dualhist;
   begingraph;
   entrytitle "Graph"; /** optional title **/
   layout overlay / xaxisopts=(label="Length");
      /** first plot: a histogram **/
      histogram y / name="Y"
          binwidth=5; 
      /** second plot: a semi-transparent histogram **/
      histogram haty / name="Predicted Y"
          binwidth=5 datatransparency=0.7
          fillattrs=(color=GraphData2:color);
      /** optional: add legend by specifying names **/
      discretelegend "y" "haty";
   endlayout;
   endgraph;
end;
run;

proc sgrender data=a template=dualhist;
run;

我需要一些关于statgraph的帮助,尤其是x轴。我希望y轴是y和haty的值。

1 个答案:

答案 0 :(得分:0)

目前尚不清楚你想要什么,但这里有一些可能有用的例子。

/* SIMPLE HISTOGRAM */
/* Overlayed histogram with counts. 
    Custom intervals can be set with the binstart= and 
    binwidth= options  */
proc sgplot data = a;
    histogram y/ scale = count;
    histogram haty / transparency = 0.7 scale = count;
run;

/* VIRTICAL HISTOGRAM */
/* Use proc univariate to bin data */
proc univariate data = a noprint;
    histogram haty y / noplot outhistogram = b;
run;
/* Get the counts for each variable in a separate column */
data c;
    merge 
        b (where = (_var_ = "y")) 
        b (where = (_var_ = "haty") rename = (_count_ = _count_hat));
    by _midpt_;
    drop _var_ _obspct_;
run;
/* Plot the histogram on its with horizontal counts */
proc sgplot data = c;
    hbar _midpt_ / 
        response = _count_ 
        legendlabel = "Y" 
        barwidth = 1;
    hbar _midpt_ / 
        response = _count_hat 
        legendlabel = "Predicted Y"  
        barwidth = 1
        transparency = 0.7;
    xaxis  integer;
run;

/* BAR CHART OF DATA */
/* This should be a scatter or series plot, not bars! */
/* Add an id variable */
data d;
    set a;
    ID = _N_;
run;
proc sgplot data = d;
    vbar ID / response = y;
    vbar ID / response = haty transparency = 0.7;
run;