search for: gasparrini

...a Poisson distribution with mean 800. y <- rpois(20,800) 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/p...

...hoice, but it's a little complicated. I've got some 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

...mer(Reaction ~ Days + (Days|Subject), sleepstudy)) (gm1 <- 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...

...ch software every time. Briefly: - Are you still planning 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/p...

...ge The procedure works, but the titles in the help pages 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

...found, using Font Substitution. Font cannot be embedded.]%% %%[ 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/ )

...to this help page doesn't work in the version installed 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

...to this help page doesn't work in the version installed 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

...solve(t(X)%*%W%*%X) %*% (t(X)%*%W%*%z) iter[i,] <- b } 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/ )

...old R was really useful. At the moment I re-installed the 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

...find a way to input simpler correlation structures (namely 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/p...

...QL' and 'glmer' expected? Which is more reliable? 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/ )

...guments are stored in list(...), but I haven't found a way to manage them, so the final call to somefunction() will include 'b' and 'c' twice if the user have included them in '...', and returns an error. Does anyone have some suggestion? Thanks for your time Antonio Gasparrini London School of Hygiene & Tropical Medicine Department of Social and Environmental Health Research 15-17 Tavistock Place, London WC1H 9SH, UK

...;s, including an detailed summary of the functions and some 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/p...

...meterized associations, among others. The package is fully documented through help pages. A package vignette will be hopefully added soon. 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 London School of Hygiene & Tropical Medicine Department of Social and Environmental Health Research 15-17 Tavistock Place, London WC1H 9SH, UK Office: 0044 (0)20 79272406 Mobile: 0044 (0)79 64925523 Skype contact: a.gasparrini http://www.lshtm.ac.uk/people/gasparrini.antonio _________________...

...rnorm(10) x2 <- rnorm(10) model <- lm(y ~ x1 + x2) 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/p...

...rted in: Peng et al., Model choice in time series studies 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...

...o search by: grep("mymatrix",names(coef(model))) 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/p...

...;s, including an detailed summary of the functions and some 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/p...

...meterized associations, among others. The package is fully documented through help pages. A package vignette will be hopefully added soon. 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 London School of Hygiene & Tropical Medicine Department of Social and Environmental Health Research 15-17 Tavistock Place, London WC1H 9SH, UK Office: 0044 (0)20 79272406 Mobile: 0044 (0)79 64925523 Skype contact: a.gasparrini http://www.lshtm.ac.uk/people/gasparrini.antonio _________________...