search for: p1t

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

Did you mean: p1
2012 Oct 05
2
problem with convergence in mle2/optim function
...))*t)*mu1* sqrt((mu1 + mu2 + p1 + p2)^2 - 4*(mu2*p1 + mu1*(mu2 + p2)))))/ exp((1/2)*(mu1 + mu2 + p1 + p2 + sqrt((mu1 + mu2 + p1 + p2)^2 - 4*(mu2*p1 + mu1*(mu2 + p2))))*t)/(2*(mu2*p1 + mu1*(mu2 + p2))* sqrt((mu1 + mu2 + p1 + p2)^2 - 4*(mu2*p1 + mu1*(mu2 + p2))))) # True parameter values p1t = 2; p2t = 2; mu1t = 0.001; mu2t = 0.001 # Function to calculate probabilities from 'true' parameter values psim <- function(x){ params <- list(p1 = p1t, p2 = p2t, mu1 = mu1t, mu2 = mu2t, t = x) eval.P1 <- eval(P1, params) eval.P2 <- eval(P2, params) P3 <- 1 - ev...
2012 Sep 27
0
problems with mle2 convergence and with writing gradient function
...))*t)*mu1* sqrt((mu1 + mu2 + p1 + p2)^2 - 4*(mu2*p1 + mu1*(mu2 + p2)))))/ exp((1/2)*(mu1 + mu2 + p1 + p2 + sqrt((mu1 + mu2 + p1 + p2)^2 - 4*(mu2*p1 + mu1*(mu2 + p2))))*t)/(2*(mu2*p1 + mu1*(mu2 + p2))* sqrt((mu1 + mu2 + p1 + p2)^2 - 4*(mu2*p1 + mu1*(mu2 + p2))))) # True parameter values p1t = 2; p2t = 2; mu1t = 0.001; mu2t = 0.001 # Function to calculate model probabilities from true parameter values psim <- function(x){ params <- list(p1 = p1t, p2 = p2t, mu1 = mu1t, mu2 = mu2t, t = x) eval.P1 <- eval(P1, params) eval.P2 <- eval(P2, params) P3 <- 1 - eval.P...