This is really a statistical issue. What do you think the Intercept term represents? See ?contrasts.
Cheers, Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Thu, Feb 15, 2018 at 5:27 PM, Marc Girondot via R-help < r-help@r-project.org> wrote: > Dear R-er, > > I try to get the standard error of fitted parameters for factors with a > glm, even the reference one: > > a <- runif(100) > b <- sample(x=c("0", "1", "2"), size=100, replace = TRUE) > df <- data.frame(A=a, B=b, stringsAsFactors = FALSE) > > g <- glm(a ~ b, data=df) > summary(g)$coefficients > > # I don't get SE for the reference factor, here 0: > > Estimate Std. Error t value Pr(>|t|) > (Intercept) 0.50384827 0.05616631 8.9706490 2.236684e-14 > b1 -0.03598386 0.07496151 -0.4800311 6.322860e-01 > b2 0.03208039 0.07063113 0.4541962 6.507023e-01 > > # Then I try to change the order of factors, for example: > > df$B[df$B=="0"] <- "3" > g <- glm(a ~ b, data=df) > summary(g)$coefficients > > By I get the same... > > Any idea ? > > Thanks > > Marc > > ______________________________________________ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posti > ng-guide.html > and provide commented, minimal, self-contained, reproducible code. > [[alternative HTML version deleted]] ______________________________________________ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.