[R] Help with maximum likelihood estimation

[avadhoot velankar]

I am analyzing animal movement pattern using levy flight pattern and want
to fit power function to observed data and estimate exponent using Maximum
Likelihood Estimation.

I am using

lff<-function(mu){1-1/mean(log(x))}
library(stats4)
mle(lff, start = list(mu = 1))

where x is the observed data.

I am getting following error

Error in solve.default(oout$hessian) :
  Lapack routine dgesv: system is exactly singular: U[1,1] = 0

This is first time i am writing function and I am as good at algebra
as brick is good at flying.

Thank you in advance for your help.

-- 
*Avadhoot D. Velankar*

	[[alternative HTML version deleted]]