set.seed(13)
temps_obs <- 5 # panel size
n_subj <- 50 # sample size
id <- rep(1:n_subj, each = temps_obs)
asym <- 0.5
panSizeVec <- unname(unlist(lapply(split(id, id), function(x) length(x))))
yvec <- c(mvtnorm::rmvnorm(n_subj, sigma = diag(rep(1,temps_obs))))
aweight <- rep(asym, temps_obs * n_subj)
aweight[!(yvec > mean(yvec))] = 1 - asym
dexpectilizeVecR(yvec, aweight, panSizeVec)
Run the code above in your browser using DataLab