[R] error while plotting

[Uwe Ligges]

Perhaps you explain us how you really generate the data that are results 
from the multcomp package. Then it would probbaly be clear how to proceed.

Uwe Ligges



Nair, Murlidharan T wrote:
> I am trying to plot the confidence limits form multiple comparison analysis. How do I need to construct the object to plot it now. 
> 
> Thanks ../Murli
> 
> 
> -----Original Message-----
> From: Uwe Ligges [mailto:ligges at statistik.tu-dortmund.de] 
> Sent: Sunday, September 27, 2009 1:17 PM
> To: Nair, Murlidharan T
> Cc: r-help at r-project.org
> Subject: Re: [R] error while plotting
> 
> 
> 
> Nair, Murlidharan T wrote:
>> I am getting the following errors when I am trying to plot the data below. I cannot figure out the error.
>> Error in plot.window(...) : need finite 'xlim' values
>> In addition: Warning messages:
>> 1: In min(x) : no non-missing arguments to min; returning Inf
>> 2: In max(x) : no non-missing arguments to max; returning -Inf
>> 3: In min(x) : no non-missing arguments to min; returning Inf
>> 4: In max(x) : no non-missing arguments to max; returning -Inf
>>
>>
>> #I am using the following code
>> #======================================================
>> library(multcomp)
>> sig.data<-structure(list(X = 1:10, Cell.lines = structure(c(2L, 5L, 8L, 
>> 9L, 3L, 6L, 10L, 1L, 4L, 7L), .Label = c("T(70%)a-N(0%)c", "T(70%)a-N(0%)f", 
>> "T(70%)a-N(0%)i", "T(70%)c-N(0%)c", "T(70%)c-N(0%)f", "T(70%)c-N(0%)i", 
>> "T(80%)a-N(0%)c", "T(80%)a-N(0%)f", "T(90%)-N(0%)f", "T(90%)-N(0%)i"
>> ), class = "factor"), estimate = c(9859.74333, -5553.64802, 6227.17947, 
>> 8063.6472, 6548.86032, -8864.53103, 4752.7642, 9057.72021, -6355.67115, 
>> 5425.15635), lower = c(5560.57875, -9852.8126, 1928.01489, 3764.48262, 
>> 2249.69575, -13163.69561, 453.59962, 4758.55563, -10654.83573, 
>> 1125.99177), upper = c(14158.90791, -1254.48344, 10526.34405, 
>> 12362.81178, 10848.0249, -4565.36645, 9051.92877, 13356.88479, 
>> -2056.50657, 9724.32092), p.val.raw = c(1.15e-08, 5.78e-05, 1.36e-05, 
>> 3.21e-07, 6.91e-06, 6.97e-08, 0.000331, 4.87e-08, 1.04e-05, 7.63e-05
>> ), p.val.bon = c(2.66e-06, 0.0133, 0.00315, 7.41e-05, 0.0016, 
>> 1.61e-05, 0.0764, 1.13e-05, 0.0024, 0.0176), p.val.adj = c(2.65e-13, 
>> 0.000592, 2.82e-05, 9.72e-08, 6.56e-05, 8.76e-09, 0.0117, 6.22e-09, 
>> 6.44e-06, 0.000334)), .Names = c("X", "Cell.lines", "estimate", 
>> "lower", "upper", "p.val.raw", "p.val.bon", "p.val.adj"), class = "data.frame", row.names = c("T(70%)a-N(0%)f", 
>> "T(70%)c-N(0%)f", "T(80%)a-N(0%)f", "T(90%)-N(0%)f", "T(70%)a-N(0%)i", 
>> "T(70%)c-N(0%)i", "T(90%)-N(0%)i", "T(70%)a-N(0%)c", "T(70%)c-N(0%)c", 
>> "T(80%)a-N(0%)c"))
>>
>> rownames(sig.data)<-sig.data[,2]
>> my.hmtest <- structure(list(
>>   estimate = t(t(structure(sig.data[,"estimate"], .Names = rownames(sig.data)))),
>>   conf.int = sig.data[,4:5],
>>   ctype = "ABCC4-2007"),
>>   class = "hmtest")
>> par(mex=0.5) #This helps to accomodate the margins when text is getting cut off
>> plot(my.hmtest, cex.axis=0.7)
> 
> 
> There is not method plot.hmtest defined anywhere. Hence plot.default is 
> used and that one does not know hoe to handle an object like the one you 
> just defined.
> 
> Uwe Ligges
> 
> 
> 
>>  
>>
>> ______________________________________________
>> R-help at r-project.org mailing list
>> https://stat.ethz.ch/mailman/listinfo/r-help
>> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
>> and provide commented, minimal, self-contained, reproducible code.
> 
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.