Displaying 1 result from an estimated 1 matches for "sqrtincom".
Did you mean:
sqrtincome
2011 Mar 13
1
troubles with logistic regression
...ying to
run logistic regression. here is my code:
data <- read.csv(file="C:/Users/fieder.data.2000.csv")
# creating subset of men
fieder.male<-subset(data,data[,8]==1)
unmarried.male<-subset(data,data[,8]==1&data[,6]==1)
# glm fit
agesq.male<-(unmarried.male[,5])^2
male.sqrtincome<-sqrt(unmarried.male[,9])
fieder.male.mar.glm<-glm(as.factor(unmarried.male[,6])~
factor(fieder.male[,7])+fieder.male[,5]+agesq.male+
male.sqrtincome,binomial(link="logit") )
par(mfrow=c(1,1))
plot(c(0,300),c(0,1),pch=" ",
xlab="sqrt income, truncated at 90000...