<<Combinatorica`
t3dist=StudentTDistribution[3]
ndist=NormalDistribution[0,1]
zup=Quantile[ndist,.975]

Mean[t3dist]
Variance[t3dist]
n=25
B=1000
tfn=(Quantile[ndist,.75]-Quantile[ndist,.25])

tf=(Quantile[t3dist,.75]-Quantile[t3dist,.25])/tfn

tfs[z_List] = (Quantile[z,.75]-Quantile[z,.25])/tfn
CIBoot[xx_List]:={
RI=RandomInteger[{1,n},n];
xxb=Table[xx[[RI[[j]]]],{j,1,n}];
tfb=tfs[xxb]}

CIJK[xx_List]:={
CG=CyclicGroup[n];
CGM=Table[Most[CG[[j]]],{j,1,Length[CG]}];
xxm=Table[Table[xx[[CGM[[k]][[j]]]],{j,1,n-1}],{k,1,n}];
Out1=Table[tfs[xxm[[j]]],{j,1,n}];
Out2=Table[n*tfSample-(n-1)*Out1[[j]],{j,1,n}];
Out3={Out1,Out2}}

TBoot0={}
TBoot1={}
TBoot2={}
For[j=1,j<501,j++,
xx=RandomReal[t3dist,n];
tfSample=tfs[xx];
LBoot=Table[CIBoot[xx],{j,1,B}];
sd=Sqrt[Variance[Flatten[LBoot]]];
ciNB={tfSample-zup*sd,tfSample+zup*sd};
ciPB={Quantile[Flatten[LBoot],.025],Quantile[Flatten[LBoot],.975]};
LJK=CIJK[xx];
LJK2=Flatten[LJK[[1]][[2]]];
JKvar=Sqrt[Variance[LJK2]];
ciJK={tfSample-zup*Sqrt[JKvar/n],tfSample+zup*Sqrt[JKvar/n]};
TBoot0=Append[TBoot0,UnitStep[ciJK[[2]]-tf]*UnitStep[tf-ciJK[[1]]]];
TBoot1=Append[TBoot1,UnitStep[ciNB[[2]]-tf]*UnitStep[tf-ciNB[[1]]]];
TBoot2=Append[TBoot2,UnitStep[ciPB[[2]]-tf]*UnitStep[tf-ciPB[[1]]]]
]

N[Total[TBoot0]/Length[TBoot0]]
N[Total[TBoot1]/Length[TBoot1]]
N[Total[TBoot2]/Length[TBoot2]]

<<Combinatorica`
t3dist=StudentTDistribution[3]
ndist=NormalDistribution[0,1]
zup=Quantile[ndist,.975]

Mean[t3dist]
Variance[t3dist]
n=25
B=1000
tfn=(Quantile[ndist,.75]-Quantile[ndist,.25])

tf=(Quantile[t3dist,.75]-Quantile[t3dist,.25])/tfn

tfs[z_List] = (Quantile[z,.75]-Quantile[z,.25])/tfn
CIBoot[xx_List]:={
RI=RandomInteger[{1,n},n];
xxb=Table[xx[[RI[[j]]]],{j,1,n}];
tfb=tfs[xxb]}

CIJK[xx_List]:={
CG=CyclicGroup[n];
CGM=Table[Most[CG[[j]]],{j,1,Length[CG]}];
xxm=Table[Table[xx[[CGM[[k]][[j]]]],{j,1,n-1}],{k,1,n}];
Out1=Table[tfs[xxm[[j]]],{j,1,n}];
Out2=Table[n*tfSample-(n-1)*Out1[[j]],{j,1,n}];
Out3={Out1,Out2}}

TBoot0={}
TBoot1={}
TBoot2={}
For[j=1,j<501,j++,
xx=RandomReal[t3dist,n];
tfSample=tfs[xx];
LBoot=Table[CIBoot[xx],{j,1,B}];
sd=Sqrt[Variance[Flatten[LBoot]]];
ciNB={tfSample-zup*sd,tfSample+zup*sd};
ciPB={Quantile[Flatten[LBoot],.025],Quantile[Flatten[LBoot],.975]};
LJK=CIJK[xx];
LJK2=Flatten[LJK[[1]][[2]]];
JKvar=Sqrt[Variance[LJK2]];
ciJK={tfSample-zup*Sqrt[JKvar/n],tfSample+zup*Sqrt[JKvar/n]};
TBoot0=Append[TBoot0,UnitStep[ciJK[[2]]-tf]*UnitStep[tf-ciJK[[1]]]];
TBoot1=Append[TBoot1,UnitStep[ciNB[[2]]-tf]*UnitStep[tf-ciNB[[1]]]];
TBoot2=Append[TBoot2,UnitStep[ciPB[[2]]-tf]*UnitStep[tf-ciPB[[1]]]]
]

N[Total[TBoot0]/Length[TBoot0]]
N[Total[TBoot1]/Length[TBoot1]]
N[Total[TBoot2]/Length[TBoot2]]

0.`