协偏度的程序如何实现？

mzyoung

4917

收藏 2014-02-25

数据结构如下code、R1、R2、DATE（日期）、MTH(连续的月份1、2、3或者2、3、4等依次类推，不过一张表只有3个月），表的内容是多个code的3个月日交易数据，且每个code的记录条数不等。数据附件：

s1003_01.xls
大小:(32 KB)

马上下载

现在需要统计协偏度，公式如下

若在matlab中就可以用以下的程序实现单个code的协偏度计算，多个code在一张表的还是不会编程

a=r1-average(r1);

b=(r2-average(r2))^2;

c=(r2-average(r2))^3;

cs=sum(a*b)/sum(c);

但是在sas中就不能这样编程，因为mean、sum是对行数据进行统计，而我需要的是对每个code的列向量进行统计。

请诸位大神多多指导，非常感谢！

扫码加我拉你入群

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

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

全部回复

ziyenano

2014-2-25 10:48:34

proc import out=ex
datafile="E:\s1003_01.xls"
replace;
run;

proc iml;
use ex;
read all var{r1 r2} into g;
a=g[,1]-g[:,1];
b=(g[,2]-g[:,2])##2;
c=(g[,2]-g[:,2])##3;
cs=(a#b)[+]/(c[+]);
print cs;
quit;

扫码加我拉你入群

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

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

yongyitian

2014-2-25 12:06:41

复制代码

扫码加我拉你入群

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

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

mzyoung

2014-2-25 13:35:48

yongyitian 发表于 2014-2-25 12:06

谢谢，你的程序计算结果与另一位朋友提供的方法的结果一致，非常感谢

扫码加我拉你入群

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

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

mzyoung

2014-2-25 13:41:55

ziyenano 发表于 2014-2-25 10:48
proc import out=ex
datafile="E:\s1003_01.xls"
replace;

这种方法只计算出一个值，谢谢你的热心帮助

扫码加我拉你入群

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

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

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

分享

扫码加好友，拉您进群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群