[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]]

__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Re: [R] R beginner

[avadhoot velankar]

Hi andrew

Stick to one tutorial before you get confident enough. I found this one
very good. informative, easy to understand and best thing you basically try
everything out on sample files provided. word of caution though instead of
using ctrl+R command in script file, manually type it in consol initially
to learn about syntax

http://www.unt.edu/rss/class/Jon/R_SC/

welcome to R

avadhoot


On Mon, Dec 3, 2012 at 10:34 AM, andrewcd andre...@gmail.com wrote:

 My $.02: re-do an analysis that you did in another software package in R,
 making sure that you get the same results.  A good way to learn any
 language.



 --
 View this message in context:
 http://r.789695.n4.nabble.com/R-beginner-tp4651719p4651759.html
 Sent from the R help mailing list archive at Nabble.com.

 __
 R-help@r-project.org mailing list
 https://stat.ethz.ch/mailman/listinfo/r-help
 PLEASE do read the posting guide
 http://www.R-project.org/posting-guide.html
 and provide commented, minimal, self-contained, reproducible code.


[[alternative HTML version deleted]]

__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

[R] Species scores on PCoA

[avadhoot velankar]

Dear members,

I am conducting principal coordinates anaysis on a set of morphometry data.
I have 6 variables and 27 samples. I used almost all the packages for PCoA,
 and plotted ordination using orditkplot.

I want to know is there any way i can add variables scores in the same
plot. I want to check which variable is contributing most to variation.

Thank you in anticipation

Avadhoot

[[alternative HTML version deleted]]

__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

[R] help interpreting dudi.pco

[avadhoot velankar]

I am working on morphometry of hairs and want to see if selected variables
are giving significantly distinct groups.

I am new to both R and principal coordinate analysis. Scatter plot is
showing distinct groups but i dont know how to refine the analysis and
interprete it.

[[alternative HTML version deleted]]

__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.