search for: ungul

Hi, I'm working on the codes below however every time I run them when they get to OpenBUGS I keep getting the error message: array index is greater than array upper bound for hab. Any help would be greatly appreciated, Suzie Codes: ungulate <- read.csv(file.choose ()) #ungulate ungulate <- as.matrix(ungulate);colnames(ungulate)<-NULL;rownames(ungulate)<-NULL habitat <- read.csv(file.choose ()) #ungulate habitat habitat <- habitat[,3] site.data <- read.csv(file.choose ())...

I have been working on the following code but keep getting an err message. My current thinking is that the problem is on the indexing but do not know how to fix it. Any help please? ungulate <- read.csv("Ungulate.csv",row.names=1) ungulate <- as.matrix(ungulate);colnames(ungulate)<-NULL;rownames(ungulate)<-NULL habitat <- read.csv("Ungulate_vegetation.csv") habitat <- habitat[,3] site.data <- read.csv("Ungulate_site.csv") site <-...

...named percod, with 11 levels) of a bird presence/absence (pres.f) on 14 artificial ponds (charca.f). The ponds are integrated in 4 greater spatial units (zepa). I have 11 explanatory variables, some of them are continuous (prof1,prof2,prof4,capac,perim_veget,div_veg,div_herp) and some are factors (ungul.f,ganado.f,turb.f,peces.f). Not all artificial ponds have been monitored in all seasons, so I have some missing values in the dependent variable. I have thought that the correct model should be, ignoring the spatial pseudorreplication: m1<-lmer(pres.f~ganado.f+ungul.f+turb.f+prof1+prof2+pro...

...urrently working on the below codes however whenever I run it in openbugs it gives an error message saying: unknown type of logical function error pos 76. Any help would be appreciated. ## bugs code library(R2OpenBUGS) sink("C:/Users/CCF/Documents/Suzie Work/PTY Project/Waterhole Correction/ungulate.txt") cat(" model{ # hyperparameters # habitat effects for each functional group g <- length(table(G)) for(i in 1:g){ # number of functional group for(j in 1:3){ # number of habitat type mu.h[i,j] ~ dnorm(0,0.0001) I(-5,5) sigma.h[i,j] ~ dunif(0,5) tau.h[i,j] <- 1/(sigma.h[i,...