Displaying 2 results from an estimated 2 matches for "lnwagevec".
2017 Nov 07
2
Using MLE on a somewhat unusual likelihood function
...temw-(1-alpha)*xhat)/alpha))/(eta+lam*survivalf(statemw)))
} else if (cpsdata2[i,1,]>statemw) {
wagevec[i,] <-
c(lam*plnorm((cpsdata2[i,1,]-(1-alpha)*xhat)/alpha,mean=mu,sd=logalpha)/(alpha*(eta+lam*survivalf(statemw))))
}
else {
wagevec[i,] <- NA
}
}
lnwagevec <- log(wagevec)
-sum(lnwagevec>-200 & ln2wagevec<200, na.rm=TRUE)
}
fit <- mle(wagefn, start=listmat, method= "L-BFGS-B",lower=
c(-Inf,0),upper=c(Inf,Inf)
I know the likelihood function is a handful but it does return a reasonable
looking vector of values. The "l...
2017 Nov 07
0
Using MLE on a somewhat unusual likelihood function
...*survivalf(statemw)))
> } else if (cpsdata2[i,1,]>statemw) {
> wagevec[i,] <-
> c(lam*plnorm((cpsdata2[i,1,]-(1-alpha)*xhat)/alpha,mean=mu,sd=logalpha)/(alpha*(eta+lam*survivalf(statemw))))
> }
> else {
> wagevec[i,] <- NA
> }
> }
> lnwagevec <- log(wagevec)
> -sum(lnwagevec>-200 & ln2wagevec<200, na.rm=TRUE)
> }
>
> fit <- mle(wagefn, start=listmat, method= "L-BFGS-B",lower=
> c(-Inf,0),upper=c(Inf,Inf)
>
>
> I know the likelihood function is a handful but it does return a reasonabl...