Description Usage Arguments Details Value Author(s) See Also Examples
The model has a correlation parameter which is estimated and theoretically bound between -1 and +1. To ensure that the estimated parameters are within the theoretical bounds a transformation is necessary. The chosen transformation is:
f(ρ): ρ = 2/(1-exp(-θ))- 1
Whereas ρ is the actual correlation coefficient and θ is the parameter we estimate in the model. This parametrization has been worked into the likelihood function and ensures that ρ will be between -1 and +1.
The variance covariance matrix thus contains entries based on θ but not ρ. Hence, this function takes the variance of the transformed correlation parameter (θ) and produces the value correct for ρ.
To create the correct measure of var(ρ) this function simulates 1,000 θ's and then transforms them to ρ's. The variance of these ρ's is then reported. Note, this means that the variance-covariance returned by StratSel
is only correct for all diagonals and off-diagonals for the parameters (β) but for the correlation coefficient only the variance is correct. Given that there is no reason to use the full variance-covariance for post-estimation commands this is not a problem.
1 | fetch.rho.v(v, b)
|
v |
Variance-covariance matrix based on the regular parameters (β) and θ. |
b |
Coefficient vector, first n elements are for the regular parameter estimates and the last element is θ. |
This function is for internal use but documented as a regular function to enable any user to assess the estimator and its functionality.
Returns the correct variance estimate for the estimate of the correlation coefficient ρ.
Lucas Leemann lleemann@gmail.com
1 2 | fetch.rho.v(matrix(c(1,0,0,1),2,2),c(0,0))
fetch.rho.v(matrix(c(1,0,0,2),2,2),c(0,0))
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.