1

我正在为以下示例数据创建交叉表：

data have1;
   input username $  betdate : datetime. stake winnings;
   dateOnly = datepart(betdate) ;
   format betdate DATETIME.;
   format dateOnly ddmmyy8.;
   datalines; 
    player1 12NOV2008:12:04:01 90 -90 
    player1 04NOV2008:09:03:44 100 40 
    player2 07NOV2008:14:03:33 120 -120 
    player1 05NOV2008:09:00:00 50 15 
    player1 05NOV2008:09:05:00 30 5 
    player1 05NOV2008:09:00:05 20 10 
    player2 09NOV2008:10:05:10 10 -10 
    player2 15NOV2008:15:05:33 35 -35 
    player1 15NOV2008:15:05:33 35 15 
    player1 15NOV2008:15:05:33 35 15 
run;
PROC PRINT; RUN;

Proc rank data=have1 ties=mean out=ranksout groups=2;
     var    stake;
     ranks  stakeRank;
run;

PROC TABULATE DATA=ranksout NOSEPS;
    VAR stake;
    class stakerank;
    TABLE stakerank, stake*N;
        TABLE stakerank, stake*(N Mean Skewness);
RUN;

我想在 PROC REPORT 中复制我在 PROC TABULATE 中所做的事情，因为我需要为均值差异检验和其他一些事情添加 p 值。但是，Skewness 似乎不是 Proc Report 中的内置功能。我该如何计算？

PROC REPORT DATA=ranksout NOWINDOWS;
  COLUMN stakerank stake, (n mean);
  DEFINE stakerank / GROUP id 'Rank for Variable Stake' ORDER=INTERNAL;
  DEFINE stake / ANALYSIS '';
  define n/format=8. ;
RUN;

感谢您在这方面的任何帮助

4

1 回答 1

0

可以按如下方式完成。

向 rankouts1 表添加一个额外的中间变量：

proc sql;
    create table withCubedDeviationsas
    select *,
    ((stake - (select avg(stake) from ranksout1 where stakeRank = main.stakeRank and  winnerRank = main.winnerRank))/(select std(stake) from ranksout1 where stakeRank = main.stakeRank and  winnerRank = main.winnerRank)) **3 format=8.2 as cubeddeviations
    from ranksout1 main;    
quit;

PROC REPORT DATA=withCubedDeviationsNOWINDOWS out=report;
    COLUMN stakerank winnerrank, ( N stake=avg cubeddeviations skewness);
    DEFINE stakerank / GROUP  ORDER=INTERNAL '';
    DEFINE winnerrank / ACROSS  ORDER=INTERNAL '';
    DEFINE cubeddeviations / analysis 'SumCD' noprint;
    DEFINE N / 'Bettors';
    DEFINE avg / analysis mean 'Avg' format=8.2;
    DEFINE skewness / computed format=8.2 'Skewness';
    COMPUTE skewness;
        _C5_ =  _C4_ * (_C2_ / ((_C2_ -1) * (_C2_ - 2)));
        _C9_ =  _C8_ * (_C6_ / ((_C6_ -1) * (_C6_ - 2)));
    ENDCOMP;
RUN;

他们为什么不直接将偏度添加到 PROC REPORT 中允许的统计信息列表中？

于 2013-09-13T19:59:21.137 回答