similar to: random effects in analysis of variance

One of our students did something like: x[a < b < c] instead of x[a < b & b < c] But why is 3 < 2 < 1 # [1] TRUE ??? Is there any reason? Or wouldn't it be better to get a warning / error? Uwe Ligges -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html

Please pardon an extremely naive question. I see related earlier posts, but no responses which answer my particular question. In general, I'm very confused about how to do variance decomposition with random and mixed effects. Pointers to good tutorials or texts would be greatly appreciated. To give a specific example, page 193 of V&R, 3d Edition, illustrates using raov assuming pure

I am working with a RCB 2x2x3 ANCOVA, and I have noticed a difference in the calculation of sum of squares in a Type III calculation. Anova output is a follows: > Anova(aov(MSOIL~Forest+Burn*Thin*Moisture+ROCK,data=env3l),type=3) Anova Table (Type III tests) Response: MSOIL Sum Sq Df F value Pr(>F) (Intercept) 22.3682 1 53.2141 3.499e-07 *** Forest

Dear all, In a greenhouse experiment we tested performance of 4 different species (B,H,P,R) under 3 different water levels in 10 replications. As response variable e.g. the number of emerging sprouts were measured on three dates. A simple Anova considering every measurement date separately shows a higly significant effect of species and moisture (and partly the interaction of both). The

I am trying to it a particular nonlinear model common in Soil Science to moisture release data from soil. I have written the function as shown below according to the logist example in Ch8 of Pinheiro & Bates. I am getting the following error (R version 2.1.1) *Error in qr(attr(rhs, "gradient")) : NA/NaN/Inf in foreign function call (arg 1)* Below is the function and data. /#

Hello all R listers, I'm struggling to select an appropriate statistical method for my data set. I have collected soil moisture measurements every hour for 2 years. There are 75 sensors taking these automated measurements, spread evenly across 4 treatments and a control. I'm not interested in being able to predict soil future soil moisture trends, but rather in knowing whether the

Hi All, I am trying to lable the y-axis on my scatterplot with the following: "Soil moisture content (m3m-3)" I am using the following coding for plotting the graph: plot(soilmoisture~gradientlevel, xlab="Levels of droughting gradient", ylab="Soil moisture content (m3m-3)", bty="l", font.main="2", pch=16, las=1, cex.lab="1.13") I have

Hi everyone, I have a question about calculating r-squared in R. I have tried searching the archives and couldn't find what I was looking for - but apologies if there is somewhere I can find this... I carried out a droughting experiment to test plant competition under limited water. I had: - 7 different levels of watering treatment (1 -7 - from most watered to least watered/) - 15

Hello, For an experiment, I selected plots of land within a forest either with honeysuckle or without honeysuckle. Thus, my main factor is fixed, with 2 levels: "honeysuckle present"(n=11) and "honeysuckle absent"(n=8). Within each plot of land, I have a "trenched" subplot and an "untrenched" subplot. Within each subplot of every plot, I measured soil

Hello, I am very new to R but would like to use the software to analyse the attached data. The experiment followed a split-split plot design There were two blocks and the whole plot is CO2 with two levels. The sub-plot is soil temperature with three levels and the sub-sub plot is soil moisture content with three levels (low, intermediate and high-similar for soil temperature). I had 7 plants per

Hi All: I am working with a dataset on Arsenic toxicity, and I am trying to calculate the 20th, 40th, 60th, 80th, and highest percentiles for a variable, dietary Moisture (variable name dMoist). The inbuilt function quantile(dMoist) would print 0, 25th, 50th, 75th, and 100th percentile. Does there exist a function that can calculate xth percentile (where x = 10th, 20th, ... etc) values?

Hi, I want to know the area under a curve, which is not given as a function, but as values in a time series. It is not a smooth curve, but switches often between positive values and zero (the values describe the moisture in the soil over a year, one entry is one day). I already tried area.between.curves, but got only 0 as result. I guess, it doesn't work because of these multiple changes

Hi, For my master thesis I have 24 micro-plots on which I did measurements during 3 months. The measurements were: - Rainfall and runoff events throughout 3monts (runoff being dependant on the rainfall, a coefficient (%) has been made per rainfall event and per 3 months) - Soil texture (3 different textures were differentiated) - Slope (3 classes of slopes) - Stoniness (one time measurement)

I have created a scatter plot that has come out okay but I am having trouble with the x axis. My data consists of 4 treatments but these treatments are days so R keeps reading them as numeric and making my x axis continuous. Here is what I have so far: plot(pair$MC~pair$Day, pch=c(19,24)[as.factor(Cookie)], main='Paired t Test', xlab='Days in Field', ylab='Moisture

I have nc files for global soil moisture,here is one file https://echange-fichiers.inra.fr/get?k=f9DDllPKdUKs5ZNQwfq from the metadata ,the projection is cylindrical and the resolution is 25 km(it is based on authalic sphere based on International 1924 ellipsoid).As I want to compare with other data, I have to make them identical. - my other data are in WGS84 with 0.25*0.25 degree resolution

Hi everyone! What's the best filesystem to use for the mail spool? I'm debating between Debian with xfs or FreeBSD with zfs. I'm not sure which way to go. I'm migrating from cyrus. I have about 50 users so it's not a large setup. Cheers! Monika -- Monika Janek Systems Administrator, Side Effects Software Toronto, Ontario Canada 416-504-9876 x207 www.sidefx.com

I would like to calculate the mean of tree leader increment growth over 5 years (I1 through I5) where each tree is a row and each row has 5 columns. So far I have achieved this using rowMeans when all columns are numeric type and used in the calculation: Data1 <- data.frame(cbind(I1 = 3, I2 = c(0,3:1, 2:5,NA), I3 =c(1:4,NA,5:2),I4=2,I5=3)) Data1 Data1$mean_5 <- rowMeans(Data1, na.rm =T)

Dear members, I have started to work with R recently and there is one thing which I could not solve so far. I don't know how to define macros in R. The problem at hand is the following: I want R to go through a list of 1:54 and create the matrices input1, input2, input3 up to input54. I have tried the following: for ( i in 1:54) { input[i] = matrix(nrow = 1, ncol = 107)

I have site data with variables that vary across sites (like wind speed, moisture content) and within sites I have experimental units (logs) with associated variables (like decay class, suspension). With normal response y one can use R to get the between and within site information using aov(y~u+v+Error(site)) where u and v are independent variables. Problem is we have a binary response (burn or

Hi, I'm new to R. I have some daily soil moisture data for the year 1979 in netcdf format such as these sm19790101.1.nc sm19790102.1.nc . . . sm19791231.1.nc I need to average a variable called "sm" to monthly resolution. I've done these days = formatC(1:31, width=2, flag="0") ncfiles = lapply(days, function(d){ filename = paste("sm197901", d,