|
21
|
1 options(show.error.messages=F,warn=-1, error=function(){cat(geterrmessage(),file=stderr());q("no",1,F)})
|
|
0
|
2 sink(stdout(), type = "message")
|
|
21
|
3 sink(stdout() ,type = "output")
|
|
0
|
4 suppressWarnings(suppressMessages(library(affy)))
|
|
|
5 suppressWarnings(suppressMessages(library(affyPLM)))
|
|
|
6 suppressWarnings(suppressMessages(library(batch)))
|
|
|
7 suppressWarnings(suppressMessages(library(annotate)))
|
|
|
8 suppressWarnings(suppressMessages(library(limma)))
|
|
|
9 suppressWarnings(suppressMessages(library(markdown)))
|
|
|
10 suppressWarnings(suppressMessages(library(knitr)))
|
|
|
11 source_local <- function(fname){
|
|
|
12 argv <- commandArgs(trailingOnly = FALSE)
|
|
|
13 base_dir <- dirname(substring(argv[grep("--file=", argv)], 8))
|
|
|
14 source(paste(base_dir, fname, sep="/"))
|
|
|
15 }
|
|
|
16 file_path <- function(fname){
|
|
|
17 argv <- commandArgs(trailingOnly = FALSE)
|
|
|
18 base_dir <- dirname(substring(argv[grep("--file=", argv)], 8))
|
|
|
19 pato <- paste(base_dir, fname, sep="/")
|
|
|
20 return(pato)
|
|
|
21 }
|
|
|
22 base_dir <- function(){
|
|
|
23 argv <- commandArgs(trailingOnly = FALSE)
|
|
|
24 base_dir <- dirname(substring(argv[grep("--file=", argv)], 8))
|
|
|
25 return(base_dir)
|
|
|
26 }
|
|
|
27 source_local("Affymetrix_Preprocessing_Functions.R")
|
|
|
28 listArguments = parseCommandArgs(evaluate=FALSE)
|
|
|
29 thefunctions=listArguments[["thefunctions"]]
|
|
|
30 listArguments[["thefunctions"]]=NULL
|
|
|
31 h=listArguments[["h"]]
|
|
|
32 listArguments[["h"]]=NULL
|
|
|
33 w=listArguments[["w"]]
|
|
|
34 listArguments[["w"]]=NULL
|
|
|
35 if (!is.null(listArguments[["image"]])){
|
|
|
36 load(listArguments[["image"]])
|
|
|
37 listArguments[["image"]]=NULL
|
|
|
38 }
|
|
|
39 listArguments[["rawdata"]]=MicroArray_Object$affy_object
|
|
|
40 if(datasetsource=="intern"){
|
|
|
41 designo<-MicroArray_Object$designo
|
|
|
42 }
|
|
|
43
|
|
|
44 if(datasetsource=="extern"){
|
|
|
45 listArguments<-append(listArguments,list(datasetsource=datasetsource,listfullnames=listfullnames))
|
|
|
46 }
|
|
|
47 Prepro_object<-do.call(thefunctions,listArguments)
|
|
|
48 if(datasetsource=="extern"){
|
|
|
49 rownames(Prepro_object$data.bg@phenoData@data)<-listfullnames
|
|
|
50 colnames(exprs(Prepro_object$data.bg))<-listfullnames
|
|
|
51 colnames(exprs(Prepro_object$data.sm))<-listfullnames
|
|
|
52 colnames(exprs(Prepro_object$data.norm))<-listfullnames
|
|
|
53 colnames(exprs(Prepro_object$data.bg))<-listfullnames
|
|
|
54 colnames(exprs(Prepro_object$data.norm))<-listfullnames
|
|
|
55 rownames(Prepro_object$data.norm@phenoData@data)<-listfullnames
|
|
|
56 rownames(Prepro_object$data.norm@protocolData@data)<-listfullnames
|
|
|
57 }
|
|
|
58 par(las=2,mar=c(15,2,1,1))
|
|
|
59 png(filename ="boxplot_before_NM.png",width = w, height = h, units = "px", pointsize = 14, bg = "white")
|
|
|
60 par(las=2,mar=c(15,2,1,1))
|
|
|
61 boxplot((na.omit(as.data.frame(exprs(Prepro_object$data.bg)))), main="Boxplot of intensities before Normalization",col="red",las=2,mar=c(15,2,1,1))
|
|
|
62
|
|
|
63 dev.off()
|
|
|
64
|
|
|
65
|
|
|
66
|
|
|
67 par(las=2,mar=c(15,2,1,1))
|
|
|
68 png(filename ="boxplot_after_NM.png",width = w, height = h, units = "px", pointsize = 14, bg = "white")
|
|
|
69 par(las=2,mar=c(15,2,1,1))
|
|
|
70 boxplot((na.omit(as.data.frame(exprs(Prepro_object$data.sm)))), main="Boxplot of intensities After Normalization",col="red",ylab="(intensities)",las=2,mar=c(15,2,1,1))
|
|
|
71
|
|
|
72 dev.off()
|
|
|
73 png(filename ="MA_plot.png",width = w, height = h)
|
|
|
74 MAplot((Prepro_object$data.norm) ,
|
|
|
75 show.statistics = F, span = 2/3, family.loess = "gaussian",
|
|
|
76 cex = 2, plot.method = as.character("smoothScatter"),
|
|
|
77 azdd.loess = TRUE, lwd = 1, lty = 1, loess.col = "red")
|
|
|
78
|
|
|
79 dev.off()
|
|
|
80 png(filename = "densities_plot_before_NM.png",width = w, height = h)
|
|
|
81 plotDensities(exprs(Prepro_object$data.bg),log=T)
|
|
|
82 dev.off()
|
|
|
83
|
|
|
84 png(filename = "densities_plot_after_NM.png",width = w, height = h)
|
|
|
85 plotDensities(exprs(Prepro_object$data.norm),log=T)
|
|
|
86 dev.off()
|
|
|
87 AffymetrixRmd=file_path("Affymetrix.Rmd")
|
|
|
88 Style=file_path("look.css")
|
|
|
89 suppressWarnings(suppressMessages(knit2html(AffymetrixRmd,output="PreprocessingPlots.html",quiet = T, stylesheet=Style)))
|
|
|
90 rm(listArguments)
|
|
|
91 save.image("MicroArray.Preprocessing.RData")
|
|
21
|
92 sink()
|
|
|
93 sink()
|
