From 507be9fc10be00f524942c5d29a432124ec1f61d Mon Sep 17 00:00:00 2001 From: Quarto GHA Workflow Runner Date: Wed, 3 Apr 2024 03:41:38 +0000 Subject: [PATCH] Built site for gh-pages --- .nojekyll | 2 +- 01-introduccion.html | 40 +- 02-flujo-basico-2.html | 10 +- 02-flujo-basico.html | 10 +- 03-modelos-genericos.html | 26 +- 05-dags.html | 134 +- .../figure-html/unnamed-chunk-10-1.png | Bin 30601 -> 30514 bytes .../figure-html/unnamed-chunk-13-1.png | Bin 520912 -> 543051 bytes .../figure-html/unnamed-chunk-14-1.png | Bin 513369 -> 534480 bytes .../figure-html/unnamed-chunk-22-1.png | Bin 19878 -> 19996 bytes .../figure-html/unnamed-chunk-25-1.png | Bin 14018 -> 16540 bytes .../figure-html/unnamed-chunk-28-1.png | Bin 12901 -> 22737 bytes .../figure-html/unnamed-chunk-29-1.png | Bin 26951 -> 26804 bytes .../figure-html/unnamed-chunk-32-1.png | Bin 100587 -> 106228 bytes .../figure-html/unnamed-chunk-33-1.png | Bin 210967 -> 229154 bytes .../figure-html/unnamed-chunk-37-1.png | Bin 27432 -> 27456 bytes .../figure-html/unnamed-chunk-38-1.png | Bin 28957 -> 28694 bytes .../figure-html/unnamed-chunk-41-1.png | Bin 79957 -> 84292 bytes .../figure-html/unnamed-chunk-42-1.png | Bin 186379 -> 190299 bytes .../figure-html/unnamed-chunk-50-1.png | Bin 99146 -> 99695 bytes .../figure-html/unnamed-chunk-52-1.png | Bin 91636 -> 85664 bytes .../figure-html/unnamed-chunk-8-1.png | Bin 28611 -> 28864 bytes 06-calculo-do.html | 68 +- 07-buenos-malos-controles.html | 34 +- 08-mcmc.html | 2503 ++++++++++++++++- .../figure-html/unnamed-chunk-37-1.png | Bin 0 -> 334802 bytes .../figure-html/unnamed-chunk-38-1.png | Bin 0 -> 23239 bytes .../figure-html/unnamed-chunk-39-1.png | Bin 0 -> 23384 bytes .../figure-html/unnamed-chunk-47-1.png | Bin 0 -> 72322 bytes .../figure-html/unnamed-chunk-51-1.png | Bin 0 -> 57296 bytes .../figure-html/unnamed-chunk-52-1.png | Bin 0 -> 67020 bytes .../figure-html/unnamed-chunk-52-2.png | Bin 0 -> 64374 bytes .../figure-html/unnamed-chunk-55-1.png | Bin 0 -> 112484 bytes .../figure-html/unnamed-chunk-56-1.png | Bin 0 -> 79095 bytes index.html | 2 +- search.json | 42 +- 36 files changed, 2632 insertions(+), 239 deletions(-) create mode 100644 08-mcmc_files/figure-html/unnamed-chunk-37-1.png create mode 100644 08-mcmc_files/figure-html/unnamed-chunk-38-1.png create mode 100644 08-mcmc_files/figure-html/unnamed-chunk-39-1.png create mode 100644 08-mcmc_files/figure-html/unnamed-chunk-47-1.png create mode 100644 08-mcmc_files/figure-html/unnamed-chunk-51-1.png create mode 100644 08-mcmc_files/figure-html/unnamed-chunk-52-1.png create mode 100644 08-mcmc_files/figure-html/unnamed-chunk-52-2.png create mode 100644 08-mcmc_files/figure-html/unnamed-chunk-55-1.png create mode 100644 08-mcmc_files/figure-html/unnamed-chunk-56-1.png diff --git a/.nojekyll b/.nojekyll index d16c61b..1ce1f84 100644 --- a/.nojekyll +++ b/.nojekyll @@ -1 +1 @@ -27796c71 \ No newline at end of file +1e1005d6 \ No newline at end of file diff --git a/01-introduccion.html b/01-introduccion.html index 0d5adec..5b1735e 100644 --- a/01-introduccion.html +++ b/01-introduccion.html @@ -2,7 +2,7 @@ - + @@ -355,8 +355,8 @@

Ejemp B -chicos -mejora +grandes +sin_mejora A @@ -364,24 +364,24 @@

Ejemp mejora -B +A chicos mejora -A -chicos +B +grandes mejora -A -chicos -mejora +B +grandes +sin_mejora B chicos -mejora +sin_mejora B @@ -389,18 +389,18 @@

Ejemp mejora -B -chicos -mejora - - A grandes +sin_mejora + + +B +chicos mejora A -chicos +grandes mejora @@ -654,8 +654,8 @@

Ejemp ", width = 200, height = 50)
-
- +
+

Es decir, el tamaño de los cálculos es una causa común de tratamiento (T) y resultado (M). Veremos más adelante que la decisión de condicionar a el tipo de cálculos proviene de un análisis relativamente simple de este diagrama causal, independientemente de los métodos que usemos para estimar las proporciones de interés (en este ejemplo, examinar las tablas cruzadas es equivalente a hacer estimaciones de máxima verosimlitud).

@@ -776,8 +776,8 @@

Eje ", width = 200, height = 50)
-
- +
+

Nótese que el análisis más apropiado no está en los datos: en ambos casos la tabla de datos es exactamente la misma. Los supuestos acerca del proceso que genera los datos sin embargo nos lleva a respuestas opuestas.

diff --git a/02-flujo-basico-2.html b/02-flujo-basico-2.html index 55dcf5d..c4ad151 100644 --- a/02-flujo-basico-2.html +++ b/02-flujo-basico-2.html @@ -2,7 +2,7 @@ - + @@ -314,8 +314,8 @@

")#, width = 200, height = 50)
-
- +
+

Donde vemos ahora que el estado real de cada persona de la prueba es desconocido, aunque el resultado de la prueba depende de ese estado, y la cantidad de positivos que observamos es ahora \(N_{obs}\), que depende también de la sensibilidad y especificidad de la prueba.

@@ -557,8 +557,8 @@

")#, width = 200, height = 50)
-
- +
+

Usando argumentos como los del modelo original, las distribuciones de esp y sens son beta y podemos incorporarlas en la simulación de la posterior. Nuestra nueva función para simular el proceso generativo es:

diff --git a/02-flujo-basico.html b/02-flujo-basico.html index 39b6e35..ad25564 100644 --- a/02-flujo-basico.html +++ b/02-flujo-basico.html @@ -2,7 +2,7 @@ - + @@ -333,8 +333,8 @@

", width = 300, height = 100)
-
- +
+

Que también podríamos simplificar (suponiendo la \(N\) fija y conocida, pues \(N_+\) y \(M\) dan \(N_{-}\)) como:

@@ -367,8 +367,8 @@

", width = 300, height = 100)
-
- +
+

Y ahora construimos el modelo generativo. Supondremos que la muestra de \(N\) personas se toma de manera aleatoria de la población (una población grande, así que podemos ignorar el efecto de muestreo). Supondremos provisionalmente, además, que la prueba es perfecta, es decir, no hay falsos positivos o negativos.

diff --git a/03-modelos-genericos.html b/03-modelos-genericos.html index 3f40f0c..6e35d3a 100644 --- a/03-modelos-genericos.html +++ b/03-modelos-genericos.html @@ -2,7 +2,7 @@ - + @@ -339,8 +339,8 @@

")#, width = 200, height = 50)
-
- +
+

Nótese que no consideramos \(W\to H\), porque podemos pensar en varias intervenciones que podrían cambiar el peso por no cambian la estatura. Por otro lado, es difícil pensar en alguna intervención que cambie la estatura pero no cambie el peso de una persona. Adicionalmente, hay otros factores desconocidos no observados \(U\) que afectan el peso de cada persona adicionalmente a su estatura.

@@ -934,8 +934,8 @@

", width = 200, height = 50)
-
- +
+

Omitiendo del diagrama las variables no observadas que también son causas únicamente de \(S\) y \(W, H\):

@@ -960,8 +960,8 @@

", width = 200, height = 50)
-
- +
+

Si queremos saber cómo influye el sexo en el peso, este diagrama indica que hay dos tipos de preguntas que podemos hacer:

@@ -1451,8 +1451,8 @@

")#, width = 200, height = 50)
-
- +
+

En este caso, el modelo causal es como sigue: conocemos la distancia \(D\) al hoyo en cada tiro. El éxito (\(Y=1\)) o fracaso (\(Y=0\)) depende de la distancia, junto con la velocidad a la que sale la pelota (muy alto o muy bajo puede dar un tiro fallido), y el ángulo \(\theta\) de salida. Adicionalmente, hay otros factors \(U\) que pueden afectar la probabilidad de éxito. Nótese que no escribiríamos, por ejemplo \(Y \leftarrow D\), porque la distancia no cambia causalmente con el resultado del tiro, aunque es cierto que si intervenimos en la distancia, esperaríamos obtener tasas de éxito diferentes. Igualmente, es necesario poner una flecha de \(V\) a \(D\) y \(V\) a \(Y\).

@@ -2201,12 +2201,12 @@ 

Warning: 236 of 4000 (6.0%) transitions hit the maximum treedepth limit of 10.
diff --git a/05-dags.html b/05-dags.html
index 2f0b996..e87dd49 100644
--- a/05-dags.html
+++ b/05-dags.html
@@ -2,7 +2,7 @@
 
 
 
-
+
 
 
 
@@ -349,8 +349,8 @@ 
", width = 150, height = 40)
-
- +
+

Nótese que no describimos exactamente cómo son las funciones que relacionan las variables, sino más bien qué variables son causas directas de qué otras. Por ejemplo, aunque en nuestro ejemplo de arriba \(Y\) puede estar correlacionado con \(Z\), no hay una causa directa a \(Y\), porque cambios en \(Z\) afectan a \(X\), y es el cambio en \(X\) que es causa directa de \(Y\).

@@ -407,8 +407,8 @@

")
-
- +
+

En este ejemplos no podemos saber \(U1\) y \(U2\), y no nos interesa modelar la física de monedas, manera de lanzarlas, etc. En este ejemplo también no consideraremos qué hace que un día sea soleado o lluvioso (no nos interesa modelar el clima). En este momento, en teoría tenemos ecuaciones determinísticas para todas las variables, y si conocemos todas las variables exógenas \(U1,U2,U3,U4\) podríamos determinar exactamente lo que va a suceder con la ganancia, por ejemplo, o cualquier otra variable del sistema.

@@ -439,8 +439,8 @@

")
-
- +
+
@@ -515,13 +515,13 @@

Ejemplo

simular_juego(5)
# A tibble: 5 × 5
-      x d           s1    s2     g
-  <dbl> <chr>    <int> <int> <int>
-1 0.618 soleado      4     2     4
-2 0.601 soleado      4     1     4
-3 0.938 soleado      5     5     5
-4 0.816 lluvioso     4     3     7
-5 0.718 soleado      3     2     3
+ x d s1 s2 g + <dbl> <chr> <int> <int> <int> +1 0.805 soleado 3 5 3 +2 0.891 soleado 5 4 5 +3 0.0458 lluvioso 0 0 0 +4 0.158 soleado 0 1 0 +5 0.226 soleado 2 1 2
@@ -583,8 +583,8 @@

", width = 200, height = 50)
-
- +
+

En este caso,

@@ -633,8 +633,8 @@

Ejemplo (si 
cor(sims_confusor |> select(x,y)) |> round(3)
       x      y
-x  1.000 -0.424
-y -0.424  1.000
+x 1.000 -0.417 +y -0.417 1.000

Sin embargo, si condicionamos a \(Z\), que puede tomar los valores 0 o 1, vemos que \(X\) y \(Y\) son independientes, o dicho de otra manera, la condicional de \(Y\) dada \(Z\) y \(X\) sólo depende de \(Z\):

@@ -660,14 +660,14 @@

Ejemplo (si 
cor(sims_confusor |> filter(z == 1) |> select(x,y)) |> round(3)
       x      y
-x  1.000 -0.003
-y -0.003  1.000
+x 1.000 -0.004 +y -0.004 1.000 
cor(sims_confusor |> filter(z == 0) |> select(x,y)) |> round(3)
-
       x      y
-x  1.000 -0.007
-y -0.007  1.000
+
      x     y
+x 1.000 0.002
+y 0.002 1.000

Un ejemplo con variables continuas podría ser como sigue:

@@ -767,8 +767,8 @@

E ", width = 200, height = 50)
-
- +
+

Por la discusión de arriba, es claro que es necesario considerar la edad al casarse si queremos estimar el efecto de tasa de matrimonio en la tasa de divorcio. Es posible que la correlación entre estas dos tasas puede ser explicada solamente por la edad al casarse, y que en realidad al flecha \(M\to D\) sea muy débil o inexistente.

@@ -860,7 +860,7 @@

E All 4 chains finished successfully. Mean chain execution time: 0.1 seconds. -Total execution time: 0.6 seconds. +Total execution time: 0.9 seconds.
@@ -919,8 +919,8 @@

", width = 200, height = 50)

-
- +
+

Es decir, borramos todas las flechas que caen en \(M\) (pues la estamos interveniendo al valor que queramos), y luego simulando \(D\).

@@ -936,7 +936,7 @@

@@ -988,8 +988,8 @@

", width = 150, height = 20)
-
- +
+

En este caso,

@@ -1056,15 +1056,15 @@ 

cor(sims_mediador |> filter(z == 1) |> select(x,y)) |> round(3)
-
       x      y
-x  1.000 -0.006
-y -0.006  1.000
+
      x     y
+x 1.000 0.001
+y 0.001 1.000
cor(sims_mediador |> filter(z == 0) |> select(x,y)) |> round(3)
       x      y
-x  1.000 -0.002
-y -0.002  1.000
+x 1.000 -0.004 +y -0.004 1.000

Podemos también hacer un ejemplo continuo:

@@ -1157,8 +1157,8 @@

Ejemplo: Burks")
-
- +
+

Como el NSE es del hogar (una medida general de estatus social), se consideró en principio como una variable pre-tratamiento a la inteligencia de los niños por la que tradicionalmente se controlaba. Burks notó que hacer esto tenía no era apropiado, pues tiene como consecuencia cortar parte del efecto total de la inteligencia sobre el la inteligencia de los hijos. En otras palabras: la inteligencia de los padres hace más probable mejor NSE, y mejor NSE presenta mejores condiciones de desarrollo para sus hijos. Estatificar por esta variable bloquea este efecto.

@@ -1185,8 +1185,8 @@

", width = 200, height = 50)
-
- +
+