similar to: Regresión Beta: más rápido?

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

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",

Hola Jorge, Puesto que trabajas en Mac, lo tienes muy inmediato el hacerlo sobre shell-script con los comandos "cron" y "crontab" si no recuerdo mal. En el cron/crontab lo que dirías es que se lanzara un script de R (este script de R incluye el generador las funciones y hasta su cierre con "q()"). De esta forma no tienes por qué tener una instancia de R abierta de

Dear all, I am trying to get an estimate of uncertainty surrounding a single predicted value from a beta regression model (this is similar to a logistic glm - in that it involves a link function and linear predictor - but it uses the beta distribution rather than discrete binomial). For example: library(betareg) data("GasolineYield")

Hola Jorge, En Windows hay un equivalente al "cron": http://stackoverflow.com/questions/132971/what-is-the-windows-version-of-cron Pero para que sea plataforma-independiente, o bien incluyes en "R" con la solución que has planteado. O bien, haces que un Windows sea como un Linux/Unix con Cygwin que también tiene cron/crontab (esto para mí es matar moscas a cañonazos), pero es

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

Buenos dias a todos, Me gustaria genear n muestras de una N(0,1) a determinados intervalos de tiempo. Asi por ejemplo, si el vector de tiempos es tiempos <- c(2, 4, 6) la primera columan debe generarse a los 2 segundos de ejecutado el script, la 2a columna a los 4 segundos y la ultima a los 6 segundos. Al final el tiempo total de ejecucion sera aproximadamente 12 segundos para completar

Dear all, we received a bug report for betareg, that in some cases the optim call in betareg.fit would hang the R session and the command cannot be interrupted by Ctrl-C? We narrowed down the problem to the dbeta function which is used for the log likelihood evaluation in betareg.fit. Particularly, the following command hangs the R session to a 100% CPU usage in all systems we tried it (OS X

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:

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

Dear useRs, version 2.2-2 of the "betareg" package has just been released on CRAN http://CRAN.R-project.org/package=betareg accompanied by an article in the Journal of Statistical Software http://www.jstatsoft.org/v34/i02/ The package provides beta regression for data in the unit interval (0, 1) such as rates and proportions. The manuscript replicates several practical

Dear useRs, version 2.2-2 of the "betareg" package has just been released on CRAN http://CRAN.R-project.org/package=betareg accompanied by an article in the Journal of Statistical Software http://www.jstatsoft.org/v34/i02/ The package provides beta regression for data in the unit interval (0, 1) such as rates and proportions. The manuscript replicates several practical

Muchas gracias, Carlos por tu respuesta. En este caso una solucion solo en R pareceria la manera de resolver el problema y hacia alla apuntare. El algoritmo que necesito ejecutar es mas complejo que el que describi en mi mensaje inicial, asi que aun hay mucho camino por recorrer. Saludos, Jorge.- 2013/10/9 Carlos Ortega <cof@qualityexcellence.es> > Hola Jorge, > > En Windows

Hola Carlos. Muchas gracias por la sugerencia. En realidad me gustaria tener la posibilidad de que el script fuera multiplataforma. Alguna otra sugerencia? Gracias! --JIV Sent from my phone. Please excuse my brevity and misspelling. On Oct 8, 2013, at 10:57 AM, Carlos Ortega <cof@qualityexcellence.es> wrote: Hola Jorge, Puesto que trabajas en Mac, lo tienes muy inmediato el hacerlo sobre

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

Stephen, thanks for your effort. The more appropriate list for this discussion is probably R-devel (as far as I understand it) so I've moved the discussion there. Related topics have already been discussed in the past. Specifically, I remember contributions by Paul Johnson ("rockchalk" package) and John Fox ("effects" and "car" package) as their packages

Hi, Im just doing some Beta-Regressions with the betareg package. My question is now, is there a possibility to calculate the marginal effects with the betareg package or is there another package which can handle marginal effects on regression output for the "beta" class? I try to calculate the marginal effects also by hand but then i have the problems with the standard errors :-( Thx

Buenas tardes a todos, Estoy interesado en realizar varios gráficos y luego guardarlos individualmente como PDF. Hasta ahí todo bien. Sin embargo, el archivo resultante tiene aproximadamente 20MB. Podría alguien por favor indicarme cómo reducir su tamaño? A continuación un ejemplo de lo que intento hacer y mi sessionInfo(). Estoy usando Mac OS X. Muchas gracias, Jorge Ivan Velez # R code

Hi, what's mean "/" in command: betareg(inf~Grupo/Sexo, data=dados) it's a effect nested? -------------------------------------- Silvano Cesar da Costa Departamento de Estat?stica Universidade Estadual de Londrina Fone: 3371-4346