Questions tagged [bioconductor]

Bioconductor provides tools for the analysis and comprehension of high-throughput genomic data in the R language.

Filter by
Sorted by
Tagged with
2
votes
1answer
37 views

DESeq2. Which is better: using the altHypothesis argument in the results function or manually filtering for your desired P-value and Log2FoldChange?

In DESeq2 you can identify your differentially expressed (DE) genes using the results function. I noticed people in my lab using the results() function with the minimum number of arguments supplied (...
2
votes
1answer
41 views

ensembldb equivalent in python

Is there a python equivalent to ensembldb? I want to get genome coordinates for a transcript (like the function transcriptToGenome) but need to do it in python. I ...
1
vote
0answers
26 views

Is New Tuxedo approach really better than the earlier ones?

I am planning to perform multiple routine NGS analyses for a hospital lab on patients' samples on an everyday basis. Before I used a new tuxedo approach with ...
2
votes
0answers
39 views

Simultaneously get data from multiple applied gates in flowCore

I'm applying two parallel and non-overlapping gates to a gatingSet directly under "root": ...
0
votes
2answers
31 views

Bioconductor Conda release delay

How long does it take for the latest Bioconductor to end up on Conda? I think I saw a blog post about the challenges with this process, but I cannot find it. For example, Bioconductor 3.13 was ...
0
votes
1answer
34 views

RDAVIDWebService (R package) manual installation, or alternative tool

Recently, I had to re-install all my R packages (R version 4.1.0 (2021-05-18) -- "Camp Pontanezen") and it turns out that "RDAVIDWebService" is ...
1
vote
0answers
14 views

Bioconductor, genefilter() returns NULL dimensions, is my filtering process wrong?

I'm using R and bioconductor in order to conduct some gene analysis on an Affymetrix dataset (GSE173360) but I'm having trouble at filtering genes. I'm trying to create a dataset called small.eset by ...
0
votes
1answer
51 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 ...
1
vote
1answer
34 views

Calling GEOquerry error while loading the library

Note: this question has also been asked on Biostars I started getting a strange error when trying to load GEOquerry with library(GEOquerry). ...
0
votes
1answer
37 views

STRINGdb: reverse map STRING identifiers to gene names?

Is there a functionality within STRINGdb package or in another environment to reverse map the STRING identifiers to Gene names? For example, the entire STRING network can be downloaded using the ...
0
votes
1answer
38 views

Retrieving all genes for a gene ontology term

I'm looking for an easy way to retrieve all the genes in a list that are associated with a certain GO term, preferably using R/Bioconductor packages. I'm not interested in under/overrepresentation or ...
0
votes
0answers
20 views

Trouble aligning next generation sequencing data to reference genomes using QuasR package in Bioconductor. Cannot import .txt

I'm trying to check the quality of my paired end read sequencing data. I am following this pipeline (https://f1000research.com/articles/4-1062#ref-21) which uses QuasR in the first step. My list of ...
1
vote
0answers
45 views

How do I include silent mutations in an oncoplot using maftools?

I have annotated files from annvar and I changed them into MAF by using annovarToMaf from ‘maftools’. Then I used read.maf to ...
4
votes
1answer
87 views

installing multiple bioconductor packages at once

I was wondering if there is more elegant way of installing and loading multiple packages in Bioconductor similar to pacman with CRAN packages. I tried: ...
1
vote
1answer
39 views

Retrieve overlap widths with Bioconductor

I want to retrieve all overlaps between two sets of genomic intervals (GRanges) using Bioconductor. I am interested in the widths of the resulting overlapping ranges. In cases with discontinuous ...
1
vote
0answers
31 views

workflow for processed RNAseq data analysis in R [closed]

I'm learning to use R in data analysis. I'm getting fluent in baseR and the tidyverse, but thus far only dealt with medium throughput plate-based experiments. I am currently trying to learn how to ...
2
votes
1answer
91 views

`TxDb.Hsapiens.UCSC.hg19.knownGene` and `EnsDb.Hsapiens.v75` can't be installed

I think there is a problem with TxDb.Hsapiens.UCSC.hg19.knownGene and EnsDb.Hsapiens.v75. I'm struggling to install all of ...
0
votes
1answer
39 views

Get list of urls of GSM data set of a GSE set

I have this GSE dataset ( GSE104279 ) (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE104279). I want to make a table with set IDs and ftp urls to use it as a table in galaxy.org I know that we ...
1
vote
1answer
61 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)...
0
votes
1answer
75 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
15 views

Methods to analyze MeDIP-seq data

I have MeDIP-seq data. I have 6 groups with 15 replicates and a total of 90 human samples. I am trying to choose the best method to analyze them. I would highly appreciate if you can give me your ...
1
vote
2answers
358 views

Python equivalent to R GRanges and IRanges

Problem: I am trying to convert some codes written in R to Python and part of that conversion process is find classes equivalent to the GRanges and IRanges from the GenomicRanges R package in Python. ...
0
votes
1answer
57 views
1
vote
0answers
17 views

Questions about EBSeq (in R)

I have made a UMAP of malignant cells and the result is split into 3 clusters, as seen below. For the sake of example, let's say I want to find the differentially expressed genes in the uppermost ...
1
vote
1answer
94 views

Power calculation for GWAS/EWAS

I want to investigate, how much sample size i needed to obtain 80% power for GWAS/EWAS studies. Phynotype trait is discrete (not case/control) for human disease. I wonder, does anyone has came across ...
0
votes
0answers
36 views

a proper Design Matrix for several drug treatments with both control negative and control positive

I have a dataset of RNA-seq samples for testing different drugs on the presence of another drug. One of my samples is the normal cells with no drugs (control negative) and another is the cells with ...
1
vote
1answer
35 views

DiffBind, diffferentially binding site

I have data for 3 histone marks (2 for silencing and 1 for activation) each mark has three replicates. when I run the diffBind package I have three contrast: ...
2
votes
3answers
365 views

How to reduce the occupied RAM when you are dealing with a very sparse matrix in a single-cell Experiment in R?

I'm dealing with a very large and sparse dataset and the first issues I met occurred when I tried to use quickCluster that reported me this error: ...
0
votes
1answer
34 views

Is there a way to customize what observations are plotted in ComplexHeatmap?

So right now my observations in the dataset I've processed is stored as sample identifier columns and gene rows: an example of this would be for sample 1 and gene TRIM21 the observation is Missense ...
1
vote
1answer
125 views

Plotting metadata from GenomicRanges object

I have a GRanges object with metadata that I want to plot. Many of the online tutorials which manipulate these GRanges for plotting show how to plot the main columns (chromosome,start/end,strand) but ...
0
votes
0answers
27 views

Grabbing all Immune related genes with databases in R

I am having trouble grabbing specific pathway info using databases in R. I have RNAseq results and I want to remove immune related genes from the current list that I have. With a vector of gene ...
0
votes
1answer
150 views

WGCNA co-expression network analysis with less than 20 samples

I have two cancer subtypes data. Subtype A is 14 samples and Subtype B is 23 samples. I'm interested in identifying the functions of some LncRNAs in the Subtype A group. For this I'm using all protein-...
0
votes
1answer
294 views

PCA plot did not work in single cell RNA-seq

I run plotPCA for single cell RNA-seq data, while I get error message (I use R 4.0). I attached the code and error message here. Did anyone know the reason and how to deal with it. Thanks! ...
0
votes
2answers
902 views

calculateQCMetrics defunct, how to calculate the quality metrics by perCellQCMetrics

I am trying to follow a tutorial from Sanger institute (from May 2019) on analysis of single cell RNA Seq data. They use calculateQCMetric function to calculate the quality metrics, but I am getting ...
0
votes
1answer
108 views

Get regions' information from DNA sequence data (bsgenome.hsapiens.ucsc.hg19)

I have a problem in R. I have the following dataSet (the first three rows shown) (the 5th number is the methylation level in its region): ...
0
votes
3answers
873 views

isSpike function in SingleCellExperiment package is deprecated?

I am trying to follow a tutorial from Sanger institute (from May 2019) on analysis of single cell RNA Seq data. They use isSpike function to filter out ERCC (control) and MT (mitochondrial RNA) reads, ...
1
vote
1answer
544 views

Using bioconductor from Python

Has anyone used bioconductor from Python? Is there any reason I would choose to use it from R instead of Python? It seems like there is a Python extension for it. Also, is there any reason to use ...
0
votes
0answers
641 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 ...
1
vote
2answers
35 views

Identify differentially covered genes only between two samples

I have a question about finding differentially covered regions (coverage represents methylation level which goes from 0 to several thousands). I'm using enrichment based method which can be summarized ...
1
vote
2answers
67 views

Is there a tool for synonymous recoding of DNA sequences?

I've got a DNA sequence that I'd like to make synonymous mutations throughout, thereby preserving the amino acid sequence. Does anyone know of a tool to achieve this
1
vote
1answer
60 views

How to calculate average minor allele frequency (MAF) difference between two populations (i.e. European, Hispanic)

The data is present in the PLINK format .bed (individual genotypes), .bim (genetic markers) and .fam (sample IDs and disease phenotype file). I have also separate text files containing information on ...
1
vote
1answer
155 views

qvalue and p.adjust functions

Between the qvalue package and the p.adjust function, which is more appropriate to use when trying to calculate the q-values of a dataset? According to the manual for the q-value package, the q-value ...
1
vote
1answer
83 views

When running ChIPQC bioconductor package after peak calling, should I use sorted BAM files or unsorted BAM files?

I have called peaks using MACS2 and now I'm trying to run ChIPQC bioconductor package. Can someone please tell me should I use sorted BAM files or unsorted BAM files? I have created sorted BAM files ...
4
votes
1answer
84 views

Can I submit a R package to Bioconductor or CRAN if I have already published it a journal?

I have written a bioinformatics package in R that I want to publish in a bioinformatics Journal. Presently, I am maintaining a local repo of that package and I want to put in the Bioconductor ...
0
votes
1answer
36 views

Error in createMAE function: non-unique values when setting 'row.names' in TCGA LIHC Data

I am working with the createMAE() function of the ELMER Bioconductor package. While executing the createMAE() function using ...
0
votes
2answers
516 views

Making a volcano plot look less cluttered

I have generated a volcano plot with a differential expression file. Code for inputing file: ...
1
vote
0answers
78 views

include a glimma interface in a shiny app

I am trying to code a shiny app for RNA-Seq data analysis. I would like to include glimma interactive plots in it. However, in my current interface, clicking the action button ...
0
votes
2answers
98 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
2answers
151 views

How can I find all the columns available in UniProt.ws package for R?

I am trying to find the subcellular localization of my 10,000 proteins using UniProt.ws package for R. However, I am unable to find all the columns available for query. I used another package named ...
0
votes
2answers
155 views

GenomicRanges: get nearest neighbor distances for random genes using a for loop

I want to compute the nearest neighbour distances using a for loop. I can do this with a random subset by doing this: ...