This is an R Markdown document. Markdown is a simple formatting syntax for authoring HTML, PDF, and MS Word documents. For more details on using R Markdown see http://rmarkdown.rstudio.com.
install ecoliLeucine and preprocess it with RMA.
source("http://bioconductor.org/biocLite.R")
## Bioconductor version 2.14 (BiocInstaller 1.14.2), ?biocLite for
## help
biocLite(c("genefilter", "ecoliLeucine"))
## BioC_mirror: http://bioconductor.org
## Using Bioconductor version 2.14 (BiocInstaller 1.14.2), R version
## 3.1.1.
## Installing package(s) 'genefilter' 'ecoliLeucine'
## package 'genefilter' successfully unpacked and MD5 sums checked
## package 'ecoliLeucine' successfully unpacked and MD5 sums checked
##
## The downloaded binary packages are in
## C:\Users\knishida\AppData\Local\Temp\RtmpU7yE5H\downloaded_packages
## Warning: installed directory not writable, cannot update packages 'class',
## 'mgcv'
library("genefilter")
##
## Attaching package: 'genefilter'
##
## The following object is masked from 'package:base':
##
## anyNA
library("ecoliLeucine")
## Loading required package: affy
## Loading required package: BiocGenerics
## Loading required package: parallel
##
## Attaching package: 'BiocGenerics'
##
## The following objects are masked from 'package:parallel':
##
## clusterApply, clusterApplyLB, clusterCall, clusterEvalQ,
## clusterExport, clusterMap, parApply, parCapply, parLapply,
## parLapplyLB, parRapply, parSapply, parSapplyLB
##
## The following object is masked from 'package:stats':
##
## xtabs
##
## The following objects are masked from 'package:base':
##
## anyDuplicated, append, as.data.frame, as.vector, cbind,
## colnames, do.call, duplicated, eval, evalq, Filter, Find, get,
## intersect, is.unsorted, lapply, Map, mapply, match, mget,
## order, paste, pmax, pmax.int, pmin, pmin.int, Position, rank,
## rbind, Reduce, rep.int, rownames, sapply, setdiff, sort,
## table, tapply, union, unique, unlist
##
## Loading required package: Biobase
## Welcome to Bioconductor
##
## Vignettes contain introductory material; view with
## 'browseVignettes()'. To cite Bioconductor, see
## 'citation("Biobase")', and for packages 'citation("pkgname")'.
##
## Loading required package: ecolicdf
## Warning: replacing previous import by 'utils::head' when loading 'ecolicdf'
## Warning: replacing previous import by 'utils::tail' when loading 'ecolicdf'
##
data("ecoliLeucine")
eset = rma(ecoliLeucine)
## Background correcting
## Normalizing
## Calculating Expression
export rowttest result as csv (p.value < 0.05).
r = rowttests(eset, eset$strain)
filtered = r[r$p.value < 0.05,]
write.csv(filtered, file="ttest.csv")
head(filtered)
## statistic dm p.value
## aceA_b4015_st -2.665 -0.2314 0.03729
## alpA_b2624_st -2.594 -0.3141 0.04097
## amn_b1982_st 2.501 0.3008 0.04648
## ampE_b0111_st 3.227 0.4919 0.01797
## aroA_b0908_st -2.628 -0.6776 0.03916
## arp_b4017_st -3.098 -0.2036 0.02118