R scripts to estimate power in bivariate LGCM