[R] help with lme random = syntax and groupedData syntax

2010-12-02 Thread Huso, Manuela
Hello, all, I have a fairly complicated experimental design and would really appreciate some help on correctly specifying my random effects. I have a CRD split plot design with one covariate measured on the whole plot unit, one on the subplot unit and 2 other covariates measured on the subsam

Re: [R] How do I set the dispersion parameter in poisson glm?

2010-10-07 Thread Huso, Manuela
-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Ben Bolker Sent: Thursday, October 07, 2010 12:04 PM To: r-h...@stat.math.ethz.ch Subject: Re: [R] How do I set the dispersion parameter in poisson glm? Huso, Manuela oregonstate.edu> writes: > > Thank you so much for your

Re: [R] How do I set the dispersion parameter in poisson glm?

2010-10-07 Thread Huso, Manuela
rdson Hall Department of Forest Ecosystems and Society Oregon State University Corvallis, OR 97331 ph: 541-737-6232 fx: 541-737-1393 -Original Message- From: Achim Zeileis [mailto:achim.zeil...@uibk.ac.at] Sent: Thursday, October 07, 2010 9:48 AM To: Huso, Manuela Cc: r-help@r-project.org

[R] How do I set the dispersion parameter in poisson glm?

2010-10-07 Thread Huso, Manuela
Dear R users, I would like to fit a glm with Poisson distribution and log link with a known dispersion parameter. I do not want to estimate the dispersion parameter. I know what it is, so I simply want to fix it at a constant for this and other models to follow. My simple, no covariate model

[R] fixing the dispersion parameter in glm

2010-09-28 Thread Huso, Manuela
I would like to fit a glm with Poisson distribution and log link with a known dispersion parameter. I do not want to estimate the dispersion parameter. I know what it is, so I simply want to fix it at a constant for this and other models to follow. My simple, no covariate model is: Tall.glm<

Re: [R] in continuation with the earlier R puzzle

2010-07-12 Thread Huso, Manuela
Using Ted Harding's example: news1o <- runif(100) s2o<- runif(100) pt1 <- proc.time() s <- numeric(length(news1o))-1 # Set all of s to -1 s[news1o>s2o] <-1# Change to 1 only those values of s # for which news1o>s2o pt2<-

Re: [R] parameterization of glm nested design

2010-07-01 Thread Huso, Manuela
onses are below: -Original Message- From: Charles C. Berry [mailto:cbe...@tajo.ucsd.edu] Sent: Thursday, July 01, 2010 9:30 AM To: Huso, Manuela Cc: 'r-help@r-project.org' Subject: Re: [R] parameterization of glm nested design On Wed, 30 Jun 2010, Huso, Manuela wrote: > Dear

[R] parameterization of glm nested design

2010-06-30 Thread Huso, Manuela
Dear R community, I am new to R, a reforming SAS user :) I am running R 2.10.1 on a Windows XP machine. I would like to write linear functions of my coefficient parameter estimates from a glm, but am having a difficult time understanding the parameterization R uses. In the toy example belo