source("http://bioconductor.org/biocLite.R")
biocLite(c("GEOquery", "affy", "mouse4302.db"))

library(GEOquery)
library(affy)
library(mouse4302.db)

mouse.probe.keys <- mouse4302ENTREZID
mapped_probes <- mappedkeys(mouse.probe.keys)
probe2gene <- as.data.frame(mouse.probe.keys[mapped_probes])
write.csv(probe2gene, "probe2gene.csv")

geo.data <- getGEO("GSE46443")
geo.exprs <- exprs(geo.data[[1]])
write.csv(geo.exprs, "expressiondata.csv")

-- match rows which the value of GSM1130181 is larger than 14
GSM1130181 > 14

-- match rows which the value of GSM1130181 is larger than the value of GSM1130182
GSM1130181 > GSM1130182

-- match rows which the value of GSM1130181 or GSM1130182 is larger than 14
max(GSM1130181, GSM1130182) > 14

-- match rows which the value of GSM1130181 is larger than 12 and GSM1130182 is smaller than 11
GSM1130181 > 12 AND GSM1130182 < 11

-- match rows which the value of V starts with 143600
V LIKE "143600%"

SELECT * FROM probe2gene INNER JOIN expressiondata ON probe2gene.probe_id = expressiondata.V

SELECT gene_id, expressiondata.* FROM probe2gene INNER JOIN expressiondata ON probe2gene.probe_id = expressiondata.V

CREATE VIEW newview AS SELECT gene_id, expressiondata.* FROM probe2gene INNER JOIN expressiondata ON probe2gene.probe_id = expressiondata.V

SELECT gi.gene_info.Symbol, gi.gene_info.description, newview.* FROM newview INNER JOIN gi.gene_info ON newview.gene_id = gi.gene_info.GeneID

SELECT gi.gene_info.Symbol, gi.gene_info.description, GSM1130181 - GSM1130188 as "diff",newview.* FROM newview INNER JOIN gi.gene_info ON newview.gene_id = gi.gene_info.GeneID ORDER BY diff DESC