Displaying 1 result from an estimated 1 matches for "ordfile".
Did you mean:
oldfile
2012 Aug 01
1
optim() for ordered logit model with parallel regression assumption
...thing along that line among those if
conditions. Thanks in advance for any suggestion.
Jun Xu, PhD
Assistant Professor
Department of Sociology
Ball State University
Muncie, IN 47306
####################################################################
library(foreign)
readin <- read.dta("ordfile.dta", convert.factors=FALSE)
myvars <- c("depvar", "x1", "x2", "x3")
mydta <- readin[myvars]
# remove all missings
mydta <- na.omit(mydta)
# theta is the parameter vector
ologit.lf <- function(theta, y, X) {
n <- nrow(X)
k <- ncol(...