Most Popular
1500 questions
3
votes
2 answers
pvalue calculation of protein-protein network with permutation test
Note: this question has also been asked on Biostars
I have performed a network recreation analysis based on the interactions of proteins from String db.
I wanted to know if the interactions have significance by reshuffling the proteins with…
user3698773
- 179
- 3
3
votes
1 answer
What is the R/Bioconductor preferred way of accessing TCGA data?
It seems there are two packages providing access to TCGA data these days:
TCGAbiolinks (older)
GenomicDataCommons (newer)
Will GenomicDataCommons supersede TCGAbiolinks? Is there anything that I might miss if I decide to use only…
Ramiro Magno
- 165
- 1
- 7
3
votes
1 answer
What is the significance of the effect size on WOMAC scores for TKA?
https://www.medicaljournals.se/jrm/content/html/10.2340/16501977-2210
I understand the difference between statistical significance and effect size, but I am having trouble interpreting the data on Table III in the article linked above. I am…
CharlieCal
- 39
- 1
3
votes
3 answers
Extract LINE-1 methylation status from WGBS data
I've been trying to figure out the correct way of extracting methylation status of particular elements from WGBS data, and found several possible approaches, but would like to discuss which one could be the best:
using software dedicated to the…
Kraken
- 405
- 2
- 9
3
votes
1 answer
genes with highest RPKM/FPKM for a human RNA-seq experiment?
Which genes/transcripts show up as having the highest RPKM/FPKM on a human RNA-seq experiment? Is it usually the same genes/transcripts?
EDIT: mRNAs (e.g. poly-A RNA-seq preps) and not including small RNAs.
719016
- 2,324
- 13
- 19
3
votes
4 answers
True bam file encoding and viewing as text
We all have been using samtools to convert between sam and bam files. SAM is tab-delimited file. My question resides on the ability to visualize bam files as text. Is there any way to visualize any bam file as text in linux?
Why not work only…
andresito
- 385
- 1
- 3
- 9
3
votes
1 answer
How to retrieve InChI key for KEGG compound?
I would like to retrieve the InChI representation for a given KEGG compound but I fail to find a direct solution for this (preferably doing it via bioservices).
One could do it via ChEBI like this:
from bioservices import *
kegg_con =…
Cleb
- 743
- 7
- 18
3
votes
1 answer
Database of copy number alterations in different cancer types
Is there a database of somatic copy number alterations (CNAs) in different cancer types? That is, annotations of which CNAs are more commonly seen in different cancer types.
719016
- 2,324
- 13
- 19
3
votes
1 answer
Greedy Motif Search Using Probability Matrices
In my greedy search with pseudocounts algorithm in my bioinformatics course, I did not follow the pseudocode since I wanted to solve the problem in my own way. Unfortunately, although my answer is correct, the answer we are supposed to give is the…
DrJessop
- 131
- 1
- 3
3
votes
4 answers
How to calculate p-values for fold changes?
I'm currently working with data from a Luminex multiplex assay. In this assay, the concentrations of 17 different analyte proteins were identified in 12 groups in triplicate. One of these 17 groups was used as the control, and the log2 fold changes…
J0HN_TIT0R
- 541
- 1
- 4
- 7
3
votes
2 answers
Can gene co-expression networks be used to help identify differentially expressed genes?
In my RNAseq dataset of differentiated stem-cell lines, some samples have far fewer significantly differentially expressed genes than others. QC shows that this is because there are way fewer reads for one sample than the others.
Can gene…
jaslibra
- 524
- 2
- 9
3
votes
2 answers
Why are some of the gene peptides returned by biomaRt missing an asterisk in the end?
Note: this question has also been asked on BioStars
I am getting peptides using biomaRt library in R for the gene 'BRCA1'. 27 different aminoacid sequences are returned, and 12 of those have an asterisk (*) in the end, while 15 do not. What does…
user5054
- 305
- 1
- 8
3
votes
1 answer
IGV faster loading remote bigWig files
I noticed remote bigWig files (http URL) take sometimes about 1 minute to load on IGV. Once they loaded, jumping from one region to another can also take several seconds.
My understanding is that bigWig files have an internal index for different…
719016
- 2,324
- 13
- 19
3
votes
2 answers
How do I summarise VCF files?
I have got a number of VCF files that all contain data on one particular variant, with every VCF being created from an individual sample. So there is one entry per file (headers…
Nereus
- 209
- 1
- 5
3
votes
1 answer
Annovar index out of date
I'm trying to use Annovar to annotate some variants with their CADD and FATHMM scores. I've downloaded the latest versions of the software and the databases but when I run it I get an error saying the index is out of date.
WARNING: Your index file…
Greg
- 831
- 6
- 12