search for: ntreat

...ces? (I haven't found anything relevant in Pinheiro and Bates, but I may not have looked in the right place ...) thanks Ben Bolker self-contained but unnecessarily complicated simulation code/demonstration: --------------- library(lme4) library(lattice) simfun <- function(reefeff,ntreat=2,nreef=12, nreefpertreat=3, t.eff=10, totvar=25,seed=NA) { if (!is.na(seed)) set.seed(seed) ntot = nreef*nreefpertreat npertreat=ntot/ntreat reef = gl(nreef,nreefpertreat) treat = gl(ntreat,npertreat) r.sd = sqrt(totvar*reefeff)...

Hi, Is there a way to stack bars in a barchart as well as "beside" bars for the same treatment? eg.... I have one barchart like this: bio<-matrix(c(10,23,9,25),nrow=2,byrow=T) ntreat<-c("n0","n96") colnames(bio)<-ntreat barplot(bio,beside=T) now i want a similar barchart but with stacked bars: ld<-matrix(c(25,25,50,10,30,60,15,35,50,30,30,40),nrow=4,byrow=T) ntreat2<-c("n0","n96","n0","n96") rownames(ld)&...

In S-Plus, I can obtain polynomial contrasts for an ordered factor with contr.poly(). The function also exists in R, however is limited to factors where the levels are equally spaced. In S-Plus, one can obtain the contrasts for a set of numeric values representing unequally spaced ordered factors. Has anyone implemented this in R? I see that the S-Plus function calls another function (poly.raw())

Dear All I would like to produce interaction boxplots and this seems to work: par(mfrow=c(2,2)) A=sample(rnorm(50,50,10)) B=sample(rnorm(50,100,10)) Test=merge(A,B,by=0)#by=0 where 0 is the row.names TreatA=(gl(2,50,100,labels=c("High","Low"))) TreatB=rep(gl(2,25,50,labels=c("High","Low")),2) Newdata=data.frame(TreatA,TreatB,Test)