Displaying 3 results from an estimated 3 matches for "raceblack".
2003 Jan 21
1
Logistic regression: At times correlation matrix of coefficients gets messed up
...quot;,
"black" and "other") in my logistic regression model, the correlation
matrix of the coefficients gets messed up. I get something like:
-----------------------------------------
Correlation of Coefficients:
( A L RACEb
AGE , 1
LWT , 1
RACEblack 1
RACEother . .
attr(,"legend")
[1] 0 ` ' 0.3 `.' 0.6 `,' 0.8 `+' 0.9 `*' 0.95 `B' 1
-------------------------------------
I couldn't figure out how to interpret it. Here is the sequence of
commands and the complete output. (I am using R 1.6....
2011 Sep 29
0
geeglm estimates and standard deviation are too large
...tr="exchangeable")
summary(fit1)$coef gives too large estimates and standard deviation:
Estimate Std.err Wald Pr(>|W|)
(Intercept) 3.07e+16 7.20e+14 1821.29 0.00000
age 6.38e+13 2.22e+13 8.24 0.00409
RACEBlack 1.48e+16 6.28e+14 555.35 0.00000
RACEOther -1.84e+16 6.17e+14 887.78 0.00000
SEXFemale 1.84e+16 5.23e+14 1235.19 0.00000
FEVERYes -4.41e+15 4.74e+14 86.73 0.00000
FEVERUnknown 1.76e+16 1.60e+15 120.55 0.00000
compare...
2003 May 03
1
expand.grid
...moke from the table:
> row.factors(ft)
race smoke
1 white FALSE
2 white TRUE
3 black FALSE
4 black TRUE
5 other FALSE
6 other TRUE
in such a way that they can be directly used in glm:
> glm(ft ~ race + smoke, family=binomial, data = row.factors(ft))
...
Coefficients:
(Intercept) raceblack raceother smokeTRUE
1.841 -1.084 -1.109 -1.116
...
Note that the reference level for race is "white" (the first row).
PS - Obviously, the same analysis is very easy from the original
dataframe (which here is supposed to be missing):
> glm(low ~ race +...