统计步骤：
先将rp和sp都是缺失值的观测删除。
然后rp不是空值的做rp*s，得到变量rtp;
sp不是空值的做sp*s，得到变量stp;
然后根据area,style,ros进行分类；ros是页分类，area是行分类，style是列分类；
在分类里面：对rtp加总，再除以s的分类加总，得到各个分类下的平均值rap；
在分类里面：对stp加总，再除以s的分类加总，得到各个分类下的平均值sap；
我只写出了部分的代码，但是分类加总不会写了~求高手解决一下！
（我的程序就是狗屎，不看也罢！）
data s1;
set ww;
if rp=. and sp=. then delete;
run;
data s2;
set s1
if rp^=. then rtp=rp*s;
if sp^=. then stp=sp*s;
run;
data s3 s4;
set s2;
if rp^=. then output s3;
if sp^=. then output s4;
run;
proc tabulate data=s3;
class ros style area;
var rtp;
table ros,area,sum=''*rtp=''*(style='' all);
run;
proc tabulate data=s4;
class ros style area;
var stp;
table ros,area,sum=''*stp=''*(style='' all);
run;
下面提供了工作表

拜求解决！！！！