%macro kde(x1);
proc kde data=work.&x1;
univar  korea00 korea08 / percentiles unistats out=a1 plots=histdensity;
run;
data b1;
set a1;
if var="korea00" then do;
  vkorea00&x1=value;dkorea08&x1=density;
end;
if var="korea08" then do;
  vkorea08&x1=value;dkorea08&x1=density;
end;

drop var value density;

goptions reset=global gunit=pct border cback=white ctext=black
   colors=(white black green red)
   ftitle=swissb ftext=swiss htitle=6 htext=4;
proc gplot data=work.a1;
plot dkorea00&x1*vkorea00&x1=1 dkorea08&x1*vkorea08&x1=2 /overlay ca=black;
symbol1 interpol=join v=none i=join c=red l=1;
symbol2 interpol=join v=none i=join c=blue l=1;

run;
%mend kde;
%kde(korea);

我是要把数据korea中的korea00和korea08两个变量的取值生成核密度数据，然后画出概率密度的图像

但是我运行之后生成的数据a1里面的desnity居然有大于1的数，因为我对核密度理解的不是特别深刻

主要是我想密度应该都是小于1的，我就觉得可能我的程序出了什么问题！~

而且我的图像也是画不出来，，，不知道出了什么问题！~


麻烦大家帮帮忙啊，，，如果能解决生成核密度数据的程序就可以了，，，画图像我还有另一种方法，可能就是麻烦点，，，
但是生成数据这一步我自己真的不知道出现了什么问题！


谢谢谢谢大家了。