变量熵计算SAS宏代码_Andy的原创帖16

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

计算变量的熵值。
输入参数： 1、数据集名 2、变量名
输出：1、宏变量：变量名_ent

%macro var_ent(ins,var)/store;

%global %sysfunc(cats(&var,_ent));

/*1 按选择变量汇总*/

proc meansdata=&ins(keep=&var) noprint;

class &var;

var &var;

output out=ent1 n=freq;

run;

/*2 根据频数统计量计算熵*/

data ent2;

retain mot;

set ent1 end=last;

if _type_=0 then mot=freq;

if _type_=1 then class_p=freq/mot;

if class_p then do;

ent=(-1)*log2(class_p);

ent_sum+ent;

end;

if last then call symput("&var"||'_ent',ent_sum);

run;

%mend;

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：Andy SAS宏 sysfunc output Global store 统计

ent=(-1)*log2(class_p); 这个公式少乘了一个概率，应该为：
ent=(-1)*class_p*log2(class_p);