*! version 1.00 94/8/1          STB-26: snp6.1
program define warpstep
  version 3.1
*This program calculates Hardle-Scott warping density estimation
*(step function)
*using a TurboPascal program which utilizes modified algorithms from
*Hardle (1991) "Smoothing Techniques with Implementations in S", 
*Springer-Verlag Series in Statistics, New York and Scott (1992) 
*"Multivariate Density Estimation. Theory, Practice, and Visualization, 
*John Wiley & Sons, New York.
!warpings
tempvar midpoi densi inter lowcut
infile `densi' `midpoi' using resfile.
gen `inter'=`midpoi'[2]-`midpoi'[1]
gen `lowcut'=`midpoi'-(`inter'/2)
label variable `densi' "Density"
label variable `lowcut' "Lower cutoff"
graph `densi' `lowcut',  c(J) s(.) 
!del resfile
end