search for: wheat2

...emove data from my dataset for which the runoff is larger than the rainfall on the same day. (Don't worry, the rainfall data is in fact the rainfall associated to the three days around the runoff event). On days that no runoff was recorded the data includes "NA" The dataset is called wheat2 I use the following: for (i in 1:nrow(wheat2)) { # Insert NA's if runoff data > rainfall data if (is.na(wheat2[i,5:7])==FALSE && any(wheat2[i,5:7]>wheat2[i,8])) { wheat2[i,5:7] <- NA } } I tried this: wheat1 <- replace(wheat2[,5:7],is.na(wheat2[i,5:7...

...ke to make the points different colors for various levels of a variable. I have tried a simple 'if' statement in the plot command, but I get an error message. Here is what I have tried and the error message I get: > plot(test, if(test[,1]<8) col="steelblue2" else col="wheat2")Error in xy.coords(x, y, xlabel, ylabel, log) : 'x' and 'y' lengths differIn addition: Warning message:In if (test[, 1] < 8) col = "steelblue2" else col = "wheat2" : the condition has length > 1 and only the first element will be usedI know 'x&...

...ed argument(s) (method ...) It looks like I might not have loaded a necessary package. I did load nls and lattice (the latter not actually being listed in the package description as a dependency). I checked the R-help archives thoroughly and also examined the library nlme scripts. The chapter 5 Wheat2 example from Pinheiro & Bates (2000) didn't work either: > library( nlme ) > library( lattice ) > data( Wheat2 ) > fm1Wheat2 <- gls( yield ~ variety - 1, Wheat2 ) > Variogram( fm1Wheat2, form = ~ latitude + longitude ) Error in as.data.frame.default(x[[i]], optional = TRU...

...of the same type of plant, are together in its own block of boxplot, because each plant has its own unit of measure. Let's say wheat is measured in ton, tomato in pound and cucumbers as counts. :-) x<-rnorm(1000,mean=0,sd=1,main="Right screen") wheat1<-rnorm(100,mean=0,sd=1) wheat2<-rnorm(150,mean=0,sd=2) tomatos3<-rnorm(200,mean=0,sd=3) tomatos4<-rnorm(250,mean=0,sd=4) cucumbers5<-rnorm(300,mean=0,sd=5) cucumbers6<-rnorm(400,mean=0,sd=6) par(mfrow=c(1,2)) hist(x, main="Left screen OK") boxplot(wheat1,wheat2,tomatos3,tomatos4,cucumbers5,cucumbers6)...

I'm looking at Pinheiro & Bates "Mixed-Effects Models in S and S-PLUS" at the moment. Several datasets are used, one of which is called "PBIB" (a partially balanced incomplete block design). All the other datasets can be found somewhere or other in R. However, I cannot locate PBIB, and it does not seem to be mentioned in the latest edition of the R Full Reference

...aize : num 5.8 4.7 6.4 4.9 6 5.2 $ Oats : num 8.3 6.1 7.8 7 5.5 7.2 > stack <- function(data) + data.frame(values = unlist(data), + lev = factor(rep(names(data), lapply(data, length)))) > stack(ex10.09) # try it on the example values lev Wheat1 5.2 Wheat Wheat2 4.5 Wheat Wheat3 6.0 Wheat Wheat4 6.1 Wheat Wheat5 6.7 Wheat Wheat6 5.8 Wheat Barley1 6.5 Barley Barley2 8.0 Barley Barley3 6.1 Barley Barley4 7.5 Barley Barley5 5.9 Barley Barley6 5.6 Barley Maize1 5.8 Maize Maize2 4.7 Maize Maize3 6.4 Maiz...

Hi, I'm fitting a relatively simple growth model to some forest plot data. Two species of trees were planted in different mixtures in 10 (nearly-adjacent) plots and measured on four occasions over 10 years. The model is constructed in terms of the diameter increments (per year; DI) in the 3 intervals, in which DI is modelled as a function of mid-interval D and DSQ. The details of the

...;, "violetred1", "#FF3E96", "violetred2", "#EE3A8C", "violetred3", "#CD3278", "violetred4", "#8B2252", "wheat", "#F5DEB3", "wheat1", "#FFE7BA", "wheat2", "#EED8AE", "wheat3", "#CDBA96", "wheat4", "#8B7E66", "whitesmoke", "#F5F5F5", "yellow", "#FFFF00", "yellow1", "#FFFF00", "yellow2", "#EEE...

  Hola,   Por favor, ¿podríais indicarme qué recursos (librerías o ideas) pueden resultar de utilidad para crear un gráfico del estilo del de la figura 3.8 del siguiente link?   http://www.tsc.uvigo.es/BIO/Bioing/ChrLDoc3.html#3.5   Actualmente estoy utilizando funciones muy básicas y la verdad es que no me encuentro muy satisfecha con el resultado.   Muchas gracias.   Eva [[alternative HTML