楼主: priss111
3297 1

[实际应用] 请问如何对多个数据集批量处理后,批量进行卡方检验?谢谢。 [推广有奖]

  • 0关注
  • 5粉丝

副教授

15%

还不是VIP/贵宾

-

威望
0
论坛币
153 个
通用积分
47.5396
学术水平
11 点
热心指数
14 点
信用等级
10 点
经验
14938 点
帖子
473
精华
0
在线时间
805 小时
注册时间
2008-3-30
最后登录
2024-4-12

+2 论坛币
k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群
赵安豆老师微信:zhaoandou666

经管之家联合CDA

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

感谢您参与论坛问题回答

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

+2 论坛币
大家好.
请问如何对多个数据集批量整理后,再批量进行卡方检验?
如何用宏实现红色的代码?
谢谢.



第1步:如何对导入的多个txt数据集(保存在逻辑库work中)批量进行整理?
             其中txt导入数据集a为例,其他导入的txt数据集格式与a同:
data a;
        format var1 best1.;
        format var2 best2.;
        input var1 var2;  *var1=阳性例数,var2=总例数;
    cards;
    7 30
    1 32
    ;
run;

    data group;
        do group = 1 to 2;
            do r = 1 to 2;
            output;
            end;
        end;
    run;

    %do i = 1 %to &nobs.;
     data work.&name.(keep=freq);
        set work.&name.;
        var1_=var2-var1;

        freq=var1; output;
        freq=var1_; output;
    run;
   
    /***将(1)与group数据集横向合并生成所需数据格式;*/

    data work.&name.;
        merge group work.&name.;
    run;

    %end;


第2步:如何批量对步骤1中的数据集逐个进行卡方检验?
   /***proc freq 生成卡方检验的结果,并将各组频数、卡方值与P值合并在一个数据集中;***/

    %do i = 1 %to &nobs.;
    proc freq data=work.&name.;
        weight freq;
        tables group*r/expected chisq out=count(keep=count) exact;
        output out=chisq(drop=df_pchi) pchi;
    run;

    data chisq&name.;
        set chisq;
            label _pchi_="卡方值"
                  p_pchi="P值";
    run;

    proc transpose data=count out=count&name.;
    run;

    data count&name.;
        set count&name.;
        keep count;
        rename  COL1=A1
                COL2=A2
                COL3=B1
                COL4=B2;
    run;

    data countchisq&name.;
        merge count&name. chisq&name.;
    run;


二维码

扫码加我 拉你入群

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

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

关键词:批量处理 卡方检验 数据集 Transpose Expected

沙发
Cherieling 发表于 2019-12-15 15:38:31 |只看作者 |坛友微信交流群
这个楼主知道咋整了吗

使用道具

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群

京ICP备16021002-2号 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问:王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2024-4-20 06:11