help fixing my R code script (copulas)

Project Description:
hello. i have some R code that estimates an index popular in the quantitative social sciences called the polychoric correlation coefficient. the difference is that instead of assuming a multivariate normal distribution an elliptical copula distribution is assumed as a model.

the R script works fine, i generate good enough estimates but i cannot, for the life of me, get appropriate standard errors (SE). there are functions in R that estimate this index via maximum likelihood so i have been using these as comparison to evaluate my results. they should be the same (within reasonable rounding error of course) but mines are far, far off.

your job would be to look at my R script (which otherwise works perfectly) and help me fix it to obtain standard errors that match the ones from the already-developed functions in the 'polycor' R package. i know it is impossible for them to be the exact same, but they should be within a reasonable range of which other. ideally, the ones from the copula model should be smaller

PS- as a little bit of a guideline, this is a little more complicated than just grabing the 'optim' function, setting hessian=T and taking the square root of the inverse of the diagonal. i'm wondering if something is wrong with my likelihood function (although i hope not because the estimates i get match the ones i would expect)

i am comfortable around statistics and R but this is definitely escaping my abilities a little bit.