*****************************************************

            Mínimos Quadrados não Linear

******************************************************


use "c:\arquivos\stataprg\naolinear_ex.dta", clear

nl ( Y= {beta}*X1 + ({beta}^2)*X2 )

/* exemplo pag. 500 e 507 Judge et al. */



**********************************************************

       outro exemplo - CES

**********************************************************

    

  use "c:\arquivos\stataprg\naolinear_CES.dta", clear

  gen lnprod = ln(produto)

     nl  ( lnprod = {b0=1} + {eta=-1}*ln({delta=0.5}*mao_obra^{rho=-1} + (1-{delta})*capital^{rho}))

 
/* ver arquivo CES.doc p/ fórmula */


***************************************

           Probite e Logite

***************************************

infile x y using c:\arquivos\stataprg\griffiths747.dat

probit y x

dprobit y x

/* para obter efeito marginal ou use */

mfx

mfx, at(x=20)

predict p

list

logit y x 

mfx

***************************************

           Logite Multinomial

***************************************


infile norte seste centro sul jovem velho branco pardo 
trabalha escola idadecri edumae edupai peso classe rendanaosalpc using c:\arquivos\stataprg\boyrural.dat

/* classe 0 = só estuda, classe 1 = só trabalha, classe 2 = não estuda nem trabalha, classe 3 = trabalha e estuda */

mlogit classe  edumae edupai idadecri jovem velho rendanaosalpc seste centro sul branco pardo

mfx, predict(outcome(0))
mfx, predict(outcome(1))
mfx, predict(outcome(2))
mfx, predict(outcome(3))

*****************************************

         Probite Bivariado

*****************************************

biprobit (trabalha =  edumae edupai rendanaosalpc idadecri norte seste centro sul jovem velho branco pardo) 
(escola = edumae edupai rendanaosalpc idadecri norte seste centro sul jovem velho branco pardo)

mfx, predict(p11) /* trabalhar e estudar */
mfx, predict(p10) 
mfx, predict(p01) 
mfx, predict(p00)
mfx, predict(pcond1)  /* trabalhar dado que estuda */
mfx, predict(pcond2)  /* estudar dado que trabalha */


*************************************

          Tobite

**************************************

use c:\arquivos\stataprg\auto.dta

gen wgt = weight/1000

regress mpg wgt

/* quanto maior o peso do automovel menor é a milhagem */

replace mpg=17 if mpg<=17

tobit mpg wgt, ll(17)

/* ll= left censoring limit (ul=right censoring limit) */

mfx 


****************************************

         HECKMAN

**************************************

set memory 250m

set matsize 200

infile idade idade2 black parda urban north south central seste industry service  union
spouse son_daug relative less5 age514 age15ml age15fm v4729 anosestu amarelo nlincpc rendhora
using c:\arquivos\stataprg\womensp.dat

gen logw = ln(rendhora)

heckman logw idade idade2 anosestu  amarelo black parda urban industry  service  union, 
select (idade idade2 anosestu amarelo black parda urban  spouse less5 age514 age15ml age15fm nlincpc)twostep mills(lambda)

mfx, predict(yexpected)

mfx, predict(ycond)