Displaying 1 result from an estimated 1 matches for "women13".
Did you mean:
women
2007 Apr 18
3
Problems in programming a simple likelihood
...l <- function(mu, y, X) {
n <- nrow(X)
zeta <- X%*%mu
llik <- 0
for (i in 1:n) {
if (y[i]==1)
llik <- llik + log(pnorm(zeta[i,], mean=0, sd=1))
else
llik <- llik + log(1-pnorm(zeta[i,], mean=0, sd=1))
}
return(-llik)
}
women <- read.table("~/R/Examples/Women13.txt", header=TRUE) # DATA
# THE DATA SET CAN BE ACCESSED HERE
# women <-
read.table("http://wps.aw.com/wps/media/objects/2228/2281678/Data_Sets/ASCII/Women13.txt", header=TRUE)
# I HAVE CHANGED THE NAMES OF THE VARIABLES
# J is changed to "work"
# M is changed to "...