similar to: LMER quasibinomial

Hi, my dependent variable is a proportion ("prob.bind"), and the independent variables are factors for group membership ("group") and a covariate ("capacity"). I am interested in the effects of group, capacity, and their interaction. Each subject is observed on all (4) levels of capacity (I use capacity as a covariate because the effect of this variable is normatively

Hi, I am writing a function to plot a pdf of a distribution, GNL.pdf.fn = function(x,mu,sigma,alpha,beta,rho) { y = x-rho*mu cf.fn = function(s){ cplex = complex(1,0,1) temp1 = alpha*beta*exp(-sigma*s^2/2) temp2 = (alpha-cplex*s)*(beta+cplex*s) out = (temp1/temp2)^rho out } temp.fn = function(s){ (Mod(cf.fn(s)))*cos(Arg(cf.fn(s))-s*y) } int.fn =

I have /etc/auto.master and /etc/auto.smb such that ls /smb/<name of system> works for the two Win98 machines on my network (the machines get automatically mounted and the directory contains "C"). When I try the same command on the Windows XP machine, I get No such file or directory The auto.smb file that came with my OS (Libranet 3) uses smbclient to get a list of

Hi, We''re trying to follow the instructions at http://reductivelabs.com/trac/puppet/wiki/UsingMongrel (Just using a single puppetmaster/mongrel instance and plain mod_proxy for now, instead of balancer one) We got past the "Server is not a class" error by modifying mongrel.rb We also kind-of got past the: /opt/bin/puppetmasterd:293: undefined method `daemonize'' for

Hi, I've got a problem that sounds a lot like this, http://r.789695.n4.nabble.com/Re-R-R-2-12-0-hangs-while-loading-RGtk2-on-FreeBSD-td3005929.html under windoze 7. but it seems to hang with this stack trace, #0? 0x77830190 in ntdll!LdrFindResource_U () ?? from /cygdrive/c/Windows/system32/ntdll.dll building goes as follows, $ ./R CMD INSTALL --no-test-load nlme_3.1-97.tar.gz *

Hi! I'm having trouble accessing Samba shares via autofs. Let me begin by saying that I can access the samba shares directly using smbclient without any trouble. For example: smbclient //Matsa/Pisi works fine. I first started out by copying a config line from an older system (redhat 7) that has in its auto.misc: (all on one line, of course) pisi

Dear R-help members, I'm trying to run LME model on some behavioral data and need confirmations about what I'm doing... Here's the story... I have some behavioral reaction time (RT) data (participants have to detect dome kind of auditory stimuli). the dependant variable is RT measured in milliseconds. 61 participants were tested separated in 4 age groups (unblanced groups,

Hello list, i have a glm with family=binomial, link=logit but there is over-dispersion. So, in order to take into account for this problem i choose to do a glm with family=quasibinomial(). I'm not an expert on this subject and i ask if someone could validate my approach (i'm not sure for the tests) : quasi_glm(myformula,quasibinomial(),start=mystart) summary(quasi) # test t for

I've originally made 48 GLM binomial models and compare the AIC values. But dispersion was very large: Example: Residual deviance: 8811.6 on 118 degrees of freedom I was suggested to do a quasibinomial afterwards but found that it did not help the dispersion factor of models and received a warning: Residual deviance: 3005.7 on 67 degrees of freedom AIC: NA Number of Fisher Scoring

Full_Name: Henric Nilsson Version: 2.3.0 Patched (2006-05-09 r38014) OS: Windows 2000 SP4 Submission from: (NULL) (83.253.9.137) When supplying an unavailable link to `quasibinomial', the error message looks strange. E.g. > quasibinomial("x") Error in quasibinomial("x") : 'x' link not available for quasibinomial family, available links are "logit",

Dear all, I've got some questions, probably due to misunderstandings on my behalf, related to fitting overdispersed binomial data using glm(). 1. I can't seem to get the correct p-values from anova.glm() for the F-tests when supplying the dispersion argument and having fitted the model using family=quasibinomial. Actually the p-values for the F-tests seems identical to the p-values for

I am looking for the correct commands to do the following things: 1. I have a binomial logistic regression model and i want to test for overdispersion. 2. If I do indeed have overdispersion i need to then run a quasi-binomial model, but I'm not sure of the command. 3. I can get the residuals of the model, but i need to then apply a shapiro wilk test to test them. Does anyone know the command

After looking at 48 glm binomial models I decided to try the quasibinomial with the top model 25 (lowest AIC). To try to account for overdispersion (residual deviance 2679.7/68 d.f.) After doing so the dispersion factor is the same for the quasibinomial and less sectors of the beach were significant by p-value. While the p-values in the binomial were more significant for each section of the

Dear all, I try to consider overdispersion in a lmer model. But using family=quasibinomial rather than family=binomial seems to change the fit but not the result of an anova test. In addition if we specify test="F" as it is recomanded for glm using quasibinomial, the test remains a Chisq test. Are all tests scaled for dispersion, or none? Why is there a difference between glm and lmer

I am unable to calculate confidence intervals for the slope estimate in a quasibinomial glm using confint(). Below is the output and the package info for MASS. Thanks in advance! R 2.9.2 MASS 7.2-48 > confint(glm.palive.0.str) Waiting for profiling to be done... Error: dims [product 37] do not match the length of object [74] > glm.palive.0.str Call: glm(formula = cbind(alive, red) ~ str,

Hello R-list. I am a "long time listener - first time caller" who has been using R in research and graduate teaching for over 5 years. I hope that my question is simple but not too foolish. I've looked through the FAQ and searched the R site mail list with some close hits but no direct answers, so... I would like to estimate QAIC (and QAICc) for a glm fit using the

Hola a tod en s, mi compañero y yo intentamos ver la correlación de nuestros datos mediante regresiones logísticas. Trabajamos con proporciones (1 variable dependiente y 1 independiente) mediante modelos mixtos (los datos están agrupados porque hay pseudoreplicación). Hemos usado el paquete "lme4" y la función "lmer". Encontramos "overdispersion" en el resultado

I am running a binomial glm with response variable the no of mites of two species y->cbind(mitea,miteb) against two continuous variables (temperature and predatory mites) - see below. My model shows overdispersion as the residual deviance is 48.81 on 5 degrees of freedom. If I use quasibinomial to account for overdispersion the dispersion parameter estimate is 2501139, which seems

Dear R-Users, I can't understand the behaviour of quasibinomial in lmer. It doesn't appear to be calculating a scaling parameter, and looks to be reducing the standard errors of fixed effects estimates when overdispersion is present (and when it is not present also)! A simple demo of what I'm seeing is given below. Comments appreciated? Thanks, Russell Millar Dept of Stat U.

Sadly you failed to set your email program to send plain text and the data is corrupted at my end. I also think you need to reduce the size of the data set... the intent here is to increase your understanding, not debug your particular analysis. I will say that I am having a very challenging time understanding what you are trying to accomplish though. What are the equations that you think need