Tutorial: From Data import to Normalization in Microarray Analysis using in R (Part I)
4
gravatar for majuang66
4.3 years ago by
majuang6690
majuang6690 wrote:

I want to share my re-analysis and understanding of published microarray data of GEO database and others. I will update consistently the results...... and I can not write English well ^.^!!! Example 1: GSE4536 (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE4536)

# Open R
# Install the bioconductor packages, if not already installed
source("http://bioconductor.org/biocLite.R")
# Affymetrix Chip
biocLite("affy")

library(affy)
# Data Download for CEL files - I use Filezillla

Connect ftp.ncbi.nlm.nih.gov

Download the GSE4536 and unzip the GSE4536 zip file using the 7-zip

# Import into R from your local storage
setwd("e:/GSE4536 microarray")

mydata_GSE4536<-ReadAffy()
# Normalize to raw data (mydata_GSE4536)
# RMA normalization
mydata_GSE4536_rma<-rma(mydata_GSE4536)
# MAS5 normalization
mydata_GSE4536_mas5<-mas5(mydata_GSE4536)
# GC-RMA normalization
biocLite("gcrma")

library(gcrma)

mydata_GSE4536_gcrma<-gcrma(mydata_GSE4536)
# Farms normalization
biocLite("farms")

library(farms)

mydata_GSE4536_farms<-qFarms(mydata_GSE4536)
# dchip normalization
mydata_GSE4536_dchip<-expresso(mydata_GSE4536,normalize.method="invariantset",bg.correct=FASE,pmcorrect.method="pmonly",summary.method="liwong")

Reference

  1. http://repository.countway.harvard.edu/xmlui/bitstream/handle/10473/4710/Bioconductor_Manual.doc?sequence=1
  2. Using FARMS for summarization Using I/NI-calls for gene filtering
tutorial R gene • 3.6k views
ADD COMMENTlink modified 4.3 years ago by Devon Ryan96k • written 4.3 years ago by majuang6690

Your English is very good - do not worry.

ADD REPLYlink written 2.1 years ago by Kevin Blighe65k
1

Thank you for your reply. I have studied english everyday!

ADD REPLYlink written 2.1 years ago by majuang6690
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 975 users visited in the last hour