这个图是一次性的。不过我经常碰到要划类似图的情况。

第一个图是stata 做的图吧

敏锐。是stata 做的图

我想到一个让两个纵坐标错开的办法，不过应该有更好的办法。

proc sgplot data=stat_qtr1 NOAUTOLEGEND;
   vline  yqtr / response=all y2axis;
   vbar yqtr / response=p_pos TRANSPARENCY=0;
   vbar yqtr / response=p_neg TRANSPARENCY=0;
   format yqtr yyqc4.;
   xaxis values=('1jan00'd to '1jul07'd by year) label = "Year";
   yaxis min=-15 max=30 label="Percentage" ;
   y2axis min =-4000 label="# of observations" ;
run;
quit;
quit;


效果图：

9# BraveMadMan

拥挤的x轴还是存在。也就是说你对x轴做value的定义是无效的。

不懂sas作图

stata作图容易一些

是呀，不知道是什么原因。

按道理讲不应该是discrete axis，因为x轴的数据是时间。但是不知为何SAS认为是discret数据。

连续的数轴是我们所愿望的。但是无论vline还是vbar用的都是discrete的数轴，因为他们都是classification 类的制图。举个例子，如果你改变date的format成year2.， 你会得到一个更宽的条形图。因为这个画法，就会造成不兼容，也就是说，你很可能不能用series一类的要求continuous 轴的画法。

我记得你曾经解决过一个问题过于scatter 和density plot 并于一张图的例子。其实，说白了，你不得不cheat sas，用univariate生成的数据来画density。用生成的对频率的估计值来替代sas在histogram里自动计算的频率。另外，或许你没有注意，你的图是依赖于endpoints的选择。这也是你为什么看到一些点的频率居然为0.也就说，你‘捏造’了一些子虚乌有的点。

一本书上的一个例子

Proc template;
define statgraph needle;
   begingraph;
    entrytitle 'Stock Trends';
    layout overlay;
      needleplot x=date y=close;
    endlayout;
   endgraph;
end;
run;

proc sgrender data=sashelp.stocks template=needle;
run;

proc sgplot data=sashelp.stocks;
title 'Stock Trends';
needle x=date y=close;
run;

16# BraveMadMan

原因是，needle plot用的是连续轴。needle 画法的特点是基于基准线做有宽度的垂线。既然如此，我猜想可以用0做基线，向上用正值，下为负值。然后在用series做vline。前提是你得把你的数据变一变。

先这样吧，虽然还能再调。

goptions reset=(symbol pattern axis);

proc gplot data=stat_qtr1;
format yqtr yyqc.;
  plot p_pos*yqtr=1 p_neg*yqtr=2/ overlay hminor=0 vaxis=axis1 haxis=axis2;
  plot2 all*yqtr=3 / overlay vaxis=axis3 ;
  symbol1 i=needle c=red w=4;
  symbol2 i=needle c=blue w=4;
  symbol3 i=join c=green w=3 v= h=2;
  title1 "Overlay bar graph with line";
  axis1 c=blue order=(-15 to 30 by 10) label=(angle=90 rotate=0 'Percentage');
  axis2 c=black label=('Year Qtr');
  axis3 c=green order=(-5000 to 5000 by 1000)  label=(angle=90 rotate=0 '# of observations');
run;
quit;

20# jingju11

多谢jinggu11。看你的图，那个-5000在图上是不显示的，不过为什么我的就不行呢？



并且proc format的log显示有错误：
214  title;
215  proc format;
216  value Yfmt -5000 = ' ' 0-5000 = [best.];
NOTE: The BEST (in)format was specified on the right-hand side of an equal sign, but without a
      length specification. PROC FORMAT will assume a default length of at least 40 for the
      format being generated. If this is an insufficient width, you can rerun PROC FORMAT with an
      explicit width for the BEST (in)format, or provide a sufficient DEFAULT= option.
NOTE: Format YFMT is already on the library.
NOTE: Format YFMT has been output.
217  value pfmt -15 - 15 = [best.] 35 = ' ';
NOTE: The BEST (in)format was specified on the right-hand side of an equal sign, but without a
      length specification. PROC FORMAT will assume a default length of at least 40 for the
      format being generated. If this is an insufficient width, you can rerun PROC FORMAT with an
      explicit width for the BEST (in)format, or provide a sufficient DEFAULT= option.
NOTE: Format PFMT is already on the library.
NOTE: Format PFMT has been output.
218  run;

NOTE: PROCEDURE FORMAT used (Total process time):
      real time           0.00 seconds
      cpu time            0.00 seconds

21# BraveMadMan

first, it should be (or the best) like that  [best7.] ..
另外，why 35 not gone?原因是用的是format 6.2 而不是你定义的。

22# jingju11
all done!

不是不能带价放配置吗?馒头短路了?再看看lz的内存。2g的才145。那里copy的配置啊?几年前的高端配置?

蛮复杂的，学习了。

SAS作图还可如此，学习了！

楼主是否可以把数据再放上来，以飨众人！

数据重新发了。

真讨厌，刚刚用随机数模拟了一把。

proc format;
picture posval low-high='000,009';
run;
proc sgplot data=stat_qtr1 NOAUTOLEGEND;
/*proc sgplot data=sss.bravemadman NOAUTOLEGEND;*/
needle x =yqtr y =p_pos/baseline =0 lineattrs =(color =blue thickness =3);
needle x =yqtr y =p_neg/baseline =0 lineattrs =(color =red thickness =3);
series x =yqtr y =all/y2axis lineattrs =(color =green thickness =2);;
  format yqtr year2. all p_neg p_pos posval. ;
xaxis label = "Year"; yaxis label="Percentage" ;
yaxis grid values =(-10 to 10 by 5 20) label="# of observations" ;
y2axis  values =(-5000 to 5000 by 1000) label="# of observations" ;
run; quit;

顺便问一下，你的数据格式问题，为什么我的logo报这个：
1397  proc sgplot data=sss.bravemadman NOAUTOLEGEND;



ps：楼主，我在博客里面用了这个数据，希望楼主能够同意，否则我就要我那丑陋的随机模拟数据了，谢谢。

29# sxlion

我也不知道为什么会出现这样的错误信息。