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
53 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
12 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
77 views

Making a volcano plot look less cluttered

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

Error in my RNAseq analysis

I was following the RNAseq analysis tutorial online(https://combine-australia.github.io/RNAseq-R/06-rnaseq-day1.html) but am not obtaining the same variance values for each row in the logcounts matrix....
0
votes
0answers
15 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
0answers
19 views

Geneoverlap library data input

GeneOverlap I am trying to use this library so in the default data sets it has different list of genes which are categorized as high low or medium expressed genes. So my question is can I run this one ...
0
votes
0answers
11 views

co-expression between LncRNA and mRNA based on P.value and correlation

I have two sets of data, each of which has two treatments with two replications. The first group consisted of 2261 LncRNA (Long non-coding RNA) with FPKM each LncRNA and the second group consisted of ...
0
votes
2answers
45 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
18 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
56 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: ...
0
votes
0answers
38 views

Row labelling in a heatmap in R

Could I have specific row labels in bold text but not all of them when creating a heatmap using R? (eg. using heatmap.2 function)
1
vote
1answer
52 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
58 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
13 views

GEOquery, Bioconductor question: Is there a more efficient way to tease out datasets that contain data I need?

We are doing a study analyzing the expression of certain genes and correlating that with response to chemotherapy. So far I have been manually going through every dataset on the NCBI website and ...
0
votes
0answers
26 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 ...
2
votes
1answer
289 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
0answers
20 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
85 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. ...
0
votes
0answers
8 views

how to create database and taxonomy for arsM amplicon analysis

Could you please guide me on how I can create the arsM database and taxonomy like 16s data for the arsM gene using the primer set. I need to analyse the arsM microbial diversity using amplicon ...
0
votes
0answers
29 views

How to order grouped elements in Gviz

Gviz is useful for visualization of gene transcripts by genomic coordinates. I'm using it to show the difference between splice variant transcripts. However, it seems that the ordering of transcript ...
2
votes
2answers
156 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
77 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
72 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
105 views

install bioconductor

I used if (!requireNamespace("BiocManager", quietly = TRUE)) install.packages("BiocManager") BiocManager::install() for installing "bioconductor" but I ...
0
votes
1answer
42 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
20 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
48 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
135 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 ...
2
votes
0answers
105 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: ...
3
votes
1answer
125 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
135 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
73 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: ...
-1
votes
1answer
78 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 ...
2
votes
1answer
908 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
178 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
243 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
157 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
75 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: ...
3
votes
1answer
108 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. ...
2
votes
1answer
852 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
78 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
61 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
27 views

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

I would like to analyze my RNA seq profiles from bulk tissue samples (Paired-End, 50M reads/sample, tumour-normal pairs) with varscan2 to detect mutations. Then I would like to use those detected ...
6
votes
2answers
313 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
43 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
95 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 ...
2
votes
1answer
95 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
0answers
133 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
121 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
35 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....