similar to: Time-series analysis with treatment effects - statistical approach

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

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

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

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)

Hello all, I'm attempting to conduct spatial analysis of trees within a plot. I want to see if the trees are spatially correlated to soil characteristics, say pH, or moisture content. I think one way to do it is with mpl, however, my soil characteristics were not taken at exactly the same locations as my trees and further, the vectors aren't the same length. I'm getting the

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

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

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

I am trying to fit a nonlinear regression to infiltration data in order to determine saturated hydraulic conductivity and matric pressure. The original equation can be found in Bagarello et al. 2004 SSSAJ (green-ampt equation for falling head including gravity). I am also VERY new to R and to nonlinear regressions. I have searched the posts, but am still unable to determine why my data come up

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,

I am performing a binary classification using a classification tree. Ironically, the data themselves are 2483 tree (real biological ones) locations as described by a suite of environmental variables (slope, soil moisture, radiation load, etc). I want to separate them from an equal number of random points. Doing eda on the data shows that there is substantial difference between the tree and random

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 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 : Factor

I'm using Samba 3.0.3 on Fedora Core 2 with OpenLDAP 2.1.29 for a backend. I'm getting to typical "The user name could not be found." error upon trying to join a Windows box. I've gone through every digest on lists.samba.org and other sites and nothing has worked yet. Any suggestions: Here's what I've done so far: 1. Installed everything via RPMS: [root@smbtest

Greetings, After pouring over the documentation for the 'colorspace' package, I have not been able to figure out how the plot() method converts colorspace coordinates to RGB values for display on the screen. I am convert between colorspaces with the various as() methods... but cannot seem to find a way to extract RGB (i.e. for displaying on a computer screen) triplets from color space

After, oh, six months of attempts here and there to read everyone's experiences with Samba/LDAP and inability for a windows 2000/XP machine to join the domain, I finally discovered what was not working properly. In my smb.conf I put: add machine script = /usr/local/sbin/smbldap-useradd -w "%u" As instructed by many How-to's and Idealx. However, I thought to myself, %m

Hi, forgive me if someone has already posted about this but I have had a look and cannot find the answer, also I am very new to R and been getting the grips with this. I have been trying to use Adonis to find out if there are significant difference between groups on data that I have analyses with NMDS, and have been struggling with getting this to work and understanding what is going on. I am

Dear all; First of all, this is probably a more conceptual question than a R-related one, but still want to give it a try. When working with the interpolation function "interp" from the package akima and the triangulation function "tri.mesh" from package tripack I've got NA's for the interpolation and "error" for the triangulation. The data is arranged in a