search for: legumes

...What I get is the following dataframe (simplified, there are more categories, not real data): trial cover_of dateofsurvey cover 1 ZU-316 Cover dead material 2004-09-16 0 2 ZU-316 Cover grasses 2004-09-16 16 3 ZU-316 Cover herbs 2004-09-16 14 4 ZU-316 Cover legumes 2004-09-16 10 5 ZU-316 Cover open soil 2004-09-16 30 6 ZU-316 Cover dead material 2005-09-18 5 7 ZU-316 Cover grasses 2005-09-18 26 8 ZU-316 Cover herbs 2005-09-18 14 9 ZU-316 Cover legumes 2005-09-18 15 10 ZU-316 Cover open soil 2005-09-1...

Dear List, I have got error message when I run the R-code. Can anyone has a suggestion? v.code <- df.bm7[,c(10:31)]; v.code[1:3,] names(v.code) CM = v.code # variable binomial code sim.sp <- function(data,CM,n,N) { C <- matrix(rep(NA,N),ncol=1) for(i in 1:N) { j <- n xx <- which(colSums(CM[j,])==1) V <- names(xx) V <- paste(V,

Dear list users, I have trouble with covariance analysis. I measured nitrate concentrations in the soil (NO3) and the percentage of legumes (LEG, continuous), affected by 2 different CO2 concentrations (CO2, discrete). I suspect that CO2 has an effect on LEG and NO3, but also that LEG has an effect on NO3, so this is the formula I wrote to test this: NO3 ~ CO2 + LEG + CO2:LEG Will LEG be considered continuous if I use aov(NO3 ~ CO2...

Dear List, I'm using gam in a multiple imputation framework -- specifying the knot locations, and saving the results of multiple models, each of which is fit with slightly different data (because some of it is predicted when missing). In MI, coefficients from multiple models are averaged, as are variance-covariance matrices. VCV's get an additional correction to account for how

Dear List, I'm getting an error in mgcv, and I can't figure out where it comes from. The setup is the following: I've got a fitted GAM object called "MI", and a vector of "prediction data" (with default values for predictors). I feed this into predict.gam(object, newdata = whatever) via the following function: makepred = function(varstochange,val){ for