help fixing my R code script (copulas)

Avg Bid (USD)
Project Budget (USD)
$250 - $750

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.

Skills required:
Algorithm, Engineering, Mathematics, Matlab & Mathematica, Statistics
About the employer:
Public Clarification Board
Bids are hidden by the project creator. Log in as the employer to view bids or to bid on this project.
You will not be able to bid on this project if you are not qualified in one of the job categories. To see your qualifications click here.

$ 750
in 5 days
$ 750
in 15 days
$ 500
in 4 days
Hire YyvTao
$ 600
in 10 days
Hire JepBo
$ 400
in 5 days