search for: grazing

Hi, Try either: res1 <- apply(mydata[,1:2],2,mean) ?res2 <- colMeans(mydata[,1:2]) ?identical(res1,res2) #[1] TRUE # Also if you need to find means for each group ("Ungrazed vs. "Grazed") by(mydata[,-3],mydata[,3],colMeans) #or if column names are "V1", "V2", "V3" aggregate(.~V3,mydata,mean) #or library(plyr)

Hi everyone! I have a simple model that i would like to plot with 95% CIs. It is like follows: m1<-lmer(Richness~Grazing+I(Grazing^2)+(1|Plot),family=poisson) By using the effects package I get two plots, one for the linear term and one for the squared term. Q1: Can I get all in one? I.e. with one line for the whole model? Q2: Can I also visualize the random effects? I would be very happy for your answers! -- Mvh...

Dear R- Community, to learn i reanalysed some data provided and analysed by Zuur et. al. in their book "Mixed effect models and Extensions in Ecology with R". When i run the last command i get a warning message i dont understand. Loyn<- read.table(file = "loyn.txt",header = TRUE) Loyn$L.AREA<- log10(Loyn$AREA) fGRAZE <-factor(Loyn$GRAZE) M0<- lm(ABUND~ L.AREA

...5 3 2 6 HE close no 1 F 0 8 1 6 HE close no 2 F 0 5 2 6 AHE close no 1 F 0 8 2 6 AHE close no 2 F 1 11 2 The response variables are "Litter", "Bareground", "Logs" (measured by visually estimating % cover in a plot in a grazing exclusion experiment) The three factors are fixed and are "Site" - 6 levels; "Treatment" - 3 levels; and "Time" - 6 levels Additional factors are "Graze" and "Dam" but are constructs on top of the sites, that is Sites 1,3,5 have dams adjacent whil...

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

I have a two part question one about statistical theory and the other about implementations in R. Thank you for all help in advance. (1) Am I correct in understanding that Heteroscedasticity is a problem for Generalized Additive Models as it is for standard linear models? I am asking particularly about the GAMs as implemented in the mgcv package. Based upon my online search it seems that some

...0.000000 Phytoplankton Zooplankton Grazing Rate Constant for Group 1 (1/day)...

...oring (over a period of a few years) a number of paired plots that were flooded and / or burnt. The plots are located in two topographical settings, some were burnt, some were flooded, some were burnt & flooded and some were not affected at all. At each location one plot was fenced to exclude grazing while another was left unfenced to include the effect of grazing. Effectively, I guess, we have a split plot design in a longitudinal study. The plots were established opportunistically, so there were some practical aspects that might not be ideal... but then when is this ever the case... I now...

Hello and thanks in advance I am running a glm in R the code is as follows with residual diagnostic code below model4<-glm(Biomass~(Treatment+Time+Site)^2, data=bobB, family=quasi(link="log", variance="mu")) par(mfrow=c(2,2)) plot(model2) to test the effect of grazing exclusion of feral horses for a Phd with following factors: Treatment - 3 levels which are grazed controls (C), ungrazed horse exclosures (H)with 3 wire strand fence to allow other herbivores but not horses and ungrazed full netting exclosures to exclude all herbivores (A) Site - 6 levels for 6...

Hallo List, I'm trying to implemement a restricted permutation scheme in permutest(). More precisely I have dependence in my data that should be allowed for in the permutation - I simulated the problem in the example of the vegan documentation p.24: library(vegan) data(varespec) ## Bray-Curtis distances between samples dis <- vegdist(varespec) ## First 16 sites grazed, remaining 8 sites

...es Hello and thanks in advance I am running a glm in R the code is as follows with residual diagnostic code below model4<-glm(Biomass~(Treatment+Time+Site)^2, data=bobB, family=quasi(link="log", variance="mu")) par(mfrow=c(2,2)) plot(model2) to test the effect of grazing exclusion of feral horses for a Phd with following factors: Treatment - 3 levels which are grazed controls (C), ungrazed horse exclosures (H)with 3 wire strand fence to allow other herbivores but not horses and ungrazed full netting exclosures to exclude all herbivores (A) Site - 6 levels for 6...

Hello, I would appreciate if somebody could help me clear my mind about the below issues. I have a factorial experiment to study the effects of Grazing and Fire on Forest biomass production. The experimental unit (to which the treatment combinations are applied) are PLOTs. The measures were made repeatedly for 13 years. I am planning to use the linear mixed effect model function lme in R for this. I know that in software like SPSS, using Repeated...

..., add herd effect as covariate, etc.) will be more then welcome, but seem to always violate some assumption (equal slopes of covariate, independence from predictors), or lose much of the analytical power. 1. Shiyomi, M., Tsuiki, M., 1999. Model for the spatial pattern formed by a small herd in grazing cattle. Ecological Modelling 119, 231-238. 2. Shiyomi, M., 2004. How are distances between individuals of grazing cows explained by a statistical model? Ecological Modelling 172, 87-94 3. Blakely, T.A., Woodward, A.J., 2000. Ecological effects in multi-level studies. J. Epidemiol. Comm...

A rather strange situation here and I've not found the source of the problem. The point is to print a ternary plot matrix of compositional data with ellipses enclosing 95% of the variance in each plot. The ellipses display on the monitor, dev = x11cairo (see attached winters-x11cairo.pdf), but not when sent directly to a file, dev = pdf (see attached winters-pdf.pdf). Here's

Heya, i got this problem with stubborn win98 machine refusing to show up on the browse list, my samba is the local/domain master browser and login server for the whole local network, using user security. the workstation has all the same protocols and settings like the other 6, only difference is the ip and the netbios name, the login used on that machine matches existing username on my linux

Dear List, I am writing a custom panel function and xyplot method to plot the results of a procrustes analysis from the vegan package. I am having trouble getting the call to panel.arrows to work as I wish when conditioning. The attached file contains the function definitions for the xyplot method and the custom panel and prepanel functions I am using. This example, using data and functions from

hi,all I want set up a puppet HA structure. but if the puppetd can set more than one server address ,the work will be simple :D if one puppetmaster have error, the client can auto try the next. I think I can use the heartbeat or dns to do this same thing,but not so good. so I want know the puppetd can direct set in the configure file ? if not,will be add this option? -- Huang

All, I'm having a problem with an Innomedia 6328-2Re (old Sunrocket Gizmo). It keeps losing registration after a period of time ranging from a few minutes to a few hours. It seems that right before it loses registration, it fails to send a second register (after the 401 unauthorized). Here's a transcript from wireshark (at the end). The last message is all that's received and

All, I've migrated a share from an older samba server (3.0.33) running on a Sun Solaris box to a newer samba server (3.6.23) running on a RedHat Linux 6 box. For the most part I've tried to keep the configuration the same (simply to make the migration easier since I'm not an expert on Samba config options). I got the authentication working (security=domain) but now I'm