Testing homogeneity with homoStat!

Hi everyone,
I conducted a two-stage metaSEM with fixed effect model, since the random effect model reported estimation problem (after more then 5 hours running the first stage it never gave an output!).

Anyway, I am trying to compute the homogeneity statistics by using the homoStat function:
**homoStat(y, v)**
where, **y**: A vector of effect size for univariate meta-analysis or a k x p matrix of effect sizes for multivariate meta-analysis where k is the number of studies and p is the number of effect sizes.
**v**: A vector of the sampling variance of the effect size for univariate meta-analysis or a k x p* matrix of the sampling covariance matrix of the effect sizes for multivariate meta-analysis where p* = p(p+1)/2. It is arranged by column major as used by vech. It is assumed that there is no missing value in v if y is complete. If there are missing values in v due to the missingness on y, the missing values in v will be removed automatically.

I managed to compute the y matrix, but I am having difficulties computing the covariance matrix of the effect sizes for multivariate meta-analysis.
The coding is: **sampcov <- vech(cov(y))**
The output is a matrix of 1 column and 78 rows, and when I put it in the homoStat function I get this error message: *Error in homoStat(y, sampcov) : The expected no. of columns in v is 78 while the observed no. of columns in v is 1.*.
I tried to switch the column and rows, but then I get this other error message: *Error in V[!miss.index, !miss.index, drop = FALSE] : (subscript) logical subscript too long*.

What am I missing?
I do not understand what I did wrong.

Thank you in advance!