limma error message
2
0
Entering edit mode
Jianping Jin ▴ 890
@jianping-jin-1212
Last seen 10.2 years ago
Dear BioConductor list, I got an error, Error: subscipt out of bounds, when I started to run a limma program normalizeWithinArrays(). What I wanted to do was to normalize eight cDNA slides. >dim(RG): 34970 8. I used >debug(normalizeWithinArrays) and then >normalizeWithinArrays(RG) to debug by either Browse[1]> n or c as well. But when I retried the program it gave me the same error message above. I will highly appreciate it if any one could give it shot. Best regards! JP- xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx x Jianping Jin Ph.D. x x Bioinformatics scientist x x Center for bioinformatics x x 3133 Bioinformatics Building x x CB# 7104 x x University of North Carolina x x Chapel Hill, NC 27599 x x Tel: (919)843-6105 x x Fax: (919)843-3103 x x E-mail: jjin at email.unc.edu x xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
• 1.0k views
ADD COMMENT
0
Entering edit mode
Naomi Altman ★ 6.0k
@naomi-altman-380
Last seen 3.6 years ago
United States
I got this message when I used the wrong layout file. --Naomi At 03:03 PM 7/28/2005, Jianping Jin wrote: >Dear BioConductor list, > >I got an error, Error: subscipt out of bounds, when I started to run a >limma program normalizeWithinArrays(). What I wanted to do was to normalize >eight cDNA slides. >dim(RG): 34970 8. > >I used >debug(normalizeWithinArrays) and then >normalizeWithinArrays(RG) to >debug by either Browse[1]> n or c as well. But when I retried the program >it gave me the same error message above. > >I will highly appreciate it if any one could give it shot. > > >Best regards! > >JP- > >xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx >x Jianping Jin Ph.D. x >x Bioinformatics scientist x >x Center for bioinformatics x >x 3133 Bioinformatics Building x >x CB# 7104 x >x University of North Carolina x >x Chapel Hill, NC 27599 x >x Tel: (919)843-6105 x >x Fax: (919)843-3103 x >x E-mail: jjin at email.unc.edu x >xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx > >_______________________________________________ >Bioconductor mailing list >Bioconductor at stat.math.ethz.ch >https://stat.ethz.ch/mailman/listinfo/bioconductor Naomi S. Altman 814-865-3791 (voice) Associate Professor Bioinformatics Consulting Center Dept. of Statistics 814-863-7114 (fax) Penn State University 814-865-1348 (Statistics) University Park, PA 16802-2111
ADD COMMENT
0
Entering edit mode
Thanks Naomi! Could you explain how to specify parameter "layout"? JP- --On Thursday, July 28, 2005 3:35 PM -0400 Naomi Altman <naomi at="" stat.psu.edu=""> wrote: > I got this message when I used the wrong layout file. > > --Naomi > > At 03:03 PM 7/28/2005, Jianping Jin wrote: >> Dear BioConductor list, >> >> I got an error, Error: subscipt out of bounds, when I started to run a >> limma program normalizeWithinArrays(). What I wanted to do was to >> normalize eight cDNA slides. >dim(RG): 34970 8. >> >> I used >debug(normalizeWithinArrays) and then >normalizeWithinArrays(RG) >> to debug by either Browse[1]> n or c as well. But when I retried the >> program it gave me the same error message above. >> >> I will highly appreciate it if any one could give it shot. >> >> >> Best regards! >> >> JP- >> >> xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx >> x Jianping Jin Ph.D. x >> x Bioinformatics scientist x >> x Center for bioinformatics x >> x 3133 Bioinformatics Building x >> x CB# 7104 x >> x University of North Carolina x >> x Chapel Hill, NC 27599 x >> x Tel: (919)843-6105 x >> x Fax: (919)843-3103 x >> x E-mail: jjin at email.unc.edu x >> xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx >> >> _______________________________________________ >> Bioconductor mailing list >> Bioconductor at stat.math.ethz.ch >> https://stat.ethz.ch/mailman/listinfo/bioconductor > > Naomi S. Altman 814-865-3791 (voice) > Associate Professor > Bioinformatics Consulting Center > Dept. of Statistics 814-863-7114 (fax) > Penn State University 814-865-1348 (Statistics) > University Park, PA 16802-2111 > xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx x Jianping Jin Ph.D. x x Bioinformatics scientist x x Center for bioinformatics x x 3133 Bioinformatics Building x x CB# 7104 x x University of North Carolina x x Chapel Hill, NC 27599 x x Tel: (919)843-6105 x x Fax: (919)843-3103 x x E-mail: jjin at email.unc.edu x xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
ADD REPLY
0
Entering edit mode
Jianping Jin ▴ 890
@jianping-jin-1212
Last seen 10.2 years ago
Dear Naomi and BioConductor list, Thanks for your kind reply! But I got the same error message after re- run the program as the following: > targets <- readTargets() > RG <- read.maimages(targets$FileName, source="genepix") Read 0341_3289_A5_052_HO34K_sample-1_HRef_scan2.gpr Read 0341_3290_A5_053_HO34K_sample-2_HRef_scan2.gpr Read 0341_3291_A5_054_HO34K_sample-3_HRef.gpr Read 0341_3292_A5_055_HO34K_sample-4_HRef.gpr Read 0341_3293_A5_056_HO34K_sample-5_HRef.gpr Read 0341_3294_A5_057_HO34K_sample-6_HRef.gpr Read 0341_3295_A5_058_HO34K_sample-7_HRef.gpr Read 0341_3296_A5_059_HO34K_sample-8_HRef.gpr > summary(RG) Length Class Mode R 279760 -none- numeric G 279760 -none- numeric Rb 279760 -none- numeric Gb 279760 -none- numeric targets 1 data.frame list genes 5 data.frame list > dim(RG) [1] 34970 8 > getLayout(RG$genes) $ngrid.r [1] 12 $ngrid.c [1] 4 $nspot.r [1] 28 $nspot.c [1] 28 attr(,"class") [1] "PrintLayout" > MA <- normalizeWithinArrays(RG, layout=list(ngrid.r=12, ngrid.c=4, nspot.r=28, nspot.c=28)) Error: subscript out of bounds I checked the data files and found nsopt.r and nspot.c in the last one third blocks of every files were changed from 28 to 24 and 23 respectively. That may be the reason for the error message. Do you have any ways to handle this? Many thanks! JP- xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx x Jianping Jin Ph.D. x x Bioinformatics scientist x x Center for bioinformatics x x 3133 Bioinformatics Building x x CB# 7104 x x University of North Carolina x x Chapel Hill, NC 27599 x x Tel: (919)843-6105 x x Fax: (919)843-3103 x x E-mail: jjin at email.unc.edu x xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx --On Thursday, July 28, 2005 5:17 PM -0400 Naomi Altman <naomi at="" stat.psu.edu=""> wrote: > layout=list(ngrid.r= ,ngrid.c= ,nspot.r= , nspot.c= ) > > Your array is laid out in blocks according to the print tip configuration > of the printer. > nspot.r is the number of rows in each block > nspot.c is the number of columns in each block > > ngrid.r is the number of rows of blocks > ngrid.c is the number of columns of blocks > > e.g. a typical 8 tip design might have 4 rows and 2 columns of blocks > > If you don't know the right numbers to put in, look at the raw data - the > block number and row and column are output to your expression file. > > --Naomi > > At 04:49 PM 7/28/2005, you wrote: >> Thanks Naomi! Could you explain how to specify parameter "layout"? >> JP- >> >> --On Thursday, July 28, 2005 3:35 PM -0400 Naomi Altman >> <naomi at="" stat.psu.edu=""> wrote: >> >>> I got this message when I used the wrong layout file. >>> >>> --Naomi >>> >>> At 03:03 PM 7/28/2005, Jianping Jin wrote: >>>> Dear BioConductor list, >>>> >>>> I got an error, Error: subscipt out of bounds, when I started to run a >>>> limma program normalizeWithinArrays(). What I wanted to do was to >>>> normalize eight cDNA slides. >dim(RG): 34970 8. >>>> >>>> I used >debug(normalizeWithinArrays) and then >>>> >normalizeWithinArrays(RG) to debug by either Browse[1]> n or c as >>>> well. But when I retried the program it gave me the same error message >>>> above. >>>> >>>> I will highly appreciate it if any one could give it shot. >>>> >>>> >>>> Best regards! >>>> >>>> JP- >>>> >>>> xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx >>>> x Jianping Jin Ph.D. x >>>> x Bioinformatics scientist x >>>> x Center for bioinformatics x >>>> x 3133 Bioinformatics Building x >>>> x CB# 7104 x >>>> x University of North Carolina x >>>> x Chapel Hill, NC 27599 x >>>> x Tel: (919)843-6105 x >>>> x Fax: (919)843-3103 x >>>> x E-mail: jjin at email.unc.edu x >>>> xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx >>>> >>>> _______________________________________________ >>>> Bioconductor mailing list >>>> Bioconductor at stat.math.ethz.ch >>>> https://stat.ethz.ch/mailman/listinfo/bioconductor >>> >>> Naomi S. Altman 814-865-3791 (voice) >>> Associate Professor >>> Bioinformatics Consulting Center >>> Dept. of Statistics 814-863-7114 (fax) >>> Penn State University 814-865-1348 (Statistics) >>> University Park, PA 16802-2111 >> >> >> >> xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx >> x Jianping Jin Ph.D. x >> x Bioinformatics scientist x >> x Center for bioinformatics x >> x 3133 Bioinformatics Building x >> x CB# 7104 x >> x University of North Carolina x >> x Chapel Hill, NC 27599 x >> x Tel: (919)843-6105 x >> x Fax: (919)843-3103 x >> x E-mail: jjin at email.unc.edu x >> xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx > > Naomi S. Altman 814-865-3791 (voice) > Associate Professor > Bioinformatics Consulting Center > Dept. of Statistics 814-863-7114 (fax) > Penn State University 814-865-1348 (Statistics) > University Park, PA 16802-2111 >
ADD COMMENT
0
Entering edit mode
28*28*12*4 [1] 37632 But you have 34970 genes. The simplest thing to do is to add the missing lines to your GPR files. You could probably do this in your favorite text manipulation language but with 8 arrays, you might be able to do it manually more quickly. --Naomi At 12:13 PM 7/29/2005, Jianping Jin wrote: >Dear Naomi and BioConductor list, > >Thanks for your kind reply! But I got the same error message after re-run >the program as the following: > >>targets <- readTargets() >>RG <- read.maimages(targets$FileName, source="genepix") >Read 0341_3289_A5_052_HO34K_sample-1_HRef_scan2.gpr >Read 0341_3290_A5_053_HO34K_sample-2_HRef_scan2.gpr >Read 0341_3291_A5_054_HO34K_sample-3_HRef.gpr >Read 0341_3292_A5_055_HO34K_sample-4_HRef.gpr >Read 0341_3293_A5_056_HO34K_sample-5_HRef.gpr >Read 0341_3294_A5_057_HO34K_sample-6_HRef.gpr >Read 0341_3295_A5_058_HO34K_sample-7_HRef.gpr >Read 0341_3296_A5_059_HO34K_sample-8_HRef.gpr >>summary(RG) > Length Class Mode >R 279760 -none- numeric >G 279760 -none- numeric >Rb 279760 -none- numeric >Gb 279760 -none- numeric >targets 1 data.frame list >genes 5 data.frame list > >>dim(RG) >[1] 34970 8 > >>getLayout(RG$genes) >$ngrid.r >[1] 12 > >$ngrid.c >[1] 4 > >$nspot.r >[1] 28 > >$nspot.c >[1] 28 > >attr(,"class") >[1] "PrintLayout" > >>MA <- normalizeWithinArrays(RG, layout=list(ngrid.r=12, ngrid.c=4, >nspot.r=28, nspot.c=28)) >Error: subscript out of bounds > >I checked the data files and found nsopt.r and nspot.c in the last one >third blocks of every files were changed from 28 to 24 and 23 >respectively. That may be the reason for the error message. Do you have >any ways to handle this? > >Many thanks! > >JP- > >xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx >x Jianping Jin Ph.D. x >x Bioinformatics scientist x >x Center for bioinformatics x >x 3133 Bioinformatics Building x >x CB# 7104 x >x University of North Carolina x >x Chapel Hill, NC 27599 x >x Tel: (919)843-6105 x >x Fax: (919)843-3103 x >x E-mail: jjin at email.unc.edu x >xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx > >--On Thursday, July 28, 2005 5:17 PM -0400 Naomi Altman ><naomi at="" stat.psu.edu=""> wrote: > >>layout=list(ngrid.r= ,ngrid.c= ,nspot.r= , nspot.c= ) >> >>Your array is laid out in blocks according to the print tip configuration >>of the printer. >>nspot.r is the number of rows in each block >>nspot.c is the number of columns in each block >> >>ngrid.r is the number of rows of blocks >>ngrid.c is the number of columns of blocks >> >>e.g. a typical 8 tip design might have 4 rows and 2 columns of blocks >> >>If you don't know the right numbers to put in, look at the raw data - the >>block number and row and column are output to your expression file. >> >>--Naomi >> >>At 04:49 PM 7/28/2005, you wrote: >>>Thanks Naomi! Could you explain how to specify parameter "layout"? >>>JP- >>> >>>--On Thursday, July 28, 2005 3:35 PM -0400 Naomi Altman >>><naomi at="" stat.psu.edu=""> wrote: >>> >>>>I got this message when I used the wrong layout file. >>>> >>>>--Naomi >>>> >>>>At 03:03 PM 7/28/2005, Jianping Jin wrote: >>>>>Dear BioConductor list, >>>>> >>>>>I got an error, Error: subscipt out of bounds, when I started to run a >>>>>limma program normalizeWithinArrays(). What I wanted to do was to >>>>>normalize eight cDNA slides. >dim(RG): 34970 8. >>>>> >>>>>I used >debug(normalizeWithinArrays) and then >>>>> >normalizeWithinArrays(RG) to debug by either Browse[1]> n or c as >>>>>well. But when I retried the program it gave me the same error message >>>>>above. >>>>> >>>>>I will highly appreciate it if any one could give it shot. >>>>> >>>>> >>>>>Best regards! >>>>> >>>>>JP- >>>>> >>>>>xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx >>>>>x Jianping Jin Ph.D. x >>>>>x Bioinformatics scientist x >>>>>x Center for bioinformatics x >>>>>x 3133 Bioinformatics Building x >>>>>x CB# 7104 x >>>>>x University of North Carolina x >>>>>x Chapel Hill, NC 27599 x >>>>>x Tel: (919)843-6105 x >>>>>x Fax: (919)843-3103 x >>>>>x E-mail: jjin at email.unc.edu x >>>>>xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx >>>>> >>>>>_______________________________________________ >>>>>Bioconductor mailing list >>>>>Bioconductor at stat.math.ethz.ch >>>>>https://stat.ethz.ch/mailman/listinfo/bioconductor >>>> >>>>Naomi S. Altman 814-865-3791 (voice) >>>>Associate Professor >>>>Bioinformatics Consulting Center >>>>Dept. of Statistics 814-863-7114 (fax) >>>>Penn State University 814-865-1348 (Statistics) >>>>University Park, PA 16802-2111 >>> >>> >>> >>>xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx >>>x Jianping Jin Ph.D. x >>>x Bioinformatics scientist x >>>x Center for bioinformatics x >>>x 3133 Bioinformatics Building x >>>x CB# 7104 x >>>x University of North Carolina x >>>x Chapel Hill, NC 27599 x >>>x Tel: (919)843-6105 x >>>x Fax: (919)843-3103 x >>>x E-mail: jjin at email.unc.edu x >>>xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx >> >>Naomi S. Altman 814-865-3791 (voice) >>Associate Professor >>Bioinformatics Consulting Center >>Dept. of Statistics 814-863-7114 (fax) >>Penn State University 814-865-1348 (Statistics) >>University Park, PA 16802-2111 > > > > Naomi S. Altman 814-865-3791 (voice) Associate Professor Bioinformatics Consulting Center Dept. of Statistics 814-863-7114 (fax) Penn State University 814-865-1348 (Statistics) University Park, PA 16802-2111
ADD REPLY

Login before adding your answer.

Traffic: 979 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