Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Does the scale of one outcome affect the bivariate correlation? #1

Open
andkov opened this issue Oct 4, 2017 · 2 comments
Open

Does the scale of one outcome affect the bivariate correlation? #1

andkov opened this issue Oct 4, 2017 · 2 comments
Assignees

Comments

@andkov
Copy link
Member

andkov commented Oct 4, 2017

Fev vs Fev100

Question: will a bivariate model yield the same inference regarding bivariate linear correlation if one of the processes is linearly transformed?

Observations so far by @ampiccinin
Based on the covariances, the fev100 had the desired effect - we can now see slope covariance values larger than 0.00, and the p values become consistently lower - but they vary by how much

@ampiccinin
Copy link
Member

While p-values are probably not the best way to compare the results, they are generally what people rely on in drawing conclusions about their data. With this in mind:

Of, for example, the 18 measures in one study:

  • 2 have identical p-values both for FEV and FEV100, and based on covariance and correlation. These also have identical (Mplus estimated) correlations, but neither the SEs nor the covariances are exactly the same.

  • 1 has n.s. association according to FEV, but "statistically significant" for FEV100 (for cov and corr). Identical estimated correlations, but much smaller SE (.53 vs .14)

  • 3 have n.s. association according to FEV and FEV100cov, but "statistically significant" for FEV100corr. Again, with identical estimated correlations, but much smaller SE (~.8 vs .2)

  • the rest (12) have slightly smaller SE and p-value for FEV100, but are essentially unchanged between FEV and FEV100.

An additional issue is that the computed correlation CIs (I don't see the SEs in the table) are a lot smaller than the Mplus esimated CIs (e.g., -.94, .96 vs -.05, .07)

Highlighted in the attached file (green, yellow, red, grey/white, and orange, respectively):
pulmonary-meta-fev-vs-fev100-MAP-2017-10-03 - AMP.xlsx

@andkov
Copy link
Member Author

andkov commented Oct 26, 2017

@ampiccinin
A few notes.

  1. The computed correlations should not have SE. due to the fact that correlation coefficient is not distributed normally, you can’t just throw a 95%CI from Gaussian onto a point estimate and call it a 95%CI of a correlation.
  2. The Mplus estimate this by just throwing a 95CI from Gaussian onto a point estimate and calling this a 95%CI of a correlation.
  3. It seems that Mplus, by sssuming that point estimate of a correlation is distributed normally, overestimate its value, compared to Fisher-Z transform employed during our “in-house” computation of the correlation point estimate.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants