[BioC] replace error reading 2color txt data into limma
Susan J. Miller
sjmiller at email.arizona.edu
Thu Oct 27 23:14:49 CEST 2005
I am hoping someone can tell me what is causing this problem. I am
using R 2.1.0, limma 2.0.2 on a Windows 2000 system. I have four .txt
files containing 2-color microarray data. If I read each individual txt
file it seems to work, but when I try to read all of the files I get the
error shown below:
> RG=read.maimages("C1_vs_O1.txt",columns=list(Rf="Spot Median
Intensity (w685)",Gf="Spot Median Intensity (w595)",Rb="Background
Median Intensity (w685)",Gb="Background Median Intensity (w595)"))
Read C1_vs_O1.txt
(seems OK)
> RG=read.maimages(ext="txt",columns=list(Rf="Spot Median Intensity
(w685)",Gf="Spot Median Intensity (w595)",Rb="Background Median
Intensity (w685)",Gb="Background Median Intensity (w595)"))
Read C1_vs_O1.txt
Error in "[<-"(`*tmp*`, , i, value = c(168, 175, 265, 313, 32495, 32475,
: number of items to replace is not a multiple of replacement length
The first four lines in the file C1_vs_O1.txt are:
Number Name Array/Plate Flagged As Control Spot Multi-Set Multi-Set Row
Multi-Set Column Set Set Row Set Column Row Column Comments Clone ID
UniGene cluster ID Gene symbol Locus ID Chromosome ID GenBank
accession number Tissue Gene map Gene name Organism Sequence Detection
Method (w595) Detection Method (w685) Detection Status (w595) Detection
Status (w685) X Centroid (pixels) (w595) X Centroid (pixels) (w685) Y
Centroid (pixels) (w595) Y Centroid (pixels) (w685) X Spot Deviation
(um) (w595) X Spot Deviation (um) (w685) Y Spot Deviation (um) (w595) Y
Spot Deviation (um) (w685) X Set Deviation (um) (w595) X Set Deviation
(um) (w685) Y Set Deviation (um) (w595) Y Set Deviation (um) (w685)
Major Extent (um) (w595) Major Extent (um) (w685) Minor Extent (um)
(w595) Minor Extent (um) (w685) Compactness (%) (w595) Compactness (%)
(w685) Threshold (w595) Threshold (w685) Spot Saturation (%) (w595) Spot
Saturation (%) (w685) Spot Mean Intensity (w595) Spot Mean Intensity
(w685) Spot Median Intensity (w595) Spot Median Intensity (w685) Spot
Standard Deviation (w595) Spot Standard Deviation (w685) Spot Number of
Pixels (w595) Spot Number of Pixels (w685) Spot Total Intensity (w595)
Spot Total Intensity (w685) Spot Peak Intensity (w595) Spot Peak
Intensity (w685) Background Mean Intensity (w595) Background Mean
Intensity (w685) Background Median Intensity (w595) Background Median
Intensity (w685) Background Standard Deviation (w595) Background
Standard Deviation (w685) Background Number of Pixels (w595) Background
Number of Pixels (w685) Spot Confidence (w595) Spot Confidence (w685)
Spot Normalized Intensity (w595) Spot Normalized Intensity (w685) Spot
Ratio (w595) Spot Ratio (w685) Spot SNR (w595) Spot SNR (w685) Spot CV
(%) (w595)
1 A1_AB IGERT_C1vsO1_11_10_04 NORMAL NO 1 1 1 A 1 1 1 1 None None
None None None None None None None None None None AUTO AUTO DETECTED
DETECTED 122.388885 122.388885 19.125 19.125 13.961801 13.961801
-7.780155 -7.780155 -9.756 -9.756 0 0 10.134623 10.134623 9.486059
9.486059 34.932919 34.932919 338 338 0 0 225.117279 117.413582 216
115 40.293694 25.080153 162 162 36469 19021 317 184 174.990311
86.96489 174 86 15.307359 12.753221 826 826 100 100 42 29 1.448276
0.690476 3.274698 2.38753 74.356911
2 A1_AB IGERT_C1vsO1_11_10_04 NORMAL NO 1 1 1 A 1 1 1 2 None None
None None None None None None None None None None AUTO AUTO DETECTED
DETECTED 150.583328 150.583328 19.666666 19.666666 17.027 17.027
-2.151189 -2.151189 -9.756 -9.756 0 0 10.197981 10.197981 9.409546
9.409546 46.076748 46.076748 377 377 0 0 252.658386 129.621124 239
124 46.958637 25.637608 161 161 40678 20869 404 229 197.486191
101.923172 197 103 19.706385 15.429781 833 833 100 100 42 21 2 0.5
2.799712 1.795097 77.255592
3 A2_AB IGERT_C1vsO1_11_10_04 NORMAL NO 1 1 1 A 1 1 1 3 None None
None None None None None None None None None None AUTO AUTO DETECTED
DETECTED 177.301987 177.301987 18.797029 18.797029 5.694523 5.694523
-10.29089 -10.29089 -9.756 -9.756 0 0 11.567095 11.567095 9.996995
9.996995 19.474066 19.474066 521 521 0 0 400.362488 219.306244 393
213 108.947014 62.532211 160 160 64058 35089 685 387 228.681213
121.866669 224 119 29.274176 18.376568 825 825 100 100 169 94
1.797872 0.556213 5.864598 5.302382 61.125092
If anyone can tell me how to get around this problem I would be grateful
for the help!
Thanks,
Susan J. Miller
Biotechnology Computing Facility
Arizona Research Laboratories
University of Arizona
Tucson, AZ 85721
(520) 626-2597
More information about the Bioconductor
mailing list