*! version 1.0.0 PR 20Feb99.       STB-49 sg81.1
program define frac_rep
* 1=descriptor e.g. FRACTIONAL POLYNOMIAL
* 2=param descriptor e.g. df
* 3=param names e.g. powers
* Assumes existence of e(fp_*)s.
version 6
args desc param paramv
local i 1
local l=length("`paramv'")
while `i'<=e(fp_nx) {
	local l=max(`l',length("`e(fp_k`i')'"))
	local i=`i'+1
}
local l=`l'+44
local title "Final multivariable `desc' model for `e(fp_depv)'"
local lt=length("`title'")
di _n in gr "`title'"
di in gr _dup(`lt') "-"
di _n in gr "Variable |" _col(15) _dup(5) "-" _col(20) "Initial" /*
 */ _col(27) _dup(5) "-" _col(42) _dup(5) "-" _col(47) "Final" /*
 */ _col(52) _dup(5) "-"
di in gr _col(10) "| `param'"/*
 */ _col(21) "Select" /*
 */ _col(30) "Alpha" /*
 */ _col(39) "Status" /*
 */ _col(47) "`param'" /*
 */ _col(55) "`paramv'"
di in gr _dup(9) "-" "+" _dup(`l') "-"
local i 1
while `i'<=e(fp_nx) {
	local pars `e(fp_k`i')'
	if "`pars'"=="" | "`pars'"=="." {
		local final 0
		local status out
		local pars
	}
	else {
		local status in
		local final=e(fp_fd`i')
	}
	local name `e(fp_x`i')'
	local skip=8-length("`name'")
	di in gr _skip(`skip') "`name' |" in ye /*
	 */ _col(15) e(fp_id`i') /*
	 */ _col(20) %7.4f e(fp_se`i') /*
	 */ _col(29) %7.4f e(fp_al`i') /*
	 */ _col(41) "`status'" /*
	 */ _col(49) "`final'" /*
	 */ _col(55) "`pars'"
	local i=`i'+1
}
di in gr _dup(9) "-" "+" _dup(`l') "-"
`e(cmd)'
di in gr "Deviance:" in ye %9.3f e(fp_dev) in gr "."
end