Estimating own weights in TWAS/FUSION

Estimating own weights in TWAS/FUSION - manual

Hi,

I am trying to make own weights with RNAseq and genotype data. I have created plink files by Gene (using GeneName - FTO). However, I did not find detailed doc/manual to see how to run Weights generation.

I am recovering step by step the various commands but I was wondering whether there is a manual somewhere which explains explicitly what tmp and other options stand for.

Best
Christian

twas-fusion twas

• 1,420 views • link

updated 3.0 years ago by Ram

• written 3.0 years ago by Genotepes

Google search for "twas fusion" leads me to a webpage that has a link to this manual: http://gusevlab.org/projects/fusion/#command-line-parameters

• 0 views • link

written 3.0 years ago by Ram

Thanks for sharing. However, I have this but still so far I observed that the tmp was not totally free format - in the example shell file, the tmp needs to have some keywords. I have advanced since my message and think I will solve this. Thanks again

• 0 views • link

written 3.0 years ago by Genotepes

Please add an answer and accept it once you solve the problem - that will be useful for future users.

• 0 views • link

written 3.0 years ago by Ram

0 answers

No answers yet.

More posts like this

GATK GenomicsDBImport update database error using snakemake
written 2.1 years ago by Peter Chung

I am new to use snakemake and now I am able to apply it in the GATK GenomicsDBImport steps combining 500 genotype vcf files. now …
GATK memory error with Java
written 3.0 years ago by anikcropscience •

Hello, I have genotype by sequencing data for 400 samples. I am trying to run a SNP calling pipeline using `GATK`. I could manage until …
"split chromosome" after sorting with bcftools
written 4.1 years ago by dec986

I am running bcftools in a series of commands First I run `bcftools sort` and then plink --bcf tmp.sort.bcf --allow-extra-chr --make-bed --out merged.eigenstrat.34string.length.vcf.gz which works, …
Filtering out samples by FORMAT attributes using bcftools
written 5.5 years ago by scottbrainard •

The `-e` and `-i` options of the `bcftools filter` command appear, by default, to only allow for including or excluding _sites_. E.g., `-e 'FMT/DP < …
Error in `rownames<-`(`*tmp*`, value = colnames(countData)): attempt to set 'rownames' on an object…
written 6.3 years ago by katarzyna.wieciorek •

Hello, I am using DESeq2 library following the manual [3.2 Starting from count matrices][1]. I have my countdata and coldata imported from CSV files. I …
Hisat2 error: --read-lengths arg must be at least 20
written 6.3 years ago by aindap

I am new to Hisat2 and I am trying to align paired-end rnaseq fastq files with hisat2 v 2.2.0. I downloaded the hg38 transcriptome index …
Error message whengenerating read coverage tables for Ballgown by stringtie
written 10.1 years ago by hu.guangzhen9 •

Background Differential expression analysis workflow by stringtie includes the following steps: 1. for each RNA-Seq sample, map the reads to the genome with HISAT2 using …
Material for calling in NGS - statistical view
written 10.4 years ago by Genotepes

Hi all, I am preparing a short course on statistical tests for exome analysis (association in NGS) and I'd like to discuss a bit the …
Job: Postodoctoral position (29 months) in bioinformatics Lab. of Plant Physiology and Molecular Ge…
written 12.0 years ago by chermans •

**Genetic and nutritional influences on the root architecture of Brassicaceae species** Nitrogen (N) is the quantitatively most important nutrient required by crops, however a considerable …
Whole-Exome Data Cleanup / Filtering
written 15.3 years ago by Khader Shameer

<p>What tools do you you use to clean up your whole-exome paired end reads: I am running a FastQC check for various parameters followed by …