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

GWAS finishing effort - 691 missing variants #18

Closed
yk-tanigawa opened this issue Jul 4, 2020 · 3 comments
Closed

GWAS finishing effort - 691 missing variants #18

yk-tanigawa opened this issue Jul 4, 2020 · 3 comments
Assignees

Comments

@yk-tanigawa
Copy link
Contributor

It turned out that the summary statistics generated in array-combined/gwas/current does not match the number of expected lines (1,080,969).

There were 407 files with 1080278 lines (meaning that 691 lines are missing in each file).

@yk-tanigawa
Copy link
Contributor Author

We identified that of those 691 variants,

  • 369 variants are on both arrays
  • 322 variants are on one array

and performed association analysis.

For some phenotypes, we saw

Skipping ... since all samples are controls

@yk-tanigawa yk-tanigawa self-assigned this Jul 4, 2020
@yk-tanigawa
Copy link
Contributor Author

Summary of counts

  • We applied association analysis for 407 (pop, GBE_ID)
    • 180 have some results
      • 176 are finished with full results (691 variants)
      • 4 have 322 variants (369 variants are missing)
        • # samples <= # predictor columns
    • 227 were skipped
      • 2: all cases
      • 196: all controls
      • 17: constant quantitative phe
      • 12: # samples <= # predictor columns

@yk-tanigawa
Copy link
Contributor Author

180 files were merged.

https://github.com/rivas-lab/ukbb-tools/tree/master/04_gwas/extras/202006-GWAS-finish/gwas691

From this analysis, we know that

  • 4 files have 1,080,599 lines
    • 4: # samples <= # predictor columns
  • 227 files have 1,080,278 lines
    • 2: all cases
    • 196: all controls
    • 17: constant quantitative phe
    • 12: # samples <= # predictor columns

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

1 participant