help fixing my R code script (copulas)

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: Algorithm, Engineering, Mathematics, Matlab and Mathematica, Statistics

See more: copula code, matlab works, match statistics, match results statistics, life works, job se, job in mines, index bit, help me fix this, help me fix it, help fix it, fix it help, fixing my, code that, code 42, c&m engineering, bit index, get code, statistics R, statistics matlab

About the Employer:
( 12 reviews ) Vancouver, Canada

Project ID: #1654093

Awarded to:


I did my statistics dissertation computation work in R. I think I can solve your interesting problem.

$600 USD in 10 days
(2 Reviews)

4 freelancers are bidding on average $563 for this job


Hi We can help you Best regards

$750 USD in 5 days
(43 Reviews)

I have very strong hands on R and matlab. I am sure you will be happy to have me working you, and might look forward to re-hire me like many others. Please see pm.

$750 USD in 15 days
(13 Reviews)

sir, please check PMB.... thanks

$500 USD in 4 days
(12 Reviews)

Hi! The first thing I'd try would be to test your functions with normal data in place of elliptical. Thus we could know if the likelihood is good.

$400 USD in 5 days
(0 Reviews)