Hi all,
    I have the following two function f1 and f2.

f1 <- function(lambda,z,p1){
lambda*(p1*exp(-3*z-9/2)+(0.2-p1)*exp(4*z-8))-(1-lambda)*0.8}

f2 <- function(p1,cl, cu){
 0.8*(pnorm(cl)+(1-pnorm(cu)))/(0.8*(pnorm(cl)+(1-pnorm(cu)))+p1*(pnorm(cl+3)+(1-pnorm(cu+3)))+(0.2-p1)*(pnorm(cl-4)+(1-pnorm(cu-4))))}-0.05

First fix p1 to be 0.15.
(i) choose a lambda value, say lamda=0.6,
(ii) there should be two solutions for the equation f1(z, p1=0.15,
lambda=1)=0, say cl and cu,
(iii) evaluate f2 at cl and cu found in (ii) and p1=0.15,

   How do we find the lambda value such that f2=0?

    Thank you very much!
                                                           Hannah

	[[alternative HTML version deleted]]

