More posts like this
-
Error in eigen(if (doDykstra) R else Y, symmetric = TRUE) : infinite or missing values in 'x' when …
written by Agatha Nabilla L •Dear All, I am doing some spatial trasncriptomic analysis and when doing it on particular sample, I got this error message as the result: > …
-
Remove unknown batch effect in microarrays data
written by antonio.dos.hermanas •Hi everyone! Currently, I am working with microarray raw data. After background correction, normalization and summarization I realized that my data had a huge batch …
-
PCA on FPKM data
written by ipalmisa •Hi I am running a PCA analysis using prcomp function on my FPKM data from Rnaseq analysis. I am scaling the FPKM data by using …
-
How to remove batch effects using surrogate variables?
written by hansong798 •My methylation data has about 800 sample and 300,000 variables and I found 26 surrogate variables using 'DoISVA' package in r. The problem is that …
-
Remove batch effects in meta-analysis
written by Morteza Hadizadeh •Hi every one, In Meta-Analysis, after combining two data (expression matrix from two microarray study) with the following command, all<- merge(data1, data2, by= "Entrez") and …
-
seqlm installation error
written by niutsterHi, I am installing "seqlm" package in R 3.2.1 in windows OS, by these commands: library(devtools) install_github("raivokolde/seqlm") but i get this error: Error in curl::curl_fetch_disk(url, …
-
Remove batch effect for differential expression
written by RonHi, I want to do differential expression between two groups(two batches).For clustering analysis of these samples from two groups,I did batch normalization using COMBAT (sva …
-
SVA Analysis Error Question
written by questionaire •Hello All: I'm currently doing SVA analysis on my data set. I have gotten back the following error: Error in svd(res) : infinite or missing …
-
Error in installation R package
written by niutsterHi, I am trying to install "minfi" package with these commamnds: source("https://bioconductor.org/biocLite.R") biocLite("minfi"), but when R starts to install RSQLite package on Ubuntu workstation, i …
-
Batch ID in TCGA database
written by Na SedHi, I have downloaded some data from TCGA and I need to remove their batch effects using 'combat' algorithm from sva package, but 'combat' function …
Please include a reproducible example or at least the commands you used until you reached the error, otherwise the only think we can do is guess. Now this is a guess: do you have missing values in your data?
I am using SVA package on a methylation data file, the phenotype of interest was Mvalue(methylation value) that faced error in eigen, then to solve this error i changed my interested phenotype from Mvalue to Betavalue. In Betavalue the code was run succsessfully, this code was run: pheno=pData(gset) betavalue=getbeta(gset) mod=model.matrix(~as.factor(sample_group), data=pheno) mode0=model.matrix(~1,data=pheno) sva.results=sva(betavalue,mod,mod0) I am new to batch effects and related packages, in normalization the used package normalizes the data but in batch effect, only estimates the bach effects related to interested phenotype for example Betavalue but it is not able to correct it, could you guide me how can i remove these batch effects?
I would suggests you to read the sva vignette where how to correct for batch effect is explained in detail. In one approach that I prefer you obtain some surrogate variables that you can add to your linear model. But in the vignette there are several approaches described.