楼主: lnlhckao123
4539 2

怎样用SAS程序画卡方分布的图 [推广有奖]

  • 0关注
  • 3粉丝

副教授

61%

还不是VIP/贵宾

-

威望
0
论坛币
14832 个
通用积分
9.1001
学术水平
7 点
热心指数
4 点
信用等级
4 点
经验
12566 点
帖子
754
精华
0
在线时间
399 小时
注册时间
2010-8-21
最后登录
2024-2-8

100论坛币
怎样用SAS程序画卡方分布的图,而且自由度可以调整?谢谢各位高手指点!

最佳答案

ssliuna 查看完整内容

data chiprobt; retain n 5; retain e 0.05; do x=0 to 10 by e*2; p=probchi(x,n,0); output; end; run; goptions reset=global gunit=pct cback=white border htitle=7 htext=2 ftext=swissb colors=(black); proc gplot data=chiprobt; plot p*x/vaxis=axis1 haxis=axis2; symbol1 i=join v=none l=1 h=6.5 w=2; axis1 label=('probability') ; axis2 labe ...
关键词:sas程序 卡方分布 高手指点 自由度 程序 SAS
已有 1 人评分经验 学术水平 收起 理由
耕耘使者 + 100 + 3 对论坛有贡献

总评分: 经验 + 100  学术水平 + 3   查看全部评分

即使在人大经济论坛这个网络世界,我仍以真诚为基础与我的好友进行交往!
沙发
ssliuna 发表于 2010-10-8 17:32:06 |只看作者 |坛友微信交流群
data chiprobt;
retain n 5;
retain e 0.05;
  do x=0 to 10 by e*2;
  p=probchi(x,n,0);  output;
  end;
run;
goptions reset=global gunit=pct cback=white border     htitle=7 htext=2 ftext=swissb colors=(black);
proc gplot data=chiprobt;
     plot p*x/vaxis=axis1 haxis=axis2;
         symbol1 i=join v=none l=1 h=6.5 w=2;
         
         axis1 label=('probability') ;
         axis2 label=('X value') ;
         title1 'Chi distribution (df=5 nc=0)';
run;
quit;

这是自由度为5时候的程序,如果想调整自由度,得自己修改n的值。   
我刚开始接触SAS,水平很有限 ……………
已有 1 人评分论坛币 热心指数 收起 理由
耕耘使者 + 20 + 3 对论坛有贡献

总评分: 论坛币 + 20  热心指数 + 3   查看全部评分

使用道具

藤椅
ssliuna 发表于 2010-11-12 15:53:57 |只看作者 |坛友微信交流群
data chiprobt_;
do x=0.01 to 20 by 0.05;
   pdf_df1=pdf('chisquare',x,1);
   pdf_df2=pdf('chisquare',x,2);
   pdf_df3=pdf('chisquare',x,3);
   pdf_df5=pdf('chisquare',x,5);
   pdf_df10=pdf('chisquare',x,10);
output;
end;
run;

proc gplot data=chiprobt_;
plot (pdf_df1 pdf_df2 pdf_df3 pdf_df5 pdf_df10)*x/overlay vaxis=(0 to 1 by 0.1) legend;
title1 'chisquare distribution';
title2 'df=1 df=2 df=3 df=5 df=10';
symbol1 i=j c=red ;
symbol2 i=j c=black;
symbol3 i=j c=green;
symbol4 i=j c=orange;
symbol5 i=j c=blue;
run;
quit;


你也可以试试这个,这个是将几个不同的自由度的图画在一起。希望对你有所帮助
已有 1 人评分学术水平 热心指数 收起 理由
lnlhckao123 + 1 + 1 好的意见建议

总评分: 学术水平 + 1  热心指数 + 1   查看全部评分

使用道具

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群

京ICP备16021002-2号 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问:王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2024-4-30 16:12