# [R] practical help ... solving a system...

From: Carsten Steinhoff <carsten.steinhoff_at_stud.uni-goettingen.de>
Date: Mon 20 Jun 2005 - 06:58:05 EST

Hello,

I want to estimate the parameters of a binomial distributed rv using MLE. Other distributions will follow.
The equation system to solve is not very complex, but I've never done such work in R and don't have any idea how to start...

The system is:

(1) n*P = X

(2) [sum {from j=0 to J-1} Y{j} /(n-j)] = -n * ln (1-X / n)

where * only X is given (empirical mean)

• J is maximum observed
• Y is the number of observations above j
• n and P are the parameters of the binomial distribution to find....

Who could help me with an example-solution for my "first" distribution ? I also need a hint how to make the sum-element.

Maybe there's another - more simple - way to estimate the parameters... first I tried via "FITDISTR" but without success.

Thanx a lot for your help.

Carsten

[[alternative HTML version deleted]]

R-help@stat.math.ethz.ch mailing list