请问如何用SAS画 scaled t distribution of means?

shiyiming

图形请见下图.

大概原理：

假设 小组均值y1, y2,...yn 服从t 分布，自由度n-1, 位置参数（mu）为mean(y), 峰度(sigma)为sqrt(Mse/n). 要在图上画出y1,...yn点的位置（x轴），以及相关t分布的density plot. 请问如何用sas实现？ 谢谢

（实在不行的话，请推荐些用sas画图的资料给我，谢~）  

[img:35c8gy93]http://down8928.pinggu.org/UploadFile_20082009/2008-9/20089251550096702.jpg[/img:35c8gy93]

shiyiming

我不太理解你的意思，不过SAS画图可用proc plot 或proc gplot. 下面是t分布的gpolt code.
[code:2q3rpjrf]data aa;
     do x=-5 to 5 by 0.2;
            v=8;
                pi=3.1415926;
            y1=gamma((v+1)/2);
        y2=(1+x**2/v)**-(v+1)/2;
                y3=sqrt(v*pi)*gamma(v/2);
                y=y1*y2/y3;
                output;
         end;
run;

goptions reset=global gunit=in htext=0.1 htitle=0.1 ftext=simplex ftitle=simplex ;

title1 j=l 'T distribution';
axis1 order=(-5 to 5 by 0.2) label=(j=c 'X value');
axis2 order=(0 to 0.2 by 0.01) label=(j=c a=90 'Y value');

symbol1 i=join value=dot color=black h=0.1;

proc gplot data=aa;
     plot y*x/haxis=axis1 vaxis=axis2 noframe;
run;
quit;[/code:2q3rpjrf]

shiyiming

不好意思，刚才T分布公式写错了[code:3uf26nk7]data aa;
     do x=-5 to 5 by 0.2;
       v=8;
      pi=3.1415926;
       y1=gamma((v+1)/2);
        y2=(1+x**2/v)**(-(v+1)/2);
      y3=sqrt(v*pi)*gamma(v/2);
      y=y1*y2/y3;
      output;
    end;
run;

data aa;
     set aa;
         if -0.1<x<1E-12 then x=0;
run;

goptions reset=global gunit=in htext=0.1 htitle=0.1 ftext=simplex ftitle=simplex ;

title1 j=l 'T distribution';
axis1 order=(-5 to 5 by 0.2) label=(j=c 'X value');
axis2 order=(0 to 0.5 by 0.05) label=(j=c a=90 'Y value');

symbol1 i=join value=dot color=black h=0.1;

proc gplot data=aa;
     plot y*x/haxis=axis1 vaxis=axis2 noframe;
run;
quit;[/code:3uf26nk7]

shiyiming

好人啊~太感谢了！！！