[Rd] weights in glm (PR#8720)
robert.pusz at wp.pl
robert.pusz at wp.pl
Tue Mar 28 17:01:54 CEST 2006
Full_Name: Robert Pusz
Version: 2.2.1
OS: Windows
Submission from: (NULL) (157.25.9.126)
Hello,
In my opinion something is wrong with 'weights' option in glm.
My code is following:
###begin of the code####
cl <- c(5012, 106, 3410, 5655, 1092, 1513, 557, 1351, 3133, 2063, 3257,
4179, 5582, 5900, 8473, 4932, 3463, 5596, 2262, 0, 2638, 1111,
4881, 4211, 6271, 5257, 6926, 6165, 0, 0, 898, 5270, 2268, 5500,
6333, 1233, 1368, 0, 0, 0, 1734, 3116, 2594, 2159, 3786, 2917,
0, 0, 0, 0, 2642, 1817, 3479, 2658, 225, 0, 0, 0, 0, 0, 1828,
0, 649, 984, 0, 0, 0, 0, 0, 0, 599, 673, 603, 0, 0, 0, 0, 0,
0, 0, 54, 535, 0, 0, 0, 0, 0, 0, 0, 0, 172, 0, 0, 0, 0, 0, 0,
0, 0, 0)
w <- c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0,
1, 1, 1, 1, 1, 1, 1, 1, 0, 0, 1, 1, 1, 1, 1, 1, 1, 0, 0, 0, 1,
1, 1, 1, 1, 1, 0, 0, 0, 0, 1, 1, 1, 1, 1, 0, 0, 0, 0, 0, 1, 0,
1, 1, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0,
0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0)
row <- c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10,
1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10,
1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10,
1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10,
1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10)
col <- c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2,
3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 5,
5, 5, 5, 5, 5, 5, 5, 5, 5, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 7, 7,
7, 7, 7, 7, 7, 7, 7, 7, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 9, 9, 9,
9, 9, 9, 9, 9, 9, 9, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10)
Row <- ordered(as.factor(row))
Col <- ordered(as.factor(col))
fit <- glm(cl ~ Row + Col, family = quasipoisson, weights = w)
###end of the code####
When I have written summary(fit) I have got:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 7.64459 NA NA NA
Row2 -0.05467 NA NA NA
Row3 0.24541 NA NA NA
Row4 0.42035 NA NA NA
Row5 0.43961 NA NA NA
Row6 0.04532 NA NA NA
Row7 -0.04881 NA NA NA
Row8 0.25370 NA NA NA
Row9 -0.14976 NA NA NA
Row10 -0.01267 NA NA NA
Col2 0.69283 NA NA NA
Col3 0.62603 NA NA NA
Col4 0.27695 NA NA NA
Col5 0.06056 NA NA NA
Col6 -0.19582 NA NA NA
Col7 -0.83044 NA NA NA
Col8 -1.27914 NA NA NA
Col9 -1.93235 NA NA NA
Col10 -2.49709 NA NA NA
When I omited 'weights=w' above table was filled in with numbers, but the
results were wrong (because of taking zeros in regression).
Could you tell me what's wrong?
Kind regards,
Robert
More information about the R-devel
mailing list