Using symput allows one to keep the desired column for each AREF&i._&j data sets and solves the post. %MACRO bootR; %do i=1% to 34; %DO j=1% to 3; DATA AREF&i._&j(keep=an1 col&j treat); set ref&i._&j; call symput ("col&j", col1); Treat='BREF'; run; ... View more

The reason the Upper and Lower values are not connecting is very simple: Series connects the values in each Series statement, not across multiple Series statements. So the UCL at time 03hr gets connected to  the UCL at 36hr. ... View more

Actually, OUTHITS should have produced the desired duplicates.   Here's a simplified example: data have; input wsub info; cards; 1 11 1 12 2 21 2 22 3 31 3 32 ; proc surveyselect data=have method=urs rep=2 rate=1 seed=3495 out=want1; cluster wsub; run; The resulting output dataset (WANT1) contains only two subjects per replicate (but this depends on the seed value): Number Obs Replicate wsub info Hits 1 1 1 11 1 2 1 1 12 1 3 1 2 21 2 4 1 2 22 2 5 2 2 21 2 6 2 2 22 2 7 2 3 31 1 8 2 3 32 1 Variable NumberHits contains the number of times each subject was selected. In this example it happened that wsub=2 was selected twice in replicate 1 and (accidentally) also twice in replicate 2. The total number of subjects (including the duplicates) in each of the two bootstrap samples is, of course, 3 (=number of subjects in dataset HAVE), as it should with rate=1.   Now, using the OUTHITS option ... proc surveyselect data=have method=urs rep=2 rate=1 seed=3495 out=want outhits; cluster wsub; run; ... the samples remain unchanged. Only their representation in the output dataset (WANT) is different: Number Obs Replicate wsub info Hits 1 1 1 11 1 2 1 1 12 1 3 1 2 21 2 4 1 2 22 2 5 1 2 21 2 6 1 2 22 2 7 2 2 21 2 8 2 2 22 2 9 2 2 21 2 10 2 2 22 2 11 2 3 31 1 12 2 3 32 1 Records with NumberHits>1 have now been copied NumberHits-1 times. Given that we used REP=2 and each of the three subjects in dataset HAVE had two observations, dataset WANT has now 2*6=12 observations (independent of the seed value). Variable NumberHits contains the multiplicities as before, but unlike dataset WANT1 the new output dataset WANT is not aggregated and thus NumberHits is actually redundant. ... View more