[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