_crcsrvc.ado

来自「是一个经济学管理应用软件 很难找的 但是经济学学生又必须用到」· ADO 代码 · 共 47 行

*! version 3.0.0  01/08/92
program define _crcsrvc
	version 3.0
	local varlist "req ex"
	local if "opt"
	local in "opt"
	parse "`*'"
	parse "`varlist'", parse(" ")
	if "`3'"!="" {
		local byc "by `3':"
	}
	cap drop _surv
	cap drop _stds
	cap drop _vlogs
	tempvar xi surv stds vlogs dead
	quietly {
		gen byte `dead' = cond(`2'==.,.,`2'!=0)
		gen `stds' = - `dead'
		sort `3' `1' `stds' /* deaths first */
		gen long `xi' = 0
		`byc' replace `xi' = 1 `if' `in'
		`byc' replace `xi' = 0 if `1'==. | `dead'==.
		`byc' replace `xi' = sum(`xi')
		`byc' replace `xi' = `xi'[_N] + 1 - `xi'
		gen `surv' = 1
		replace `surv' = (`xi' - `dead') / `xi' `if' `in'
		replace `surv' = 1 if `1'==. | `dead'==.
		`byc' replace `surv' = `surv' * `surv'[_n-1] if _n>1
		by `3' `1': replace `surv' = `surv'[_N]
		label var `surv' "Survival Probability"
		replace `stds' = `dead' / `xi' / (`xi' - 1)
		replace `xi'=-1 `if' `in' /* -1 ==> in the sample */
		replace `xi'=0 if `1'==. | `dead'==.
		replace `stds'=0 if `xi'!=-1
		`byc' gen `vlogs' = sum(`stds')
		by `3' `1': replace `vlogs' = `vlogs'[_N]
		replace `surv'=. if `xi'!=-1
		replace `vlogs'=. if `xi'!=-1
		replace `stds' = `surv'*sqrt(`vlogs')
		label var `stds' "Greenwood Survival S.D."
		label var `vlogs' "Var(log(_surv))"
	}
	rename `surv' _surv
	rename `stds' _stds
	rename `vlogs' _vlogs
end