search for: treatm

...at in a 4x2 plot array On top of each plot I would like to have two lines. The 1st line will indicate what is plotted The 2nd line will indicate a a reference Graphically speaking I would like a result similar to this (simplified) one: par(mfrow=c(2,4)) for(aux in 1:4){ plot(1,1) mtext("1st treatm", cex=1, line=1.7, adj=0.1, col=1) mtext("2nd treatm", cex=1, line=1.7, adj=0.5, col=2) mtext("3rd treatm", cex=1, line=1.7, adj=0.9, col=3) mtext("1st treatm", line=0.6, adj=0.1, col=1) mtext("1st treatm", line=0.6, adj=0.5, col=1) mtext("1st treat...

...ated post on this list, but a full answer was never provided. I understand that the routine lmer (or lmer2) in the lme4 package is often recommended in such a case, but actually implementing it is where I've hit a wall. In a nutshell, the experiment involved presenting females from two groups (treatment, control) with an opportunity to mate with a virgin male every 6 hours for 48 hours. Every female was presented this opportunity at every time step (i.e., whether or not she mated at 6 hr, she was again presented with a male at 12 hr, and so on). In addition to which group a female belongs to, w...

...ulations). The response variables in the numbe of seconds invested in each probe. Further, I have plant floral display as a covariate, as it may influence visitation rates. I think I have to analyse this desing considering population, pollinator group and their interaction as fixed effects, and treatment nested within plant, and its interaction with population and pollinator group, as random factors. In SAS terminology, the model looks like this: proc mixed data=flwfunc.visitflower covtest method=reml; class site pollclass treatm plantid; model time = site|pollclass flwinflor / chisq; ran...

Dear list, I am new to this list and I am new to the world of R. Additionally I am not very firm in statistics either but have to deal. So here is my problem: I have a dataset (which I attach at the end of the post) with a binomial response variable (alive or not) and three fixed factors (trapping,treat,sex). I do have repeated measures and would like to include one (enclosure) random factor. I

...xperience with this test and with R,so I have some problems to use the leveneTest function. I read some previous posts on this topic but I cannot understand what I have to do. I have a five-columns data frame organised in the following way: number of observation, number of subject, score, type of treatment (I have 4 of them), number of replication (3 replications). My idea was: leveneTest(score~treatm*rep,data.ex,center=median) where: treatm is the treatments' column rep is the replications' column data.ex is the data frame. thank you giovanna -- View this message in context: http...

...ease of the lme4 package? Using simulated data for a quite standard mixed-model anova (a balanced two-way design; see code for the function SimMixed pasted below), I compared the output of lmer, for three slightly different models, with the output of aov. For an example where there is no fixed treatment effect (null hypothesis is true), with 4 blocks, 2 treatments, and 40 observations per treatment-block combination, I find that lmer gives more statistical significances than it should, whereas aov does not have this problem. An example of output I generated by calling > SimMixed(1000)...