However, if package multcomp or multtest (in Bioconductor) suggested by Dieter Menne will work in your application, they would probably be better than Bonferroni. Bonferroni is like using a sledge hammer to open a can. The sledge hammer would work, but a can opener would be preferable. Bonferroni will work when other things won't, so I would use something like multcomp or multtest when I could and Bonferroni otherwise.
hope this helps. spencer graves
John Fox wrote:
Dear Cara,
You could use the confint() function, setting the level argument to 1 - alpha/length(coefficients(model)), where model is the linear model that you've fit and alpha is the complement of the level of confidence. If you're interested only in a subset of say p coefficients, then use 1 - alpha/p.
I hope this helps, John
--------------------------------
John Fox
Department of Sociology
McMaster University
Hamilton, Ontario
Canada L8S 4M4
905-525-9140x23604
http://socserv.mcmaster.ca/jfox --------------------------------
-----Original Message-----
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Cara Gormally
Sent: Thursday, March 10, 2005 8:41 PM
To: [email protected]
Subject: [R] Bonferroni simultaneous confidence intervals for multipleregression
Hi,
I'm having no luck figuring out how to find Bonferroni simultaneous confidence intervals to obtain a family of estimates in R. Does anyone know how to do this?
Thank you!
______________________________________________
[email protected] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
______________________________________________
[email protected] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
______________________________________________ [email protected] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
