Hello SAS community,
I am looking for resources on how to compute jacknife estimate of a parameter using BY processing. For example, at the end of the execution, I hope to have the jacknife estimate of treatment effect (treat) for N=4 and effects=0.15, jacknife estimate for N=4 and effects=0.20, jacknife estimate for N=4 and effects=0.3, etc. These will then be used to compute the corresponding Acceleration factor for each of these BY processing. Example, toy data are below for N=4,6, and effects=0.15, 0.20, 0.30. I tried the following SAS code, but it appears not to be giving me what I want. The %jack() macro is from the SAS jackboot macro. If there is a way to generate the jacknife samples with SURVEYSELECT Procedure, that will simplify analysis. Any help is appreciated.
Thanks.
%ODSOff ; ods output estimates=est; proc glm data=chg; class treat; model post= b treat ; estimate "difference" treat 1 -1; by N effects; run;
%let by=N effects;
%ODSOff ; %macro analyze(data=,out=); ods output estimates=&out(where=(parameter eq 'difference') drop=stderr tvalue probt); proc glm data=&data; class treat; model post=b treat ; estimate "difference" treat 1 -1; %bystmt; run; %mend;
%jack(data=chg, id=parameter);
N
b
treat
effects
post
4
6.36910
0
0.15
6.97049
4
6.52892
0
0.15
7.40859
4
6.26562
1
0.15
7.56820
4
7.05193
1
0.15
7.28983
4
6.36910
0
0.20
6.97049
4
6.52892
0
0.20
7.40859
4
6.26562
1
0.20
7.56820
4
7.05193
1
0.20
7.28983
4
6.36910
0
0.30
6.97049
4
6.52892
0
0.30
7.40859
4
6.26562
1
0.30
7.56820
4
7.05193
1
0.30
7.28983
6
6.23206
0
0.15
6.61017
6
6.56842
0
0.15
6.79354
6
6.85177
0
0.15
7.60327
6
6.87010
1
0.15
6.53018
6
6.48628
1
0.15
7.16911
6
7.74127
1
0.15
8.13158
6
6.23206
0
0.20
6.61017
6
6.56842
0
0.20
6.79354
6
6.85177
0
0.20
7.60327
6
6.87010
1
0.20
6.53018
6
6.48628
1
0.20
7.16911
6
7.74127
1
0.20
8.13158
6
6.23206
0
0.30
6.61017
6
6.56842
0
0.30
6.79354
6
6.85177
0
0.30
7.60327
6
6.87010
1
0.30
6.53018
6
6.48628
1
0.30
7.16911
6
7.74127
1
0.30
8.13158
... View more
