Stats: mle with two parameters

Click For Summary
SUMMARY

The discussion centers on finding the maximum likelihood estimates (MLEs) for two parameters, x1 and x2, in a genetics context where both variables are constrained (x1 > 0, x2 > 0, and x1 + x2 < 1). The provided MLE formulas are x1(hat) = (2N1 + N4 + N5)/2n and a similar expression for x2(hat), where n is the sum of all Ni values. The challenge arises when differentiating the likelihood function, as the term for p2 does not include x1, leading to the loss of the necessary N2 term. This indicates a potential misunderstanding of MLEs with multiple variables, which requires careful consideration of constraints and relationships between parameters.

PREREQUISITES
  • Understanding of maximum likelihood estimation (MLE) principles
  • Familiarity with probability distributions in genetics
  • Knowledge of calculus, particularly differentiation
  • Basic concepts of genotype probabilities and allele inheritance
NEXT STEPS
  • Study the derivation of MLEs for multiple parameters in statistics
  • Learn about constrained optimization techniques in MLE problems
  • Explore the application of the multinomial distribution in genetics
  • Review examples of MLE in genetic studies involving multiple alleles
USEFUL FOR

Statisticians, geneticists, and students studying statistical methods in genetics who are working with maximum likelihood estimation in multi-parameter scenarios.

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.
 

Similar threads

  • · Replies 10 ·
Replies
10
Views
49K
Replies
3
Views
5K
  • · Replies 3 ·
Replies
3
Views
9K
  • · Replies 177 ·
6
Replies
177
Views
31K