Skip to contents

Build the wide longitudinal file from per-wave surveys

Source: R/panel-design.R
panel_merge.Rd

Joins a named list of per-wave data.frames on a stable unit key into the wide, one-row-per-unit file that the longitudinal recipe consumes. Each wave's non-key columns are suffixed with the wave name, and per-wave presence indicators .wf_in_<wave> (and, if responded is given, response indicators .wf_resp_<wave>) are added so a later step_attrition() can model the response pattern.

Usage

panel_merge(
  waves,
  by,
  responded = NULL,
  require = c("any", "all"),
  suffix = "_"
)

Arguments

waves

a named list of per-wave data.frames; the names become the wave labels and column suffixes.

by

one or more column names (a character vector) forming the unit key, present in every wave – e.g. "ID" for a household or c("ID", "nper") for a person.

responded

optional string: the name of a response indicator present in every wave, used to build .wf_resp_<wave>.

require

one of "any" (union of units, default) or "all" (intersection).

suffix

string inserted before the wave label when renaming columns (default "_"), e.g. edad in wave T1 becomes edad_T1.

Value

a wide data.frame, one row per unit, ready for panel_design() / weighting_spec().

Details

require = "any" (the default) keeps every unit seen in at least one wave – which the attrition model needs, because it has to be fitted over responders and non-responders. require = "all" keeps only the intersection (CEPAL's s(2) = s1 cap s2 cap ...); prefer to reach it with step_attrition() after estimating, not by dropping units before the model can see them.

See also

panel_design()

Examples

# reshape two waves of the long panel into one wide, one-row-per-unit file
wide <- panel_merge(
  list(T1 = subset(panel_ine, ola == 1), T2 = subset(panel_ine, ola == 2)),
  by = c("id_hogar", "nper"), require = "any")
names(wide)                       # per-wave columns are suffixed _T1 / _T2
#>  [1] "id_hogar"          "nper"              "id_persona_T1"    
#>  [4] "estrato_T1"        "psu_T1"            "region_T1"        
#>  [7] "sexo_T1"           "edad_T1"           "ola_T1"           
#> [10] "grupo_rotacion_T1" "mes_en_muestra_T1" "w_base_T1"        
#> [13] "disp_T1"           "condicion_T1"      "ocupado_T1"       
#> [16] "desocupado_T1"     "ingreso_T1"        ".wf_in_T1"        
#> [19] "id_persona_T2"     "estrato_T2"        "psu_T2"           
#> [22] "region_T2"         "sexo_T2"           "edad_T2"          
#> [25] "ola_T2"            "grupo_rotacion_T2" "mes_en_muestra_T2"
#> [28] "w_base_T2"         "disp_T2"           "condicion_T2"     
#> [31] "ocupado_T2"        "desocupado_T2"     "ingreso_T2"       
#> [34] ".wf_in_T2"