Displaying 2 results from an estimated 2 matches for "survivalf".
Did you mean:
survival
2017 Nov 07
2
Using MLE on a somewhat unusual likelihood function
...ular: U[1,1] = 0
This error sometimes indicates that the list of starting values is too far
from optimum but this is unlikely since I picked values close 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...
2017 Nov 07
0
Using MLE on a somewhat unusual likelihood function
...times indicates that the list of starting values is too far
> from optimum but this is unlikely since I picked values close 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*l...