Wednesday, January 27, 2010

One example on how to upload and select data

## upload data
gefdata = read.table(file="gef-combo-commoneffects.txt",sep = ' ',header=TRUE)
docdata = read.table(file="doc-combo-commoneffects.txt",sep = ' ',header=TRUE)
d1= gefdata[(gefdata$Ecombo>gefdata$Egef&gefdata$Ecombo>2000)==TRUE,]
d2= d1[(d1$DoseGef=="20um")==TRUE,]
d3= d2[(d2$DoseDoc=="1.2um")==TRUE,]

> d3
Gene Egef Ecombo DoseGef DoseDoc
2074 STMN1 268.9287 2806.164 20um 1.2um
2094 AREG 2252.8792 2300.873 20um 1.2um
2175 ACTB 545.9648 3802.323 20um 1.2um
2229 ARHGEF6 2394.9913 2396.652 20um 1.2um
2263 UBB 1610.2985 5845.116 20um 1.2um
2265 ACTG1 261.1798 4023.568 20um 1.2um
2287 XRCC6 815.7438 2749.311 20um 1.2um
2289 KRT8 2966.2987 10506.174 20um 1.2um
2290 KRT18 2470.1614 10265.253 20um 1.2um
2291 NCOR2 681.8869 2255.428 20um 1.2um
2296 STATIP1 1418.3960 4557.809 20um 1.2um
2310 CALM2 481.1912 2865.798 20um 1.2um
2403 GAPDH 2970.3675 6717.153 20um 1.2um
2421 ATP1A1 725.7277 2798.023 20um 1.2um
2502 SLC9A3R1 1838.9419 6886.804 20um 1.2um
2521 HMGN2 117.6828 3845.859 20um 1.2um

No comments:

Post a Comment