Hi all, I need to combine 10-fold cross-validation and Bootstrapping in the same macro. But i have the impression that I am missing the Bootstrapping loop. I was wondering if you would be able to advise me what is wrong on my code below please? Is there a better way to combine 10-fold cross-validation and Bootstrapping in the same macro? data kyphosis;
* infile cards dlm='09'x;
input y1 y2 d;
cards;
13 100 0
64 437 0
12 334 0
618 285 0
104 150 0
65 136 0
1573 523 0
291 927 0
84 62 0
13 54 0
338 248 1
758 917 1
189 305 1
260 88 1
223 257 1
604 231 1
366 1106 1
1094 658 1
176 65 1
1499 147 1
69 319 1
,
run;
%macro bootstrap (bootnum); *bootnum=12;
*******************************************************************************
BOOTSTRAPPING
*******************************************************************************;
proc surveyselect data=kyphosis NOPRINT seed=1234
out=kyphosis1(rename=(Replicate=bootsample))
method=urs
samprate=100
outhits
reps=&bootnum;
run;
*******************************************************************************
10-FOLD CROSS VALIDATION
*******************************************************************************;
data kyphosis1;
set kyphosis1;
theRandom = ranuni(0);
*by bootsample;
run;
proc rank data = kyphosis1 out=kRanked groups=10;
*by bootsample;
var theRandom;
run;
%do x = 0 %to 9;
******************************************************************************
Create Training Data
*******************************************************************************;
data training&x;
set kRanked;
where theRandom ne &x;
*by bootsample;
run;
******************************************************************************
Generate Binary variable beta1 for Biomarker 1
*******************************************************************************;
proc sql;
create table temp&x as
select
a.y1 as compare_y1,
b.*,
a.y1 < b.y1 as beta
from training&x as a, training&x as b;
*order by compare_y1, y1;
quit;
******************************************************************************
Generate Binary variable beta1 for Biomarker 2
*******************************************************************************;
proc sql;
create table temp2&x as
select
a.y2 as compare_y2,
b.*,
a.y2 < b.y2 as beta1
from training&x as a, training&x as b;
*order by compare_y2, y2;
quit;
******************************************************************************
Merge Temp and Temp2
*******************************************************************************;
data want&x;
merge temp&x temp2&x;
run;
******************************************************************************
Generate Combined Binary variable beta1_beta2 and Append
*******************************************************************************;
data final&x;
set want&x;
if beta=1 & beta1=1 then beta1_beta2=1;
if beta=1 & beta1=0 then beta1_beta2=1;
if beta=0 & beta1=1 then beta1_beta2=1;
if beta=0 & beta1=0 then beta1_beta2=0;
run;
******************************************************************************
Calculate se, spec, ect...
*******************************************************************************;
proc sort data=final&x;
by compare_y1 compare_y2;
run;
proc freq data = final&x order=data noprint ;
by compare_y1 compare_y2;
tables beta1_beta2*d / out=freq_results&x OUTPCT sparse ;
*output out=mnocol;
run;
******************************************************************************
Create Test Data
*******************************************************************************;
data test&x; * The test dataset, called test0 (when x=0), test1 (when x=1), etc.;
set kRanked;
where theRandom = &x;
run;
******************************************************************************
APPEND THE RESULTS OF THE SELECTED CUTOFFS IN THE TRAING SET
*******************************************************************************;
proc append base = Base force data = freq_results&x;
run;
%end; *end 10-FOLD CROSS VALIDATION loop;
%mend;
%bootstrap(12);
... View more
