search for: ln2wagevec

Displaying 2 results from an estimated 2 matches for "ln2wagevec".

2017 Nov 07
2
Using MLE on a somewhat unusual likelihood function
...} 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 "lnwagevec>-200" etc is an inelegant way of preventing...
2017 Nov 07
0
Using MLE on a somewhat unusual likelihood function
...mw) { > 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 "lnwagevec>-200&quot...