1. I generate a random sample x=(x1,…,xn) of size n from dist A with parameters (a1,a2) in R.

2. I have a likelihood function which involves another 2 dist B & C with parameters (b1,b2,b3) & (c1,c2) respectively. So, I want to maximize this likelihood to obtain mle of (b1,b2,b3,c1,c2).

3. I use R to solve this. But, I have a problem in getting the initial/starting values of (b1,b2,b3,c1,c2).

a) I’m thinking of using method of moment (henceforth MOM) to get the initial values but I’m not sure on how to do this.

b) Should I apply MOM to dist B & C respectively and then plug in the xi values into those moments function? Does this make sense since xi are coming from A and not B or C?

c) If I apply MOM to dist A, what is the benefit I get from it? I already know and don’t care about (a1,a2) and I believe this is not going to help me getting the hints for (b1,b2,b3,c1,c2) values.

d) As far as I know, we usually apply MOM onto the dist we are interested in estimating its parameter only and usually, this dist is the dist we use to generate the random sample. But, this is not the case we considered here. The dist of interest are B & C but the generated dist are A. So, I think b) or c) are not the right thing to do. Am I right?

e) Or should I apply MOM onto the likelihood function in 2? Can we really find the moments of likelihood? It sounds complicated to me.

f) Anyway, does anyone know how to connect MOM and likelihood function we considered here?

g) Or do you have any ideas on getting the initial values of (b1,b2,b3,c1,c2)?

Correct me if any of my belief above is wrong.I really need advice/opinion regarding this matter and any help will be really appreciated. Thanks in advanced.