svyrich missingness handling

[ajdamico]

library(survey)
library(convey)
library(laeken)

data(eusilc) ; names( eusilc ) <- tolower( names( eusilc ) )

# add a missing
eusilc[ 1 , 'eqincome' ] <- NA

# linearized design

des_eusilc <- svydesign( ids = ~rb030 , strata = ~db040 ,  weights = ~rb050 , data = eusilc )
des_eusilc <- convey_prep( des_eusilc )

# works
svyrich(~eqincome, des_eusilc, type_measure = "FGTT1" , g=0,  abs_thresh=30000,na.rm=T)
svyrich(~eqincome, des_eusilc, type_measure = "FGTT1" , g=1,  abs_thresh=30000,na.rm=T)

# incorrectly errors
svyrich(~eqincome, des_eusilc, type_measure = "FGTT1" , g=0, type_thresh= "relq" ,na.rm=T)
svyrich(~eqincome, des_eusilc, type_measure = "FGTT1" , g=1, type_thresh= "relq" ,na.rm=T)

# incorrectly returns NAs
svyrich(~eqincome, des_eusilc, type_measure = "FGTT1" , g=0, type_thresh= "relm" ,na.rm=T)
svyrich(~eqincome, des_eusilc, type_measure = "FGTT1" , g=1, type_thresh= "relm" ,na.rm=T)