[R] simulate() and glm fits

] base =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-= Dr William Valdar ++44 (0)1865 287 589 Wellcome Trust Centre [EMAIL PROTECTED] for Human Genetics, Oxford www.well.ox.ac.uk/~valdar __ R-help@stat.math.ethz.ch

Re: [R] simulate() and glm fits

2007, William Valdar wrote: Dear All, I have been trying to simulate data from a fitted glm using the simulate() function (version details at the bottom). This works for lm() fits and even for lmer() fits (in lme4). However, for glm() fits its output does not make sense to me -- am I missing

[R] get ylim of current plot()

- xyplot(rnorm(100), rnorm(100))) print(obj$y.limits) Thanks in advance, Will =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-= Dr William Valdar ++44 (0)1865 287 589 Wellcome Trust Centre [EMAIL PROTECTED] for Human Genetics, Oxford www.well.ox.ac.uk

[R] survreg anova: problem with indirect invocation

package version 2.18. Many thanks, William =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-= Dr William Valdar ++44 (0)1865 287 717 Wellcome Trust Centre [EMAIL PROTECTED] for Human Genetics, Oxford www.well.ox.ac.uk/~valdar

[R] setting value arg of pdSymm() in nlme

. Is this functionality deprecated in lme4 and excluded from lmer? Many thanks, William Version details: running R 2.1.0 on windows XP, using nlme 3.1-57. =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-= Dr William Valdar ++44 (0)1865 287 717 Wellcome Trust Centre [EMAIL

Re: [R] setting value arg of pdSymm() in nlme

not be possible. William =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-= Dr William Valdar ++44 (0)1865 287 717 Wellcome Trust Centre [EMAIL PROTECTED] for Human Genetics, Oxford www.well.ox.ac.uk/~valdar __ R-help

Re: [R] refitting lm() with same x, different y

RSS's are on a level footing and map trivially to the p-value. I can therefore take the RSS at, say, the 5th percentile, turn it into a p-value and report that as my 5% threshold. Thanks again, William =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-= Dr William Valdar

[R] refitting lm() with same x, different y

thanks, William =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-= Dr William Valdar ++44 (0)1865 287 717 Wellcome Trust Centre [EMAIL PROTECTED] for Human Genetics, Oxford www.well.ox.ac.uk/~valdar __ R-help

[R] anova grouping of factors in lme4 / lmer

19.5 9.5 4242.0 12.7213 0.0003655 *** Experimenter 68.5 1.4 4242.0 1.9067 0.0759496 . Does anyone know how to make anova() group properly with lmer objects? Many thanks, William =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-= Dr William Valdar ++44 (0

[R] Re: anova grouping of factors in lme4 / lmer

Hi. Yesterday a revised version of lme4 came out that makes my question unnecessary. I was using the latest version available from my mirror (which was 6 days old rather than 1 day old). Please disregard my earlier posting. Thanks. William On Fri, 25 Feb 2005, William Valdar wrote: Hi. I'm