Questions tagged [r]

Please supplement your question with a minimal reproducible example. Use dput() for data and specify all non-base packages with library calls. For statistical questions please use http://stats.stackexchange.com.

138 questions with no upvoted or accepted answers
Filter by
Sorted by
Tagged with
0
votes
0answers
28 views

Differential miRNA expression using RPM

I have microRNA (miRNA) expression data in RPM. I would like to do differential gene expression between two groups. How can I do this? I have considered edgeR and DESeq2 in R, but it looks like they ...
0
votes
2answers
125 views

Subsetting a Seurat object based on colnames

I'm trying to subset my seurat object based on colnames. I have gone ahead and labeled each cluster and now I want to subset all the colnames that are in Cancer_human for human_colnames and all the ...
0
votes
0answers
51 views

functions in seurat to calculate the gene count per cluster

I'm running FindAllMarkers ...
0
votes
0answers
203 views

How to group rarefaction curves by metadata

I'm trying to generate a rarefaction curve plot from a phyloseq object. I'm using the function ggrare using rarecurve function from the vegan package: The function rarecurve allow the ploting of all ...
0
votes
0answers
44 views

co-occurrence analysis and visualization for amplicon microbial data

I am trying to learn co-occurrence analysis for microbial data. I got relative abundance at the genera level for 50 genera in 19 samples. I understand that for network analysis first I need to do ...
0
votes
0answers
21 views

Is SNPRelate the goto package for PCA of SNP data in R?

I'm trying to evaluate a method that I'm not familiar with. They used SNPRelate. I'm wondering if this is the best / only choice... what other tools would you recommend? Many thanks,
0
votes
0answers
45 views

seqff script got trouble after samtools treatment

I'm trying to study fetal fraction by using a seqff script from here. The instruction says that I need a headless sam file to work so after alignment with bwa, I process my sam file by using samtools: ...
0
votes
2answers
84 views

Determining what RNAseq data is filtered on volcano plot

I am using RNA seq data to analyze genes via a volcano plot comparing differential gene expression of bacteria with and without antibiotic in R. After having created my plot, I am unsure why some of ...
0
votes
0answers
24 views

How to retreive all the genes belonging to one GO term

I'd like to retreive all the genes belonging to a GO term. I tried (in R): library(org.Hs.eg.db) xx <- as.list(org.Hs.egGO2ALLEGS) ...
0
votes
0answers
24 views

How can I retreive specific qualifiers from genbank files in R?

I'm trying to extract the host field from the features table using the biofiles package. From what I understand, the function filter (from the same package) could be of help but I can't extract the ...
0
votes
0answers
56 views

Adding matrix as a confounder in EnrichmentBrowser

Brain data, particualrly human brain data have very heterogeneus cell types, so it's important to normalise by cell type/add proportoins of different cell types to the formula when performing ...
0
votes
0answers
41 views

FastqCleaner in R, not recognizing paired-end fastq file?

I've recently downloaded the FastqCLeaner to process paired-end fastq files locally using R (4.0.3). I'm having trouble with the shiny app. It is not uploading my second fastq file, and is repeated ...
0
votes
0answers
49 views

tissue cross talk using gene expression data

I am studying about tissue cross-talk using gene expression data. My first goal is to identify the genes responsible for communication signaling(cross talk). Let say the tissues are fat and muscle. ...
0
votes
1answer
95 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. ...
0
votes
0answers
70 views

Help importing a biom file into R

How can I import a biom file? When I run the following code: ...
0
votes
0answers
84 views

workaround for error: Bad CPU type in executable c

I use clustal() from library ape every now and then, but now it doesnt work... code: sylvia.clus = clustal(sylvia.seq) Error: ...
0
votes
0answers
26 views

About getting rs id from chromosome and position

Hi I have a dataset with chromosome and position. VARIANT_ID chromosome position 17_26797415_147499 17 26797415 17_26797556_147500 17 26797556 How can I ...
0
votes
0answers
107 views

Calculate Entropy for DNA Multiple Sequence Alignment in R

I am pretty new to R. So I apologize for asking maybe a very basic question. Let's say I have a fasta file with sequence below: ...
0
votes
0answers
40 views

How can I obtain gene lengths programmatically for genes in Wormbase?

I would like to programmatically, in R, obtain gene lengths (excluding introns) of genes of c. elegans given their Wormbase identifier (for example: WBGene00000004). How can I do so?
0
votes
1answer
125 views

Which module to select for Pathway analysis based on Module trait correlation and pvalue?

I have a total of 35 tumor samples classified into 4 subtypes. Subtype A, B, C, and D. I have RNAseq data. I'm interested in identifying modules related to each subtype with co-expression network ...
0
votes
0answers
55 views

Step-by-Step Construction of Gene Co-expression Networks from High-Throughput Arabidopsis RNA Sequencing Data

I am following the tutorial here to learn how to process raw RNA-seq data and get gene read counts. I've had no issues until I got to step 3.9, "Assignment of RNA-Seq Reads to Genes." Here, the ...
0
votes
0answers
122 views

about the scaling after normalization and transformation in RNA-seq data

When we use count data in RNA-seq analysis, we usually use normalization and sometimes vst, rlog transformation (DESeq2)or log2(CMP+4) transformation (edgeR) to perform K-means clustering. Can we use ...
0
votes
0answers
772 views

Error in as.data.frame.default(x[[i]], optional = TRUE, stringsAsFactors = stringsAsFactors) : cannot coerce class ‘"SeqFastaAA"’ to a data.frame

The example that I am trying to follow is this PGA tutorial. I want to use the information in the analysis of the mzML raw files from the proteomics data analysis. When I try to load the fasta file, I ...
0
votes
1answer
60 views

Finding difference between groups

I have mutational catalogues of 4 samples like below ...
0
votes
2answers
2k views

Percentage of each cluster in Seurat

I am using Seurat to analyze my single cell data. I have 2 conditions, treated and untreated. I am trying to create a stacked bar graph in order to show the differences in cell types for each ...
0
votes
0answers
209 views

Create a Heatmap that uses modules of genes rather than individual genes in Seurat

Is it possible to create a Heatmap in Seurat that takes in and displays modules of genes as its features rather than individual genes? So if I have a module of genes associated with a trait or ...
0
votes
0answers
169 views

module differential connectivity (MDC) analysis

Paper in this paper I am citing the lines below . ) Intra-modular connectivity of resilient and susceptible network modules was compared to that of corresponding genes in control mice to identify ...
0
votes
0answers
196 views

Module preservation analysis WGCNA

I have done a WGCNA (weighted gene co-expression network analysis) analysis for various brain disorders. So along with control, I have 3 different disorders such as BPD, MDD and SCZ. Now I have to ...
0
votes
0answers
39 views

How can extract the list of genes name from the raw data in GEO?

The format of the raw data that upload in GEO is different, like TXT, CEL or other form. How can I extract the list of genes name from these formats? And then how can I normalize them with R? I would ...
0
votes
0answers
60 views

Mapping GPL Platform Number to Microarray Name

I am trying to download a dataset from GEO using the GEOquery R package. Using this package, I can retrieve the GPL platform number(s) for the dataset, and I hope to use a CDF file from BrainArray on ...
0
votes
0answers
34 views

How to find essential genes by differential expression from TCGA 3rd level CNV data?

I have two types of data one have normal patients copy number variation (CNV) data and another have tumor patients CNV data. In each file there are three columns one contain patient ID, second ...
0
votes
0answers
69 views

Any way to filter out highly correlated genes with limma linear model?

I have Affymetrix gene level expression matrix (genes in the rows and sample ID on the columns), and I have annotation data of microarray experiment observation where sample ID in the rows and ...
0
votes
0answers
64 views

Normalization of data with rpkm

I'm very i difficult with normalization of my data. I was searching for transposable elements in my genome, and after this step, I made counts of reads in some transcripts. I produced something like ...
0
votes
0answers
155 views

maf file format: Variant_Classification value for non frameshift substitution

I have a vcf containing SNV and MNV (short indells or substitutions), I annotate that file whit annovar and in the resulting txt file I have this row ...
0
votes
1answer
114 views

Too slow issue of BioMart

I am using library(biomaRt) in R to annotate SNP rsid. It takes too much time for only 1000 SNPs....... I would like to know how I can run this code faster. ...
0
votes
1answer
28 views

multi--platform GEO analysis based in expression of single gene

This question has also been asked on Biostars I have a challenge in analyzing GSE34941( It has 4 platforms) according to the expression of a specific gene in R. First, which is the best way to group ...
-1
votes
0answers
24 views

how to get the info if a given gene expressed in single cell sequencing or not

I have been trying to understand if TP53 is expressed in various organs of Tabula Muris database. I tried few packages but I cannot get to figure out how two check each genes and visualize it. for ...
-2
votes
1answer
96 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? ...

1 2
3