Displaying 1 result from an estimated 1 matches for "mydatanew".
Did you mean:
mydat_new
2010 Oct 31
1
Questions about Probit Analysis
...uot;))
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")
mydatanew$successes <-as.numeric(md[ ,2])
mydatanew$failures <-as.numeric(md[ ,1])
########################################################################
where first I carry out a logit regression of...