MuscleGain.median=function(gamma=.9){
   weight.diff=c(6, 7, 5, 10.5, 8.5, 3.5, 6.1, 4, 4.6, 4.5, 5.9, 6.5)
   N=length(weight.diff)
   Z.vec=sort(Z)
   i0=qbinom((1-gamma)/2,N,.5)
   if(pbinom(i0,N,.5)==(1-gamma)/2) i0=i0+1
   index=c(i0,N-i0+1); names(index)=c("i0","N-i0+1")
   Interval=c(Z.vec[i0],Z.vec[N-i0+1]); names(Interval)=c("Lower","Upper")
   estimate=median(Z); names(estimate)="median"
   list(estimate=estimate,index=index,Interval=Interval,
       achieved.confidence=1-2*pbinom(i0-1,N,.5))
}