help fixing my R code script (copulas)

This project was successfully completed by YyvTao for $600 USD in 10 days.

Get free quotes for a project like this
Project Budget
$250 - $750 USD
Completed In
10 days
Total Bids
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.

Completed by:

Looking to make some money?

  • Set your budget and the timeframe
  • Outline your proposal
  • Get paid for your work

Hire Freelancers who also bid on this project

    • Forbes
    • The New York Times
    • Time
    • Wall Street Journal
    • Times Online