大神看看我的代码问题出现在哪里

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

总是提示where 子句运算符要求兼容的变量，我做的是一个10-fold cross- validation 交叉验证，因变量是字符型，自变量是数字，

1. %macro xval(dsn=,outcome=,covars=,k=10,sel=stepwise,outdsn=_xval_,outdsn2=comparison);
   
2. 
   
3. data _modif;
   
4.   set &dsn;
   
5.   unif=&k*ranuni(20052905);
   
6.   xv=ceil(unif);
   
7. run;
   
8. 
   
9. %do i=1 %to &k;
   
10.   proc logistic data=_modif(where=(xv ne &i)) outmodel=_mod&i;
    
11.   model &outcome (event="1") =&covars / selection=&sel;
    
12.   ods output association=assoc&i;
    
13.   run;
    
14. %if print^=0 %then %do;proc printto file='junk.txt';%end;
    
15.   proc logistic inmodel=_mod&i;
    
16.   score data=_modif(where=(xv=&i)) out=out&i;
    
17.   run;
    
18. 
    
19. ods select none;
    
20. ods output KolSmir2Stats=KS&i;
    
21. proc npar1way data= out&i edf;
    
22. where &outcome^=.;
    
23. class &outcome;
    
24. var  P_1;
    
25. run;
    
26. ods select all;
    
27. 
    
28. ods select none;
    
29. ods output WilcoxonScores=Wil&i;
    
30. proc npar1way wilcoxon data= out&i;
    
31. where &outcome^=.;
    
32. class &outcome;
    
33. var  P_1;
    
34. run;
    
35. ods select all;
    
36. 
    
37. data AUC&i;
    
38. set Wil&i end=eof;
    
39. retain v1 v2 1;
    
40. if _n_=1 then v1=abs(ExpectedSum - SumOfScores);
    
41. v2=N*v2;
    
42. if eof then do;
    
43. d=v1/v2;
    
44. /*Gini=d * 2; */
    
45. Scoring_AUC= d + 0.5;
    
46. put Scoring_AUC=;
    
47. put "****Open work.results dataset to see results of training datasets....";
    
48. keep Scoring_AUC;
    
49. output;
    
50. end;
    
51. run;
    
52. 
    
53. %if print^=0 %then %do;proc printto;run;%end;
    
54. %end;
    
55. data &outdsn;
    
56. set %do j=1 %to &k;out&j %end;;
    
57. run;
    
58. 
    
59. data training (keep =label2 nvalue2 rename= (nvalue2=Training_AUC));
    
60. set %do j=1 %to &k;assoc&j %end;;
    
61. where label2= 'c';
    
62. if label2='c' then label2 ='AUC';
    
63. run;
    
64. 
    
65. data ks (keep =label2 nvalue2 rename= (nvalue2=Scoring_KS));
    
66. set %do j=1 %to &k;ks&j %end;;
    
67. where label2= 'D';
    
68. if label2='D' then label2 ='KS';
    
69. run;
    
70. 
    
71. data validation;
    
72. set %do j=1 %to &k;AUC&j %end;;
    
73. run;
    
74. 
    
75. data &outdsn2 (drop = label2);
    
76. merge training validation ks;
    
77. run;
    
78. 
    
79. ods select none;
    
80. ods output WilcoxonScores=WilcoxonScore;
    
81. proc npar1way wilcoxon data= &outdsn;
    
82. where &outcome^=.;
    
83. class &outcome;
    
84. var  P_1;
    
85. run;
    
86. ods select all;
    
87. 
    
88. data AUC;
    
89. set WilcoxonScore end=eof;
    
90. retain v1 v2 1;
    
91. if _n_=1 then v1=abs(ExpectedSum - SumOfScores);
    
92. v2=N*v2;
    
93. if eof then do;
    
94. d=v1/v2;
    
95. Gini=d * 2;
    
96. AUC= d + 0.5;
    
97. put AUC=  GINI=;
    
98. put "****Open work.results dataset to see results of training datasets....";
    
99. keep AUC Gini;
    
100. output;
     
101. end;
     
102. run;
     
103. 
     
104. %mend;
     
105. 
     
106. %xval(dsn=data,outcome=var_13,covars= t569 t868 t725 t291 t323,k=10,sel=stepwise, outdsn=kfold, outdsn2=comparison);

复制代码

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：运算符

日志报错的地方贴出来

learsaas 发表于 2018-6-29 14:06
日志报错的地方贴出来

我已经解决了，谢谢。

我自己已经解决了

非常好，good！

ransuoqu 发表于 2018-6-29 22:08
我自己已经解决了

请问怎么解决的

ransuoqu 发表于 2018-6-29 22:08
我自己已经解决了

请问怎么解决的