Using MAANOVA functions
2
0
Entering edit mode
@heather-maughan-1545
Last seen 10.2 years ago
Dear Bioconductor-users: I am using the package MAANOVA to analyze microarray data and have encountered problems when trying to plot data. I have tried emailing a MAANOVA discussion group, as well as the author of the package, and have not yet received a response so I am hoping that someone on this listserv can be of assistance. There are several functions in MAANOVA (riplot, resiplot) which call the "plot" function. For some unknown reason, when I use these functions, it appears that "xlim" has not been specified and I get an error (see below). However, when reading the code for each function, a command for how to calculate "xlim" has been specified but for some reason it does not get done. Does anyone have experience using MAANOVA commands and is willing to help? Here is an example of the error I get. It is identical also for the "resiplot" command: > riplot(spore) Error in plot.window(xlim, ylim, log, asp, ...) : need finite 'xlim' values Many thanks, Heather --
Microarray maanova Microarray maanova • 1.2k views
ADD COMMENT
0
Entering edit mode
@james-w-macdonald-5106
Last seen 2 hours ago
United States
Heather Maughan wrote: > Dear Bioconductor-users: > > I am using the package MAANOVA to analyze microarray data and have > encountered problems when trying to plot data. I have tried emailing a > MAANOVA discussion group, as well as the author of the package, and have not > yet received a response so I am hoping that someone on this listserv can be > of assistance. > > There are several functions in MAANOVA (riplot, resiplot) which call the > "plot" function. For some unknown reason, when I use these functions, it > appears that "xlim" has not been specified and I get an error (see below). > However, when reading the code for each function, a command for how to > calculate "xlim" has been specified but for some reason it does not get > done. Does anyone have experience using MAANOVA commands and is willing to > help? > > Here is an example of the error I get. It is identical also for the > "resiplot" command: > > >>riplot(spore) > > Error in plot.window(xlim, ylim, log, asp, ...) : > need finite 'xlim' values Looks like the calculation of xlim results in either an Inf or -Inf. Your best bet is to run riplot() under the debugger and see what is being passed to plot(). debug(riplot) riplot(spore) You will be able to step through the function and see what you get for max(I) and min(I), which may help figure out what is happening. HTH, Jim > > Many thanks, > Heather > -- > > _______________________________________________ > Bioconductor mailing list > Bioconductor at stat.math.ethz.ch > https://stat.ethz.ch/mailman/listinfo/bioconductor -- James W. MacDonald Affymetrix and cDNA Microarray Core University of Michigan Cancer Center 1500 E. Medical Center Drive 7410 CCGC Ann Arbor MI 48109 734-647-5623
ADD COMMENT
0
Entering edit mode
Hi Heather, >>Here is an example of the error I get. It is identical also for the >>"resiplot" command: >> >>>riplot(spore) >> >>Error in plot.window(xlim, ylim, log, asp, ...) : >> need finite 'xlim' values > > > Looks like the calculation of xlim results in either an Inf or -Inf. > Your best bet is to run riplot() under the debugger and see what is > being passed to plot(). > > debug(riplot) > riplot(spore) > > You will be able to step through the function and see what you get for > max(I) and min(I), which may help figure out what is happening. Another, quite powerful, possibility is to set options(error=recover) before running the command. When an error occurs, you will be presented with a menu of evaluation frames (this is R-speak for different subroutine levels), you can choose any of them and study the values of the variables just before the error occured. With "c" you can get back to the menu, and with "Q" you can exit. Best regards Wolfgang ------------------------------------- Wolfgang Huber European Bioinformatics Institute European Molecular Biology Laboratory Cambridge CB10 1SD England Phone: +44 1223 494642 Fax: +44 1223 494486 Http: www.ebi.ac.uk/huber
ADD REPLY
0
Entering edit mode
Xiangqin Cui ▴ 20
@xiangqin-cui-1550
Last seen 10.2 years ago
Hi Heather, Are you trying the riplot function on one-color data? It only works for two-color as far as I know from previous versions. I have not tried the new version yet. Otherwise, you may have some missing or strange entries in some cells. Good luck, Xiangqin -----Original Message----- From: bioconductor-bounces@stat.math.ethz.ch [mailto:bioconductor-bounces at stat.math.ethz.ch] On Behalf Of Heather Maughan Sent: Tuesday, December 20, 2005 5:19 PM To: bioconductor at stat.math.ethz.ch Subject: [BioC] Using MAANOVA functions Dear Bioconductor-users: I am using the package MAANOVA to analyze microarray data and have encountered problems when trying to plot data. I have tried emailing a MAANOVA discussion group, as well as the author of the package, and have not yet received a response so I am hoping that someone on this listserv can be of assistance. There are several functions in MAANOVA (riplot, resiplot) which call the "plot" function. For some unknown reason, when I use these functions, it appears that "xlim" has not been specified and I get an error (see below). However, when reading the code for each function, a command for how to calculate "xlim" has been specified but for some reason it does not get done. Does anyone have experience using MAANOVA commands and is willing to help? Here is an example of the error I get. It is identical also for the "resiplot" command: > riplot(spore) Error in plot.window(xlim, ylim, log, asp, ...) : need finite 'xlim' values Many thanks, Heather -- _______________________________________________ Bioconductor mailing list Bioconductor at stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/bioconductor
ADD COMMENT
0
Entering edit mode
Dear list, I have two sets of genes from the same experiment, > PinC Expression Set (exprSet) with 1310 genes 8 samples phenoData object with 2 variables and 8 cases varLabels FileName: read from file Target: read from file > PinS Expression Set (exprSet) with 2891 genes 8 samples phenoData object with 2 variables and 8 cases varLabels FileName: read from file Target: read from file How can I merge these two sets? I tried union() on two vectors created from the probe IDs but failed. Any hints? thanks in advance, David
ADD REPLY

Login before adding your answer.

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