*! v1.0.0 PR 25-Feb-94.         (sg26.3: STB-25)
program define _jprfpdt	/* "`devthr'" `dev' `pj' `pk' `yl' `yh' `stdp' */
	version 3.1
	if "`1'"=="" { exit }
	if `2'>`1' { exit }
	local pj `3'
	local pk `4'
	local yl `5'
	local yh `6'
	local stdp `7'
	global S_P1 "$S_P1 `pj'"
	if `pk'!=. { global S_P2 "$S_P2 `pk'" }
	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 `stdp' = max(`stdp',`s')
	}
end