Untitled
user_1300604
plain_text
10 months ago
3.9 kB
4
Indexable
### Opció A, B i C per a Catalunya (estat de salut afectat) ```{r} # Fer el datasets de Catalunya per a opcions A, B i C AB_cat <- AB_esp |> filter(CCAA == 9) C_cat <- C_ESP |> filter(CCAA == 9) # Ponderació de Catalunya utilitzant pes CCAA disseny_cat_opcio_a_b <- svydesign(ids = ~1, data = AB_cat, weights = ~PESOCCAA) # Càlcul d'A Catalunya ## Canvi de nom de la var AB_cat <- AB_cat |> mutate(vio_a_cat = vio_a_esp) preval_vio_a_cat <- svytotal(~vio_a_cat, disseny_cat_opcio_a_b) (preval_vio_a_cat / 808)* 20736963 # Càlcul de B Catalunya ## Canvi de nom de la var AB_cat <- AB_cat |> mutate(vio_b_cat = vio_b_esp) preval_vio_b_cat <- svytotal(~vio_b_cat, disseny_cat_opcio_a_b) (preval_vio_b_cat / 808)* 20736963 # Càlcul de C Catalunya # Variable usada: "vio_c_estatafectat" per a opció C disseny_cat_opcio_c <- svydesign(ids = ~1, data = C_cat, weights = ~PESOCCAA) resultatprova <- svytotal(~C_cat$vio_c_estatafectat, disseny_cat_opcio_c) (resultatprova / 808)* 20736963 ``` #### Opció A, Opció B i Opció C: activitat remunerada (Espanya) ```{r} #Opció A Espanya si realitza act. remunerada AB_esp <- AB_esp |> mutate(vio_a_esp_tact_remu = if_else(vio_a_esp == 1 & (M0P0 == 1 | M0P0 == 2), 1, 0)) |> mutate(vio_a_esp_tact_remu = ifelse(is.na(vio_a_esp_tact_remu), 0, vio_a_esp_tact_remu)) # Càlcul de dones per a l'opció A per a total d'Espanya vio_a_esp_tact_remu <- svytotal(~vio_a_esp_tact_remu, disseny_ab_esp) (vio_a_esp_tact_remu / 9568) * 20736963 #Opció B Espanya si realitza act. remunerada AB_esp <- AB_esp |> mutate(vio_b_esp_tact_remu = if_else(vio_b_esp == 1 & (M0P0 == 1 | M0P0 == 2), 1, 0)) |> mutate(vio_b_esp_tact_remu = ifelse(is.na(vio_b_esp_tact_remu), 0, vio_b_esp_tact_remu)) # Càlcul de dones per a l'opció B per a total d'Espanya vio_b_esp_tact_remu <- svytotal(~AB_esp$vio_b_esp_tact_remu, disseny_ab_esp) (vio_b_esp_tact_remu / 9568) * 20736963 #Opció C Espanya si realitza act. remunerada C_ESP <- C_ESP |> mutate(vio_c_esp_tact_remu = if_else(vio_c_esp == 1 & (M0P0 == 1 | M0P0 == 2), 1, 0)) |> mutate(vio_c_esp_tact_remu = ifelse(is.na(vio_c_esp_tact_remu), 0, vio_c_esp_tact_remu)) # Càlcul de dones per a l'opció C per a total d'Espanya vio_c_esp_tact_remu <- svytotal(~C_ESP$vio_c_esp_tact_remu, disseny_c_esp) (vio_c_esp_tact_remu / 9568) * 20736963 # Els results estan sobrestimats respecte els nombres que surten a la taula (Opció A, B i C per a la fila realitza activitat remunerada) ``` #### Opció A, Opció B i Opció C: activitat remunerada (Catalunya) ```{r} #Opció A Catalunya si realitza act. remunerada AB_cat <- AB_cat |> mutate(vio_a_cat_act_remu = if_else(vio_a_cat == 1 & (M0P0 == 1 | M0P0 == 2), 1, 0)) |> mutate(vio_a_cat_act_remu = ifelse(is.na(vio_a_cat_act_remu), 0, vio_a_cat_act_remu)) # Càlcul de dones per a l'opció A per a total de Catalunya vio_a_cat_act_remu <- svytotal(~AB_cat$vio_a_cat_act_remu, disseny_cat_opcio_a_b) (vio_a_cat_act_remu / 808) * 3358100 #Opció B Catalunya si realitza act. remunerada AB_cat <- AB_cat |> mutate(vio_b_cat_act_remu = if_else(vio_b_cat == 1 & (M0P0 == 1 | M0P0 == 2), 1, 0)) |> mutate(vio_b_cat_act_remu = ifelse(is.na(vio_b_cat_act_remu), 0, vio_b_cat_act_remu)) # Càlcul de dones per a l'opció B per a total de Catalunya vio_b_cat_act_remu <- svytotal(~AB_cat$vio_b_cat_act_remu, disseny_cat_opcio_a_b) (vio_b_cat_act_remu / 808) * 3358100 #Opció C Catalunya si realitza act. remunerada C_cat <- C_cat |> mutate(vio_c_cat_tact_remu = if_else(vio_c_estatafectat == 1 & (M0P0 == 1 | M0P0 == 2), 1, 0)) |> mutate(vio_c_cat_tact_remu = ifelse(is.na(vio_c_cat_tact_remu), 0, vio_c_cat_tact_remu)) # Càlcul de dones per a l'opció C per a total de Catalunya vio_c_cat_tact_remu <- svytotal(~C_cat$vio_c_cat_tact_remu, disseny_cat_opcio_c) (vio_c_cat_tact_remu / 808) * 3358100 ```
Editor is loading...
Leave a Comment