Giter Site home page Giter Site logo

gwaspipeline's Introduction

gwaspipeline

This repositor contains a set of bash and python scripts used for GWAS analyses mostly based on plink.

If you are looking for scripts related to the Sweden exome dataset, you are mostly interested in the following master scripts:

install.sh
makegenes.sh
genovese2016.sh
genovese2016.R

Current version was updated on October 3rd 2016.

csv2xlsx.py

Basic script to convert text tables into XLSX tables (supports multiple sheets)

genome2ped.py

Identifies trios from a plink genome file and optionally generates PI_HAT x Z0 plot

install.sh

Installs several resources useful for genetic analyses in a ~/res/ directory in your system

markerqc.sh

Perform basic quality control for markers in a plink dataset

plink2lite.sh

Extracts a core set of highly confident non-rare 1000 Genomes project variants to analyze a given dataset befor marker QC

plink2score.py

Compute several polygenic scores from either binary or dosage plink files and optionally adjusts them for covariates (e.g. principal components)

makegenes.sh

This scripts installs and generates several gene tables and lists useful for enrichment analyses

sampleqc.sh

Performs basic quality control for samples in a plink dataset

skat.py

Wrapper for performing burden analysis with SKAT (see http://www.hsph.harvard.edu/skat/) with several plink-like options

vcf2annot.sh

This script annotates a given VCF using the snpEff/SnpSift software and the dbNSFP database

vcf2plink.py

This scripts converts VCF to plink following best practices (see http://bit.ly/1sgOcuP)

gwaspipeline's People

Contributors

freeseek avatar

Stargazers

 avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar

Watchers

 avatar  avatar  avatar  avatar  avatar

gwaspipeline's Issues

Error: Missing headers in the plink assoc file

Hi!
I get this error but I think my input is OK.
Three first lines of the file pasted:

CHR SNP BP A1 F_A F_U A2 P OR
1 1:565433:C:T 565433 T 0 0 C 1 NA
1 1:567667:C:T 567667 T 0 0 C 1 NA
1 1:568208:T:C 568208 C 0 0.5 T 1 0

The assoc file cames from kgp2anc.
What is the problem here in the code?

if (args[1]=='-') { df <- read.table(file('stdin'), header=TRUE) } else { df <- read.table(args[1], header=TRUE) } if(!all(names(df) %in% c('CHR', 'BP', 'P'))) { stop('Missing headers in the plink assoc file') }

Recommend Projects

  • React photo React

    A declarative, efficient, and flexible JavaScript library for building user interfaces.

  • Vue.js photo Vue.js

    ๐Ÿ–– Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.

  • Typescript photo Typescript

    TypeScript is a superset of JavaScript that compiles to clean JavaScript output.

  • TensorFlow photo TensorFlow

    An Open Source Machine Learning Framework for Everyone

  • Django photo Django

    The Web framework for perfectionists with deadlines.

  • D3 photo D3

    Bring data to life with SVG, Canvas and HTML. ๐Ÿ“Š๐Ÿ“ˆ๐ŸŽ‰

Recommend Topics

  • javascript

    JavaScript (JS) is a lightweight interpreted programming language with first-class functions.

  • web

    Some thing interesting about web. New door for the world.

  • server

    A server is a program made to process requests and deliver data to clients.

  • Machine learning

    Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.

  • Game

    Some thing interesting about game, make everyone happy.

Recommend Org

  • Facebook photo Facebook

    We are working to build community through open source technology. NB: members must have two-factor auth.

  • Microsoft photo Microsoft

    Open source projects and samples from Microsoft.

  • Google photo Google

    Google โค๏ธ Open Source for everyone.

  • D3 photo D3

    Data-Driven Documents codes.