SAS GRAPH:对齐2个图的轴与格子行

时间:2015-10-13 07:27:14

标签: graph plot sas

我真的在努力解决这个问题,我真的找不到解决方案,我希望在这里张贴,有人会帮忙。

问题是,使用这个2格子行模板,我无法对齐轴。

BIG PROBLEM

现在,我详细介绍了我用代码重现每一步的所有内容:

我有这个情节First_plot

这里是图表和样本数据集的代码(它是这个的简单版本:http://support.sas.com/kb/39/092.html)。

示例数据集:

data immune;
   format sival f6.3;
   input trt $ 1-8 cyc $ 10-18 pt lbparm $ xval sival;
   cards;
Drug A   Cycle 1   1   C3  1   1.120
Drug A   Cycle 1   1   C4  1   0.147
Drug A   Cycle 1   1   C3  2   1.080
Drug A   Cycle 1   1   C4  2   0.131
Drug A   Cycle 1   1   C3  3   0.887
Drug A   Cycle 1   1   C4  3   0.113
Drug A   Cycle 2   1   C3  4   1.440
Drug A   Cycle 2   1   C4  4   0.278
Drug A   Cycle 2   1   C3  5   1.180
Drug A   Cycle 2   1   C4  5   0.234
Drug A   Cycle 2   1   C3  6   1.360
Drug A   Cycle 2   1   C4  6   0.281
Drug A   Cycle 3   1   C3  7   1.190
Drug A   Cycle 3   1   C4  7   0.282
Drug A   Cycle 3   1   C3  8   1.000
Drug A   Cycle 3   1   C4  8   0.228
Drug A   Cycle 3   1   C3  9   1.040
Drug A   Cycle 3   1   C4  9   0.228
Drug A   Cycle 4   1   C3 10   0.917
Drug A   Cycle 4   1   C4 10   0.163
Drug A   Cycle 4   1   C3 11   0.789
Drug A   Cycle 4   1   C4 11   0.136
Drug A   Cycle 4   1   C3 12   0.861
Drug A   Cycle 4   1   C4 12   0.148
Drug A   Cycle 1   2   C3  1   1.180
Drug A   Cycle 1   2   C4  1   0.264
Drug A   Cycle 1   2   C3  2   0.942
Drug A   Cycle 1   2   C4  2   0.184
Drug A   Cycle 1   2   C3  3   1.010
Drug A   Cycle 1   2   C4  3   0.160
Drug A   Cycle 2   2   C3  4   1.050
Drug A   Cycle 2   2   C4  4   0.134
Drug A   Cycle 2   2   C3  5   0.980
Drug A   Cycle 2   2   C4  5   0.120
Drug A   Cycle 2   2   C3  6   1.020
Drug A   Cycle 2   2   C4  6   0.126
Drug A   Cycle 3   2   C3  7   0.961
Drug A   Cycle 3   2   C4  7   0.110
Drug A   Cycle 3   2   C3  8   0.859
Drug A   Cycle 3   2   C4  8   0.091
Drug A   Cycle 3   2   C3  9   0.928
Drug A   Cycle 3   2   C4  9   0.097
Drug A   Cycle 4   2   C3 10   1.380
Drug A   Cycle 4   2   C4 10   0.330
Drug A   Cycle 4   2   C3 11   1.210
Drug A   Cycle 4   2   C4 11   0.281
Drug A   Cycle 4   2   C3 12   1.180
Drug A   Cycle 4   2   C4 12   0.278
Drug A   Cycle 1   3   C3  1   1.180
Drug A   Cycle 1   3   C4  1   0.269
Drug A   Cycle 1   3   C3  2   1.010
Drug A   Cycle 1   3   C4  2   0.213
Drug A   Cycle 1   3   C3  3   1.040
Drug A   Cycle 1   3   C4  3   0.200
Drug A   Cycle 2   3   C3  4   1.200
Drug A   Cycle 2   3   C4  4   0.332
Drug A   Cycle 2   3   C4  5   0.371
Drug A   Cycle 2   3   C4  6   0.316
Drug A   Cycle 3   3   C4  7   0.271
Drug A   Cycle 3   3   C3  8   1.050
Drug A   Cycle 3   3   C4  8   0.246
Drug A   Cycle 3   3   C3  9   1.100
Drug A   Cycle 3   3   C4  9   0.248
Drug A   Cycle 4   3   C3 10   1.090
Drug A   Cycle 4   3   C4 10   0.234
Drug A   Cycle 4   3   C3 11   0.937
Drug A   Cycle 4   3   C3 12   0.980
Drug A   Cycle 1   4   C3  1   1.220
Drug A   Cycle 1   4   C4  1   0.182
Drug A   Cycle 1   4   C3  2   0.983
Drug A   Cycle 1   4   C4  2   0.132
Drug A   Cycle 1   4   C3  3   0.979
Drug A   Cycle 1   4   C4  3   0.128
Drug A   Cycle 2   4   C3  4   1.190
Drug A   Cycle 2   4   C4  4   0.134
Drug A   Cycle 2   4   C3  5   1.010
Drug A   Cycle 2   4   C4  5   0.076
Drug A   Cycle 2   4   C3  6   1.100
Drug A   Cycle 2   4   C4  6   0.083
Drug A   Cycle 3   4   C3  7   1.140
Drug A   Cycle 3   4   C4  7   0.108
Drug A   Cycle 3   4   C3  8   1.140
Drug A   Cycle 3   4   C4  8   0.104
Drug A   Cycle 3   4   C3  9   1.120
Drug A   Cycle 3   4   C4  9   0.080
;
run;

proc sort data=immune out=sasuser.immune;
   by pt;
run;

现在是图表的模板代码:

proc template;
  define statgraph gsv00251;
    begingraph / designwidth=7in designheight=4.5in;
      entrytitle 'Immunology Profile by Treatment';
      layout gridded;
        layout datalattice rowvar=lbparm /
               headerlabeldisplay=value includemissingclass=false
               columnaxisopts=(label="Cycle Day" griddisplay=on type=discrete
               discreteopts=(tickdisplaylist=
               ("0" "15" "30" "0" "15" "30" "0" "15" "30" "0" "15" "30"))  
               offsetmin = .03 offsetmax = .03)
               rowaxisopts=(offsetmax=.1 label="Values Converted to SI Units "
               griddisplay=on) rowdatarange=union;
        layout prototype;
          blockplot x=xval block=cyc / datatransparency = .75
                                       display=(outline fill ) 
                                       name="block" filltype=alternate;
          seriesplot x=xval y=sival / group=pt index=pt name='a'
                                      display=all lineattrs=(pattern=1 thickness=2) 
                                      markerattrs=(symbol=circlefilled size=8);    
        endlayout;
        endlayout;
        entry ' ';
        discretelegend 'a' / title='Patient' across=8;
      endlayout; 
    endgraph;
  end;
run;

现在我想在每个块的顶部添加一些信息,所以我在网上进行了一些搜索,我发现实现这个结果的最好方法是更改​​proc模板,添加2个格子行:

1)包含具有这些信息的块图的第一个格子行,如此图https://support.sas.com/documentation/cdl/en/grstatug/62464/HTML/default/viewer.htm#n07sfa0fbx6yh5n1dyutt7boc3kp.htm(具有箱线图和统计数据的那个,第二个示例)。

2)第二个格子行包含先前发布的图。

所以我创建了这个原始样本数据集:

data firstrow;
length ROW $20 TITLE $20 REASON $50 CHAR $20 VAR $20;
infile datalines dlm=',';
input ROW TITLE REASON CHAR VAR;
datalines;
ROW 1,ROW 1,FIRST REASON,FIRST ROW,FIRST VAR,
ROW 2,ROW 2,SECOND REASON,SECOND ROW,SECOND VAR,
ROW 3,ROW 3,THIRD REASON,THIRD ROW,THIRD VAR,
ROW 4,ROW 4,FOURTH REASON,FOURTH ROW,FOURTH VAR,
;
proc transpose data=firstrow out=firstrow_;
by row;
var title reason char var;
run;

并将其与之前的数据集合并(否):

proc sort data=immune out=immune;
   by pt;
run;

data immune;
merge immune firstrow_;
run;

然后修改了proc模板添加layout lattice和blockplot。

proc template;
  define statgraph gsv00251;
    begingraph / designwidth=7in designheight=4.5in;
      entrytitle 'Immunology Profile by Treatment';

     layout lattice / columns=1 rowweights=(.15 .85);
              blockplot x=row block=col1  / class=_name_
                        datatransparency = .75
                        display=(outline fill values label)
                        name="stats" filltype=alternate 
                        valuehalign=center
                        includemissingclass=false
                        labelattrs=GraphDataText valueattrs=GraphDataText ;
      layout gridded;
        layout datalattice rowvar=lbparm /
               headerlabeldisplay=value includemissingclass=false
               columnaxisopts=(label="Cycle Day" griddisplay=on type=discrete
               discreteopts=(tickdisplaylist=
               ("0" "15" "30" "0" "15" "30" "0" "15" "30" "0" "15" "30"))  
               offsetmin = .03 offsetmax = .03)
               rowaxisopts=(offsetmax=.1 label="Values Converted to SI Units "
               griddisplay=on) rowdatarange=union;
        layout prototype;
          blockplot x=xval block=cyc / datatransparency = .75
                                       display=(outline fill ) 
                                       name="block" filltype=alternate;
          seriesplot x=xval y=sival / group=pt index=pt name='a'
                                      display=all lineattrs=(pattern=1 thickness=2) 
                                      markerattrs=(symbol=circlefilled size=8);    
        endlayout;
        endlayout;
        entry ' ';
        discretelegend 'a' / title='Patient' across=8;
      endlayout; 

    endlayout;

    endgraph;
  end;
run;

结果很好,但并非最不重要,我的意思是,柱轴未与两个图对齐,结果是第一个贴图。

基本上有两个问题:右侧的标题标签和左侧的轴标签。我可以改变这个情节,以获得正确的结果吗? (轴对齐)

由于我发布了大量代码,因此您只能找到所有这些代码的最终版本

 data immune;
   format sival f6.3;
   input trt $ 1-8 cyc $ 10-18 pt lbparm $ xval sival;
   cards;
Drug A   Cycle 1   1   C3  1   1.120
Drug A   Cycle 1   1   C4  1   0.147
Drug A   Cycle 1   1   C3  2   1.080
Drug A   Cycle 1   1   C4  2   0.131
Drug A   Cycle 1   1   C3  3   0.887
Drug A   Cycle 1   1   C4  3   0.113
Drug A   Cycle 2   1   C3  4   1.440
Drug A   Cycle 2   1   C4  4   0.278
Drug A   Cycle 2   1   C3  5   1.180
Drug A   Cycle 2   1   C4  5   0.234
Drug A   Cycle 2   1   C3  6   1.360
Drug A   Cycle 2   1   C4  6   0.281
Drug A   Cycle 3   1   C3  7   1.190
Drug A   Cycle 3   1   C4  7   0.282
Drug A   Cycle 3   1   C3  8   1.000
Drug A   Cycle 3   1   C4  8   0.228
Drug A   Cycle 3   1   C3  9   1.040
Drug A   Cycle 3   1   C4  9   0.228
Drug A   Cycle 4   1   C3 10   0.917
Drug A   Cycle 4   1   C4 10   0.163
Drug A   Cycle 4   1   C3 11   0.789
Drug A   Cycle 4   1   C4 11   0.136
Drug A   Cycle 4   1   C3 12   0.861
Drug A   Cycle 4   1   C4 12   0.148
Drug A   Cycle 1   2   C3  1   1.180
Drug A   Cycle 1   2   C4  1   0.264
Drug A   Cycle 1   2   C3  2   0.942
Drug A   Cycle 1   2   C4  2   0.184
Drug A   Cycle 1   2   C3  3   1.010
Drug A   Cycle 1   2   C4  3   0.160
Drug A   Cycle 2   2   C3  4   1.050
Drug A   Cycle 2   2   C4  4   0.134
Drug A   Cycle 2   2   C3  5   0.980
Drug A   Cycle 2   2   C4  5   0.120
Drug A   Cycle 2   2   C3  6   1.020
Drug A   Cycle 2   2   C4  6   0.126
Drug A   Cycle 3   2   C3  7   0.961
Drug A   Cycle 3   2   C4  7   0.110
Drug A   Cycle 3   2   C3  8   0.859
Drug A   Cycle 3   2   C4  8   0.091
Drug A   Cycle 3   2   C3  9   0.928
Drug A   Cycle 3   2   C4  9   0.097
Drug A   Cycle 4   2   C3 10   1.380
Drug A   Cycle 4   2   C4 10   0.330
Drug A   Cycle 4   2   C3 11   1.210
Drug A   Cycle 4   2   C4 11   0.281
Drug A   Cycle 4   2   C3 12   1.180
Drug A   Cycle 4   2   C4 12   0.278
Drug A   Cycle 1   3   C3  1   1.180
Drug A   Cycle 1   3   C4  1   0.269
Drug A   Cycle 1   3   C3  2   1.010
Drug A   Cycle 1   3   C4  2   0.213
Drug A   Cycle 1   3   C3  3   1.040
Drug A   Cycle 1   3   C4  3   0.200
Drug A   Cycle 2   3   C3  4   1.200
Drug A   Cycle 2   3   C4  4   0.332
Drug A   Cycle 2   3   C4  5   0.371
Drug A   Cycle 2   3   C4  6   0.316
Drug A   Cycle 3   3   C4  7   0.271
Drug A   Cycle 3   3   C3  8   1.050
Drug A   Cycle 3   3   C4  8   0.246
Drug A   Cycle 3   3   C3  9   1.100
Drug A   Cycle 3   3   C4  9   0.248
Drug A   Cycle 4   3   C3 10   1.090
Drug A   Cycle 4   3   C4 10   0.234
Drug A   Cycle 4   3   C3 11   0.937
Drug A   Cycle 4   3   C3 12   0.980
Drug A   Cycle 1   4   C3  1   1.220
Drug A   Cycle 1   4   C4  1   0.182
Drug A   Cycle 1   4   C3  2   0.983
Drug A   Cycle 1   4   C4  2   0.132
Drug A   Cycle 1   4   C3  3   0.979
Drug A   Cycle 1   4   C4  3   0.128
Drug A   Cycle 2   4   C3  4   1.190
Drug A   Cycle 2   4   C4  4   0.134
Drug A   Cycle 2   4   C3  5   1.010
Drug A   Cycle 2   4   C4  5   0.076
Drug A   Cycle 2   4   C3  6   1.100
Drug A   Cycle 2   4   C4  6   0.083
Drug A   Cycle 3   4   C3  7   1.140
Drug A   Cycle 3   4   C4  7   0.108
Drug A   Cycle 3   4   C3  8   1.140
Drug A   Cycle 3   4   C4  8   0.104
Drug A   Cycle 3   4   C3  9   1.120
Drug A   Cycle 3   4   C4  9   0.080
;
run;


data firstrow;
length ROW $20 TITLE $20 REASON $50 CHAR $20 VAR $20;
infile datalines dlm=',';
input ROW TITLE REASON CHAR VAR;
datalines;
ROW 1,ROW 1,FIRST REASON,FIRST ROW,FIRST VAR,
ROW 2,ROW 2,SECOND REASON,SECOND ROW,SECOND VAR,
ROW 3,ROW 3,THIRD REASON,THIRD ROW,THIRD VAR,
ROW 4,ROW 4,FOURTH REASON,FOURTH ROW,FOURTH VAR,
;

proc transpose data=firstrow out=firstrow_;
by row;
var title reason char var;
run;

proc sort data=immune out=immune;
   by pt;
run;

data immune;
merge immune firstrow_;
run;

proc template;
  define statgraph gsv00251;
    begingraph / designwidth=7in designheight=4.5in;
      entrytitle 'Immunology Profile by Treatment';

     layout lattice / columns=1 rowweights=(.15 .85);
                columnaxes;
                    columnaxis / display=(ticks tickvalues);
                endcolumnaxes;



              blockplot x=row block=col1  / class=_name_
                        datatransparency = .75
                        display=(outline fill values label)
                        name="stats" filltype=alternate 
                        valuehalign=center
                        includemissingclass=false
                        labelattrs=GraphDataText valueattrs=GraphDataText ;
      layout gridded;
        layout datalattice rowvar=lbparm /
               headerlabeldisplay=value includemissingclass=false
               columnaxisopts=(label="Cycle Day" griddisplay=on type=discrete
               discreteopts=(tickdisplaylist=
               ("0" "15" "30" "0" "15" "30" "0" "15" "30" "0" "15" "30"))  
               offsetmin = .03 offsetmax = .03)
               rowaxisopts=(offsetmax=.1 label="Values Converted to SI Units "
               griddisplay=on) rowdatarange=union;
        layout prototype;
          blockplot x=xval block=cyc / datatransparency = .75
                                       display=(outline fill ) 
                                       name="block" filltype=alternate;
          seriesplot x=xval y=sival / group=pt index=pt name='a'
                                      display=all lineattrs=(pattern=1 thickness=2) 
                                      markerattrs=(symbol=circlefilled size=8);    
        endlayout;
        endlayout;
        entry ' ';
        discretelegend 'a' / title='Patient' across=8;
      endlayout; 

    endlayout;

    endgraph;
  end;
run;
ods listing close;
ods html image_dpi=100 file='Immunology.html' path='.';
ods graphics / reset noborder width=600px height=400px 
  imagename='ClinicalHandout_Immunology' imagefmt=gif noscale;

proc sgrender data=immune template=gsv00251; 
run;
ods html close;
ods listing;

1 个答案:

答案 0 :(得分:0)

来自sas社区的交叉线程(谢谢,Sanjay)。 https://communities.sas.com/t5/SAS-GRAPH-and-ODS-Graphics/SAS-Graph-Align-axis-with-lattice-rows/m-p/229696

当您使用带有单独单元格的外部LAYOUT LATTICE容器作为块图和DATAPANEL图时,所包含图的轴之间没有对齐。这些情节彼此不了解。

面板中的BLOCKPLOT具有X = XVAL,而顶部的BLOCKPLOT具有X = ROW。这些是自变量。此外,顶部图表是独立的,不需要为下部的行标题保留空间。

执行此操作的最佳方法是避免使用数据面板。使用一个LAYOUT LATTICE三行。底部的两行将包含“C3”和“C4”数据的系列图。您可以使用表达式仅提取LBPARM的一个值的数据。使用INSET或ROWHWADER放置类值。顶行将使用块图或AxisTable进行统计。确保x变量相同或相似。

了解您正在使用的SAS版本非常有用。使用SAS 9.4M3会更容易。