ldsc | LD Score Regression | Analytics library
kandi X-RAY | ldsc Summary
kandi X-RAY | ldsc Summary
ldsc is a command line tool for estimating heritability and genetic correlation from GWAS summary statistics. ldsc also computes LD Scores.
Support
Quality
Security
License
Reuse
Top functions reviewed by kandi - BETA
- Annotate a list of files
- Parse annotation file
- Return a list of chromosomes present in fh
- Parse MF header file
- Estimate H2
- Print the covariance matrix
- Checks the LD score matrix condition
- Calculate the overlap output
- Estimate genotype parameters
- Check if arguments are the same length
- Make annotation files
- Combine two joint objective functions
- Concatenate ld score from a list of files
- Make a bedtool for a gene set
- Read cts file
- Compute the weight matrix for the fitting function
- Calculate the weights for the model
- Convert an array of block values to the largest block
- Print a summary of the SNPs
- Delete values from x and y
- Summarize the correlation coefficient
- Compute the delete values for each block
- Compute the block - wise quadrature values
- Return a summary of the model
- Compute cell type specific
- Updates the weights for the model
ldsc Key Features
ldsc Examples and Code Snippets
Community Discussions
Trending Discussions on ldsc
QUESTION
I have a functioning Snakefile for Partitioned Heritability using multiple bed files. This produces a perfect list of jobs using snakemake -np
so this file only needs a minor tweak (I hope!).
My issue occurs in the merge_peaks
rule below.
At this stage, I have 25 bed files and need to run 5 calls of merge_peaks
, one call for each extension ext=[100,200,300,400,500]
, so I need only the 5 bed files containing the relevant extension to be called each time.
For example for the following merge_peaks
output file peak_files/Fullard2018_peaks.mrgd.blrm.100.bed
, I only need the following 5 bed files with ext=100
to be used as input:
ANSWER
Answered 2018-Dec-19 at 17:17I have a solution for this using double curly brackets as an escape character. This works and gives me what I need, but there is still something going on that I don't understand.
So for the input line to the merge_peaks
rule instead of:
Community Discussions, Code Snippets contain sources that include Stack Exchange Network
Vulnerabilities
No vulnerabilities reported
Install ldsc
Support
Reuse Trending Solutions
Find, review, and download reusable Libraries, Code Snippets, Cloud APIs from over 650 million Knowledge Items
Find more librariesStay Updated
Subscribe to our newsletter for trending solutions and developer bootcamps
Share this Page