SAS 求哈温平衡

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

我完全按上述输入(论坛无法显示group后及f前的美元符号）：
data a1;
input group$(a1-a2)($)f;
datalines;
U TT 34
U TG 26
U GG 5
A TT 36
A TG 33
A GG 9
;
run;
data two;
set a1;
do i=1 to f;
output two;
end;
run;
Ods html;
proc sort data=two;
by group;
proc allele data=two prefix=marker
perms=10 000 boot=1000 seed=123;
var a1 a2;
by group;
run;
Ods html close;


    结果提示出错，请教各位老师，这是什么问题呢？
  1    data a1;
2    input group$ (a1-a2)($)f;
3    datalines;

NOTE: 在第 5 行、第 1-1 列中有对“f”无效的数据。
RULE:      ----+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8----+----9----+----0
5          U TG 26
group=U a1=TT a2=34 f=. _ERROR_=1 _N_=1
NOTE: 在第 7 行、第 1-1 列中有对“f”无效的数据。
7          A TT 36
group=U a1=GG a2=5 f=. _ERROR_=1 _N_=2
NOTE: 在第 9 行、第 1-1 列中有对“f”无效的数据。
9          A GG 9
group=A a1=TG a2=33 f=. _ERROR_=1 _N_=3
NOTE: INPUT 语句到达一行的末尾，SAS 已转到新的一行。
NOTE: 数据集 WORK.A1 有 3 个观测和 4 个变量。
NOTE: “DATA 语句”所用时间（总处理时间）:
      实际时间          0.95 秒
      CPU 时间          0.04 秒


10   ;
11   run;
12   data two;
13   set a1;
14   do i=1 to f;
15   output two;
16   end;
17   run;

ERROR: 无效的 DO 循环控制信息，INITIAL 或 TO 表达式缺失，或 BY 表达式缺失、为 0 或无效。
group=U a1=TT a2=34 f=. i=1 _ERROR_=1 _N_=1
NOTE: 由于出错，SAS 系统停止处理该步。
NOTE: 从数据集 WORK.A1. 读取了 1 个观测
WARNING: 数据集 WORK.TWO 可能不完整。该步停止时，共有 0 个观测和 5 个变量。
NOTE: “DATA 语句”所用时间（总处理时间）:
      实际时间          0.15 秒
      CPU 时间          0.01 秒


18   Ods html;
NOTE: 正在写入 HTML Body（主体）文件: sashtml.htm
19   proc sort data=two;
20   by group;

NOTE: 输入数据集为空。
NOTE: 数据集 WORK.TWO 有 0 个观测和 5 个变量。
NOTE: “PROCEDURE SORT”所用时间（总处理时间）:
      实际时间          0.09 秒
      CPU 时间          0.01 秒


21   proc allele data=two prefix=marker
22   perms=10 000 boot=1000 seed=123;
              ---
              22
              200
ERROR 22-322: 语法错误，期望下列之一: ;, ALLELEMIN, ALPHA, AMIN, BOOT, BOOTSTRAP, CORRCOEFF, DATA, DELIMITER, DELTA, DPRIME,
              GENOCOL, GENOMIN, GMIN, HAPLO, HAPLOMIN, HMIN, INDIV, INDIVIDUAL, LOGNOTE, MARKER, MAXDIST, NDATA, NOFREQ, NOPRINT,
              OUTSTAT, PERMS, PREFIX, PROPDIFF, RHO, SEED, TALL, YULESQ.
ERROR 200-322: 该符号不可识别，将被忽略。
23   var a1 a2;
24   by group;
25   run;

NOTE: 由于出错，SAS 系统停止处理该步。
NOTE: “PROCEDURE ALLELE”所用时间（总处理时间）:
      实际时间          0.10 秒
      CPU 时间          0.03 秒

26   Ods html close;

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：Individual Procedure delimiter individua Bootstrap 平衡

i. 建立第一个dataset的时候datalines下面出了错，系统说了是到行尾只有有3 个观测，而frame有 4 个变量，最后一个变量f没有值读入；
ii. suggestion: 把datalines里面的‘TT TG GG ... ...’ 用空格分开，可保证有四个数值读入。

1. data a1;
   
2. input group $ (a1-a2)($) f;
   
3. datalines;
   
4. U T T 34
   
5. U T G 26
   
6. U G G 5
   
7. A T T 36
   
8. A T G 33
   
9. A G G 9
   
10. ;
    
11. run;
    
12. proc print data=a1 (obs=2);Run;
    
13. data two;
    
14. set a1;
    
15. do i=1 to f;
    
16. output two;
    
17. end;
    
18. run;
    
19. data _null_;
    
20. ans=34+26+5+36+33+9;
    
21. put ans;
    
22. Run;

复制代码

shangzhixi 发表于 2014-9-16 09:46
i. 建立第一个dataset的时候datalines下面出了错，系统说了是到行尾只有有3 个观测，而frame有 4 个变量，最 ...

多谢shangzhixi兄指点，按您的建议修改后，再把“perms=10 000” 改为“perms=10000”，然后把SAS关了重启就可以了。

20   proc allele data=two; prefix=marker perms=10000 boot=1000 seed=123;
ERROR: PROCEDURE ALLELE 没有找到。
21   var a1 a2;
22   by group;
23   run;
请教一下，同样的程序，上面是SAS日志的一部分，我的怎么procedure allele没有找到呀