Unifying an inconsistent coded categorical variable in a panel/longtitudal dataset.
There is offered the novel cat2cat procedure to map a categorical variable according to a mapping (transition) table between two different time points.
The mapping (transition) table should to have a candidate for each category from the targeted for an update period. The main rule is to replicate the observation if it could be assigned to a few categories, then using simple frequencies or modern statistical methods to approximate probabilities of being assigned to each of them.
This algorithm was invented and implemented in the paper by (Nasinski, Majchrowska and Broniatowska (2020)).
For more details please read the paper by (Nasinski, Gajowniczek (2023)).
Please visit the cat2cat webpage for more information
# install.packages("remotes")
remotes::install_github("polkas/cat2cat")
# or
install.packages("cat2cat")occup dataset is an example of unbalance panel dataset.
This is a simulated data although there are applied a real world characteristics from national statistical office survey.
The original survey is anonymous and take place every two years.
trans dataset containing mappings (transitions) between old (2008) and new (2010) occupational codes.
This table could be used to map encodings in both directions.
Panel dataset without the unique identifiers and only two periods, backward and simple frequencies:
library("cat2cat")
data("occup", package = "cat2cat")
data("trans", package = "cat2cat")
occup_old <- occup[occup$year == 2008, ]
occup_new <- occup[occup$year == 2010, ]
occup_simple <- cat2cat(
data = list(
old = occup_old, new = occup_new,
cat_var_old = "code", cat_var_new = "code", time_var = "year"
),
mappings = list(trans = trans, direction = "backward")
)Panel dataset without the unique identifiers and four periods, backward direction and ml models:
library("cat2cat")
data("occup", package = "cat2cat")
data("trans", package = "cat2cat")
occup_2006 <- occup[occup$year == 2006,]
occup_2008 <- occup[occup$year == 2008,]
occup_2010 <- occup[occup$year == 2010,]
occup_2012 <- occup[occup$year == 2012,]
library("caret")
ml_setup <- list(
data = occup_2010,
cat_var = "code",
method = c("knn"),
features = c("age", "sex", "edu", "exp", "parttime", "salary"),
args = list(k = 10, ntree = 50)
)
mappings <- list(trans = trans, direction = "backward")
# ml model performance check
print(cat2cat_ml_run(mappings, ml_setup))
# from 2010 to 2008
occup_back_2008_2010 <- cat2cat(
data = list(
old = occup_2008, new = occup_2010,
cat_var_old = "code", cat_var_new = "code", time_var = "year"
),
mappings = mappings,
ml = ml_setup
)
# from 2008 to 2006
occup_back_2006_2008 <- cat2cat(
data = list(
old = occup_2006, new = occup_back_2008_2010$old,
cat_var_new = "g_new_c2c", cat_var_old = "code", time_var = "year"
),
mappings = mappings,
ml = ml_setup
)
o_2006_new <- occup_back_2006_2008$old
o_2008_new <- occup_back_2008_2010$old # or occup_back_2006_2008$new
o_2010_new <- occup_back_2008_2010$new
o_2012_new <- dummy_c2c(
occup_2012, cat_var = "code", ml = c("knn")
)
final_data_back <- do.call(
rbind,
list(o_2006_new, o_2008_new, o_2010_new, o_2012_new)
)
# possible processing, leaving only one obs per subject and period
# still it is recommended to leave all replications and use the weights in the statistical models
library(magrittr)
ff <- final_data_back %>%
split(.$year) %>%
lapply(function(x) cross_c2c(x)) %>%
lapply(function(x)
prune_c2c(x, column = "wei_cross_c2c", method = "highest1")
) %>%
do.call(rbind, .)
all.equal(nrow(ff), sum(ff$wei_cross_c2c))
all.equal(nrow(ff), sum(final_data_back$wei_freq_c2c))More complex examples are presented in the "Get Started" vignette.