*! v1.0.0 PR 15-Mar-94.
program define _jprfpdx	/* "`devthr'" `dev' `yl' `yh' `yfl' `yfh' `stdp' */
	version 3.1
	if "`1'"=="" { exit }
	if `2'>`1' { exit }
	local yl `3'
	local yh `4'
	local yfl `5'
	local yfh `6'
	local stdp `7'
	tempvar s f l h
	local z = -invnorm((100-$S_level)/200)
	quietly {
		predict `f', xb
		predict `s', stdp
		gen `l' = `f'-`z'*`s'
		gen `h' = `f'+`z'*`s'
		replace `yl' = min(`yl',`l',`h')
		replace `yh' = max(`yh',`l',`h')
		replace `yfl' = min(`yfl',`f')
		replace `yfh' = max(`yfh',`f')
		replace `stdp' = max(`stdp',`s')
	}
end