[R] betareg help

Ben Bolker bbolker at gmail.com
Sat Mar 12 16:28:24 CET 2011


Vlatka Matkovic Puljic <v.matkovic.puljic <at> gmail.com> writes:

> 
> Maybe I should include data:
> 
> results$cond
> [1] 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001
[snip]
> [205] 0.950 0.960 0.980 0.980 0.999
> 
> 2011/3/12 Vlatka Matkovic Puljic <v.matkovic.puljic <at> gmail.com>
> 

> > I'm trying to do betareg on my dataset.
> > Dependent variable is not normally distributed and is proportion (of condom
> > use (0,1)).
> > But I'm having problems:
> >
> >  gyl<-betareg(cond ~ alcoh + drug, data=results)
> > Error in optim(par = start, fn = loglikfun, gr = gradfun, method = method,
> > :
> > initial value in 'vmmin' is not finite
> >
> > Why is R returning me error in optim()?
> > What should I do?
> >

  My first guess would have been that you had zeros and ones in
your data, but you seem to have taken care of that.
   Are you willing to post the whole data set (use dput() to get it
in a useful form)?

  Ben Bolker



More information about the R-help mailing list