Questions tagged [statistics]

For statistical questions related to bioinformatics. If the nature of question is purely statistical, consider stats.stackexchange.com instead.

Filter by
Sorted by
Tagged with
0 votes
0 answers
24 views

statistical test and visualization of genes present in an experiment

I have a list of genes (n = 120) that are involved in breast cancer. And I have a list of differentially expressed genes (n = 70) that are present in the breast cancer gene list. My question is: Are ...
user avatar
  • 75
1 vote
0 answers
17 views

Calculating the effect of differences between two groups of nucleotide sequences at every position

I'm investigating two different groups of short nucleotide sequences. The behaviour of a certain DNA-binding protein is likely influenced by common differences between these two groups leading high or ...
user avatar
1 vote
1 answer
24 views

Why all values become 1 after dcast?

I am reshaping my dataframe as follows. I want to set the Wavelength column as the column names with the code as follows. ...
user avatar
  • 81
1 vote
2 answers
56 views

Why is the p-value significance threshold for HLA association tests $5*10^{-8}$? [closed]

Typically, a p-value of $p<=5*10^{-8}$ is used for genome-wide association testing, as there are roughly $1/p$ independent LD blocks in the human genome, so it correct for multiple testing across ...
user avatar
  • 1,194
1 vote
1 answer
35 views

What is the best way to perform batch ajustment between two experiments?

I have duplicate drug efficacy experiments performed six months apart. I would like to know what is the right way to perform statistical analysis to understand fold changes after treatment with a drug....
user avatar
  • 11
1 vote
1 answer
28 views

Association analysis for methylation beta value and phenotypes in R

I am trying to learn how to do association analysis of methylation array data with numerical trait like body weight, eGFR. As I have seen some of the publications, methylation data association with ...
user avatar
0 votes
0 answers
19 views

How to do enrichment analysis with scipy.stats.hypergeom.sf function

I want to do pathway enrichment analysis in python using scipy.stats.hypergeom.sf function to select pathways that are significantly enriched in at least one tissue....
user avatar
  • 131
0 votes
1 answer
33 views

IBM SPSS- 2-years survival of cohort

Over the past few months I have been retrospectively analyzing patient files for my master's thesis in medicine. As the median overall survival of my cohort is low, I have decided that I want to ...
user avatar
  • 1
0 votes
0 answers
36 views

What is the appropriate test to study the association between the symptom and disease?

I have a stress test for coronary artery disease and I have the following results: D=0 D=1 Y=0 32 21 Y=1 11 82 Y=0 corresponds to a failure of a stress ...
user avatar
  • 103
4 votes
1 answer
103 views

LD Score Regression Derivation hard to follow

I am trying to understand the derivations from Sullivan et al. (2015) in the Supplementary Material. There, it is mentioned in the first page that the least squares estimate of the j-th SNP effect, ...
user avatar
0 votes
1 answer
43 views

Comparing differential expression across samples - is batch effect correction needed?

I have a bulk RNA-seq dataset made up of control and treatment conditions for a range of cell lines. This dataset was generated in two batches, such that the cell lines are split between batches but ...
user avatar
0 votes
2 answers
27 views

How Does Cox Proportional Hazard Model Describe Interaction between Covariates

I have read STDHA and the documentation for the coxph() function in R. I understand that the cox PH model explores how different covariates simultaneously impact survival. Based on the Lung data ...
user avatar
0 votes
0 answers
60 views

How can a stat major get into bioinformatics?

I'm in my first semester of phd in information science and I want to work with bioinformatics professors. As a stat major, what should I prepare for before contacting some of them? I thought about ...
user avatar
  • 1
2 votes
1 answer
35 views

What statistics can I use on qPCR data composed of several runs if all normality tests fail on the dataset?

I have >100 qPCR experiments that I'd like to analyze together, each containing the same set of genes (10 genes of interest and 2 reference genes). I have four different samples (untreated & 3 ...
user avatar
  • 21
0 votes
0 answers
41 views

Linking UK biobank patient IDs with their corresponding genetic data

I am analyzing ukbb data and I have patient IDs for both cases and controls. Now I want the genetic data of all patients based on their IDs. I am trying to test with chromosome 5 as I have both bgen ...
user avatar
0 votes
1 answer
41 views

Understanding the different designs in DESeq2

I am using DEseq2 and trying to understand the results obtained using different models. I have a data design with 2 genotypes and 2 time points. ...
user avatar
5 votes
1 answer
46 views

UPGMA computation

The unweighted pair group method with arithmetic mean (UPGMA) is a hierarchical clustering method, for example used in phylogeny. In a phylogenetic tree it would result in a global molecular clock ...
user avatar
1 vote
1 answer
41 views

Statistical or performance test for genomic intervals between 2 samples

Is there an appropriate statistical or performance test to perform when I want to compare genomic intervals between 2 samples. Sample A is the the benchmark of comparing the intervals and sample B is ...
user avatar
-1 votes
1 answer
608 views

Adjusted P-value from DESeq2

What does the padj really mean? I believe it's a BH method, where it attempts to quantify the likelihood of a "false positive" but I am not sure. Also, what does it mean if the resulting ...
user avatar
  • 9
1 vote
0 answers
16 views

General Question: Statistical Alternatives to Oncoprint

I am a computational genomics student who has recently joined a research lab as an RA. We often work with oncoprints (comuts/heatmaps), and we visualize molecular signatures by reorganizing the comut. ...
user avatar
0 votes
2 answers
65 views

Question about the practical relevance of chromosome position and p-value

Article of interest: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3023908/#SD1 Referring to Figure 2b of that article, if I understand correctly, the x-axis refers to the position of each variant in a ...
user avatar
0 votes
1 answer
20 views

Where can I find a basic example for testing association between a phenotype and each of the variants of a chromosome?

I am a statistics student trying to find data suitable for a particular method. From what I have read so far, it is of interest to determine if variants within a chromosome are associated with a ...
user avatar
1 vote
2 answers
31 views

What is the tP statistic from ANGSD?

We were looking for genetic diversity in the exon level of an MHC transcript and analysing the data under ANGSD (Analysis of next generation Sequencing Data) written in C++ We used ANGSD to calculate ...
user avatar
  • 85
1 vote
2 answers
24 views

P value recalculation from multiple analysis (Transcriptomics)

I've done a differential expression analysis with transcriptomic information (microarrays and RNAseq) from different datasets by applying for some of them t-test and for the rest Wald-test. I've got ...
user avatar
  • 13
0 votes
1 answer
147 views

How to create Beta Diversity distance boxplot

I'm trying to generate boxplot from Beta diversity boxplot between two conditions from a phyloseq object. I have a Morisita-horn distance generated by vegdist fonction : ...
user avatar
  • 21
1 vote
0 answers
89 views

How to compare Log2 Fold Change values

I have 5 sets of Log2 Fold Change values pulled off DualSeqDB. Each of the 5 groups corresponds to a different microbial taxon. Within the 5 groups, there are many Log2 Fold Change values, each ...
user avatar
0 votes
0 answers
30 views

How to decide between standardized mean difference (SMD) and log2fold change (log2FC) before doing meta-analysis?

I want to do meta-analysis of alpha diversity between test and control groups across 10 different populations. I want to see what is the summed effect (pool effect) of intervention between the two ...
user avatar
  • 83
2 votes
0 answers
53 views

What are the best practices for comparing .seg copy number outputs from iChorCNA to VisCapCancer?

I have shallow whole genome sequence data that I used the tool iChorCNA to find the mean copy number segmentation per region. The tool produces a seg file and can be seen here: https://github.com/...
user avatar
3 votes
1 answer
26 views

Tajima-Nei Distance estimate with BioPerl

I have been trying to estimate Tajima-Nei distance for my data (link below). The output is s phylogenetics distance (pairwise) matrix which accommodates a detailed transition transversion and the ...
user avatar
  • 175
0 votes
1 answer
18 views

How do you find differences when your metagenome sample size is small?

I received a matrix of genomics features (KEGG annotated metagenomes) abundance from 6 samples belonging to 2 groups. My collaborator is interested in finding differences between the two groups. The ...
user avatar
  • 3
0 votes
1 answer
27 views

How is the t-statistic value calculated in GEO2R or Bioconductor?

I was trying to calculate t-statistics using Python's scipy and numpy as np- ...
user avatar
-1 votes
1 answer
98 views

Multiple box plot

I am using below code to plot multiple box plot but I am still struggling to understand how I can add whisker to box plot and how to change sample order like in this figure SampleType? ...
user avatar
0 votes
0 answers
69 views

How to read E-value annotation on NCBI BLAST?

i am a bit confused how to read the current E value annotation on NCBI BLAST results. I looked into the matter but could not find an factual answer . What i understood so far: E -values are the ...
user avatar
0 votes
1 answer
312 views

Regarding canonical analysis of principal coordinates (CAP) for linking microbial data with elements

I used the below code to do CAP analysis to link microbial data with elements, but not sure how to do statistical analysis to show which element variation is responsible most for the difference in ...
user avatar
1 vote
1 answer
423 views

Confirmation about PERMANOVA and Variance partitioning (permutest)

Could you please help me to understand what is the difference between Permutational multivariate analysis of variance (PERMANOVA,adonis) and Variance partitioning (permutest). I understood that ...
user avatar
0 votes
1 answer
428 views

How to select specific lines and specific columns

Please I'm a beginner in using R and I work on database. I have a database that contain lot of columns and lines, like this for example ...
user avatar
  • 101
1 vote
1 answer
105 views

How to calculate RNA copies in qPCR using R?

I have these Ct data from qPCR for a series of samples A-H, I want to quantify the RNA copies of my samples using R: ...
user avatar
1 vote
0 answers
160 views

Discrepany between AutoDock Vina and AutoDock Tools

I have been calculating the binding Affinity between a Ligand and a Protein using both AutoDock Vina and AutoDock tools but the binding affinity estimates obtained from the 2 molecular docking ...
user avatar
0 votes
1 answer
67 views

Comparing AUCs: Discrimination of same Control from different Test Group - paired or unpaired? miRNA

I want to compare two AUCs using a bootstraping method from pROC package (roc.test). I have done a classification using Logistic Regression and Serum miRNAs as prediction factors. The Controls are in ...
user avatar
  • 13
0 votes
1 answer
270 views

Power law distribution alpha values

I did a Pearson co-expression analysis for generating networks for my tissue-specific (chondrocytes) RNA-seq data. I used R package poweRlaw to check for the power-law distribution. We got the ...
user avatar
  • 125
1 vote
1 answer
71 views

error in random forest analysis

I am now struggling to do random forest analysis, I will be thankful if you could help with code for random forest analysis. I got samples from the root, soil, and leaf from two regions (bau & mau)...
user avatar
2 votes
1 answer
124 views

alpha diversity wilcox.test

I am trying to do Wilcox test to detect the significant difference in alpha diversity but it is showing an error? physeqN2 is a phyloseq object and Season is a metadata column. ...
user avatar
1 vote
0 answers
28 views

Truncated Chemical Names

[Note: this question has also been posted on KNIME] Good Morning, I am currently researching treatment options for the Corona Virus from the CAS COVID-19 Antiviral Candidate Compounds Dataset ...
user avatar
0 votes
1 answer
353 views

Choosing Fisher's Exact or Binomial test for overrepresentation in PANTHER

The PANTHER website offers a tool to obtain the GO-based overrepresentation of a gene list (the analyzed list) versus a reference gene list. After entering both these lists, the tool asks to choose ...
user avatar
0 votes
1 answer
82 views

Odds ratio and enrichment of SNPs in gene regions?

I did a QTL analysis with a panel of 7M SNPs, and want to analyze the enrichment of the significant qtl-SNPs in different genic regions (promoters, gene bodies, TFBS, etc.). A straightforward way to ...
user avatar
2 votes
1 answer
117 views

Compute the significance of the overlap between 2 or more gene sets

I was able to compute the significance of the overlap between 2 gene sets using the cdf function of scipy hypergeometric distribution. I wish to be able to perform the same calculation for more than 2 ...
user avatar
0 votes
1 answer
68 views

B-factor analysis in nanobodies

I am doing my research on nanobodies. Currently, I am trying to do a comparison study of the flexibility of CDR loops using B-factor values. I would like to know if there is any statistical analysis ...
user avatar
1 vote
1 answer
57 views

Bootstrapping cell type identifications (scRNA-seq)

I have scRNA-seq data held in a SEURAT object, and cell-types were identified with SingleR (essentially a vector of strings - ...
user avatar
1 vote
1 answer
66 views

Bonferroni count = 1 in GO-Slim overrepresentation test (PANTHER) is it possible?

I encountered a problem while performing analysis using the PANTHER website. I run the statistical overrepresentation test for a genelist of 99 genes (GO-Slim BP, Binomial, Bonferroni) in February ...
user avatar
1 vote
1 answer
260 views

How to generate Beta diversity boxplot from phyloseq object?

I am working with profiled metagenomic taxonomy abundance data. I want to generate beta diversity (bray-curtis) boxplot from a phyloseq object where two groups (control and test) will be shown. ...
user avatar
  • 83