Traduciendo Stata a R y teniendo problemas con is.na.

Estoy traduciendo código de Stata a código de R, pero ahora tengo algunos problemas de novato como este.

Este es mi código de Stata:

gen     aposentadofam=1 if proprendaposent  > 0 &amp; proprendaposent ~=.;

replace aposentadofam=0 if proprendaposent == 0 | proprendaposent ==.;

Y esto es lo que intenté hacer en R:

# pemg &lt;- mutate(pemg, aposentadofam = NA_real_)

  # pemg &lt;- mutate(pemg, aposentadofam = case_when(proprendaposent >0 &amp; !is.na(proprendaposent) ~ 1, TRUE ~ aposentadofam))

  # pemg &lt;- mutate(pemg, aposentadofam = case_when(proprendaposent==0 | is.na(proprendaposent) ~ 0, TRUE ~ aposentadofam))

La línea con is.na() parece estar funcionando correctamente, pero la que tiene !is.na() no. Me da este mensaje de error:

LHS of case 1 (`proprendaposent > 0 &amp; !is.na(proprendaposent) ~ 1`) must be a logical vector, not a `formula` object.

¿Qué debo hacer?

davy.ai

June 18, 2023 at 8:08 am

El problema con el código R es que la función case_when() espera un vector lógico como el lado izquierdo de cada caso, pero la expresión proprendaposent > 0 & !is.na(proprendaposent) es un objeto formula.

Para solucionar esto, puede usar la función if_else(), que toma una condición lógica como primer argumento y evalúa a un valor si la condición es “VERDADERA” y otro valor si es “FALSO”. Así es como podría traducir el código Stata a R usando if_else():

pemg <- mutate(pemg, aposentadofam = if_else(proprendaposent > 0 & !is.na(proprendaposent), 1, 
                                              if_else(proprendaposent == 0 | is.na(proprendaposent), 0, NA_real_)))

Este código establece la variable aposentadofam en 1 si proprendaposent es mayor que 0 y no está perdido (!is.na()), 0 si es igual a 0 o está perdido (is.na()), y NA de lo contrario.

Answer

davy.ai

June 18, 2023 at 8:08 am
El problema con el código R es que la función case_when() espera un vector lógico como el lado izquierdo de cada caso, pero la expresión proprendaposent > 0 & !is.na(proprendaposent) es un objeto formula.

Para solucionar esto, puede usar la función if_else(), que toma una condición lógica como primer argumento y evalúa a un valor si la condición es “VERDADERA” y otro valor si es “FALSO”. Así es como podría traducir el código Stata a R usando if_else():
```
pemg <- mutate(pemg, aposentadofam = if_else(proprendaposent > 0 & !is.na(proprendaposent), 1, 
                                              if_else(proprendaposent == 0 | is.na(proprendaposent), 0, NA_real_)))
```
Este código establece la variable aposentadofam en 1 si proprendaposent es mayor que 0 y no está perdido (!is.na()), 0 si es igual a 0 o está perdido (is.na()), y NA de lo contrario.

Comments are closed.

es.davy.ai

¿Tienes una pregunta?

Traduciendo Stata a R y teniendo problemas con is.na.

Answer

davy.ai