search for: pehru

...e problems with the offset and I don't understand some of the functions. Do you know where can I find some explicit tutorial about a poisson regression with MLE or, maybe better, about maximum likelyhood estimation? Best wishes Antonio Gasparrini Public and Environmental Health Research Unit (PEHRU) London School of Hygiene & Tropical Medicine Keppel Street, London WC1E 7HT, UK Office: 0044 (0)20 79272406 Mobile: 0044 (0)79 64925523 www.lshtm.ac.uk/pehru/ antonio.gasparrini at lshtm.ac.uk

...00) model <- glm(y ~ 1, family=poisson()) logLik(model) I would like to know what's the exact formula the function logLik uses. I looked at ?extractAIC but I cannot sort it out. Can you please give me an help? Thanks a lot Antonio Gasparrini Public and Environmental Health Research Unit (PEHRU) London School of Hygiene & Tropical Medicine Keppel Street, London WC1E 7HT, UK Office: 0044 (0)20 79272406 - Mobile: 0044 (0)79 64925523 Skype contact: a.gasparrini http://www.lshtm.ac.uk/people/gasparrini.antonio ( http://www.lshtm.ac.uk/pehru/ )

...t;- glmer(cbind(incidence, size - incidence) ~ period + (1 | herd), family = binomial, data = cbpp)) Taking into account that this package is used to perform similar analyses, this could be a problem. All the best Antonio Gasparrini Public and Environmental Health Research Unit (PEHRU) London School of Hygiene & Tropical Medicine Keppel Street, London WC1E 7HT, UK Office: 0044 (0)20 79272406 - Mobile: 0044 (0)79 64925523 http://www.lshtm.ac.uk/people/gasparrini.antonio ( http://www.lshtm.ac.uk/pehru/ ) -------------- next part -------------- An embedded message was scrubbed....

...% %%[ Warning: Helvetica-Bold not found, using Font Substitution. Font cannot be embedded.]%% If I open it, I can see that it has changed the font, and often the text doesn't match the boundaries. Can you help me? Thanks a lot Antonio Gasparrini Public and Environmental Health Research Unit (PEHRU) London School of Hygiene & Tropical Medicine Keppel Street, London WC1E 7HT, UK Office: 0044 (0)20 79272406 - Mobile: 0044 (0)79 64925523 http://www.lshtm.ac.uk/people/gasparrini.antonio ( http://www.lshtm.ac.uk/pehru/ )

...g this extension? And in this case, do you have a planned date? - What are likely to be the advantages and limitations of a potential 'rma.multi' if compared to Stata's 'mvmeta'? Thank you for your time Regards, Antonio Gasparrini Public and Environmental Health Research Unit (PEHRU) London School of Hygiene & Tropical Medicine Keppel Street, London WC1E 7HT, UK Office: 0044 (0)20 79272406 - Mobile: 0044 (0)79 64925523 Skype contact: a.gasparrini http://www.lshtm.ac.uk/people/gasparrini.antonio ( http://www.lshtm.ac.uk/pehru/ )

...} summary(model1,dispersion=1)$coef summary(model2) iter[nrow(temp),] sqrt(diag(solve(t(X)%*%W%*%X))) ####################################################### Can you explain if this difference is due to an error? Thanks in advance Antonio Gasparrini Public and Environmental Health Research Unit (PEHRU) London School of Hygiene & Tropical Medicine Keppel Street, London WC1E 7HT, UK Office: 0044 (0)20 79272406 - Mobile: 0044 (0)79 64925523 http://www.lshtm.ac.uk/people/gasparrini.antonio ( http://www.lshtm.ac.uk/pehru/ )

...le? 2) Is there a way to set the parameter 'nAGO' to 7 or perform the same analysis in another way? 3) Is there a tutorial on the use of lme4, especially to handle the results (summary, coef, etc.)? Thank you for your time Antonio Gasparrini Public and Environmental Health Research Unit (PEHRU) London School of Hygiene & Tropical Medicine Keppel Street, London WC1E 7HT, UK Office: 0044 (0)20 79272406 - Mobile: 0044 (0)79 64925523 http://www.lshtm.ac.uk/people/gasparrini.antonio ( http://www.lshtm.ac.uk/pehru/ )

...he old R 2.9.2, but I hope this will change in later versions. Or at least, that the user could have the opportunity to choose the help pages he deems more suitable for his purposes. I hope someone could give some details Regards, Antonio Gasparrini Public and Environmental Health Research Unit (PEHRU) London School of Hygiene & Tropical Medicine Keppel Street, London WC1E 7HT, UK Office: 0044 (0)20 79272406 - Mobile: 0044 (0)79 64925523

...e examples of application to real data, in order to get new users started easily. I hope that this package will be useful to your work. Any kind of feedback (questions, suggestions, bug-reports, etc.) is appreciated. Sincerely, Antonio Gasparrini Public and Environmental Health Research Unit (PEHRU) London School of Hygiene & Tropical Medicine Keppel Street, London WC1E 7HT, UK Office: 0044 (0)20 79272406 - Mobile: 0044 (0)79 64925523 Skype contact: a.gasparrini http://www.lshtm.ac.uk/people/gasparrini.antonio ( http://www.lshtm.ac.uk/pehru/ ) [[alternative HTML version deleted]] _____...

...of the created package lose the usual colours and font. Conversely, the result is perfect using R 2.8.1 (changing the path) and the same version of Rtools (.29). Not a big issue, but maybe it's worth considering it. Regards, Antonio Gasparrini Public and Environmental Health Research Unit (PEHRU) London School of Hygiene & Tropical Medicine Keppel Street, London WC1E 7HT, UK Office: 0044 (0)20 79272406 - Mobile: 0044 (0)79 64925523 Skype contact: a.gasparrini

...predict(model,type="terms")[1,] cbind(x1,x2)[1,] * t(coef(model)[2:3]) I guess the true result should be the second one. Maybe the command include some additional transformation? Any suggestion? Thanks for your time Antonio Gasparrini Public and Environmental Health Research Unit (PEHRU) London School of Hygiene & Tropical Medicine Keppel Street, London WC1E 7HT, UK Office: 0044 (0)20 79272406 - Mobile: 0044 (0)79 64925523 Skype contact: a.gasparrini http://www.lshtm.ac.uk/people/gasparrini.antonio ( http://www.lshtm.ac.uk/pehru/ )

...s os air pollution and mortality. J R Stat Soc A, 2006; 162: pag 190. Unfortunately, the function logLik doesn't work for a quasipoisson link. Do you know a fast method to extract the AIC for these models? Thanks in advance Antonio Gasparrini Public and Environmental Health Research Unit (PEHRU) London School of Hygiene & Tropical Medicine Keppel Street, London WC1E 7HT, UK Office: 0044 (0)20 79272406 - Mobile: 0044 (0)79 64925523 http://www.lshtm.ac.uk/people/gasparrini.antonio ( http://www.lshtm.ac.uk/pehru/ ) [[alternative HTML version deleted]]

...)) but I don't know how to extract the name "mymatrix" from the object mymatrix. Is there a specific function, something like func(mymatrix) returning the string "mymatrix"? Thanks so much for your time. Antonio Gasparrini Public and Environmental Health Research Unit (PEHRU) London School of Hygiene & Tropical Medicine Keppel Street, London WC1E 7HT, UK Office: 0044 (0)20 79272406 - Mobile: 0044 (0)79 64925523 Skype contact: a.gasparrini http://www.lshtm.ac.uk/people/gasparrini.antonio ( http://www.lshtm.ac.uk/pehru/ ) [[alternative HTML version deleted]]

...e examples of application to real data, in order to get new users started easily. I hope that this package will be useful to your work. Any kind of feedback (questions, suggestions, bug-reports, etc.) is appreciated. Sincerely, Antonio Gasparrini Public and Environmental Health Research Unit (PEHRU) London School of Hygiene & Tropical Medicine Keppel Street, London WC1E 7HT, UK Office: 0044 (0)20 79272406 - Mobile: 0044 (0)79 64925523 Skype contact: a.gasparrini http://www.lshtm.ac.uk/people/gasparrini.antonio ( http://www.lshtm.ac.uk/pehru/ ) [[alternative HTML version deleted]] _____...

...y diagonal or identity), in the same way as the using the functions pdDiag or pdIdent with 'glmmPQL'. I wonder if this option is still to be implemented in lme4. In this case, any suggestion/comment? Thanks for your time Antonio Gasparrini Public and Environmental Health Research Unit (PEHRU) London School of Hygiene & Tropical Medicine Keppel Street, London WC1E 7HT, UK Office: 0044 (0)20 79272406 - Mobile: 0044 (0)79 64925523 Skype contact: a.gasparrini http://www.lshtm.ac.uk/people/gasparrini.antonio ( http://www.lshtm.ac.uk/pehru/ )

...directory and modifying the help files, I cannot understand where I should run the commands R CMD INSTALL, check and build. Sorry for the stupid question, but I've been stuck for weeks with this problem. Thanks for your help Antonio Gasparrini Public and Environmental Health Research Unit (PEHRU) London School of Hygiene & Tropical Medicine Keppel Street, London WC1E 7HT, UK Office: 0044 (0)20 79272406 - Mobile: 0044 (0)79 64925523

Dear all, I need to fit a gee model with an auto-regressive correlation structure and I faced some problems. I attach a simple example: ####################################################### library(gee) library(geepack) # I SIMULATE DATA FROM POISSON DISTRIBUTION, 10 OBS FOR EACH OF 50 GROUPS set.seed(1) y <- rpois(500,50) x <- rnorm(500) id <- rep(1:50,each=10) # EXAMPLES FOR

...d by the CRAN (eg. try ?mkbasis). I checked the help pages for internal functions in other packages: some of them work (e.g. ?geepack::anova.geeglm), others don't (e.g. ?gam::.First.lib). Any suggestions? Thanks for your time Antonio Gasparrini Public and Environmental Health Research Unit (PEHRU) London School of Hygiene & Tropical Medicine Keppel Street, London WC1E 7HT, UK Office: 0044 (0)20 79272406 - Mobile: 0044 (0)79 64925523 Skype contact: a.gasparrini

...d by the CRAN (eg. try ?mkbasis). I checked the help pages for internal functions in other packages: some of them work (e.g. ?geepack::anova.geeglm), others don't (e.g. ?gam::.First.lib). Any suggestions? Thanks for your time Antonio Gasparrini Public and Environmental Health Research Unit (PEHRU) London School of Hygiene & Tropical Medicine Keppel Street, London WC1E 7HT, UK Office: 0044 (0)20 79272406 - Mobile: 0044 (0)79 64925523 Skype contact: a.gasparrini