similar to: use aov or lme for split plot design?

Windows 2000 . 5.00.2195 with Service Pack 1. R 1.5.1 Output from my split-split plot aov "alerted" me that I have done something wrong. I designed an experiment with all combinations of all levels of each treatment, but lost a little data (3 out of 192 plots). With the following data, I run the following model: > collim[c(1:6,187:192),c(1,3:6,9)] plot Litter Fert

Hi, I've data from an incomplete fatorial design. One level of a factor doesn't has the levels of the other. When I use lm(), the summary() return NA for that non estimable parameters. Ok, I understant it. But I use contrast::contrast(), gmodels::estimable(), multcomp::glht() and all these fail when model has NA estimates. This is becouse vcov() and coef() has different dimensions. Is

Hi R users, I sent a message yesterday about NA in model estimates ( http://r.789695.n4.nabble.com/How-set-lm-to-don-t-return-NA-in-summary-td3722587.html). If I use aov() instead of lm() I get no NA in model estimates and I use gmodels::estimable() without problems. Ok! Now I'm performing a lot of contrasts and I need correcting for multiplicity. So, I can use multcomp::glht() for this.

Windows 2000 (v.5.00.2195), R 1.5.1 I have an lme object, fm0, which I examine with anova() and intervals(). The anova output indicates one of the interaction terms is significant, but the intervals output shows that the single parameter for that term includes 0.0 in its 95% CI. I believe that the anova is a conditional (sequential) test; is intervals marginal or approximate? Which should I

Dear R community, I am trying to fit an nlme model where I want to estimate the fixed effects of two treatments on the parameters on the following equation Photo~(a*(1-exp(-c*PARi/a)))-b I was able to fit a simple model without covariates following the method described in Mixed-Effects Methods and Classes for S and S-PLUS, version 3.0, but when I add the covariates, I get the error "

Dear R sages, I used the function rbind to combine a matrix (M) and a vector (Fert) to get a new matrix (A). This was fine. The issue is however, that the new matrix A has as its row names the name of the vector Fert, even though I set teh new vector A to have dimnames=NULL. See short code below fyi. *Fert<-c(0,1,5) * *M <- matrix(0, 2, 3) diag(M) <- c(0.3,0.5) * *A<-

Hi Rneters, I'm running a GLM model with a full factorial design in blocks and binomial error distribution. I would like to have the F values for this model but I got a message that "using F test with a binomial family is inappropriate in: anova.glm(model, test = "F")". Should I not report F statistics on this kind of analysis? I would appreciate any comment on this.

Hi, I'm trying to fit a logit model to a set of data that was collected under a split plot scheme. The structure of the data is Whole plot factors: Watering Frequency (2 levels: Hi/Lo) and Fertilizer type (3 factors A/B/C) Subplot factor: slope type (2 factors up/down) Response: Proportion of infected leaves(Infected leaves/Total leaves) of the plant (2 plants recorded from each plot)

Hi R-Helpers: I?m trying to use the function GLMM from lme4 package, (R-1.8.1, Windows 98),and I get the following error: > pd5 = GLMM(nplant~sitio+ + fert+ + remo+ + sitio:fert+ + remo:sitio+ + remo:fert+ + remo:fert:sitio + data=datos, + family=binomial, +

Using R v. 1.7.0 on Windows 2000 I would like to plot the fitted values of a model as a function of a continuous covariate, augmented with data (e.g., augPred) grouping by combinations of fixed effects. I have not been able to use augPred effectively, and am wondering if it does not handle unbalanced data (3 out of 192 missing). I include below the model and an xyplot that almost does the

Hello, I am a Masters student and I am working on my thesis modelling smallholder farms using a program in R. I have modified the original code and I am having some issues with lists that I cannot figure out. Originally, I had list file defining lists such as: Param, Crop1, Crop1, Soil, etc. (ex. Param <- list() ). Their subsets were listed as Crop1$CContent for example and there was quite

Hi I have a file saved in R, named agrexp.Rdata, shown below > agrdata fert yield 1 25 84 2 50 80 3 75 90 4 100 154 5 125 148 If I double clicked on this file, the data is displayed without problem. However if I tried to import using: > agrdata<-read.table("agrexp.Rdata") or >

I am trying to run the following model in R > lmer(leaves.eaten~Geocytotype+(1|TEST/ PLANT),data=cyphoplantfeeding,family=poisson) My experimental setup is 41 replicates (TEST) of an experiment in which there are three Geocytotypes of a plant species in each TEST, and two plant pseudoreplicates per Geocytotype in each test (i.e. 3*2=6 plants per test). So my random factors are trying

Hello, I am trying to help someone who has carried out an experiment and I'm finding it quite difficult to understand the appropriate model to use & code it. The response is a measurement - the amount of DNA extracted during the experiment. There were 2 factors to be tested - one is the condition under which the experiment took place and the other is the type of DNA to be

Dear R-help-list, I have a problem in which the explanatory variables are categorical, the response variable is a proportion, and experiment contains technical replicates (pseudoreplicates) as well as biological replicated. I am new to both generalized linear models and mixed- effects models and would greatly appreciate the advice of experienced analysts in this matter. I analyzed the

As an agronomist, sometimes i need to perform statistical analisys on "weird" field data, like 3 or 4 level split plot design. I've read both the help of the program, some books and manual about S and "Venables/Ripley, Modern Applied Statistics with S", but i still have problems in (i suppose) formulae sintax. I'm trying to perform AOV, as explained showed in

As an agronomist, sometimes i need to perform statistical analisys on "weird" field data, like 3 or 4 level split plot design. I've read both the help of the program, some books and manual about S and "Venables/Ripley, Modern Applied Statistics with S", but i still have problems in (i suppose) formulae sintax. I'm trying to perform AOV, as explained in "Random

Dear list, I was wondering if there is a function to plot league tables, sometimes also known as "caterpillar plots"? A league table is conceptually very similar to a box plot. One difference is that the inter-quartile ranges are not shown. If there isn't such a function a first attempt for a "selfmade" plot would be to tell boxplot not to plot boxes (sounds silly

________________________________ Principal Post-Doctoral Fellow PKPD modelling platform (full time) Leiden University, Leiden, The Netherlands <javascript:doredirect();> Job description ________________________________ TOP Institute Pharma (TI Pharma) has granted our proposal to set up a mechanism-based PK-PD modelling platform. This platform focuses on the transfer of

Post-Doctoral Fellowship Modelling of Cardiovascular Safety (full time) Leiden University, Leiden, The Netherlands <javascript:doredirect();> Job description TOP Institute Pharma (TI Pharma) has granted our proposal to set up a mechanism-based PKPD modelling platform. This platform focuses on the transfer of knowledge from academia to the pharmaceutical industry and is a