From rpeng at stat.ucla.edu  Sun Mar  9 19:03:10 2003
From: rpeng at stat.ucla.edu (Roger Peng)
Date: Mon Mar 10 04:03:27 2003
Subject: [Rd] Problem with example(medpolish) in R-devel
Message-ID: <Pine.GSO.4.10.10303091858040.5748-100000@quetelet.stat.ucla.edu>

Today's R-devel has a problem in `example(medpolish)'.  It appears that
the `plot.medpolish' method is not available.  Here is the output and
traceback() after starting `R --vanilla':

> library(eda)
> example(medpolish)

mdplsh> deaths <- rbind(c(14, 15, 14), c(7, 4, 7), c(8, 2, 
    10), c(15, 9, 10), c(0, 2, 0))

mdplsh> dimnames(deaths) <- list(c("1-24", "25-74", "75-199", 
    "200++", "NA"), paste(1973:1975))

mdplsh> deaths
       1973 1974 1975
1-24     14   15   14
25-74     7    4    7
75-199    8    2   10
200++    15    9   10
NA        0    2    0

mdplsh> (med.d <- medpolish(deaths))
1 : 19 
Final: 19 

Median Polish Results (Dataset: "deaths")

Overall: 8 

Row Effects:
  1-24  25-74 75-199  200++     NA 
     6     -1      0      2     -8 

Column Effects:
1973 1974 1975 
   0   -1    0 

Residuals:
       1973 1974 1975
1-24      0    2    0
25-74     0   -2    0
75-199    0   -5    2
200++     5    0    0
NA        0    3    0


mdplsh> plot(med.d)
Error in plot.window(xlim, ylim, log, asp, ...) : 
        need finite xlim values
In addition: Warning messages: 
1: no finite arguments to min; returning Inf 
2: no finite arguments to max; returning -Inf 
3: no finite arguments to min; returning Inf 
4: no finite arguments to max; returning -Inf 
> traceback()
7: plot.window(xlim, ylim, log, asp, ...)
6: plot.default(med.d)
5: plot(med.d)
4: eval.with.vis(expr, envir, enclos)
3: eval.with.vis(ei, envir)
2: source(zfile, local, echo = echo, prompt.echo = prompt.echo, 
       verbose = verbose, max.deparse.length = 250)
1: example(medpolish)
> methods(plot)
 [1] "plot.POSIXct"        "plot.POSIXlt"        "plot.TukeyHSD"      
 [4] "plot.data.frame"     "plot.default"        "plot.density"       
 [7] "plot.design"         "plot.factor"         "plot.formula"       
[10] "plot.function"       "plot.hclust"         "plot.histogram"     
[13] "plot.isoreg"         "plot.lm"             "plot.mlm"           
[16] "plot.mts"            "plot.new"            "plot.ppr"           
[19] "plot.spec"           "plot.spec.coherency" "plot.spec.phase"    
[22] "plot.table"          "plot.ts"             "plot.window"        
[25] "plot.xy"            
> 
> version
         _                           
platform i686-pc-linux-gnu           
arch     i686                        
os       linux-gnu                   
system   i686, linux-gnu             
status   Under development (unstable)
major    1                           
minor    7.0                         
year     2003                        
month    03                          
day      09                          
language R                           


-roger
_______________________________
UCLA Department of Statistics
rpeng@stat.ucla.edu
http://www.stat.ucla.edu/~rpeng


[Rd] Problem with example(medpolish) in R-devel
Roger Peng