mean vs median of pixel values within spot mask wasRe: something weird with read.maimages and .gpr files
0
0
Entering edit mode
Pita ▴ 120
@pita-1011
Last seen 9.6 years ago
It seems that I answered my own question .... There was a typo that was entered into the Genepix software so the output would not parse. Limmas read.maimages looks specifically for the F532 Mean for the R channel, instead F594 Mean was written Since it is the case that the Mean of the pixels in the pot mask are used in limma as default: are they any arguments for using the mean values over the median values? Would the assumption here be that the spot morphology must be impeccable? Peter At 03:48 PM 2/25/2005, Pita wrote: >what can trigger the following error? I have imported many genepix >files without any problem > > >> RG.new.gpr <- read.maimages(files.new,"genepix",wt.fun=wtflags(0.1)) > >Error in "[.data.frame"(obj, , columns$Gf) : > > undefined columns selected > >This is with limma version 1.8.22 > > >Peter > >P.S. Here is a snippet of the first .gpr file: > >ATF 1.0 >29 56 >"Type=GenePix Results 3" >"DateTime=2005/01/26 14:13:02" >"Settings=" >"GalFile=C:\Axon\Params\ImmunArray2.gal" >"PixelSize=10" >"Wavelengths=635 594" >"ImageFiles=Y:\one\two\IMMUNOARRAY\26.01.05\13012006_SEKA0137-01_Cy5_ 26.01.05.tif.tif >0 Y:\one\two\IMMUNOARRAY\26.01.05\13012006_HURR_Cy3_26.01.05.tif 0" >"NormalizationMethod=None" >"NormalizationFactors=1 1" >"JpegImage=C:\Axon\Data\13012006_SEKA0137-01_Cy5_0026.jpg" >"StdDev=Type 1" >"RatioFormulations=W1/W2 (635/594)" >"FeatureType=Circular" >"Barcode=" >"BackgroundSubtraction=LocalFeature" >"ImageOrigin=0, 0" >"JpegOrigin=1190, 940" >"Creator=GenePix Pro 6.0.0.76" >"Scanner=4000 680056" >"FocusPosition=0" >"Temperature=0" >"LinesAveraged=1" >"Comment=" >"PMTGain=0 0" >"ScanPower=100 100" >"LaserPower=0 0" >"Filters= " >"ScanRegion=0,0,2022,3665" >"Supplier=" >"Block" >"Column" "Row" "Name" "ID" "X" "Y" "Dia." "F635 >Median" "F635 Mean" "F635 SD" "F635 CV" "B635" "B635 >Median" "B635 Mean" "B635 SD" "B635 CV" "% > >B635+1SD" "% > B635+2SD" "F635 % Sat." "F594 Median" "F594 >Mean" "F594 SD" "F594 CV" "B594" "B594 Median" "B594 >Mean" "B594 SD" "B594 CV" "% > B594+1SD" "% > >B594+2SD" "F594 % Sat." "Ratio of Medians (635/594)" "Ratio of Means >(635/594)" "Median of Ratios (635/594)" "Mean of Ratios (635/594)" >"Ratios SD (635/594)" "Rgn Ratio (635/594)" "Rgn R2 >(635/594)" "F Pixels" "B Pixels" "Circularity" "Sum of Medians >(635/594)" "Sum of Means (635/594)" "Log Ratio >(635/594)" "F635 Median - B635" "F594 Median - B594" "F635 Mean - >B635" "F594 Mean - B594" "F635 Total Intensity" "F594 Total >Intensity" "SNR 635" "SNR 594" "Flags" "Normalize" "Autoflag" >1 1 1 "DUSP2:Hs.1183:AA759046" "1321598" >1540 1160 100 2528 2522 1254 49 0 >0 90 154 171 100 97 0 1971 1990 1077 > 54 0 0 59 123 208 100 97 0 > 1.283 1.267 >1.363 1.339 2.275 1.272 0.575 80 546 100 4499 >4512 0.359 2528 1971 2522 >1990 201729 159237 15.792 15.699 0 0 0 >1 2 1 "DUSP2:Hs.1183:AA759046" "1321598" >1750 1170 120 2686 2625 1525 58 0 >0 84 174 207 97 95 0 1841 2014 1266 > 62 0 0 47 106 225 96 95 0 > 1.459 1.303 >1.352 1.417 2.660 1.326 0.630 120 693 100 4527 >4639 0.545 2686 1841 2625 >2014 314976 241629 14.603 18.557 0 0 0 >1 3 1 "MYL4:Hs.356717:AA705225" "461425" >1940 1160 110 7604 7657 2219 28 0 >0 77 142 184 100 100 0 4055 4284 1390 > 32 0 0 65 133 204 100 100 0 > 1.875 >1.787 1.684 1.811 1.524 2.029 0.741 80 572 100 >11659 11941 0.907 7604 4055 >7657 4284 612532 342716 53.380 31.722 0 0 0 >1 4 1 "MYL4:Hs.356717:AA705225" "461425" >2140 1160 110 6585 7117 2955 41 0 >0 84 167 198 100 100 0 3261 3319 1518 > 45 0 0 73 164 224 96 95 0 > 2.019 2.144 >2.060 2.423 2.458 2.355 0.736 80 620 100 9846 >10436 1.014 6585 3261 7117 >3319 569339 265548 42.114 19.793 0 0 0 >1 5 1 "PSMB8:Hs.180062:AA181300" "624360" >2350 1160 110 2433 2535 1160 45 0 >0 79 146 184 100 98 0 2484 2641 1292 > 48 0 0 83 219 263 97 97 0 > 0.979 0.960 >0.951 0.989 2.062 0.908 0.637 80 579 100 4917 >5176 -0.030 2433 2484 2535 >2641 202824 211288 16.822 11.680 0 0 0 >1 6 1 "PSMB8:Hs.180062:AA181300" "624360" >2540 1170 120 2155 2327 1336 57 0 >0 72 139 193 97 95 0 2260 2256 1207 > 53 0 0 53 129 243 95 92 0 > 0.954 1.031 >0.961 0.998 2.537 1.101 0.660 120 692 100 4415 >4583 -0.069 2155 2260 2327 >2256 279193 270688 16.223 17.078 0 0 0 >1 7 1 "RAD51C:Hs.412587:R37145" "26997" >2750 1160 90 1743 1788 796 44 0 0 >76 143 188 100 98 0 2250 2460 1181 48 > 0 0 66 161 243 100 98 0 0.775 > 0.727 >0.719 0.730 2.097 0.617 0.472 52 420 100 3993 >4248 -0.368 1743 2250 1788 >2460 92995 127944 11.972 14.870 0 0 0 >1 8 1 "RAD51C:Hs.412587:R37145" "26997" >2950 1160 100 1184 1237 725 58 0 0 >74 134 181 100 95 0 1938 1806 998 55 > 0 0 66 147 222 95 92 0 0.611 > 0.685 >0.579 0.692 2.521 0.599 0.484 80 540 100 3122 >3043 -0.711 1184 1938 1237 >1806 98992 144478 8.679 11.837 0 0 0 >1 9 1 "SHC1:Hs.433795:R52961" "40303" >3140 1160 90 2265 2434 1153 47 0 0 83 >160 192 100 100 0 2064 2021 1040 51 0 > 0 73 140 191 98 96 0 1.097 1.204 >1.232 1.289 2.207 1.187 0.553 52 414 100 4329 >4455 0.134 2265 2064 2434 >2021 126582 105071 14.694 13.914 0 0 0 >1 10 1 "SHC1:Hs.433795:R52961" "40303" >3350 1170 100 2040 2032 1088 53 0 0 >75 122 162 100 100 0 1771 1847 1021 55 > 0 0 61 119 195 98 96 0 1.152 > 1.100 >1.104 1.141 2.092 1.121 0.569 80 532 100 3811 >3879 0.204 2040 1771 2032 >1847 162549 147753 16.041 15.008 0 0 0 >1 11 1 "TYROBP:Hs.9963:H12338" >"148469" 3550 1170 110 15761 14143 5861 41 >21 21 106 172 162 100 100 0 744 884 > 592 66 0 0 42 101 240 91 86 >0 21.156 >15.975 17.541 19.010 2.619 18.419 0.594 80 590 100 >16484 15006 4.403 15740 744 14122 884 1131446 >70752 81.610 8.337 0 0 0 >1 12 1 "TYROBP:Hs.9963:H12338" >"148469" 3750 1170 110 15557 14415 6677 46 0 >0 103 178 172 100 100 0 865 916 602 > 65 0 0 47 131 278 90 80 0 > 17.985 15.737 >16.305 18.014 2.349 18.037 0.619 80 590 100 16422 >15331 4.169 15557 865 14415 916 1153181 >73243 80.404 6.634 0 0 0 > >_______________________________________________ >Bioconductor mailing list >Bioconductor@stat.math.ethz.ch >https://stat.ethz.ch/mailman/listinfo/bioconductor >
limma trigger limma trigger • 1.0k views
ADD COMMENT

Login before adding your answer.

Traffic: 1040 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6