search for: wagevec

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

Did you mean: pagevec
2017 Nov 07
2
Using MLE on a somewhat unusual likelihood function
...to where the parameters usually end up. I have also tried switching these around a bit. Here is the code: xhat = c(statemw-(1-alpha)*rval) survivalf <- function(x) {(1-plnorm(statemw,mean=mu,sd=logalpha))} wagefn <- function(lam, eta, alpha, xhat, mu, logalpha) { n=nrow(cpsdata2) wagevec = matrix(nrow=n,ncol=1) for (i in 1:n) { if (cpsdata2[i,2] > 0){ wagevec[i,] <- c(eta*lam*survivalf(statemw)*exp(-lam*survivalf(statemw)*cpsdata2[i,2,])) } else if (cpsdata2[i,1,]==statemw) { wagevec[i,] <- c(lam*(survivalf(statemw)-survivalf((statemw-(1-alpha)*xhat...
2017 Nov 07
0
Using MLE on a somewhat unusual likelihood function
...I have also tried switching these around a bit. > > Here is the code: > > xhat = c(statemw-(1-alpha)*rval) > survivalf <- function(x) {(1-plnorm(statemw,mean=mu,sd=logalpha))} > > wagefn <- function(lam, eta, alpha, xhat, mu, logalpha) { > n=nrow(cpsdata2) > wagevec = matrix(nrow=n,ncol=1) > for (i in 1:n) { > > if (cpsdata2[i,2] > 0){ > wagevec[i,] <- > c(eta*lam*survivalf(statemw)*exp(-lam*survivalf(statemw)*cpsdata2[i,2,])) > } else if (cpsdata2[i,1,]==statemw) { > wagevec[i,] <- > c(lam*(survivalf(statemw...