Questions tagged [clustering]

Clustering is the task of grouping a set of objects in such a way that objects in the same group (called a cluster) are more similar (in some sense) to each other than to those in other groups (clusters).

Filter by
Sorted by
Tagged with
0
votes
3answers
94 views

RNAseq biological replicates not clustering in PCA plots

I have RNAseq data from 4 samples with 3 biological replicates per sample. I am currently trying to do the differential expression analysis with DESeq2 but the biological replicates will not cluster ...
0
votes
1answer
152 views

Clustering information saved in Seurat object

Seurat V2 had a option to find clustering information saved in object: PrintFindClustersParams(object = pbmc). How can I get the same clustering parameters from objects in Seurat3? Do I need to ...
2
votes
3answers
89 views

Getting hierarchy of cell populations with Drop-seq data

I plan on collecting some drop-seq data from brain tissue (from mice and humans). Using only the 75-85% of genes that have 1:1 orthology between mice and humans, I'd like to cluster the cells by ...
0
votes
0answers
23 views

which marker identifications method is recommended when single-nuclei RNA seq datasets are integrated across regions

I'm analyzing Sn-RNA seq datasets which include different brain regions and spinal cord in patients and normal controls. I would have different comparison within and between individuals. I've got a ...
0
votes
1answer
62 views

How to add cluster name to metadata in Seurat?

I'm working on a Seurat object and want to name the clusters according to 2 values alone (yes/no). So I want to add a new column to metadata and annotate the clusters (UMAP) with it. ...
0
votes
2answers
75 views

What kind of analysis is practically done on GSE data files?

I have a GSE data file in csv file format containing fields such as: ID, adj.P.Val, P.Value, t, B, logFC, Gene.symbol, Gene.title. In which adj.P.Val, P.Value, t, B, logFC fields being numeric. What ...
1
vote
0answers
18 views

What is the best approach to classify a patient cohort by the expression (low, int, high) of a single gene of interest?

I am working with a dataset of nearly 100 patients. I performed Salmon quantification with genecodeV34 and imported the results with tximeta. I normalised the TPM salmon output with TPM (using edgeR, ...
1
vote
1answer
21 views

classifying samples by TCGA signature

I have some RNA-seq samples from multiple glioblastoma tumours that I'm now trying to classify according to a specific gene signature (from Verhaak et al., 2010) using R. The gene signature is ...
6
votes
2answers
211 views

Are phylogenetic tree construction algorithms any different than general clustering algorithms?

Are phylogenetic tree construction algorithms any different from clustering algorhithms? I suspect the answer is no. Of course phylogenetic tree construction uses biological knowledge, e.g special ...
0
votes
1answer
33 views

Changing active.ident in Seurat

Im trying to change the active.ident to another column in metadata but this error keeps popping up! I recently upgraded to R version 4.0.2 from 3.6.1 The older version was working but the new one isn'...
3
votes
0answers
44 views

How are the values of prop.part() and prop.clades() calculated?

Consider the following dataset: ...
2
votes
1answer
121 views

sva for RNA-Seq data without known phenotype

I have been working on RNA-Seq data from two different cohorts, and they show very strong batch effect (~35% variance explained by 1st component in PCA). Since I am trying to do a class discovery from ...
2
votes
0answers
57 views

Clustering issue about overlapping CD4, CD8, and Foxp3 genes in mouse tumor models

bioinformatics scientists. I hope your work goes smoothly well and be safe. I'm stuck in trouble at the step of clustering. Briefly, I downloaded fastq files from a published study (https://www.ncbi....
0
votes
1answer
40 views

about the normalization of RNAseq data for calculating distance for unsupervised learning

I have been working on clustering using RNAseq data. To compute distance, what kinds of normalization is optimal? Can we use normalization using relative log normalization? I understand this should ...
1
vote
0answers
243 views

How to create a dendrogram of clusters and reconstruct phylogenetic relationship

I'm working with a Daphnia data set looking at the 16s gene pulled from the BOLD database. So far I did a multiple sequence alignment and attempted to cluster and plot a dendrogram: ...
2
votes
0answers
32 views

vst() from DESeq2 vs voom() from limma

I have used both to transform my leukaemia RNA-Seq data for subsequent hierarchical clustering. The result is quite different. Some subtypes of leukaemia only form a cluster (or at least sit closer) ...
12
votes
2answers
1k views

How to decide number of neighbors and resolution for Louvain clustering

I am using Louvain clustering (1,2) to cluster cells in scRNAseq data, as implemented by scanpy. One of the parameter required for this kind of clustering is the number of neighbors used to construct ...
1
vote
1answer
44 views

How to identify latent variables in single-cell RNA-Seq data

I have a single-cell RNASeq sample, in which I'd like to identify latent variables (e.g. response to stress) that I think might be affecting the clustering. The approach I was planning to use is to ...
1
vote
0answers
162 views

Understanding BuildClusterTree of Seurat

I am trying to understand how to use BuildClusterTree of Seurat to merge clusters. Being no informatics or statistical background, I am not understating how to use it. I am having the following ...
2
votes
2answers
49 views

Tools or softwares for clustering evaluation

I am using a few clustering algorithms as well as my own tool for protein sequences. Also, I have benchmarked clustering results (i.e. ground truth) to compare with the results. Is there any tool/...
1
vote
1answer
39 views

choosing a good representative genome subset

I'm trying to build a genomic database for DNA alignments. I started with NCBI accessions, but the data is very multiplicative, so I want to use subset of [max] N different strains for each specie. my ...
1
vote
3answers
80 views

K means clustering, would PCA be a better option?

I have the data below. I need to use a clustering method to classify them and into categories of "Heterozygotote, Allele 1, Allele 2 and No Call. The values in RFU1 and RFU2 are used to determine the ...
0
votes
0answers
89 views

Error in my RNAseq analysis

I was following the RNAseq analysis tutorial online(https://combine-australia.github.io/RNAseq-R/06-rnaseq-day1.html) but am not obtaining the same variance values for each row in the logcounts matrix....
0
votes
1answer
255 views

how to interpret of “pvclust” dendrogram and finding height for cutting dendrogram?

I study on RNA-seq expression dataset about one cancer in TCGA. I downloaded FPKM dataset and removed batch effect by ComBat() function. Now, I used pvclust for ...
1
vote
0answers
100 views

Error using bseqsc

I will be very grateful for any hint on how to overcome the error. I wish to deconvolve my bulk RNA seq data obtained from the lungs of mice using single cell RNA seq data. For practice, I am ...
2
votes
1answer
85 views

How I know which gene is a good predictor in this neural network or not?

I have 25 highly differentially expressed genes among and patients to chemotherapy. I have made a neural network of these genes. Accuracy of model is 0.73 but I don't know from 25 genes how I could ...
4
votes
2answers
1k views

Seurat for clustering bulk RNA-seq?

Is it ever ok to use Seurat for clustering bulk samples? I am looking at FPKM data from ~750 bulk RNA-seq samples generated using Cufflinks. As suggested for FPKM ...
2
votes
2answers
73 views

Find most similar genes to a set of genes of interest in time series RNA-seq data

I have gene expression data (RNA-seq) for 30 different time points (from 0 to 60 min each 2 min). I have a set of 8 genes that behave similarly (although not identically) and I want to find the top X ...
0
votes
1answer
271 views

Cluster is split in 2-3 locations on tsne plot - Suerat

I am running a single cell dataset (count data - exon) through Seurat. After running tsne I see a cluster (13) split in 3 different locations on the plot. Here are the commands I am running: ...
0
votes
1answer
43 views

Clustering using cellranger

Is it possible to run cell ranger clustering using only a gene x cell matrix? Cellranger reanalyze performs clustering but requires a hd5f file. I don't want to run the whole pipeline (including ...
1
vote
1answer
280 views

Hierarchical clustering for outlier detection - single cell RNASeq & WGCNA

I wish to conduct WGCNA on a single cell RNASeq dataset and, when choosing the optimal beta parameter after running pickSoftThreshold am presented with a rather ...
2
votes
2answers
368 views

Interpretation of hierarchical clustering with bootstrapping

I have data that includes 'cases' and 'controls' and have carried out hierarchical clustering. I used the pvclust package in R to bootstrap the results and ...
1
vote
0answers
54 views

How to write the subclusters in file?

I have constructed a gene co-expression network from RNA-seq data. The network file is in edge list format of memory around 1gb which was created by calculating Pearson correlation of each gene pairs ...
1
vote
0answers
119 views

GBS: clustering of PCoA axes with R package mclust to describe identity by missingness

I have a SNP dataset with 56 populations and 430 samples that was produced by GBS sequencing. I need to identify groups of samples that have similar patterns of genotype missingness, since there is ...
3
votes
1answer
872 views

Find Patterns in Cluster

I have a heatmap and I would like to find some rectangles. I have already used clustermap. But here, I can not calculate these rectangles. The order of the data should not be changed. This Code is ...
2
votes
2answers
82 views

Using multidimensional scaling to visualize protein sequences by functionality

I have a multiple sequence alignment (MSA) of protein sequences, with which I have performed multidimensional scaling to visualize their clustering. Using R, I used the ...
2
votes
2answers
169 views

Can I get the graph generated by cellranger

I’ve run the cellranger analysis pipeline on single cell RNASeq datasets. I can import the matrix and graph-based clusters into R. Doing this I can optimise the dimension reduction and plot cells with ...
0
votes
0answers
21 views

Comparing between different set of clusters produced from the same network

I have a set of eight gene clusters which I obtained from clustering a network. However, I had to update the genes in the network at a later stage which lead to certain deletions (genes not present in ...
0
votes
1answer
221 views

Using external list of PCs for clustering

I am going to use principal components (PCs) comes from calcPCA function in URD program for clustering my cells in Seurat; So ...
2
votes
2answers
44 views

How to combine multiple kernels of large sample datasets?

I have multiple large sample datasets in matrix format (each has 15000 rows and 5-50 columns) corresponding to different experiments. Each matrix contains the same number of samples(rows) but the ...
0
votes
1answer
30 views

Check if a set of genes fall more on a specific cluster than other clusters

I have a network which I clustered into 11 smaller clusters. I have a set of genes and I would like to check if statistically significant number of these genes fall on any specific cluster. One method ...
0
votes
0answers
47 views

Non-redundant set of indices for AAindex?

AAindex has 566 indices with real values for each of the 20 amino acids mostly from wetlab measurements. Many of these values are highly correlated, and thus do not make good features for ...
4
votes
1answer
353 views

RNA-Seq: clustering/treatment of genes with low expression

I have some RNA-Seq data from leukaemia patients. I want to do unsupervised clustering on them with some other published leukaemia RNA-Seq data and see how they cluster. There are a few problems I ...
0
votes
0answers
32 views

Quality checking benchmarks for clustering

I'm currently working with several RNAseq datasets from the ENCODE database. My approach to assessing the quality of the data, i.e, reproducibility is by automating the assessment of expression ...
4
votes
2answers
748 views

Omics data: How to interpret heatmap and dendrogram output?

How to interpret heat map and dendrogram output for biological data (omics) in words (when writing results and discussion)? What should I consider (statistics behind?) and what is the best approach? ...
1
vote
0answers
50 views

case/control phenotype prediction - population stratification with R

I'm struggling these last days on a case/control phenotype prediction problem based on genotype (individual x SNPs with more than 10k individuals and approx 18k SNPs encoded as 0,1,2). Thanks to PCA ...
3
votes
2answers
18k views

How to scale the size of heat map and row names font size?

I have an expression data matrix (120X15; 15 samples and 120 genes), my heatmap looks blurred and raw names (gene names) looks very small and can not read. How can I improve my scripts? Here is the ...
5
votes
1answer
148 views

Evaluate clusters of individuals by using their sequence data

For a dataset of several hundred individuals, I applied a hierarchical clustering to generate clusters based on a functional trait that sets them apart. My task is now to evaluate if these clusters ...
4
votes
1answer
107 views

Is removing samples based on clustering for downstream analysis a right choice?

I'm using TCGA Lung cancer data. I'm interested in doing differential analysis between Lung vs Normal. Before DEA, to check the distance between each pair of samples I plotted an MDS plot: In this, I ...
1
vote
2answers
332 views

What could be the reason for the samples not clustering?

I'm performing RNA-seq analysis. I have used Hisat2 for aligning reads to the genome and stringtie for quantification and extracted read count information directly from the files generated by ...