clear 

/* Se asume que se guardó el archivo nlsw88es.xlsx en tu directorio
de trabajo actual. Para saber cual es el directorio activo, ejecuta el 
comando -pwd-. */
import excel nlsw88es.xlsx, firstrow

generate ln_salario = ln(salario)

label variable ln_salario "Logaritmo del salario por hora"
label variable salario "Salario por hora"
label variable graduniv "Graduada de la universidad"
label variable unión "Afiliada a la unión"
label variable experiencia "Experiencia de trabajo (años)"

label define sino 0 "no" 1 "si"
label values unión sino
label values graduniv sino

encode zona, generate(nzona)

describe

summarize ln_salario

tabulate unión graduniv, column

table unión, contents(mean salario sd salario)

histogram ln_salario, normal

twoway scatter salario experiencia

twoway scatter salario experiencia, by(unión)

regress ln_salario experiencia

regress ln_salario i.graduniv i.nzona i.unión##c.experiencia

margins unión, at(experiencia=(5(5)25))
marginsplot

search multinomial logistic