search for: moisture

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 1.0954 2 1.3029 0.29282 Burn 2.6926 1 6.4058 0.01943 * Thin 0.0494...

...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 mixed-effects model with species and moisture shows a highly significant effect of species and moisture as well. However, when I included the interaction the t-values of the species dropped strongly and the SE increase and the results for the species are no...

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 func...

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 hav...

I have a data.frame paint (below-mentioned). In SPlus I used the command raov( MOISTURE ~ BATCH / PROBE ) (you could do raov( MOISTURE ~ BATCH + PROBE%in%BATCH) as well) so that the factors are taken as random. In R this function raov doesn't exist. How can I calculate the same? Maybe with lme, but how? str(paint) ---------- `data.frame': 60 obs. of 5 variables: $ BATCH...

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 treatment affected soil moisture respon...

...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 replicates at each level. Soil moisture readings were taken 4 times throughout the experiment (so I have 105 readings for each of the 4 times) and I now want to check that there was a significant decrease in soil moisture as I decreased the watering frequency, i.e. watering level 7 showed lower soil moisture units than level 1. I have c...

Dear All, I need to separate one column which is in a format like, "can_region", into two columns like, "can", and "region". Do we have a function in R to do this? Thanks. Regards, Jin [[alternative HTML version deleted]]

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 treatment combination and 252 plants in total. CO2, soil temperature and moisture are fixed factors. Blocks are always random. What package do I need to perform a split-split plot analysis? Can s...

...hout 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 moisture. Now I need to do a nested Anova to compare the soil moisture values between treatments. I don't really want to discard some data to make the sample sizes balanced. Does anyone know how to do a nested, unbalanced Anova in R? Many thanks for your help! -Steven Pfeiffer [[alternative HTML ve...

Good morning gentlemen! How using a weighted model in nls2? Values with the nls are logical since values with nls2 are not. I believe that this discrepancy is due to I did not include the weights argument in nls2. Here's an example: MOISTURE <- c(28.41640, 28.47340, 29.05821, 28.52201, 30.92055, 31.07901, 31.35840, 31.69617, 32.07168, 31.87296, 31.35525, 32.66118, 33.23385, 32.72256, 32.57929, 32.12674, 52.35225, 52.77275, 64.90770, 64.90770, 85.23800, 84.43300, 68.96560, 68.41395, 70.82880, 71.18400, 96.13240...

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? I looked for such functions in the documentation, but couldn't find on...

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 between 0 and positive values (most of the time, the values are 0 and in the certain case I tested, the positive values were...

...t textures were differentiated) - Slope (3 classes of slopes) - Stoniness (one time measurement) - Random roughness (throughout 3 months) - Land use (crop land or grazing land) - Vegetation cover (throughout 3 months) - Vegetation height (throughout 3 months, only measured on cropland) - Antecedent moisture content (throughout 3 months) Now I would like to investigate the effect of all these variables on the rainfall/runoff. For example does a steeper slope have a larger effect on the runoff than the soil texture? What are the possibilities in R? Thank you for any feedback, Francesca [[alternativ...

I am getting the following erro rmessage in ordistep. I have a number of similarly structured datasets using ordistep in a loop, and the message only occurs for some of the datasets. I cannot include a reproducible sample - the specific datasets where this is occur ing are fairly large and there are several pcnm's in the rhs of the formula. thanks for any pointers that may allow me to

...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 Content (Percent)') I have also done: as.factor(Day) to change it to 4 levels and when I do is.factor(Day) it returns with FALSE. I'm assuming this is the reason that I get the message that my x and y lengths differ when trying to convert my x axis using: axis(side=1, at=c(0,15,30,45) W...

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...

...ean_5 <- rowMeans(Data1, na.rm =T) Data1 My real dataset has many columns including several Factor type. Is it possible to specify a range of columns within a data frame using rowMeans dims= (say where there is a one Factor column called Species leading I1 to I5, and one Factor column called Moisture following, so 7 columns total) , or do I either need to extract those columns to a new data frame, calculate means, and reattach to the original data frame, or use a different function such as apply? Unfortunately I am fairly new to R and have a difficult time with some terminology and concepts in...

Hello everybody, I Would like to get some help to plot together, Principal Components Analysis (PCA) and clusters. I am handling environmental data from 25 locations spread across 5 different ecosystems.When grouped into 5 clusters, locations from different ecosystems are arranged in the same group. So, I want to plot together PCA and Clusters, in a such way that locations belonging to the same

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...