Questions tagged [statistics]
For statistical questions related to bioinformatics. If the nature of question is purely statistical, consider stats.stackexchange.com instead.
82
questions
0
votes
3answers
65 views
What is the difference between fixed effects and random effects in the context of Linear-mixed models?
The terms in LME i.e. fixed and random create confusion? What is the genesis that can distinguish between the two ?
1
vote
1answer
30 views
How to compute Shanon entropy
I want to analyze some viral sequences. I want to use shanon entropy as a measure to analyze the sequences. However I do not understand the concept of shannon entropy. Would anyone explain please how ...
3
votes
1answer
79 views
Comparing phylogenetic models with different datasets
I'm a linguist interested in phylogenetic tree inference using language data. I'm posting here because I'm using Bayesian phylogenetic methods in my work (probably using BEAST and/or RevBayes). For ...
2
votes
1answer
96 views
Prediction model for spread of coronavirus
A previous thread on Bioinformatics StackExchange citing an article by the Economist which references a cartoon of predictive outcomes by the CDC (Economist article on coronavirus) suggests a model in ...
0
votes
1answer
49 views
a package with randomization/permutation tests for genomic analyses
please would you let me know, is there a package in R/python/C that has implemented randomization/permutation tests for joint genomic analyses
(an example of joint genomic analyses -- when jointly ...
4
votes
1answer
147 views
Economist article on coronavirus
I am wondering about an article in the Economist here:
https://www.economist.com/briefing/2020/02/29/covid-19-is-now-in-50-countries-and-things-will-get-worse
There is a graph there
The explanation ...
2
votes
0answers
43 views
How to deal with zeros in the denominator when calculating percent enrichment?
I have a data.frame with results from a CRISPR screen.
I want to calculate the % enrichment of a variable of interest var from ...
1
vote
3answers
368 views
Calculating the mortality rate of a pandemic, e.g. coronavirus 2019-nCov?
When calculating how lethal a virus outbreak is, I've noticed that most sources use
CONFIRMED_DEATHS/CONFIRMED_INFECTIONS
However, given that there are only 2 ...
1
vote
2answers
43 views
How to exclude the repetition of gene-gene correlation calculation in python?
I am trying to make a paired matrix of gene-gene correlation. Considering that I have a huge matrix (13000 genes and 900 samples) and for some reasons I don't want to decrease the number of my genes, ...
-3
votes
3answers
200 views
Is it possible to do this in R?
I have a big clinical data file with 96 columns like age, gender, BMI, etc
I want to see which of these clinical characteristics respond to chemotherapy. Response to chemotherapy divides patients to ...
1
vote
1answer
28 views
Most scalable pseudotime ordering algorithm
What algorithms for linear Pseudotime trajectory construction (diffusion-based) are the most scalable to large datasets?
I'm currently using Slingshot based on the recommendation in this manuscript: ...
0
votes
0answers
10 views
Quartert and Statistical Consistency
I found that triplet is statistically consistent, meaning, Given large amount data, most occurring Gene tree containing three taxa will likely be a true species tree I understand the proof. Is ...
0
votes
1answer
26 views
Use of Electronic Phenotype in EHR
May I know what's the use of Electronic Phenotyping using EHR data?
I did refer this link but have few questions
I understand that ...
3
votes
1answer
68 views
protein secondary structure matching/comparison scoring
I want to compare two secondary structures of aligned proteins. I do not know exactly how to do it well.
Example:
...
0
votes
0answers
9 views
Recalculating fold change in microarray from a subset of samples
I have some bulk microarray data that I have pulled from online sources which describes the spatial distribution of transcription. The tissue in question has a few different subregions. I am ...
1
vote
0answers
43 views
Detect differentially expressing cells
after comments on my original post, I will ask my question again here
I have data (RNA expression values, obtained with multi-channel in situ hybridization) collected from 1mio human cells. For each ...
2
votes
0answers
26 views
Trouble optimizing Five Parameter Logistic (5PL) Standard Curve for ELISA data using Python
I am attempting to write a tool to be used by my lab to automatically generate a 5PL Standard Curve for ELISA data, using an XLSX file template. I am using pandas dataframes to hold the 96 wells as ...
3
votes
1answer
59 views
What is the state of the art for GWAS in terms of the statistical algorithm for either Case/control and Quantitative traits?
This question was also asked on Biostars
I'm trying to understand what is the best algorithm for GWAS nowadays.
I know we have many tools available like Plink and Hail, but currently, what is the ...
0
votes
1answer
38 views
Interpreting a kernel density function plot of gene expression
I need help interpreting a kernal density plot and a histogram of this data. I have provided some context below let me know if you need more information
A co-worker that I am collaborating with has a ...
1
vote
1answer
37 views
Normal score transformation
This question is related to a normal score transformation which is performed in the paper by
RUAN, Quansong, et al. Local similarity analysis reveals unique associations among marine bacterioplankton ...
1
vote
0answers
37 views
Should I use log2-CPM values (voom-limma) as input for my model?
We have created a model to integrate several OMICs data, but we realized that the maximum TPM values of RNA-Seq data were so big that had unexpected effects on our results. We hypothesized that this ...
2
votes
1answer
50 views
How to interpret Mendelian randomization results?
I did Mendelian randomization using this software:
https://cran.r-project.org/web/packages/MendelianRandomization/vignettes/Vignette_MR.pdf
...
3
votes
0answers
23 views
Find covariation of different DNA-binding protein binding between two conditions
I have ChIP-seq data of RNAPII and two other factors which we think follow RNAPII during transcription in two different conditions, and I'd like to show that the genes that lose RNAPII also lose the ...
2
votes
1answer
95 views
Confidence Interval with Wilcoxon Test in Python for log-normal Distribution
I have a couple arrays of dN/dS scores, and I would like to calculate the confidence interval for each array of data. dN/dS scores are not normally distributed but are log-normally distrbuted, so I ...
0
votes
1answer
71 views
Biplot scores from canonical correspondence analysis
I'm using the R package vegan to perform canonical correspondence analysis (CCA). As input we have two matrices, one being (sites)x(species) and the other being (sites)x(conditions).
Sample data (...
0
votes
0answers
54 views
length of 'dimnames' [1686] must match that of 'dims' [3]
Please if anyone has experience with the use of the BSEQ-SC package for the deconvolution of bulk RNA sequencing data with single cell RNA sequencing data I will be very grateful for your suggestion.
...
2
votes
1answer
188 views
Create a Grouped Barchart
I’m trying to create a Group Chart like this:
I have 7 classes for each of 4 species, not 5 like in the diagram. I want to build i diagram to show the distributions of the classes in each species
...
1
vote
1answer
32 views
Total read count of shotgun affected by few highly abundant sequences
I am building statistical models to analyse output from Illumina shotgun sequencing (HiSeq 4000) on stool samples (but RNA-seq data should behave similarily). The raw counts are a statistical sample ...
2
votes
1answer
65 views
How are residuals used as new phenotype?
I'm a novice trying to wrap my head around with fancy jargons. If I search residual, I understand that it represents the difference between the observed value of the dependent variable (y) and the ...
0
votes
0answers
14 views
Statistics in EGSEA
I am not sure if this question should be posted here or cross validated so please move if it is not appropriate.
I am planning to use ensemble of gene set enrichment analysis and trying to understand ...
2
votes
0answers
39 views
Number of reactions per metabolic pathway
What is the quickest way to calculate the number of chemical reactions present in each pathway for a given organism and which online database (kegg, metacyc,reactome..) is more helpful/reliable for ...
0
votes
2answers
62 views
Statistical test design question
I have condition as such mouse liver and spleen data.
So in liver I have WT (Wild Type) and KO (knock out) as well as in spleen WT (Wild Type) and KO (knock out).
I have to compare the parasite load ...
0
votes
0answers
58 views
Which statistical test to use for comparing two groups?
I have data for 396 transcriptional factors (TF), where I have calculated the number of Transcriptional targets for each transcriptional factor present between the two groups:
Group-I(G-I): where co-...
1
vote
0answers
42 views
What are the mathematical assumptions and algorithmic procedure behind MACS2?
I wonder if there's some mathematical details that outline what the algorithms does.
0
votes
2answers
88 views
residual Squared Coefficient of Variation (rCV²) vs Distance to Median (DM)
I am working with single cell RNA-seq data. I obtained the squared coefficient of variation (CV²) as a measure of gene expression variability:
I want a metric to express gene expression noise that ...
2
votes
1answer
90 views
Significant gene set testing - limma
While I have used my title, it is more of a general statistics question regarding how these significance tests are carried out, these include decideTests and ...
3
votes
1answer
318 views
Understanding PCHeatmap outputs
I am currently trying to understand the purpose of these PCHeatmaps - part of the seurat package in R:
All the online documentation I have searched for has only ...
0
votes
2answers
192 views
1
vote
1answer
156 views
Extracting genes differentially expressed by Wilcox test
I have 2 cell types (NOF and CAF) cultures individually and also have been co-cultured with tumour like this picture.
https://ibb.co/x8Ty1Wz
If this is the log cpm of these 4 samples
...
2
votes
2answers
116 views
Genomics Statistics Problem
I have used a python script to identify target sequences in a DNA sequence file. There are two classes of sequence: coding and non-coding.
I have identified 728 sequences of interest. 597 of these ...
1
vote
1answer
50 views
Continuing analysis based on a log-likelihood ratio? (Comparing methylation between two PoolSeq phenotypes)
I'm trying to get some pointers on where to start on the analysis of my dataset. First, some basic info:
Background:
We performed Nanopore sequencing of two poolseq samples called UU and DD referring ...
3
votes
0answers
35 views
2
votes
1answer
1k views
3D PCA group labelling
I would like to make a 3D PCA but not sure how to label group wise which i can do for 2D PCA
...
4
votes
1answer
67 views
Why are my Chi-squared test results different from those in a published table?
I recently read the paper “A novel long non-coding RNA linc-ZNF469-3 promotes lung metastasis through miR-574-5p-ZEB1 axis in triple negative breast cancer”. In this I see Table1 showing correlation ...
1
vote
1answer
33 views
What does it mean for an amino acid to match with itself 99% of the time in a PAM1 matrix?
In this video, the lecturer says:
You see that phenylalanine is probably still phenylalanine because only one percent of amino acids have changed so roughly 99 percent of the time, phenylalanine is ...
2
votes
0answers
54 views
Block wise protein imputation
I am currently working on a dataset that contains 50 samples (10 samples * 5 blocks). The data is perfectly balanced between blocks, with equal treatment representation in each block. Each block ...
2
votes
1answer
77 views
How can the Autocovariances, autocorrelations, and autocorrelation coefficients be calculated from a Protein Amino Acid Sequence?
Given a normal protein sequence with the 20 standard amino acids, how can the 'Autocovariances', 'autocorrelations', and 'autocorrelation coefficients' of the sequence be calculated?
What is meant by ...
3
votes
0answers
44 views
What is the correct way of dealing with the analysis of data from flow cytometry?
I would like to detect a change in expression of a molecule present on a cell type by flow cytometry.
Assuming I am able to detect, using an antibody, a signal that represents the amount of the ...
5
votes
1answer
224 views
Should PCA be standardized for gene expression?
This is a theory/good practice question more than a technical one. If samples are being plotted on a PCA projection of gene expression data, I'm wondering whether it is standard (and if so, why) to ...
1
vote
1answer
70 views
A basic question about phylogenetic inference with Bayes-MCMC
I am new to the world world of Bayesian-MCMC inference and have been working my way through various texts. I am having some difficulty understanding how MCMC works. I understand that it is a way of ...
