R scripts to estimate power in a bivariate LGCM