That worked, thanks!

On Fri, Oct 28, 2011 at 12:27 PM, bbolker [via R] <
ml-node+s789695n3948463h12@n4.nabble.com> wrote:

> ethan.shepherd <someguy235 <at> gmail.com> writes:
>
> >
> > I'm getting errors when running what seems to be a simple Weibull
> > distribution function:
> >
>
>   [snip]
>
> > If I change the data to this:
> >
>
>  [snip]
>
> > I get the error "Error in fitdistr(x, "weibull"): optimization failed"
> >
> > I can run a Weibull distribution in SAS with this same data, and it
> gives me
> > what looks like a reasonable answer.
> >
>
>
>  How about:
>
> x <- c(4,22,26,27,44,46,83,83,122,125,129,151,153,157,171,
> 173,186,199,201,202,205,219,239,242,252,315,326,449,771)
>
> hist(x)
> fit2 <-fitdistr(x, 'weibull',lower=c(0.01,0.01))
>
> hist(x,freq=FALSE,col="gray",ylim=c(0,0.004))
> with(as.list(coef(fit2)),curve(dweibull(x,shape=shape,scale=scale),
>      add=TRUE,col=2))
>
> ______________________________________________
> [hidden email] <http://user/SendEmail.jtp?type=node&node=3948463&i=0>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.
>
>
> ------------------------------
>  If you reply to this email, your message will be added to the discussion
> below:
>
> http://r.789695.n4.nabble.com/weibull-fitdistr-problem-optimization-failed-tp3947997p3948463.html
>  To unsubscribe from weibull fitdistr problem: optimization failed, click
> here<http://r.789695.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscribe_by_code&node=3947997&code=c29tZWd1eTIzNUBnbWFpbC5jb218Mzk0Nzk5N3wyNDkwMTc5NjY=>.
>
>



-- 
someguy235@gmail.com
(828) 242-0036

"You just lost The Game."


--
View this message in context: http://r.789695.n4.nabble.com/weibull-fitdistr-problem-optimization-failed-tp3947997p3963303.html
Sent from the R help mailing list archive at Nabble.com.
	[[alternative HTML version deleted]]

