全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 SAS专版
5849 8
2011-07-15
有一堆数据,画出直方图后,推测可能是正太,或伽马分布等等,然后我想知道用SAS的那个过程估计这些概率密度函数呢?请高手赐教
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

全部回复
2011-7-15 11:30:33
你可以先判断是否服从正态分布,然后计算出均值和方差,再代入公式求解
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2011-7-15 12:41:20
honghejing 发表于 2011-7-15 11:30
你可以先判断是否服从正态分布,然后计算出均值和方差,再代入公式求解
正态分布是最简单的情况,如果是伽马分布,K分布,对数正态等呢?
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2011-7-15 13:36:59
貌似SAS没有现成的proc,你要自己编程实现
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2011-7-15 14:36:27
4# honghejing 那也太悲剧了
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2011-7-16 11:33:27
橡树下的往事 发表于 2011-7-15 11:08
有一堆数据,画出直方图后,推测可能是正太,或伽马分布等等,然后我想知道用SAS的那个过程估计这些概率密度函数呢?请高手赐教
Here is an example to estimate ditribution parameters by maximizing log-likehood. The sample data is Laplace ditribution.

data t1;
    mu=10; b=1;
   do i=1 to 1500;
      U=ranuni(998)-0.5;
      x=mu - b*sign(u)*log(1-2*abs(u));
   output;
end;
KEEP X;
run;
proc univariate data=t1 noprint;
histogram x/normal;
inset n mean std skewness kurtosis;
run;
proc univariate data=t1 noprint;
histogram x/gamma;
inset n mean std skewness kurtosis;
run;
proc nlmixed data=t1;
parms a=1 b=2;
loglik=log(pdf('laplace', x, a, b));
MODEL X ~ GENERAL (loglik);
RUN;
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群