[R] Gaussian glm for grouped data with unequal variances

2007-01-16 Thread Dawn Ashcourt
Hello - I am fairly new to R, (i.e., ability to create functions/write programs insignificant) and was wondering if there might be a convenient way to model the following: I want to fit a gaussian glm to grouped data, while allowing for unequal variances in each of the groups. More specifically,

Re: [R] Gaussian glm for grouped data with unequal variances

2007-01-16 Thread Chuck Cleland
Dawn Ashcourt wrote: Hello - I am fairly new to R, (i.e., ability to create functions/write programs insignificant) and was wondering if there might be a convenient way to model the following: I want to fit a gaussian glm to grouped data, while allowing for unequal variances in each of the

Re: [R] Gaussian glm for grouped data with unequal variances

2007-01-16 Thread Thilo Kellermann
Hi, I would suggest to use lme for doing that: library(nlme) vf1Ident - varIdent( c(m = 0.5), form = ~ 1 | group ) fm.lme - lme(data ~ group, weights = vf1Ident. data = DATA) anova(fm.lme) Hope this helps, Thilo On Tuesday 16 January 2007 10:30, Dawn Ashcourt wrote: Hello - I am fairly new to