Stats: mle with two parameters

bennyska
Messages
110
Reaction score
0

Homework Statement


in a genetics situation, we have two variables, x1 and x2, such that both x1 and x2 >0, and x1+x2<1.
we have:
p1 = x12
p2 = x22
p3 = (1-x1-x2)2
p4 = 2x1x2
p5 = 2x1(1-x1-x2)
p6 = 2x2(1-x1-x2)

find the mles for x1 and x2.

Homework Equations



the answer (from the book): x1(hat) = (2N1 + N4 + N5)/2n, where n = sum (Ni). x2(hat) is similar.

The Attempt at a Solution


doing the usual mle stuff, i have sum(Ni*ln(pi)), take the derivative, set to zero, and solve for my parameter. when i do this work, i get an answer similar (but incorrect) to that of the book, with a little problem: in order to get n, i need all the Nis. when i take the derivative with respect to x1, there is one term (p2) that doesn't have x1 in it, so it drops out, and i lose that term (N2) that i think i need. i can't define x1 by x2 (i think, since x1 + x2 < 1 doesn't tell me too much). my teacher assigned this one to us, without having done the problem herself, and then she found herself unable to do it. we figure there's some little trick involved, but I'm not seeing it. also, i haven't really done mles with more than one variable, so maybe that's what I'm missing. any help would be greatly appreciated.
 
Physics news on Phys.org
bennyska said:

Homework Statement


in a genetics situation, we have two variables, x1 and x2, such that both x1 and x2 >0, and x1+x2<1.
we have:
p1 = x12
p2 = x22
p3 = (1-x1-x2)2
p4 = 2x1x2
p5 = 2x1(1-x1-x2)
p6 = 2x2(1-x1-x2)

.

We have them, but what are they? Are the p's functions?
 
sorry for my vagueness. the actual problem is something like this: a certain gene has 3 alleles, resulting in six possible genotypes. each p is the probability of being a certain genotype. each parent has an x1 or x2 probability of passing that allele to the child.
 
There are two things I don't understand about this problem. First, when finding the nth root of a number, there should in theory be n solutions. However, the formula produces n+1 roots. Here is how. The first root is simply ##\left(r\right)^{\left(\frac{1}{n}\right)}##. Then you multiply this first root by n additional expressions given by the formula, as you go through k=0,1,...n-1. So you end up with n+1 roots, which cannot be correct. Let me illustrate what I mean. For this...
Back
Top