Displaying 2 results from an estimated 2 matches for "rad09".
Did you mean:
ad09
2011 Aug 11
1
Subsampling data
...and females separately*
* *
library(foreign)
Datatemp <- read.spss("H:/Skjol/Data/HL/t1and2b.sav", use.value.labels = F)
> table(Datatemp$sex)
1 2
3049 3702
> attributes(Datatemp)
…
$names
[1] "nomiss" "Bin" "rad09" "year"
"sex"
"income" "adults"
[8] "children" "student" "retired" "disabled"
"homemaker" "unemployed" "employed"
[15] "occupation"...
2012 May 04
0
Converting code from gee() to geeglm()
...sample weights
One suggestion was to use geeglm()
I hope someone can help me to solve a problem that arises when converting a
code from gee to geeglm.
*Here is a code that I wrote with the original data, not weighted: *
> m1 <- gee( Bin ~ educ+agemean+ residencysize + yearx , id = rad09 ,
data = Males, subset = marp1 == 1 , + family = binomial, corstr
="unstructured" )
(Intercept) educ agemean residencysize yearx
-0.23875 -0.17931 -0.01470 -0.07418 -0.15200
> se <- summary(m1)$coefficients["yearx", &...