similar to: Negative binomial theta

Hello All, Is there a way of producing an ANOVA table split into contrasts, thus showing the contrasts sums of squares and associated p-values? Thanks, Martin. Martin Hoyle, School of Life and Environmental Sciences, University of Nottingham, University Park, Nottingham, NG7 2RD, UK Webpage: http://myprofile.cos.com/martinhoyle

Dear R Users, I suppose this is a school boy question, but here it is anyway. I'm trying to re-create the residuals for a poisson GLM with simulated data; x<-rpois(1000,5) model<-glm(x~1,poisson) my.resids<-(log(x)- summary(model)$coefficients[1]) plot(my.resids,residuals(model)) This shows that my calculated residuals (my.resids) are not the same as residuals(model). p 65 of

Hello R Users, I am investigating the basic use of the LME function, using the following example; Response is Weight, covariate is Age, random factor is Genotype model.lme <- lme (Weight~Age, random=~ 1|Genotype) After summary(model.lme), I find that the estimate of Age is 0.098 with p=0.758. I am comparing the above model with the AOV function; model.aov <- aov (Weight~Age + Genotype)

Hi, Please can someone tell me if this is correct for significance tests on count data; 1: If overdispersed, use quasipoisson, drop1(model, test="Chi"), 2: If not overdispersed, use (poisson or quasipoisson), drop1(model, test="Chi"). Thanks for your time, Martin. Martin Hoyle, School of Life and Environmental Sciences, University of Nottingham, University Park, Nottingham,

Dear Greg, Thanks for the new release. The decomposition of the SSQ is just what I need! Regards, Martin. Martin Hoyle, School of Life and Environmental Sciences, University of Nottingham, University Park, Nottingham, NG7 2RD, UK Webpage: http://myprofile.cos.com/martinhoyle >>> gregory_r_warnes at groton.pfizer.com 10/30/02 07:16PM >>> Version 0.7.3 of the gregmisc package

Dear All, I'm not sure of the interpretation of interactions with contrasts. Can anyone help? I do an ANCOVA, dryweight is covariate, block and treatment are factors, c4 the response variable. model<-aov(log(c4+1)~dryweight+treatment+block+treatment:block) summary(model); Df Sum Sq Mean Sq F value Pr(>F) dryweight 1 3.947 3.947 6.6268 0.01076 *

Dear All, I'm rather a beginner on nested ANOVAs, so here goes with my 2 questions; Qu 1: I'm modelling the number of galls on a leaf (the response variable) as a function of; the tree on which I find the leaf, the branch on which I find the leaf. Then, the tree and the branch are both random factors, and I'm quite happy that I should write; aov(galls~tree/branch +

Hi, Can someone tell me how to normalize a data set so that the mean of the set is 0 and the variance is 1. As I understand, when you calculate the principle components of a data set through correlation as < princomp( dataset, cor=T ) > then a similar calculation is performed. I would like to know how I can perform such a calulation directly. Any help would be greatly appreciated. Many

Hi all, I'm having some troubles with the Quantreg package. I am using R version 2.10.0, and have downloaded the most recent version of Quantreg (4.44) and SparseM (0.83 - required package). However, when I try to run an analysis (e.g. fit1<-rq(y~x, tau=0.5)) I get an error message saying that the function "rq" could not be found. I get the same message when I try to search

Hi all, Has anyone tested for FA in R? I need to seperate out the variance due to measurement error from variation between individuals (following Palmer & Strobeck 1986). Andy Higginson Animal Behaviour and Ecology Research Group School of Biology University of Nottingham NG7 2RD U.K. This message has been checked for viruses but the contents of an attachment may still contain

Hello all. I have a simple barplot with sixteen different segments. When I plot my data, only five or six of the labels are showing in the x-axis. How do go get them all to show? Can I set them at a 45.degree angle? Thank you. Jess -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send

I want to form a 3x3 crosstabulation for the signs of two vectors (i.e. Negative, Zero, Positive). The problem is that I am simulating the data so for some iterations one of the categories is absent. Thus the resulting table shrinks to 3x2. I want it to be 3x3 with zero column corresponding to the missing category. Moreover, I have tried but failed to give the dimension names. -- Sohail Chand

Dear All, I wonder if anyone can advise me as to whether there is a consensus as to how the effect size should be calculated from GLIM models in R for any specified significant main effect or interaction. In investigating the causes of variation in infection in wild animals, we have fitted 4-way GLIM models in R with negative binomial errors. These are then simplified using the STEP procedure,

Dear R users, I've just upgraded from R 0.64.1 to 1.1.0 on a PC running Red Hat Linux 6.0. When I ran demo(graphics), after the command: title(main = "January Pie Sales", cex.main = 1.8, font.main = 1) the plot did not materialise and I received this error message: Error in title(main = "January Pie Sales", cex.main = 1.8, font.main = 1) : X11 font at size 22 could

Hi guys, I am updating a package because of data format in data folder. So I just change an extension of a file to .txt ... nothing more. I get this error on the R CMD check ** help *** installing help indices ** building package indices ... Error in read.table(zfile, header = TRUE, as.is = FALSE) : more columns than column names ERREUR : installing package indices failed Note that the R CMD

I am new to R. I am going to by one of the following book: 1. William N. Venables and Brian D. Ripley. Modern Applied Statistics with S-Plus. Third Edition. Springer, 1999. ISBN 0-387-98825-4. 2. The Fourth Edition of the book from point 1. 3. `S Programming' by W. N. Venables and B. D. Ripley Springer. ISBN 0-387-98966-8, 2000. I can only by one of the above books. Q1. I have found

Removing rows with NAs, using na.omit(), doesn't seem to be working for me. Dataset: > str ( ex10s ) 'data.frame': 2189576 obs. of 5 variables: $ LOPNR : int 58 58 58 58 64 64 64 64 64 64 ... $ DIAGNOS: Factor w/ 173 levels "F20","F200","F2000",..: 128 128 128 128 105 105 105 160 105 105 ... $ X_DATE : int 20060821 20061207 20080102 20090904

Hi, I have access to two providers. On one of them the authuser is the same as the username, so outgoing works. On the other one I can only get incoming - what ever combination I try for outgoing I get an error. The register command has the ability to specify both usernames (which is why incoming works) but outgoing doesn't seem to, and without that I'm stuck. They are defined as:

It worries me how many threads talk about low performance. I'm about to build out a replica 3 setup and run Ovirt with a bunch of Windows VMs. Are the issues Tony is experiencing "normal" for Gluster? Does anyone here have a system with windows VMs and have good performance? *Vincent Royer* *778-825-1057* <http://www.epicenergy.ca/> *SUSTAINABLE MOBILE ENERGY SOLUTIONS*

The asterisk wiki states that it needs SIP, IAX2, IAX and RTP open to the world to work. Is this necessarily true, or does it only need some of these outgoing? I'm concerned as anyone that could guess an extension number&password could use my server to make outgoing calls. It would help if the extensions had a netmask/allowable IP setting like the iax.conf file uses, but there