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
4
votes
1answer
97 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
40 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
601 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
90 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
118 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
200 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
175 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: ...
1
vote
1answer
136 views

I have a list of protein targets. How to get structural data (atomic coordinate) from PDB?

I have a list of around 500 protein targets (human). I only have their protein symbols. As an example: EGFR AR FOXA1 CXCR3 ... I want to get structural data (atomic coordinate) from PDB (Protein ...
0
votes
1answer
92 views

After running a DRIMSeq pipeline, how do I know which genes are upregulated in the different conditions?

After running a DRIMSeq pipeline and obtaining the genes that are differently used between the null and full models, how do I select the genes that are differently used in the different conditions? ...
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 ...
3
votes
1answer
2k views

Using `install.packages` with conda-managed R

I have R installed and managed by conda (miniconda) on my MacBook Pro. The version of R I use most frequently (3.5.1) is installed on the base environment and I have other version-specific ...
1
vote
1answer
33 views

Is there a computational method/method/way to predict if a cell population in a tissue immigrated or is enriched locally?

I have single cell data which I have analysed for differential expression. In my experiment, I subjected two groups of mice (a control and a treatment group) to treatment that will lead to the ...
0
votes
1answer
385 views

Error: package or namespace load failed for 'ChIPseeker'

How do I solve the problem of loading the bioconductor Chipseeker package? I have installed in the following way and I having an error while loading it and I do not know why. ...
2
votes
2answers
286 views

Issue with visualising cladogram/phylogenetic tree with multiple sequence alignment data in R?

I would like to visualize tree with multiple sequence alignment. My din.newick file is shown below, ...
2
votes
3answers
149 views

Getting a stretch of genomic ranges from a dataframe/granges object based on metadata column

I have a "test" data.frame object in R, which is basically a small subset of a 66000 row dataframe, which looks as follows: ...
2
votes
1answer
89 views

Discordance of Gene Ontology data between Mouse Genome Informatics and Ensembl sources

In my first approach to the use of Gene Ontology, I have found a discordance if we consult either Mouse Genome Informatics (MGI) or Ensembl for these data. Concretely, for the mouse genes Zic1 (NCBI ...
1
vote
0answers
222 views

install bioconductor

I used if (!requireNamespace("BiocManager", quietly = TRUE)) install.packages("BiocManager") BiocManager::install() for installing "bioconductor" but I ...
0
votes
1answer
55 views

How can I predict if the binding site of a DNA-binding protein is located on the promoter sequence of a gene?

Is there any computational tool or any other method to predict if the binding site of DNA binding protein is actually located on a promoter sequence?
2
votes
0answers
53 views

How to create custom txDb using GenomeFeatures?

I have orf virus sequence in GFF3 format. I am trying to create a txDb using GenomicFeatures::makeTxDbFromGFF() function. But, it is not able to capture any CDS: <...
1
vote
0answers
78 views

Queries regarding MetagenomeSeq code

I am trying to understand metagenomeSeq, I will be very much thankful if you could help me to understand this: I do have three treatment condition (T1, T2, and T3) and want to see how these ...
1
vote
1answer
499 views

Extracting expression data in r from GEO and got a S4 object

I want to open the expression data GSE9838 dataset from GEO using GEOquery R package. gset <- getGEO(filename="GSE9838_family.soft.gz") The result of this ...
3
votes
1answer
339 views

load the phenotype data for ballgown

I am trying to reproduce the work of this paper [1], and I have run StringTie successfully, but after that I have to run Ballgown but could not understand this command: ...
4
votes
1answer
293 views

How to extract gene expression tables from this GEO dataset?

I've downloaded this GSE43013 dataset using GEOquery in R. My understanding is that it contains expression data from liver, kidney, and brain for multiple species. I would like to produce gene ...
3
votes
0answers
212 views

Differentially methylated position analysis in a related sample?

I'm trying to figure out how to do a DMP analysis (using minfi dmpFinder) on a related sample (if it's even possible). Right now the code (not written by me) is: ...
1
vote
1answer
257 views

cn.MOPS fails with 'missing value where TRUE/FALSE needed'

I am trying to use cn.MOPS to call CNVs on a set of whole exome sequencing data bam files. My script (up to the point where it fails) is: ...
0
votes
1answer
120 views

How do I perform a pathway related grouping of genes?

I have a list of gene names from the analysis of my mouse single seq data. I want to perform a simple grouping of the genes based on the pathways to which they belong or pathways through which they ...
3
votes
1answer
6k views

Error in as.vector(x) : no method for coercing this S4 class to a vector

I tried to run the following code in R studio. Everything worked fine, except at the last step [write.table(mdat, "recount_mdat.csv")] when I tried to export the 'mdat', I got the following error: <...
2
votes
3answers
235 views

Question about Co-expression analysis and finding targets for lncRNAs

I have a dataset with 88 tumor and 113 normal samples. Among the 50k genes after filtering there are a total of 28k genes (both mRNAs and lncRNAs) I wanted to do co-expression analysis between ...
1
vote
3answers
571 views

How can I install bioconductor-gviz and use it in jupyter notebook?

I tried to install gviz in a conda environment, but that library seems to be incompatible to python and r. I tried to setup a clean environment using ...
2
votes
0answers
488 views

R package CMD check errors

I have written some codes in R for my personal use. The package successfully installed on my personal laptop (where I have written those commands) and working correctly. But when I distributed it to ...
4
votes
1answer
154 views

How to find genomic overlaps faster?

I was trying to find a way to find overlaps between two genomic ranges. I found a post at Biostars but it couldn't consider the chromosome information. For example: ...
4
votes
1answer
212 views

How should I address batch effects in my experiment?

Let's say I have an RNA-Seq experiment, where I'm interested in the significantly differentiated genes between pre-treatment and post-treatment conditions. "rep" == biological replicate. ...
3
votes
1answer
3k views

How to transform a DNAStringSet from the Bioconductor package Biostrings to a data frame?

I am working on Mac OS X. I am using R version 3.5.1. I have important a FASTA file into R using Biostrings::readDNAStringSet. This creates a DNAStringSet object ...
1
vote
0answers
89 views

Bioconductor, how to retrieve the genes from a GEO gene expression dataset?

I'm quite new to Bioconductor and GEO. Here's my task: I've to develop a simple R script able to read a gene expression dataset on GEO and extract the genes related to it. But... I don't know the ...
3
votes
1answer
79 views

RSVSim insertions from chr1 to chr2

edit explaining python tag I would still rather have a solution based on the RSVSim package in R, but while waiting for someone who might have an answer I wanted to look at other solutions as well. ...
2
votes
0answers
36 views

Mutation detection using Varscan2 on RNA sequencing for estimating tumour clones with pyclone or other package

I would like to analyze my RNAseq profiles from bulk tissue samples (Paired-End, 50M reads/sample, tumour-normal pairs) with varscan2 to detect mutations. Then I ...
6
votes
2answers
676 views

Transcript Coordinate Ranges to Genomic Coordinates

I have 2 GFF3 files: Features using transcript IDs as the landmarks. i.e. "CDS" feature types using coordinates from transcript space. Features using chromosome IDs as the landmarks. i.e. "exon" ...
1
vote
1answer
55 views

Error in annotating row names as data points in a PCA plot

I made a PCA plot and was trying to plot the eigenvenctors in R so that each data point is actually the sample name. All the samples are in column 1. The following are the R commands I am using: <...
0
votes
1answer
138 views

Curl error when installing Bioconductor packages with Microsoft R Open 3.5.0

I'm trying to install Bioconductor packages with Microsoft R Open 3.5.0 seems to fail not sure what is the issue. This is what I'm trying to do ...
3
votes
1answer
232 views

Aligning sequence and comparing it against primer

I am looking to show how a primer is consistent among some genomic data. I have a primer of about 23bp and looking to compare it to about 5000 genomic sequences of 10kb. I am unsure how to format it ...
3
votes
1answer
168 views

Get genomic coordinates using GenomicFeatures by HGNC gene names

I want to get coordinates of human genes from my list (consisting of hgnc genes id) using GenomicFeatures and TxDb.Hsapiens.UCSC.hg19.knownGene R packages from Bioconductor. ...
2
votes
1answer
166 views

How to plot circular chromosomes with Gviz?

I am using Gviz to plot several tracks. I create a GenomeAxisTrack() and then several AnnotationTrack()s that I feed to ...
2
votes
1answer
39 views

Filtering the expression values of a gene from GEO

I'm filtering the gene expression value of a particular gene using Bioconductor in R. Following the instruction given here, I could filter the gene expression value for a gene from the dataset GDS3782....
5
votes
3answers
117 views

Sampling haplotypes

I am trying to simulate different genome of peoples, I have data (VCF files) of various genes from the 1000K Gene project. I want to simulate different whole genomes i.e generate a new population by ...
3
votes
1answer
247 views

Error in seq.default in chromPlot

I am using chromPlot to visualise the genome of C. elegans. library(chromPlot) I have created the following data frame with the lengths of C. elegans chromosomes. ...
1
vote
1answer
125 views

How to superimpose barplots in karyotype graphic using the chromPlot Bioconductor package?

I have a data frame containing 4 columns: Chrom: chromosome name of the feature Start: starting base of the feature End: end base of the feature Name: name of the feature, e.g. 5S ribosomal RNA The ...
1
vote
1answer
2k views

How to calculate logCPM across all samples?

Using edgeR for differential analysis between Tumor and Normal gave me differential expressed genes with logFC, logCPM, PValue and FDR. From the details of glmTreat function I see that logCPM is ...
3
votes
2answers
540 views

Detecting differentially expressed genes with foldchange >= 2 and FDR < 0.05

I'm using edgeR for differential analysis. Using glmTreat function I'm detecting differentially expressed genes between Tumor and Normal. I have set the arguments like below: ...
3
votes
2answers
1k views

Input normalization in ChIP-seq

If I subtract input counts from ChIP counts (for every gene, since I have one peak per gene) I get negative values for most genes. This makes sense to me, because (as can be seen in the figure) input ...
2
votes
2answers
192 views

CNSeg no method coercing S4 class to vector

I'm attempting to analyze some genomic data with the help of a package called CNSeg. The package is located here and the documentation is here. My input data looks like: ...