Displaying 1 result from an estimated 1 matches for "myprobit2".
Did you mean:
myprobit
2010 Oct 31
1
Questions about Probit Analysis
...ry.csv"))
names(mydata) <- c("outcome","x1","x2","x3")
myprobit <- glm(mydata$outcome~mydata$x1+mydata$x2+as.factor(mydata$x3),
family=binomial(link="probit"))
print(summary(myprobit))
#Now assume I can make a regression only on x1
myprobit2 <- glm(mydata$outcome~mydata$x1, family=binomial(link="probit"))
print(summary(myprobit2))
#express in terms of counts
md <- t(table(mydata$outcome, mydata$x1))
# create new dataframe
mydatanew <- data.frame(as.numeric(row.names(md)))
names(mydatanew) <- c("x1"...