On Fri, 2008-11-14 at 08:35 -0800, Robin Clark wrote: > Hi, > I am attempting to set the confidence limits on a ls means plot as follows: > > mult<-glht(lm(effectModel, data=statdata, na.action = na.omit), > linfct=mcp(mainEffect="Means")) > meanPlot <- sub(".html", "meanplot.jpg", htmlFile) > jpeg(meanPlot) > > plot(mult, main=NA, xlab=unlist(strsplit(Args[4],"~"))[1]) > > This produces 95% CIs by default but I would like to produce 99% CIs - How > do I do this?
plot(confint(mult, level = 0.99)) should do it if I understand both you and the help for multcomp correctly. HTH G > > Thanks, > Robin > -- %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% Dr. Gavin Simpson [t] +44 (0)20 7679 0522 ECRC, UCL Geography, [f] +44 (0)20 7679 0565 Pearson Building, [e] gavin.simpsonATNOSPAMucl.ac.uk Gower Street, London [w] http://www.ucl.ac.uk/~ucfagls/ UK. WC1E 6BT. [w] http://www.freshwaters.org.uk %~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~%~% ______________________________________________ R-help@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.