User: lamia_203

gravatar for lamia_203
lamia_20370
Reputation:
70
Status:
Trusted
Location:
Last seen:
1 week, 3 days ago
Joined:
1 year, 8 months ago
Email:
l********@hotmail.com

Posts by lamia_203

<prev • 23 results • page 1 of 3 • next >
0
votes
0
answers
114
views
0
answers
Spaces in chromosomes when plotting Manhattan plots
... Hi all. I am using qqman manhattan function for plotting. I have used this function several times before and never had these issues. I have three columns, CHR, BP and P and used the function as ; manhattan(Results, genomewide = -log10(threshold), suggestiveline = FALSE, ...
R manhattan qqman written 11 weeks ago by lamia_20370 • updated 11 weeks ago by i.sudbery9.2k
0
votes
1
answer
212
views
1
answers
Comment: C: Meta-Analysis effect size is not in the same direction as GWAS effect sizes.
... Thank your for the response and confirming the estimate values in R. I've looked at the meta analysis alleles and they are flipped. Not sure how that has occurred but it makes sense why your estimate value is the same as mine but the sign is the opposite. Many thanks for your help! ...
written 6 months ago by lamia_20370
1
vote
1
answer
212
views
1
answers
Comment: C: Meta-Analysis effect size is not in the same direction as GWAS effect sizes.
... Yes, I've pulled the most significant SNP from the meta-analysis then used that SNP to select from the other studies. I've managed to sort it out - the allele direction was flipped hence the correct absolute effect size but incorrect direction. I'm not sure how that has happened to some files and no ...
written 6 months ago by lamia_20370
4
votes
1
answer
212
views
1
answer
Meta-Analysis effect size is not in the same direction as GWAS effect sizes.
... Hi, I am conducting a meta-analysis and the output doesn't make sense. I ran gwas on several different treatments and some treatments work well, other do not. The effect size for the 3 gwas results are 0.016,0.012 and 0.011 with respective SE of 0.005,0.005,0.004. However the effect size for the m ...
R metal meta-analysis gwas written 6 months ago by lamia_20370 • updated 6 months ago by Lemire590
0
votes
0
answers
155
views
0
answers
Comment: C: Multiple commands in a for loop
... I'll take note thank you. It worked with the semi-colon but still produces outputs. I wanted to carry the amended file each command carries out onto the next command and avoid using the origin file each time as well as producing those files in the shell but it looks like the latter occurred. ...
written 8 months ago by lamia_20370
0
votes
0
answers
155
views
0
answers
Comment: C: Multiple commands in a for loop
... Thank you that has worked, but awk still appears with 'plink.awk' so I still have to edit commans once pasted on the shell. All my file names start with plink. but I'm not sure why the attach to the awk command. ...
written 8 months ago by lamia_20370
0
votes
0
answers
155
views
0
answers
Multiple commands in a for loop
... Hi, I am trying to make multiple commands within a for loop but it is not working as awk appears as plink.awk once pasted on the shell. I've written the command below. I have attempted to pipe is as I do not need intermediate files. Any help would be much appreciated, thank you! for f *.linear; ...
genome loops sequence unix bash written 8 months ago by lamia_20370
6
votes
1
answer
831
views
5 follow
1
answer
How to combine count matrix with mouse gene information from Ensemble using Biomart in R
... Hi, I am trying to add a row of gene id next to my rows of gene symbols. normlaised_mouse3 is my count matrix for rna-seq analysis. I've used to code: library('biomaRt') normalised_mouse_biomart <- read.delim("normalised_mouse3.txt") mart <- useDataset("mmusculus_gene_ensembl", u ...
ensembl biomart rna-seq R written 19 months ago by lamia_20370 • updated 18 months ago by zx87549.6k
0
votes
1
answer
600
views
1
answers
Comment: C: Aligning using hg18
... The genome was used before and aligned fine and gave reads at expected genes for another study I carried out (beta cells), however I am also unsure as to why it didn't map the reads I provided (they are human as well). The parameters were kept the same as before as well: for i in *.fq.gz; do ST ...
written 19 months ago by lamia_20370
0
votes
1
answer
600
views
1
answers
Comment: C: Aligning using hg18
... Thank you. I downloaded it from [UCSC][1] I believe, I couldn't find the gtf files clearly on there. This is my first time downloading genomes hence why I am unclear. [1]: http://hgdownload.cse.ucsc.edu/goldenPath/hg18/bigZips/ ...
written 19 months ago by lamia_20370 • updated 19 months ago by RamRS30k

Latest awards to lamia_203

Popular Question 8 months ago, created a question with more than 1,000 views. For How to speed up trimming using trim galorer
Popular Question 16 months ago, created a question with more than 1,000 views. For Creating Count Matrix
Popular Question 16 months ago, created a question with more than 1,000 views. For Alignment loop using STAR

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1258 users visited in the last hour