Displaying 1 result from an estimated 1 matches for "oprobit3".
Did you mean:
oprobit
2012 Nov 12
1
Invalid 'times' argument three-category ordered probit with maximum likelihood
...Here is the code I used:
setwd("C:/Users/Terry/Desktop/Terry/School/Fall 2012/ML/HW")
#Load Libraries
library(MASS)
library(tile)
library(simcf)
#Load Data
econrate <- read.csv("hw4econ3.csv", header=TRUE, sep=",")
attach(econrate)
#Ordered Probit Liklihood
llk.oprobit3 <- function(param, x, y) {
os <- rep(1, nrow(x))
x <- cbind(os, x)
b <- param[1:ncol(x)]
t2 <- param[(ncol(x)+1)]
xb <- x%*%b
p1 <- log(pnorm(-xb))
if (t2<=0) p2 <- -(abs(t2)*10000)
else p2 <- log(pnorm(t2-xb)-pnorm(-xb))
p3 <- log(1-pnorm(t2-xb))...