GWAMA(1)

✨ NAME ✨

Okay bestie, meet GWAMA 🧚 It's short for Genome-Wide Association Meta Analysis, and it's the tool you reach for when you're combining GWA study results πŸ’–

πŸ“œ SYNOPSIS πŸ“œ

GWAMA [options]

πŸ¦‹ DESCRIPTION πŸ¦‹

Okay so GWAMA (Genome-Wide Association Meta Analysis) does meta-analysis on the results of GWA studies, whether your phenotype is binary or quantitative πŸ’œ It runs both fixed-effect and random-effect meta-analyses, and it handles directly genotyped SNPs and imputed ones πŸ’…

For binary traits, it works with the allelic odds ratio and its 95% confidence interval πŸŽ€ For quantitative phenotypes, you give it the allelic effect size and standard error instead πŸ“

It's giving versatile, because you can analyse results from every genetic model: multiplicative, additive, dominant, and recessive 🌈

It also has error trapping, so it can catch strand alignment errors and allele flipping 🚨 And it runs tests of heterogeneity of effects between studies, so you can see how much the studies actually agree πŸ” No cap, it's a whole detective πŸ•΅οΈβ€β™€οΈ

πŸŽ€ OPTIONS πŸŽ€

  • -i, --filelist=filename β€” Point it at your studies' result files πŸ“‚ Default = gwama.in
  • -o, --output=fileroot β€” Set the file root for your analysis output 🌸 Default = gwama (gwama.out, gwama.gc.out)
  • -r, --random β€” Turns on random effect correction πŸ¦‹ Default = disabled
  • -gc, --genomic_control β€” Uses genomic control to adjust your studies' result files 🧹 Default = disabled
  • -gco, --genomic_control_output β€” Applies genomic control to the meta-analysis summary, so the results of the meta-analysis get corrected for gc ✨ Default = disabled
  • -qt, --quantitative β€” Switches to the quantitative trait version, which uses the BETA and SE columns πŸ“ Default = binary trait
  • -m, --map=filename β€” Picks the file name for your marker map πŸ—ΊοΈ
  • -t, --threshold=0-1 β€” The p-value threshold for showing direction in the summary effect directions πŸ”₯ Default = 1
  • --no_alleles β€” Tells it no allele information was given, so it expects the EA to always be the same πŸ™ˆ
  • --indel_alleles β€” Allele labels can be more than one letter long 🧬 No strand checks happen here, so no peeking at strands πŸ™…β€β™€οΈ
  • --sex β€” Runs the gender-differentiated and gender-heterogeneity analysis, using the method from the Magi, Lindgren & Morris 2010 paper πŸ‘‘ Your gender info has to live in the filelist file, in the second column after the file names, and it's either M or F πŸ’–
  • --filter β€” Filters based on a column header 🎯 You need 3 arguments: column name, equation [>,<,>=,<=,==,!=], and filter value. You can stack multiple filters, periodt πŸ’…
  • --name_marker β€” Alternative header for the marker name column 🏷️
  • --name_strand β€” Alternative header for the strand column 🧭
  • --name_n β€” Alternative header for the sample size column πŸ”’
  • --name_ea β€” Alternative header for the effect allele column πŸ’œ
  • --name_nea β€” Alternative header for the non-effect allele column 🩷
  • --name_eaf β€” Alternative header for the effect allele frequency column πŸ“Š
  • --name_beta β€” Alternative header for the beta column 🌈
  • --name_se β€” Alternative header for the std. err. column πŸ“
  • --name_or β€” Alternative header for the OR column ⭐
  • --name_or_95l β€” Alternative header for the OR 95L column ⬇️
  • --name_or_95u β€” Alternative header for the OR 95U column ⬆️
  • -h, --help β€” Shows you this help πŸ†˜
  • -v, --version β€” Prints the GWAMA version number πŸ”’ Honestly, bestie, this manual page is probably not enough for you 😭 For a quick reminder of gwama's options, use the help option, or head to its home page for the online documentation or the downloadable manual 🌸

πŸ”— SEE ALSO πŸ”—

http://www.geenivaramu.ee/en/tools/gwama πŸ¦‹

🫢 AUTHOR 🫢

GWAMA was written by Reedik Magi and Andrew P Morris πŸ‘‘


Version 2.2.2, August 2016, GWAMA(1) 🌈 That's the whole tea, bestie, now go run your meta-analysis πŸ’…βœ¨