Questions tagged [r]

Please supplement your question with a minimal reproducible example. Use dput() for data and specify all non-base packages with library calls. For statistical questions please use http://stats.stackexchange.com.

Filter by
Sorted by
Tagged with
0 votes
0 answers
15 views

BiocManager::install("DiffBind") mergeOne.c:341:3: note: use option -std=c99 or -std=gnu99 to compile your code

I am using BiocManager::install("DiffBind") to install the package "DiffBind" on CentOS_7, I got the error: ...
user avatar
  • 11
1 vote
1 answer
16 views

Find enriched processes in modules obtained by WGCNA

I am doing a gene enrichment of different clusters obtained by WGCNA, but there are some of them in which I don't get enriched processes, is there any way to find processes to these modules, for ...
user avatar
1 vote
1 answer
26 views

Can I Incorporate svaseq() into GSEA/GSVA analysis?

I understand GSEA/GSVA can take microarray-like expression matrix such as the output of voom() or vst(). However, I have a question on how we can also use svaseq() variables, correct the output from ...
user avatar
  • 13
1 vote
1 answer
27 views

Generate a single column from two different columns in the same format

I have multiple data frames in R in the following format: Dataframe1: ...
user avatar
  • 83
2 votes
0 answers
14 views

DeconRNASeq: Extract gene names from returned mixing proportions

This question was also asked on Biostars I am using the Bioconductor package "DeconRNASeq" to perform tissue deconvolution. Let's say I run the following code (this is from the manual): <...
user avatar
  • 21
1 vote
0 answers
22 views

Converting genetic matrix to STRUCTURE

I want to convert a genetic matrix (rows are samples, columns are alleles) into a STRUCTURE file (output .str) within R. Is there a preferred way to achieve this? I've found genind_to_structure in the ...
user avatar
0 votes
1 answer
22 views

Heatmap of motifs present in protein sequences

I have a file containing multiple protein amino acid sequences. I am interested in visually presenting the location of amino acid motifs PQG, QQG and AQG in the sequence in the form of a heatmap in R. ...
user avatar
  • 1
0 votes
0 answers
24 views

statistical test and visualization of genes present in an experiment

I have a list of genes (n = 120) that are involved in breast cancer. And I have a list of differentially expressed genes (n = 70) that are present in the breast cancer gene list. My question is: Are ...
user avatar
  • 75
1 vote
3 answers
58 views

beginner RNA-seq Replicate papers

have a good R and statistical analysis background (also with machine learning). in addition, I'm a fresh biotechnology grad. I would like to try to replicate some Rna-seq analysis with R papers (with ...
user avatar
  • 21
2 votes
1 answer
52 views

how to calculate the expression values per gene for all cells

I'm currently working with a seurat object and I'd like to calculate the expression values per gene for all cells within a particular cluster. I've gone ahead and subsetted the cluster of interest. <...
user avatar
  • 319
1 vote
1 answer
59 views

deseq2 full and reduced model interpretation

I would like to know how to interpret the output of the formula although i been thorough loads of literature but I'm not confident yet. So this is my full model I'm running ...
user avatar
  • 1,642
1 vote
0 answers
17 views

Is there an R package to calculate global FST, with confidence intervals, for tetraploid populations

I'm interested in estimating global FST among several populations of a tetraploid plant species. I'm looking for an R package that can do this, and can preferably also produce an upper and lower bound ...
user avatar
  • 11
1 vote
1 answer
25 views

Why all values become 1 after dcast?

I am reshaping my dataframe as follows. I want to set the Wavelength column as the column names with the code as follows. ...
user avatar
  • 81
1 vote
1 answer
26 views

how to create "sample file" for the qAlign() function after trimming the reads in R

I'm an absolute beginner trying to solve this question "Align the trimmed and untrimmed reads using QuasR and plot alignment statistics, did the trimming improve alignments?" I did trim the ...
user avatar
  • 21
1 vote
1 answer
28 views

Association analysis for methylation beta value and phenotypes in R

I am trying to learn how to do association analysis of methylation array data with numerical trait like body weight, eGFR. As I have seen some of the publications, methylation data association with ...
user avatar
1 vote
0 answers
13 views

What is the difference between AnnotationDbi/org.Mm.eg.db and biomaRt/Mus.musculus for converting to gene symbols?

I am interested in the differences between AnnotationDbi/org.Mm.eg.db and biomaRt/Mus.musculus. They yield the same results as seen in the code below. ...
user avatar
  • 153
1 vote
0 answers
23 views

Cytoscape for coloring the nodes in the network

I need to color different categories of the proteins in the interaction network automatically. Is there any way to color each category with a different color in cytoscape? or R? My input is 2 column ...
user avatar
1 vote
1 answer
39 views

How to build advanced boxplots from datasets which columns contain only numbers?

I have got a fold-change dataset and I simply would like to plot something like this: The problem is that I have only numbers in my columns and when I use a boxplot() function I can plot only boxes ...
user avatar
  • 103
2 votes
0 answers
130 views

How to manage memory contraints when analyzing a large number of gene count matrices? I keep running out of RAM with my current pipeline

I have several hundred scRNA-seq count matrices, each from a different sample. For my other dataset containg a few dozen samples, I simply merged everything together into one Seurat object, but that ...
user avatar
0 votes
2 answers
45 views

"Error in col2rgb(colors, alpha = TRUE) : invalid RGB specification" when running HeatmapAnnotation() (from ComplexHeatmap)

This question was also asked on github I'm trying to run the function HeatmapAnnotation (a function in the ComplexHeatmap package). When I run it, I get the error "Error in col2rgb(colors, alpha =...
user avatar
  • 329
2 votes
1 answer
45 views

baseMean threshold

I have an RNA-seq dataset and I am using DEseq2 to find differentially expressed genes between the two groups. I used pre-filtering to remove any genes that have ...
user avatar
1 vote
0 answers
94 views

Missing value imputation method for gene expression data

I am new to working with gene expression data sets and am wondering what is the most standard or best way to impute missing values in a gene expression data? I got mine from the GEO database and the ...
user avatar
  • 13
2 votes
0 answers
15 views

Subsetting with harmony

Currently merged two Seurat objects together and then ran Harmony for batch correction. Now I want to subset out a cell type of interest do I re run harmony or just do the standard PCA? ...
user avatar
  • 319
2 votes
0 answers
36 views

Subsetting out prefix from genes and merging

I'm working with PTX samples (human and mouse), they're aligned to both transcriptomes and so gene names contain hg38- and mm10- at the beginning of each gene. I have eliminated teh prefixes because I ...
user avatar
  • 319
5 votes
0 answers
41 views

Where to find asymmetric nucleotide substitution matrix with IUPAC encodings?

I am using the pairwiseAlignment function from the Biostrings package to calculate the distance between a consensus sequence and ...
user avatar
  • 443
0 votes
0 answers
32 views

Problem of writing scilab code on a project - Morris-Lecar model

This is a computational neuroscience project. We are asked to plot the Morris-Lescar model (question 3 below). The model should be plotted with scilab or R. however, I don't know where to start? For ...
user avatar
3 votes
1 answer
55 views

Grouping x-axis in scatter plot

I am trying to plot one scatter plot for each group, where each group has continous values in x-axis. Reference Data Gives continuous x-axis but all groups are combined. ...
user avatar
  • 345
1 vote
0 answers
18 views

open source tools for CNV analysis from illumina SNP array

What are the best open source tools to analyse copy number variation using Illumina SNP array data...preferably R based tools, but open to other languages as well. I'm working with leukemia sample and ...
user avatar
0 votes
0 answers
325 views

DimPlot - How to highlight cells with identity colors?

I made this wrapper for DimPlot: ...
user avatar
0 votes
0 answers
21 views

Error loading DiffBind sample sheet

I am trying to analyse a ChIP-Seq data set with DiffBind. These are the contents of my sample sheet ("diffbind_samples.csv"): ...
user avatar
3 votes
0 answers
24 views

Is there a generally accepted method used to impute missing DNA methylation data (probe beta-values)?

I'm looking at DNA methylation (DNAm) data such as TCGA (e.g., BRCA, KIRP, KIRC, etc.). Currently trying to build use my model to predict DNAm age on test sets, but many of the data sets are missing ...
user avatar
2 votes
0 answers
55 views

Pulling out genes in a scatterplot

I'm comparing gene expression among 2 different datasets (in vivo and in vitro) I have made a heatmap showing the correlation for each and then plotted the data frame to create a scatterplot. Now I ...
user avatar
  • 319
3 votes
0 answers
28 views

LRT test in deseq2 on multiple condition

I have some doubts when i'm running LRT This is my condition data-frame small subset ...
user avatar
  • 1,642
3 votes
1 answer
66 views

Constructing Two Triangle Heatmaps in One Square

I'm looking to make a two, two range colour scheme heatmap where each range corresponds to a separate triangle in the heatmap square. I am looking for a relevant guide for this heatmap. Any ...
user avatar
  • 121
3 votes
2 answers
93 views

Wrapping R scripts in WDL and run in docker

I am trying to run an R script as a task in WDL in a docker image. Since this is my first time working with these, I am stuck in something very simple (and in the initial stage) and am hoping somebody ...
user avatar
  • 71
1 vote
1 answer
75 views

A two color bar plot in R

I have raw read counts of miRNAs in two conditions like below ...
user avatar
  • 1,617
0 votes
2 answers
28 views

subsetting more than one ccle sample

I'm working with the CCLE dataset and I'm trying to subset the data to just 3 cell lines of interest but not sure how to it. Currently, I have been subsetting one by one but is there an easier way? <...
user avatar
  • 319
-1 votes
1 answer
81 views

Selecting a part of heatmap

I have a heatmap I see a block of genes for a group of samples How I can select a part of heatmap? I mean is there any way to export a list of genes from rows of a heatmap in R? I have tried which ...
user avatar
  • 1,617
3 votes
2 answers
111 views

How to subset a GRanges object based on a specific genomic window of interest?

Goal: To keep only rows of a GRanges object if their start coordinate is within 1 Mb of the start coordinate of a gene of interest. I am aware of some functions that come close: ...
user avatar
  • 115
3 votes
2 answers
232 views

Plotting a gene in Seurat

I saw in the extensive Seurat documentation for Dimplot (dimensional reduction plot), here, you can plot a gene by specifying it with group.by = "gene" but this does not work in practice. <...
user avatar
2 votes
0 answers
169 views

Featureplot by co-expression of some genes

I want to identify MDSCs cells which don't not express HLA-DRs (HLA-DR negative) but express ...
user avatar
  • 1,617
2 votes
1 answer
46 views

Displaying BIC and AICs from MSA nucleotide multisequence fasta file

I am initiating a project reproducing MEGA style trees using R to find the best DNA protein models under AIC/BIC. I have read the R AICs and BIC documentation but I'm struggling to implement this. ...
user avatar
4 votes
1 answer
51 views

Co-occurrence networks in Metagenomics studies

I have recently acquired some 16S metagenomics data, and was wondering if anyone can speak of the potential limitations, challenges as well as advantages to conducting a network-based study on ...
user avatar
  • 776
1 vote
0 answers
28 views

Separation of dihedral angles of the desired sequence in the protein

have a list of pdbid and I want to separate the dihedral angles of some sequences in protein . ...
user avatar
1 vote
0 answers
13 views

Different output on same command and same input file, different library version

I am running a variant calling pipeline to detect solid tumours and haematological alterations. The variant calling is done with Mutect2, whereas I am using an R script to perform some filtering on ...
user avatar
  • 213
3 votes
1 answer
225 views

In a RNA-Seq heatmap should you do Z-score standardisation before clustering the rows/columns or after?

I have made a heatmap using RPKM values from a RNA-Seq dataset using the pheatmap() function in R. I have log2-transformed the data before performing Z-score standardisation of the data. I have also ...
user avatar
  • 159
0 votes
1 answer
85 views

How to remedy a DESeq2 collapsing technical replicates error?

Goal: To ensure "the sum of the counts for [my samples] is the same as the counts in the [samples] columns in ddsColl" after collapsing replicates using ...
user avatar
  • 115
1 vote
0 answers
32 views

How can I get a column with real gene names in my ballgown analysis?

I am doing RNAseq analysis and I am using ballgown procedure for it. There is an option in R to calculate differentially expressed genes and use FPKM in calculating differential gene expression, which ...
user avatar
  • 103
0 votes
0 answers
26 views

Should I remove days_to_death = 0 before TCGA survival analysis

Some patients in TCGA data have days_to_death = 0 and days_to_follow_up = NA vice versa. Should I remove these samples before survival analysis? What does days_to_death = 0 mean? Does it mean the ...
user avatar
0 votes
0 answers
34 views

Combine different FASTQ belonging to different hypervariable 16s regions

I'm really newbie in 16s microbiome analyses, so sorry for this question. I would like to perform a meta-analysis using only Illumina Amplicons sequences from different studies that have different ...
user avatar
  • 101

1
2 3 4 5
16