similar to: betareg 2.2-2: Beta regression

Dear R-helpers: We have had fun using betareg to fit models with proportions as dependent variables. However, in the analysis of these models we found some wrinkles and don't know where is the best place to start looking for a fix. The problems we see (so far) are that 1. predict ignores newdata 2. anova does not work Here is the small working example: ---------------------------- x

The Journal of Statistical Software was founded by Jan de Leeuw in 1996. Currently, there are 26 volumes with 260 aricles and the journal publishes about 4 new volumes per year. The web address is www.jstatsoft.org. All articles are pdf files, since 2005 all produced from JSS LaTeX templates. JSS is edited by Jan de Leeuw, since 2005 jointly with Achim Zeileis of the Wirtschaftsuniversität Wien.

Hello, I have a dataset with proportions that vary around a fixed mean, is it possible to use betareg to look at variance in the dispersion parameter while keeping the mean fixed? I am very new to R but have tried the following: svec<-c(qlogis(mean(data1$scaled)),0,0,0) f<-betareg(scaled~-1 | expt_label + grouped_hpi, data=data1, link.phi="log",

Over the last couple of years I have written quite a few R programs for various "psychometric" techniques, and I am regularly updating and expanding what is there. I now have (wholly or partially), or have planned -- gifi package (update to homals on CRAN). Code for multiple correspondence analysis, nonlinear principal component analysis, nonlinear multiset canonical

Folks, The betareg package appears to be polished and works well. But I would like to look at the exact formulas for the underlying model being estimated, the likelihood function, etc. E.g. if one has to compute \frac{\partial E(y)}{\partial x_i}, this requires careful calculations through these formulas. I read "Regression analysis of variates observed on (0,1): percentages, proportions and

Hi I'm trying to fit a model in betareg and I'm getting errors, but have no idea what they mean or how to solve them. Does anyone have experience with this? > model <- betareg(ACT ~ ST*SoilT, data = actDL_F) Warning messages: 1: In sqrt(W) : NaNs produced 2: In sqrt(W) : NaNs produced 3: In sqrt(1 + phihat) : NaNs produced data summaries don't give any na's or problems I

At www.jstatsoft.org we now have three special volumes Volume 22, Ecology and Ecological Modelling in R (Thomas Kneib and Thomas Petzoldt, eds.) Volume 20, Psychometrics in R (Jan de Leeuw and Patrick Mair, eds.) Volume 18, Spectroscopy and Chemometrics in R (Kate Mullen and Ivo van Stokkum, eds.) Coming soon -- Econometrics in R -- Political Methodology in R Gleam in my eye --

Dear all, I get an error using betrag on this data set :http://dl.dropbox.com/u/1866110/dump.csv. I run it like this regression f2.1=betareg(Y~X1+X2,data=dump) summary(f2.1) I get : Call: betareg(formula = Y ~ X1 + X2, data = dump) Standardized weighted residuals 2: Error in quantile.default(x$residuals) : missing values and NaN's not allowed if 'na.rm' is FALSE In addition:

Version 1.0-0 of the package "AER" for "Applied Econometrics with R" has been released to CRAN (http://CRAN.R-project.org/package=AER) a few weeks ago. It accompanies Applied Econometrics with R Christian Kleiber, Achim Zeileis http://www.springer.com/978-0-387-77316-2 http://www.amazon.com/Applied-Econometrics-R-Use/dp/0387773169/ from Springer's useR!

Version 1.0-0 of the package "AER" for "Applied Econometrics with R" has been released to CRAN (http://CRAN.R-project.org/package=AER) a few weeks ago. It accompanies Applied Econometrics with R Christian Kleiber, Achim Zeileis http://www.springer.com/978-0-387-77316-2 http://www.amazon.com/Applied-Econometrics-R-Use/dp/0387773169/ from Springer's useR!

Dear R users, 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

Hi, In using the betareg package, I encounter the following error message: Error in lm.wfit(x, linkfun(y), weights, offset = offset) : NA/NaN/Inf in foreign function call (arg 4) Any help will be most appreciated. Thanks in advance. Alex

Hi, I want to fit 3 beta distributions to my data which ranges between 0 and 1. What are the functions that I can easily call and specify that 3 beta distributions should be fitted? I have already looked at normalmixEM and fitdistr but they dont seem to be applicable (normalmixEM is only for fitting normal dist and fitdistr will only fit 1 distribution, not 3). Is that right? Also, my data has 26

Hello list, I'm following the paper (http://www.jstatsoft.org/v27/i02/paper) on how to use "plm" to run panel regressions, and am having trouble with what I believe should be something very basic. When I run the command (p.9 in the paper): R> dynformula(emp~wage+capital,log=list(capital=FALSE,TRUE),lag=list(emp=2,c(2,3)),diff=list(FALSE,capital=TRUE)) I see: emp ~ wage +

Hello, dear experts. I don't have much experience in building regression models, so sorry if this is too simple and not very interesting question. Currently I'm working on the model that have to predict proportion of the debt returned by the debtor in some period of time. So the dependent variable can be any number between 0 and 1 with very high probability of 0 (if there are no payment)

Dear All, Is there any function in R for fitting linear probability model, as my response variable is a uniformly distributed. Regards, Adam [[alternative HTML version deleted]]

Dear useRs, version 1.0-0 of the "Formula" package has just been released on CRAN http://CRAN.R-project.org/package=Formula accompanied by an article in the Journal of Statistical Software http://www.jstatsoft.org/v34/i01/ It provides simple infrastructre for processing formulas like y ~ x1 + x2 | z1 + z2 where there are two groups of explanatory variables, e.g.,

Dear useRs, version 1.0-0 of the "Formula" package has just been released on CRAN http://CRAN.R-project.org/package=Formula accompanied by an article in the Journal of Statistical Software http://www.jstatsoft.org/v34/i01/ It provides simple infrastructre for processing formulas like y ~ x1 + x2 | z1 + z2 where there are two groups of explanatory variables, e.g.,

As the author of 3 of the 5 methods in optim, I think you may be wasting your time if this is for performance. My reasons are given in http://www.jstatsoft.org/v60/i02 Note that most of the speed benefits of compilation are found in the objective and gradient function, with generally more minor improvements in the method. JN On 14-10-22 06:00 AM, r-devel-request at r-project.org wrote: >

The final registration deadline for the R User Conference is June 20, 2010, one week away. Later registration will not be possible on site! Conference webpage: http://www.R-project.org/useR-2010 Conference program: http://www.R-project.org/useR-2010/program.html Registration: http://www.R-project.org/useR-2010/registration/registration.html The conference is scheduled for July 21-23, 2010,