Statalist The Stata Listserver

[Date Prev][Date Next][Thread Prev][Thread Next][Date index][Thread index]

st: Klein-Spady semi-parametric estimator

From   "D'Exelle Ben" <>
To   <>
Subject   st: Klein-Spady semi-parametric estimator
Date   Thu, 9 Mar 2006 10:40:34 +0100


Has anyone ever programmed the Klein-Spady semiparametric estimator? I
have tried to do so by creating my own ml program and including a kernel
function; the ml passes all tests but when maximizing the likelihood
function it shows "(non-concave)". Any ideas how to proceed now?

Reference: Klein, R. & R.Spady (1993) "An efficient semiparametric
estimator for discrete choice models", Econometrica, 61, 387-421. Or:
Green (2000: p.847) for the likelihood function.

Here is my code:

program klein_spady_lf

	version 9.1
	args lnf xb
	quietly {
		tempvar number_yes number_no tmp1 tmp2 Gn

		summarize $ML_y1, meanonly
		gen double `number_yes' = r(mean)
		gen `number_no' = 1 - `number_yes'
		kdensity `xb' if $ML_y1==1, nogr gaussian width(.5)
gen(`tmp1') at(`xb')
		kdensity `xb' if $ML_y1==0, nogr gaussian width(.5)
gen(`tmp2') at(`xb')
		replace `tmp1' = `tmp1'*`number_yes'
		replace `tmp2' = `tmp2'*`number_no'

		gen double `Gn' = `tmp1'/(`tmp1'+`tmp2')

		replace `lnf' = log(`Gn')/_N   if $ML_y1==1
		replace `lnf' = log(1-`Gn')/_N if $ML_y1==0



Ben D'Exelle
Institute of Development Policy and Management 
University of Antwerp

*   For searches and help try:

© Copyright 1996–2021 StataCorp LLC   |   Terms of use   |   Privacy   |   Contact us   |   What's new   |   Site index