Help redarding Limma
Hello everyone
I use the following script to normalize the Affy data of GSE23343.
library(affy)
rawdata<-ReadAffy()
PM<-probes(rawdata,which="pm")
AffyInfo<-dimnames(PM)[[1]]
cutpos<-regexpr("\\d+$",AffyInfo,perl=T)
AffyID<-substr(AffyInfo,1,cutpos-1)
probe<-as.numeric(substr(AffyInfo,cutpos,nchar(AffyInfo)))
data.bgc<-bg.correct(rawdata,method="rma")
data.bgc.q<-normalize.AffyBatch.quantiles(data.bgc,type="pmonly")
pm.bgc.q<- probes(data.bgc.q,which="pm")
normalized<-cbind(AffyID,probe,pm.bgc.q)
write.table(normalized,file="NormalR.csv",sep=",",row.names=FALSE, quote=FALSE)
The data is normalized with this method
Can anybody help me to perform the Limma on this normalized data?
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Meta-Analysis of Microarray data (Afyymetrix)
written by reza
-
Reading in plain text gene expression and phenotypic data files to limma
written by StephanieK
-
Tutorial: sharing some naive codes for microarray normalization in R with whom are too new in R ali…
written by zizigolu
-
Tutorial: Machine Learning For Cancer Classification - Part 1 - Preparing The Data Sets
written by Obi Griffith
-
Extract Pre-Normalisation Expression Data From Probe Set Values
written by munch •
-
Tutorial: Analysing Microarray Data In Bioconductor
written by Obi Griffith
Please read the limma manual which contains extensive explanation and example code. Feel free then to come back with more specific questions.
Can you please check if the script I wrote is right for expression data?
Please read the manual, it covers everything you need.