* GRAFICANDO probabilidades predichas con CLARIFY

* El do file original de este ejercicio es de David Leblang: ologit_graph_example.do

*
* Este ejercicio de CLARIFY simula la probabilidad predicha de que los respondents
* se autodenominen de izquierda (1), centro (2) o derecha (3) conforme aumenta su edad
* y manteniendo las demás variables fijas en sus medias
* Después traza en una gráfica la probabilidad de pertenecer a cada grupo ideológico,
* y su intervalo de confianza, conforme aumenta el año de nacimiento (ie, menor edad)
*

use "D:\Mis documentos\Docencia\Metodos\Leblang\class_sweden.dta", clear

desc         // Datos de una exit poll de Suecia
label list   // Muchos de estos datos necesitan recodificarse como dummies
summ


* Este comando recodifica la variable leftright de 5 a 3 categorías para 
* simplificar el análisis:

recode leftright (1=1) (2=1) (3=2) (4=3) (5=3), gen(lr)
recode gender (1=1) (2=0), gen(female)


** Estimación del modelo ORDERED LOGIT base para cada categoría ideológica (1 a 3)

estsimp ologit lr birthyear female union

* Creando variables para almacenar los valores predichos e intervalos de confianza:

gen one=.    // prob(lr=1)--left of center
gen two=.    // prob(lr=2)
gen three=.  // prob(lr=3)--right of center
gen Uone=.   // Group one upper limit
gen Lone=.   // Group one lower limit
gen Utwo=.
gen Ltwo=.
gen Uthree=.
gen Lthree=.

* Creando un contador para los años de nacimiento 

gen year=_n+1910 in 1/74   // Since birth-years go from 1911 to 1985

* Ciclo de simulaciones variando birthyear de 1911 a 1985 

* PREDICCIONES PARA HOMBRES
forvalues i=1911(1)1985 {	// Inicio del ciclo
	setx mean		// fija todo en sus medias
      setx female 0    // simulation for males
	setx birthyear `i'		// fija el año de nacimiento

	simqi, genpr(p1 p2 p3)  // simula las probabilidades

	* Calculando la media e IC al 95% de cada probabilidad

	ci p1, level(95)
	replace one=r(mean) if year==`i'
	replace Uone=r(ub) if year==`i'
	replace Lone=r(lb) if year==`i'

	ci p2, level(95)
	replace two=r(mean) if year==`i'
	replace Utwo=r(ub) if year==`i'
	replace Ltwo=r(lb) if year==`i'

	ci p3, level(95)
	replace three=r(mean) if year==`i'
	replace Uthree=r(ub) if year==`i'
	replace Lthree=r(lb) if year==`i'

	*Eliminando los valores de este turno
	drop p1 p2 p3

}     // Denota el fin del ciclo

** Graficando las probabilidades de cada grupo y el año de nacimiento

twoway (rarea Lone Uone year) (rarea Ltwo Utwo year) (rarea Lthree Uthree year)

* PREDICCIONES PARA MUJERES
forvalues i=1911(1)1985 {	// Inicio del ciclo
	setx mean		// fija todo en sus medias
      setx female 1    // simulation for females (set to 0 for males)
	setx birthyear `i'		// fija el año de nacimiento

	simqi, genpr(p1 p2 p3)  // simula las probabilidades

	* Calculando la media e IC al 95% de cada probabilidad

	ci p1, level(95)
	replace one=r(mean) if year==`i'
	replace Uone=r(ub) if year==`i'
	replace Lone=r(lb) if year==`i'

	ci p2, level(95)
	replace two=r(mean) if year==`i'
	replace Utwo=r(ub) if year==`i'
	replace Ltwo=r(lb) if year==`i'

	ci p3, level(95)
	replace three=r(mean) if year==`i'
	replace Uthree=r(ub) if year==`i'
	replace Lthree=r(lb) if year==`i'

	*Eliminando los valores de este turno
	drop p1 p2 p3

}     // Denota el fin del ciclo

** Graficando las probabilidades de cada grupo y el año de nacimiento

twoway (rarea Lone Uone year) (rarea Ltwo Utwo year) (rarea Lthree Uthree year)