Questions tagged [bioconductor]
Bioconductor provides tools for the analysis and comprehension of high-throughput genomic data in the R language.
745
questions
0
votes
0
answers
39
views
specific Bioconductor packages fail when deployed to shinyapps.io
I am trying to deploy a shiny app to shinyapps.io in RStudio and keep getting an error with a Bioconductor package.
I've tried the following in console before deploying: options(repos = BiocManager::...
0
votes
0
answers
35
views
How to calculate the frequency of every pathogenic germline variant in every disease cohort?
I am trying to calculate the frequency of every pathogenic germline variant in every disease cohort. For ex: for variant 1:17588689, there are 20 het.variants (column H), I need to report what ...
0
votes
1
answer
52
views
How to install edgeR package using R4.2.0, BiocManager?
I tried to install edgeR with the following:
if (!require("BiocManager", quietly = TRUE)); install.packages("BiocManager"; BiocManager::install("edgeR")
This is the ...
0
votes
1
answer
84
views
Why does R show a package error immediately when I open the program?
When I open R, before I type in any commands, the console immediately shows an error indicating that the curatedMetagenomicData package is not installed.
R version 4.2.0 (2022-04-22 ucrt) -- "...
0
votes
1
answer
21
views
How can I start manipulating an scRNA seq .txt matrix in R?
I'm really new to this, so please forgive my general lack of understanding.
I've been trying to take the dataset that is provided here (I download the .txt file at the bottom). It's scRNA-seq data and ...
0
votes
0
answers
10
views
Error in file(con, "r") : cannot open the connection to website
I am using this tutorial to perform some analyses:
https://bioconductor.org/packages/release/bioc/vignettes/rWikiPathways/inst/doc/Pathway-Analysis.html
However, I get an error message when running ...
0
votes
0
answers
22
views
R: Unable to install FlowSOM. Getting Warnings
I wish to install the FlowSOM package in R using Bioconductor (https://bioconductor.org/packages/release/bioc/html/FlowSOM.html)
RStudio 2022.02.2+485 "Prairie Trillium" Release (2022-04-19)...
0
votes
0
answers
35
views
Install R package without biocmanager
I am using bioconductor image in order to install R packages. The problem I am facing is that I can't install specific version of the package.
I have the following Dockerfile:
FROM bioconductor/...
1
vote
0
answers
45
views
Mutation plot splitting cohort
I am approaching some R packages to create oncoplot representing different mutation types.
I am using maftools (oncoplot function) and genvisR (waterfall function) packages.
I would like to split ...
0
votes
0
answers
33
views
Loop for subsetting SingleCellExperiment object in R
I have a SingleCellExperiment object (sce) which I would like to subset by cluster name (19 clusters) to return a list of 19 SingleCellExperiments which I can then use as input to plot expression of 9 ...
0
votes
1
answer
49
views
How to re-adjust figure size of trackViewer Vignette lollipopPlot and avoid x-axis cut off
I am new to trackViewer Vignette lollipopPlot, I am experimenting with the plot. However, my plot looked very wide and have white space for some reason that I would like to readjust the width and ...
4
votes
0
answers
93
views
package from Bioconductor has a dependency from CRAN that requires version of R that is not compatible
A particular situation we've encountered is when we try to install a specific R package version from Bioconductor. Dependencies of this package may not be pinned, which means that the latest and ...
-1
votes
1
answer
33
views
How to calculate length distribution and coverage distribution bed file
I have dataset in bed file and I want to calculate and plot the length and coverage distribution of file. How can I calculate length distribution in R.
df1:
chr21 2800 3270
chr21 3600 4152
chr2 3719 ...
0
votes
0
answers
46
views
Unable to install GEOquery
I just installed Bioconductor in RStudio 4.1.1:
if (!requireNamespace("BiocManager", quietly = TRUE)) install.packages("BiocManager")
I am currently trying to install GEOquery ...
0
votes
1
answer
40
views
Error in EnsDb.Hsapiens.v86 package to convert ensembl to symbol
I installed package EnsDb.Hsapiens.v86 to convert genes ensembl to symbol,
after trying many packages that's the best one to conversion with less NAs.
I tried to write the code but I got unsual error, ...
1
vote
1
answer
52
views
How do I use biomaRT to get the corresponding gene iDs
I have a txt file and it looks like this. I need to use biomaRT in R to get the corresponding gene IDs of a whole list of different Refseq and peptides. Along with that, I need to keep the peptide ...
0
votes
0
answers
12
views
How to map RefSeq IDs to gene IDs Bioconductor with R
I have Refseqs and peptides and I am trying to map them with Gene IDS, I got everything loaded and I can read the documents but I can't get the Gene IDS to map. How would I do that?
0
votes
0
answers
56
views
Error with C Stack Usage when installing packages
I haven't touched R in a while, and just tried running some old Bioconductor code on my new Mac (M1, 2020, running OS Monterey 12.2.1).
I am able to run the beginning of my script (where I install ...
0
votes
1
answer
42
views
Substitute string sequences for BLOSUM62 values in R
I have a df with several peptide sequences. For each sequence, I need to replace each AA for the BLOSUM62 values stored in another matrix.
|Sequences|
|:--------|
|KIFREIIHL|
|LALTYRGLI|
|DIKGLIVYR|
|....
1
vote
1
answer
59
views
GRanges to subset by overlap within 200kb of start or end of a gene
I have two GRanges datas and I would like to subset them by overlaps such that the overlap could also be present within 200kb of either start or end of the gene.
I was using the following command
...
0
votes
1
answer
67
views
AnnotationHub installation ‘Rcpp’ 1.0.6 is already loaded, but >= 1.0.7 is required [duplicate]
I am currently trying to install AnnotationHub but I am running into problems. I use R 4.0.4, renv (version 0.15.4) and BiocVersion 3.12.
When I run BiocManager::install("AnnotationHub"), I ...
0
votes
0
answers
29
views
How to add analysis of variance to a Bioconductor plot that combines jitter and boxplots?
I made a plot that shows median expression of some markers across different cell populations.
I made it using a Bioconductor package that is based on ggplot and the final result is a ggplot object. ...
0
votes
1
answer
63
views
data.table::fread error when converting MAF files to data table
I want to merge the 50 MAF files with the sample information so that I can read it as a data.table and subset it.
library(maftools)
# Load MAF files
maf = system.file("extdata", list.files(...
0
votes
0
answers
39
views
All values turned up NaN and NA after data correction and normalization in R
I have a set of 21500 data, I performed background correction and normalization using bioconductor packages MBCB.main however the result is only NaN values for one column and NA values for the other. ...
2
votes
1
answer
3k
views
R message "Can't subset columns that don't exist. x Locations 2, 3, 4, 5, 6, etc. don't exist. i There are only 1 column" MBCB packge
I am trying to do a background correction on my data sets that has only 3 columns, when submitting my code I get the error "Can't subset columns that don't exist. x Locations 2, 3, 4, 5, 6, etc. ...
0
votes
1
answer
123
views
Problems installing older version of Bioconductor's mixOmics packages in R
I've spent the day trying to load the appropriate package versions in R that I saved in a renv lockfile.
I used the package RVAideMemoire which is tied in with mixOmics in bioconductor, which can't be ...
0
votes
0
answers
32
views
Empty host volume in docker container
I am using Docker for Windows, and wanted to mount a host directory with files I would want to use in RStudio in a container for a bioconductor image. To mount the host directory I have used
docker ...
0
votes
0
answers
79
views
Contrast for Limma - Voom
I'm doing a differential expression analysis for RNA-seq data with limma - voom. My data is about a cancer drug, 49 samples in total, some of them are responders some of them are not. I need some help ...
4
votes
1
answer
60
views
R: How do I reduce the ranges of genes into a single vector?
I created var_nt dataframe by subsetting tx_df columns based on row variant - "J3", "J10", "J11", "J13".
Then, I converted the var_nt dataframe to a GRanges ...
1
vote
0
answers
18
views
does bioconductor/bioconductor_docker:RELEASE_X_Y install CRAN packages from a fixed snapshot?
Question
Does a container instantiated from bioconductor_docker:RELEASE_X_Y also pull CRAN packages from reproducible image (e.g. via MRAN) and if not how do I accomplish this?
Motivation
I'm building ...
0
votes
0
answers
220
views
Trouble installing Rhtslib in R/R studio
I'm using RStudio on Ubuntu 18 and I'm trying to install the htslib package from the Bioconductor repo, but I'm stuck now.
This is what I get:
* installing *source* package ‘Rhtslib’ ...
** using non-...
0
votes
1
answer
35
views
Convert R list to Pythonic list and output as a txt file
I'm trying to convert these lists like Python's list. I've used these codes
library(GenomicRanges)
library(data.table)
library(Repitools)
pcs_by_tile<-lapply(as.list(1:length(tiled_chr)) , ...
0
votes
0
answers
40
views
error when reading multiple samples in Seurat
I'm trying to read multiple samples using Read10X function from the Seurat package.
seurat_list <- lapply(samples, function(sample){
seurat_data <- Read10X(paste0(gsub("MNX.*", "&...
2
votes
1
answer
133
views
Error in as.data.frame.default(x[[i]], optional = TRUE) : cannot coerce class ‘structure("mmquery_bioc", package = "multiMiR")’ to a data.frame
I'm having trouble getting a table of results.
Command:
answer <- get_multimir(url = NULL, org = "hsa", mirna = "MIMAT0000450", target = NULL, disease.drug = "cancer", ...
3
votes
0
answers
59
views
Debugging unexpected S4 method dispatch
I am trying to pass an S4 object (sparse matrix of class dgCMatrix from the Matrix package) to a function in another package (Wrench). I am getting some very strange behavior that I am struggling to ...
0
votes
0
answers
62
views
is there any way to download Deseq1 package in R?
I'm trying to load some RData files with the function load, but I get this message:
Warning: namespace ‘DESeq’ is not available and has been replaced
by .GlobalEnv when processing object ‘cds3’
...
1
vote
0
answers
35
views
Why do I have a different amount of estimated gene counts?
I am analyzing bulk RNA seq data and I used Kallisto to align my data to the transcriptome. Then, I used tximport to assign the gene names from ensembl to my counts. I am comparing the results I ...
0
votes
0
answers
26
views
installation fails for package that imports CRAN packages that import Bioconductor packages
I have a CRAN package say "myPack", that imports two CRAN packages, say "cranA" and "cranB". These two packages import Bioconductor packages, say "biocA" and &...
-3
votes
1
answer
157
views
How can I plot multiple sequence alignment in R as a boxplot and scatterplot:
I've fasta file contains 16 protein and I aligned them with package ("msa").
now I need to Visualize my results of similar sequences using boxplots and scatterplot
I read about "ggmsa&...
0
votes
0
answers
68
views
Using R package Biostrings in Python
I am trying to use the R package Biostrings in Python and I am having trouble applying the function "DNAStrings" from Biostring to every row in a dataframe. Here is my code so far:
pip ...
0
votes
1
answer
24
views
readAxt() converts lowercase to uppercase
Is there any method where I can read AXT files without converting the sequences to uppercase automatically?
The link to the code is: https://github.com/alexander-nash/kurtosis_conservation/blob/master/...
0
votes
1
answer
45
views
string alignment inspired by global DNA sequence allignment
I would like to do something similar to this:
library(Biostrings)
s1 <-DNAString("...
0
votes
0
answers
21
views
GEOQuery - Cannot load a large GSE file
I'm trying to access a large GSE dataset from the GEO. I am using GEOQuery to try and do this. I have tried to access the file directly through the GEO:
gse <- getGEO("GSE13159",GSEMatrix=...
0
votes
0
answers
63
views
When downloading 'TxDb.Mmusculus.UCSC.mm10.knownGene' package, lazy loading error appears, what do?
Code:
if (!requireNamespace("BiocManager", quietly = TRUE))
install.packages("BiocManager")
BiocManager::install("TxDb.Mmusculus.UCSC.mm10.knownGene")
Error:package ...
0
votes
0
answers
98
views
adding labels to genes or segments with ggbio and geom_alignment
I'm trying to use geom_alignment in the ggbio package to plot a stretch of genes, and I want to show labels under each gene. The docs for geom_alignment say that this is done by specifying label = ...
3
votes
1
answer
96
views
Can you silently install Bioconductor packages?
I'm writing a pipeline in Snakemake that calls an R script. This R script has its own environment, with r-base, r-ggplot2 and r-biocmanager in it. I also need the package ggbio that can be installed ...
0
votes
0
answers
9
views
Will Kallisto show a decrease in the TPM value if I have knocked out 2 of the 8 exons from a gene?
I was wondering if Kallisto will show a decrease in the TPM value if I have knocked out 2 of the 8 exons from a gene? I have confirmed the knock out using qPCR however, when I align with Kallisto the ...
1
vote
2
answers
79
views
Overlapping and sorting counts from small, equally sized, continuous bins into irregular, uneven bins
I really appreciate if you are reading this and taking your precious time to help me with a problem I have.
In R, I would like to sort data from small, continuous bins from one dataframe to the (non-...
0
votes
0
answers
49
views
BiocParalell in R on SGE
I have a R-script which accepts BPPARAM to run the job in parallel.
I would like to run it on SGE cluster.
The bash script I use has this setting:
#!/bin/bash
#$ -pe smp 25
#$ -N my.job
#$ -e /home/...
0
votes
0
answers
23
views
Error in gzfile - analysis using compcodeR
I'm doing differential expression analysis using the compcodeR package.
I'm trying to insert the raw counts data file with this code:
runDiffExp(data.file = 'SRP029880.raw_counts.tsv',
...