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
1
vote
0answers
34 views

I can't launch FastqCleaner I always get a warning message and the application never starts

I tried to install all the needed and related packages but I still did not know what the problem is, Can anyone please help if anything else I can do?? I always get this over and over: ...
1
vote
1answer
31 views

Accessing expression data in an ExpressionSet

I used rma from the oligo package to normalize CEL files, returning an object of class ExpressionSet. How do I access the actual ...
0
votes
0answers
33 views

How to get a consensus sequence from a nanopore fastq files?

I am new in bioinformatic field. I would like to know a way to generate a consensus sequence from nanopore fastq files (fastq files demultiplexed). I usually generate a consensus sequence with "...
0
votes
0answers
23 views

Different results in Differential Expression analysis of Microarray Data

I'm new to R, so be patient with me. I am performing differential gene expression analysis to microarray data for type 2 diabetes donors and nondiabetic donors. when I run the code I get some ...
2
votes
1answer
47 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
46 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
40 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
34 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
78 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
57 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
35 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
45 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
51 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 ...
2
votes
1answer
75 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
153 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
41 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
33 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
106 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
50 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
81 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. ...
1
vote
2answers
451 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
58 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
109 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
37 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
368 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
36 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
143 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
28 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
170 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
356 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
1k 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
118 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
946 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, ...
2
votes
2answers
698 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
693 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
68 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
2
votes
1answer
73 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
184 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
106 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
92 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
544 views

Making a volcano plot look less cluttered

I have generated a volcano plot with a differential expression file. Code for inputing file: ...