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.

Compétences : Algorithme, Ingénierie, Mathématiques, Matlab and Mathematica, Statistiques

en voir plus : 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, sciences, r], R&R, r package, quantitative c

Concernant l'employeur :
( 12 commentaires ) Vancouver, Canada

Nº du projet : #1654093

Décerné à:


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

%selectedBids___i_sum_sub_7% %project_currencyDetails_sign_sub_8% USD en 10 jours
(2 Commentaires)

4 freelance font une offre moyenne de $563 pour ce travail


Hi We can help you Best regards

%bids___i_sum_sub_35% %project_currencyDetails_sign_sub_36% USD en 5 jours
(43 Commentaires)

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.

%bids___i_sum_sub_35% %project_currencyDetails_sign_sub_36% USD en 15 jours
(13 Commentaires)

sir, please check PMB.... thanks

%bids___i_sum_sub_35% %project_currencyDetails_sign_sub_36% USD en 4 jours
(12 Commentaires)

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.

%bids___i_sum_sub_35% %project_currencyDetails_sign_sub_36% USD en 5 jours
(0 Commentaires)