Introduction

Table 3 adds the option of nonresponse since the composition of the sampling frames differ across incentive groups.

Education

# Table with absolute frequencies for response
edu_rr <- WaveOne %>%
  filter(is.na(PUMA1)==FALSE) %>% 
  filter(valid==1) %>%
  with(table(edu,Pumavers)) 


# Table with absolute frequencies for nonresponse
edu_nr <- WaveOne %>%
  filter(is.na(PUMA1)==FALSE) %>% 
  filter(valid==0) %>%
  with(table(edu,Pumavers))

# Combine tables
for (i in 1:ncol(edu_rr)){
  sum_perc <- cbind(edu_nr[,i],edu_rr[,i]) %>%
      prop.table(1)
  if (i==1){
    edu_tab <- data.frame(
     bro_nr = paste0(edu_nr[,i],
                         " (", 
                         scales::percent(sum_perc[,1], accuracy = 0.1), 
                         ")"),
     bro_rr = paste0(edu_rr[,i],
                         " (", 
                         scales::percent(sum_perc[,2], accuracy = 0.1), 
                         ")")
    )
  } else {
    j <- (i*2)-1
   edu_tab[,j] <- paste0(edu_nr[,i],
                         " (", 
                         scales::percent(sum_perc[,1], accuracy = 0.1), 
                         ")")
   edu_tab[,j+1] <- paste0(edu_rr[,i],
                         " (", 
                         scales::percent(sum_perc[,2], accuracy = 0.1), 
                         ")")
  }
}

# Print dataframe
cbind(levels(WaveOne$edu),
      edu_tab) %>%
  setNames(c("Labels",
             rep(c("Nonresponse","Response"),4))) %>%
  kable("html") %>%
  kable_styling("striped") %>%
  column_spec(1, width = "30em") %>%
  add_header_above(c(" " = 1, 
                     "Brochure" = 2, 
                     "2€ token" = 2, 
                     "5€ token" = 2,
                     "Voucher" = 2))

Occupational position

# Table with absolute frequencies for response
job_rr <- WaveOne %>%
  filter(is.na(PUMA1)==FALSE) %>% 
  filter(valid==1) %>%
  with(table(job,Pumavers)) 


# Table with absolute frequencies for nonresponse
job_nr <- WaveOne %>%
  filter(is.na(PUMA1)==FALSE) %>% 
  filter(valid==0) %>%
  with(table(job,Pumavers))

# Combine tables
for (i in 1:ncol(job_rr)){
  sum_perc <- cbind(job_nr[,i],job_rr[,i]) %>%
      prop.table(1)
  if (i==1){
    job_tab <- data.frame(
     bro_nr = paste0(job_nr[,i],
                         " (", 
                         scales::percent(sum_perc[,1], accuracy = 0.1), 
                         ")"),
     bro_rr = paste0(job_rr[,i],
                         " (", 
                         scales::percent(sum_perc[,2], accuracy = 0.1), 
                         ")")
    )
  } else {
    j <- (i*2)-1
   job_tab[,j] <- paste0(job_nr[,i],
                         " (", 
                         scales::percent(sum_perc[,1], accuracy = 0.1), 
                         ")")
   job_tab[,j+1] <- paste0(job_rr[,i],
                         " (", 
                         scales::percent(sum_perc[,2], accuracy = 0.1), 
                         ")")
  }
}

# Print dataframe
cbind(levels(WaveOne$job),
      job_tab) %>%
  setNames(c("Labels",
             rep(c("Nonresponse","Response"),4))) %>%
  kable("html") %>%
  kable_styling("striped") %>%
  column_spec(1, width = "30em") %>%
  add_header_above(c(" " = 1, 
                     "Brochure" = 2, 
                     "2€ token" = 2, 
                     "5€ token" = 2,
                     "Voucher" = 2))