similar to: using hglm to fit a gamma GLMM with nested random effects?

Addendum: I tried a gamma fit in glmmPQL and got the same errors. *Ben Caldwell* PhD Candidate University of California, Berkeley On Tue, May 17, 2011 at 3:51 PM, Benjamin Caldwell <btcaldwell@berkeley.edu>wrote: > Hello > After seeing this ( > https://stat.ethz.ch/pipermail/r-sig-mixed-models/2011q1/005213.html) email > I thought I would check the issue with a gamma family

Sub: using glmer to fit a mixed-effects model with gamma-distributed response variable Hello, I'm currently trying to fit a mixed effects model , i.e.: > burnedmodel1.2<-glmer(gpost.f.crwn.length~lg.shigo.av+dbh+leaf.area+ bark.thick.bh+ht.any+ht.alive+(1|site/transect/plot), family=gaussian, na.action=na.omit, data=rws30.BL) If I run this code, I get the error below: Error:

Hi all, I'm trying to fit models for data with three levels of nested random effects: site/transect/plot. For example, modelincrBS<-glmer(l.ru.ba.incr~shigo.av+pre.f.crwn.length+bark.thick.bh+Date+slope.pos.num+dens.T+dbh+leaf.area+can.pos.num+(1|site/transect/plot), data=rws30.UL, family=gaussian, na.action=na.omit) but I get the following error: Error: length(f1) == length(f2) is not

Hello, I'm trying to compare the fit of two distributions, normal and gamma, to a histogram of my response variable. rate<-mean(na.omit(rwb$post.f.crwn.length))/var(na.omit(rwb$post.f.crwn.length)) shape<-rate*mean(na.omit(rwb$post.f.crwn.length)) hist((rwb$post.f.crwn.length), main="rwb$post.f.crwn.length")

I wonder if someone has tried to implement the hierarchical generalized linear model (HGLM) approach of Lee and Nelder (JRSSB, 1996, 58: 619-56) in R. Thanks in advance. Emmanuel Paradis At 17:18 01/04/02 +0100, ripley at stats.ox.ac.uk wrote: >On Mon, 1 Apr 2002, Jason Liao wrote: > >> Happy new month, everyone! >> >> I am planning to write a NIH grant proposal to study

Hi folks, I'm dealing with a dataset with a lot of NAs, and want to use subset on the data without removing the NAs from the the dataframe; e.g. rws50 <- subset(rw.fire.RW,shigo.av<50) This removes instances with NA for the value in addition to anything <50. Any suggestions much appreciated. *Ben Caldwell* [[alternative HTML version deleted]]

Does anyone have advice on how to instruct R CMD check to use a non-standard set of libraries? Here's the situation: I'm trying to do some automated checking on package dependencies of a package I maintain. In order to do that I've written code that takes the list of the dependent packages and for each package (1) downloads the most recent/available .tar.gz file; (2) installs the

R helpers, I am working with the R program performing multivariate regression trees (MRT). I have a matrix with species and environmental variables saved as a CSV file (sprot_matrix.csv), I have 42 species and 8 environmental variables (SECCHI+PH+TA+PTOT+NTOT+CHLA+AREA+ MEANDEP) for 104 samples Title SpA SpB SpC SpD Varible1 Variable2 Variable3 Sample1 Sample 2

Hi I am trying to convert SAS codes to R, but some of the result are quite different from SAS. When I ran proc mixed, I have an option ddfm=bw followed by the model. How can I show this method in R?(I am thinking that this maybe the reason that I can't get the similar results) below is my SAS codes: proc mixed data=test covtest empirical; class pair grade team school; model score = trt

Craig A Faulhaber wrote: >I am interested in using a generalized linear mixed model with data > that best fits a beta distribution (i.e., the data is bounded between > 0 and 1 but is not binomial). .. >For clarification, here's what I'm trying to model: >I have a beta-distributed response variable (y). I have a fixed-effect >explanatory variable (treatment),

CRAN (and crantastic) updates this week New packages ------------ * apcluster (1.0.1) Ulrich Bodenhofer http://crantastic.org/packages/apcluster The apcluster package implements Frey's and Dueck's Affinity Propagation clustering in R. The algorithms are analogous to the Matlab code published by Frey and Dueck. * BioPhysConnectoR (1.6-1) Franziska Hoffgaard

If I switch from fcgi to the cgi dispatcher the stylesheets work? Anyone have any ideas? Thanks in advance, Scott Thanks, Scott Caldwell Database Administrator City of Fayetteville 113 W Mountain St Fayetteville, AR 72701 (479) 575-8329 phone (479) 575-8316 fax http://www.accessfayetteville.org/ http://www.faygis.org/ >>>

Taken from http://www.twice.com/html/pagebeta.cfm?InputKey=2853 ------------------------------------------------------------------------- Thomson, Fraunhofer Developing Enhanced MP3 Codec Called mp3PRO Jan. 7, 2001 By Joseph Palenchar Las Vegas -- Thomson and the Fraunhofer Institute, the co-developers of MP3, are working with a third company to develop an improved MP3 codec intended to deliver

Hi I have used fa() to perform a factor analysis of a psychological battery which is thought to have 11 factors. I can identify which factors the loadings relate to easily enough because I can see which items are loading onto each of the columns in the $loading output. However, how can I identify which items or loadings are being used to create each of the columns in the $scores output? I have

PJ, I didn't have virsh installed (libvirt-client).? After having installed it, it still didn't see the process: virsh autostart testvm error: failed to connect to the hypervisor error: no valid connection error: Failed to connect socket to '/var/run/libvirt/libvirt-sock': No such file or directory Then, after having installed libvirt, it gives essentially the same error.

*Dear R-help: My name is Rodrigo and I have a question with nlme package in R to fit a mixed beta regression model. The details of the model are: Suppose that:* *j in {1, ..., J}* *(level 1)* *i in {1, ..., n_j}* *(level 2)* *y_{ij} ~ Beta(mu_{ij} * phi_{ij}; (1 - mu_{ij}) * phi_{ij}) y_{ij} = mu_{ij} + w_{ij} * *with* *logit(mu_{ij}) = Beta_{0i} + Beta_{1i} * x1_{ij} + b2 * x2_{ij}

I have a problem in which the Windows client isn't allowed to change the password in a direct manner. Using the Windows control panel, Passwords, I am only allowed to change the password by first changing the existing password to a blank password and then changing the blank password to the desired password. If I try to change the password directly I generate the error

Evening folks, I'm trying to print a series of graphs to .jpeg using a variable as the title, but run into the difficultly that I can't find a way to append the file extension to the .jpeg (in this case extensionless!) files. Example: ---- species.name="CussoniaHolstii" dia<-10:100 biomass = -21.4863 + 0.5797 * (dia ^ 2) biomass jpeg(species.name) plot (biomass,

I have some elements in a vector with extraneous information (e.g. file name and sample IDs) that I'd like to strip from every element. For example, I would like "SPI1.S1.str1.P3.sample.tif" "SPI1.S1.STR2.P1.sample.tif" to read "SPI1.S1.str1.P3" "SPI1.S1.STR2.P1". Will someone help me with the syntax in gsub? It needs to be something like

i haven't heard of problems with configure/lastlog on 10.20. i'll copy the development list as there are folks there who use 10.20. On Thu, 6 Sep 2001, Jesse Caldwell wrote: :one last problem with openssh -- building openssh on HPUX 10.20 fails :on loginrec.c (i don't have the errors handy but i would be more than :happy to recreate them). i get the same error with 2.9p2 and 9/6/2001