Hi All (first timer here, be gentle!!) I'm trying to extract the residual variance out of a glm with poisson data to conduct power and sample size to detect a change over years. I have equal reps within year (for now). Is the residual variance just the grand mean of the original count data, or something trickier? glm.Male_glm(M~year*strata,family=poisson,offset=dur) levels(year) [1] "2000" "2001"> levels(strata)[1] "0-4" "4.1-8" "12.1-16" "8.1-12" Norm Good Fisheries Biologist (Mathematician) Southern Fisheries Centre PO Box 76, Deception Bay Q 4508 telephone 3817 9588, fascimile 3817 9555 ********************************DISCLAIMER**************************** The information contained in the above e-mail message or messages (which includes any attachments) is confidential and may be legally privileged. It is intended only for the use of the person or entity to which it is addressed. If you are not the addressee any form of disclosure, copying, modification, distribution or any action taken or omitted in reliance on the information is unauthorised. Opinions contained in the message(s) do not necessarily reflect the opinions of the Queensland Government and its authorities. If you received this communication in error, please notify the sender immediately and delete it from your computer system network. -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request at stat.math.ethz.ch _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._