User: Sam
Sam • 2.2k
- Reputation:
- 2,160
- Status:
- Trusted
- Location:
- London
- Last seen:
- 5 days, 22 hours ago
- Joined:
- 7 years, 1 month ago
- Email:
- c***********@gmail.com
Posts by Sam
<prev
• 340 results •
page 1 of 34 •
next >
0
votes
1
answer
186
views
1
answers
Answer:
A: Base file of PRSice software
... The base file is a GWAS association file. Usually you can obtain it from external data base (e.g. PGC). For Type 2 Diabetes, you might be able to get the summary statistics from WTCCC. Alternatively, if you have an independent data set (independent to your target data), you can also perform the GWAS ...
written 5 weeks ago by
Sam • 2.2k
1
vote
1
answer
203
views
1
answers
... The thing is, we only perform clumping on SNPs that were found in both the target and base file. By using the imputed data, you increase the coverage and might result in more post-clumped SNPs, therefore increase amount of information and possibly increase the performance of your PRS model ...
written 5 weeks ago by
Sam • 2.2k
0
votes
1
answer
530
views
1
answers
... I'd suggest you to perform a partitioning of heritability and see if majority of the heritability lies within these regulatory elements. [LDSC][1] is one of the tool that can perform that
Though the problem of this approach is that it doesn't provide you with a SNP based probability
[1]: https:/ ...
written 13 months ago by
Sam • 2.2k
0
votes
3
answers
1.2k
views
3
answers
Comment:
C: 1000 genomes PCA for whole genome
... You can read [here][1]. Though I do suggest you to follow Kevin's tutorial
[1]: https://www.cog-genomics.org/plink/1.9/input#vcf_filter ...
written 13 months ago by
Sam • 2.2k
1
vote
3
answers
1.2k
views
3
answers
Comment:
C: 1000 genomes PCA for whole genome
... As mentioned by GabrielMontenegro, you will need to do the following:
1. clean the multi-allelic sites from your files
2. remove or rename duplicated SNPs
3. perform pruning using either `--indep` or `--indep-pairwise`
4. Extract pruned SNPs with `--extract`
5. perform PCA using `--pca`
(You can ...
written 13 months ago by
Sam • 2.2k
1
vote
3
answers
1.2k
views
3
answers
... If your genotype is in vcf format, you can convert it to plink format then use the --pca function from PLINK.
If you have large number of samples (which shouldn't be the case if you are using 1000G), then you can consider [flashPCA][1]
Note: If you are directly converting the vcf to plink format, ...
written 13 months ago by
Sam • 2.2k
2
votes
1
answer
1.2k
views
1
answers
Answer:
A: Undrestanding BGEN files
... You can find the detail information of the bgen format here:
http://www.well.ox.ac.uk/~gav/bgen_format/
Basically, this is a type of dosage data.
Unfortunately, the bgen format is in binary. So you will either have to write your own programme or use some external programe to handle them. If all yo ...
written 17 months ago by
Sam • 2.2k
0
votes
0
answers
1.3k
views
0
answers
... Strange. PLINK does not produce empty lines filtered. Maybe you should try and add --geno 0.1 in your plink command to remove SNPs with high missingness? ...
written 20 months ago by
Sam • 2.2k
0
votes
0
answers
1.3k
views
0
answers
... Try to do the following code
sed 's/^chrM\s/25\t/g; s/^chrX\s/23\t/g; s/^chrY\s/24\t/g; s/^chr//g' your.bim > fixed.bim
Should help you to modify the bim file. As for your error, it is likely that your replace has changed the number of line of bim, leading to the problem ...
written 20 months ago by
Sam • 2.2k
1
vote
0
answers
1.3k
views
0
answers
... You can change that in the plink bim file ...
written 20 months ago by
Sam • 2.2k
Latest awards to Sam
Great Question
27 days ago,
created a question with more than 5,000 views.
For Piping Markduplicates
Popular Question
8 months ago,
created a question with more than 1,000 views.
For High Kmer Content in middle of the read
Good Answer
10 months ago,
created an answer that was upvoted at least 5 times.
For A: High Sequence Duplication levels in FastQC
Commentator
11 months ago,
created a comment with at least 3 up-votes.
For C: difference between assembly hg19 and GRCh37
Appreciated
16 months ago,
created a post with more than 5 votes.
For A: High Sequence Duplication levels in FastQC
Popular Question
17 months ago,
created a question with more than 1,000 views.
For High Kmer Content in middle of the read
Popular Question
21 months ago,
created a question with more than 1,000 views.
For RNA Sequencing Pooling design and statistic analysis
Teacher
23 months ago,
created an answer with at least 3 up-votes.
For A: Gene regulatory networks from RNA-seq data
Teacher
2.0 years ago,
created an answer with at least 3 up-votes.
For A: Gene regulatory networks from RNA-seq data
Teacher
2.0 years ago,
created an answer with at least 3 up-votes.
For A: Tutorial (How to analyze) on Whole Exome sequencing. Common Errors. Best Practic
Great Question
2.1 years ago,
created a question with more than 5,000 views.
For Normalization Of Rna Sequencing Counts (By Ercc / Gene Length)
Teacher
2.4 years ago,
created an answer with at least 3 up-votes.
For A: Gene regulatory networks from RNA-seq data
Teacher
2.8 years ago,
created an answer with at least 3 up-votes.
For A: Gene regulatory networks from RNA-seq data
Scholar
2.8 years ago,
created an answer that has been accepted.
For A: Gene regulatory networks from RNA-seq data
Popular Question
3.0 years ago,
created a question with more than 1,000 views.
For Piping Markduplicates
Teacher
3.2 years ago,
created an answer with at least 3 up-votes.
For A: Gene regulatory networks from RNA-seq data
Scholar
3.2 years ago,
created an answer that has been accepted.
For A: Gene regulatory networks from RNA-seq data
Teacher
3.2 years ago,
created an answer with at least 3 up-votes.
For A: Gene regulatory networks from RNA-seq data
Teacher
3.3 years ago,
created an answer with at least 3 up-votes.
For A: Gene regulatory networks from RNA-seq data
Scholar
3.3 years ago,
created an answer that has been accepted.
For A: Gene regulatory networks from RNA-seq data
Scholar
3.3 years ago,
created an answer that has been accepted.
For A: Gene regulatory networks from RNA-seq data
Guru
3.3 years ago,
received more than 100 upvotes.
Use of this site constitutes acceptance of our User
Agreement
and Privacy
Policy.
Powered by Biostar
version 2.3.0
Traffic: 1645 users visited in the last hour