[BioC] variation in output -using gcrma: bg correction
Sucheta Tripathy
sutripa at vbi.vt.edu
Mon Dec 11 18:02:24 CET 2006
Dear Group
when I run background(using GCRMA) correction for chips:
"10_1331.1I.R208L.CEL" "9_C101.1F.R9F2L.CEL" and for chips
"10_1195.1E.R119U.CEL" "9_C101.1F.R9F2L.CEL" keeping
9_C101.1F.R9F2L.CEL as constant, For both the runs I get slightly
different outputs for chip 9_C101.1F.R9F2L.CEL (Data below). This does not
happen when GSB.adjust is set to FALSE. Also we have noticed the variation
seem to get lesser when we try more chips together.Any reason why it
should be happening this way?
Many thanks
Sucheta
My code:
-------
obj<-ReadAffy(filenames=file)
tmp1<-bg.adjust.gcrma(obj,optical.correct=T,affinity.info=affinfo)
outfile<-paste(outdir, celfile[i,1], "_bg.txt", sep="")
write.table(pm(tmp1), file=outfile, row.names=FALSE, col.names=TRUE,
sep="\t")
----------------------------------------------------------------------
"10_1331.1I.R208L.CEL" "9_C101.1F.R9F2L.CEL"
---------------------------------------------
12.5854077145271 4.27009199115734
135.577788495271 115.177321956904
113.171763014501 124.781182045557
2.69508615100678 13.3527712251500
14.5462397159524 36.7662738574386
94.6480456343775 48.972931999473
230.750229599019 212.534925776574
242.058765502072 337.251490520153
233.871977591185 339.413283090951
134.355125363486 108.047618775561
136.685235380017 126.239703801863
76.2625488684543 45.3623158856175
148.251139760089 120.315265011794
328.121529159736 319.366509787047
291.865380794849 361.992073209429
3.97218796139000 4.20902848830535
14.6419514021341 31.6450267847893
28.8699339528603 67.3237575282198
110.496037483661 111.806067290232
22.2876523362755 41.5481196151883
205.194483921088 276.995029856902
84.8559677803056 96.0077121491932
778.222646971625 559.863168099135
4.83202857118797 11.3149613102092
3.66166154528273 39.6717762417948
67.7761672656053 94.6394189279102
112.500673213605 190.985739398352
171.007645212454 200.035151512867
492.236603432973 352.584270678829
130.800353580651 155.324195356517
146.303381286542 160.442554666605
4.34069757710996 61.4777011118527
88.9169964732945 84.2320613854983
69.1622655457769 60.1255829409642
382.758342400691 429.238548268896
3.90087434250685 4.02923006445781
245.047319536046 186.199916603738
4.60503013970163 16.7179032017834
45.8974167676779 58.6091870023334
152.737507677336 87.4646339988696
===================================================
"10_1195.1E.R119U.CEL" "9_C101.1F.R9F2L.CEL"
---------------------------------------------
34.9453663026134 4.35753552341506
120.903734762231 116.962904914229
111.293345773408 133.421207928372
2.96910233809292 14.4050686356057
34.7269027797637 35.8835240852081
85.3109066158579 45.2284068764267
222.081465596656 201.346914814921
378.991271676870 336.282180468649
324.355847430672 337.586327159885
148.252162519851 107.763579952789
146.760939437243 127.027138121574
68.9112298899552 45.5584187250637
178.300028533109 125.831740523771
457.092497175424 331.253094567131
349.374485133996 357.721065463285
8.39630366380937 4.30526177658911
50.2501294702745 31.1539961383794
51.2772385586404 66.8937884059372
169.473350033569 117.993461506227
47.5601861370703 41.0301951390256
238.568680512591 271.103448196271
122.668905036686 92.7991203760259
710.506330695352 568.473304793075
14.1193442910318 11.2664752571536
29.6858453280301 41.0915584209232
78.9811024872745 98.557906977124
143.137824860583 200.312672914717
228.952312559228 208.952150636992
414.93711559597 364.157958616779
193.573470074735 156.743939369768
151.710594739511 165.182915986475
16.3790389361525 62.1220611327726
63.6344501007354 81.4281142367824
52.6876908746722 58.4074520261634
356.890313613082 426.263965724419
4.0847500724139 4.13161532597561
233.32353554667 185.840713238473
4.71348339142599 16.7565881676473
89.4105877815801 55.7427731671095
162.930238504639 91.1574660327713
27.6720869896949 30.8456197258235
107.041950691333 17.5365838616139
132.491703246086 94.8472761626705
165.348305450273 118.806052283819
428.453296492502 438.259964898252
38.5832188166168 20.0643402347992
--
Sucheta Tripathy, Ph.D.
Virginia Bioinformatics Institute Phase-I
Washington street.
Virginia Tech.
Blacksburg,VA 24061-0447
phone:(540)231-8138
Fax: (540) 231-2606
web page: http://staff.vbi.vt.edu/sutripa
blog : http://genomics-array.blogspot.com/
quick ref:
http://faculty.ucr.edu/~tgirke/Documents/R_BioCond/R_BioCondManual.html
More information about the Bioconductor
mailing list