Displaying 4 results from an estimated 4 matches for "gender1".
Did you mean:
gender
2012 May 04
2
Binomial GLM, chisq.test, or?
...viance Residuals:
Min 1Q Median 3Q Max
-1.84648 -0.96954 -0.00036 1.11269 2.03933
Coefficients: (12 not defined because of singularities)
Estimate Std. Error z value
(Intercept) -1.657e+01 2.400e+03 -0.007
gender1 -2.231e-01 9.220e-01 -0.242
colonyother 9.531e-02 8.006e-01 0.119
cohort2002 1.717e-08 3.393e+03 0.000
cohort2003 1.766e+01 2.400e+03 0.007
cohort2004 1.807e+01 2.400e+03 0.008
co...
2004 Dec 16
1
help with multiple imputation using imp.mix
...o understand, but
there's just this one step I don't get...perhaps someone knows the answer?
Thanks!
Jens
My code runs:
data<-read.table('http://www.courses.fas.harvard.edu/~gov2001/Data/immigrati
on.dat',header=TRUE)
library(mix)
rngseed(12345678)
# Preare data for imputation
gender1<-c()
gender1<-as.integer(data$gender)
gender1[gender1==1]<-2
gender1[gender1==0]<-1
data$gender<-gender1
x<-cbind(data$gender,data$ipip,data$ideol,data$prtyid, data$wage1992)
colnames(x)<-c("gender","ipip", "ideol", "prtyid","wag...
2008 May 08
2
poisson regression with robust error variance ('eyestudy
...a bit off.
### Paul Johnson 2008-05-08
### sandwichGLM.R
system("wget http://www.ats.ucla.edu/stat/stata/faq/eyestudy.dta")
library(foreign)
dat <- read.dta("eyestudy.dta")
### Ach, stata codes factor contrasts backwards
dat$carrot0 <- ifelse(dat$carrot==0,1,0)
dat$gender1 <- ifelse(dat$gender==1,1,0)
glm1 <- glm(lenses~carrot0, data=dat, family=poisson(link=log))
summary(glm1)
library(sandwich)
vcovHC(glm1)
sqrt(diag(vcovHC(glm1)))
sqrt(diag(vcovHC(glm1, type="HC0")))
### Result:
# > summary(glm1)
# Call:
# glm(formula = lenses ~ carrot0...
2004 Oct 11
0
lmList - strange problem
...which takes 'Data' partitioned according to the
levels of a grouping factor (gf) and individual 'lm' fits are obtained
for each 'data'partition, using a model defined as in "lm". So my call
to lmList looks something like:
> mg.lis <- lmList(strmbf ~ age1c + gender1 + sysbp.clinic +
+ diabp.clinic + ldl1 + mets.total + bmi1c | gf,
+ data = Data, subset= avblock=="normal");
but I get a strange error message:
"Error in model.frame.default(formula = form, data = dat, na.action =
na.action, : object is not a matrix"
If...