search for: vegtyp

...centrations at sample point locations. Ultimately I'd like to use the probabilities for each response category in a simulation model such that these probabilities are used to define a realized ecological niche. When I try the following it works for a single response. Typha.glm <- glm(fwc$VegType == "Cattail" ~ fwc$TP + fwc$TC + fwc$TN + fwc$BD + fwc$LOI + fwc$Total_Mg, family = poisson) But if I try this without specifying a specific VegType it fails. plants.glm <- glm(fwc$VegType ~ fwc$TP + fwc$TC + fwc$TN + fwc$BD + fwc$LOI + fwc$Total_Mg, family = poisson) Error in y...

...aving issues with a the mvabund package (http://cran.r-project.org/web/packages/mvabund/index.html), and would be great to get some help! Here is the code (and files are attached): library(mvabund) ##visualizing data florabund <- read.csv("CPL_floristics_abund_v1d.csv", header=TRUE) vegtype <- read.csv("CPL_vegtype_all_v1b.csv", header=TRUE) transect <- read.csv("CPL_transect_all_v1b.csv", header=TRUE) cplflorabund <-mvabund(florabund, check.rows=FALSE, check.names=TRUE) vegtype=as.factor(vegtype$vegtype) plot(cplflorabund~vegtype) #there is a column cal...

...6-pc-linux-gnu (Suse 8.2) Submission from: (NULL) (134.84.86.22) Two bugs (perhaps related, perhaps independent) revealed by the same Poisson regression with offset mydata <- read.table(url("http://www.stat.umn.edu/geyer/5931/mle/seeds.txt")) out.fubar <- glm(seedlings ~ burn01 + vegtype * burn02 + offset(log(totalseeds)), data = mydata, family = poisson) summary(out.fubar) out.barfu <- glm(seedlings ~ burn01 + vegtype * burn02, offset = log(totalseeds), data = mydata, family = poisson) summary(out.barfu) out.ok <- glm(seedlings ~ vegtype * burn02 + burn01, offse...

I am sure this question has come up, but searching hasn't given me any results. So I need to enter this line: mx1 <- randomForest(X1 ~ elevation + slope + vegtype, data = moths.train) But the problem is that X1 is currently hard coded. I would instead like to be able to put in the value of X1 through a list. For example: list <- list("X1", "X2", "X3") #Then, instead of X1, I want to put list[1] mx1 <- randomForest(list...

...omparison to be based on probabilities. An example of the graph I'm seeking can be seen here: http://www.statmethods.net/graphs/density.html (scroll to the bottom and examine the Comparing Groups via Kernel Density). A graph (attached) produced using my data: > attach(SoilVegHydro) > vegtype.f <- factor(Physiogomy) > sm.density.compare(meanAnnualDepthAve, vegtype.f, xlab="Mean Annual Depth") # create color fill patterns >colfill <- c(2:(2+length(levels(vegtype.f)))) > legend(locator(1), levels(vegtype.f), fill=colfill) (See attached file: MeanAnnualD...

...summary(mlogit(mode ~ pr + ca -1, data = Fish)) Call: Error in deparse(object at call) : trying to get slot "call" from an object (class "mlogit") that is not an S4 object Using my data I've done the following: xx <- fwc # copies and protects original data.frame xx$VegType <- as.factor(xx$VegType) xx.dat <- mlogit.data(xx, varying = "NULL", choice = "VegType", shape = "wide") xx.mlogit <- mlogit(VegType ~ TP, data = xx.dat,, reflevel = "1") Error in relevel.factor(attr(mf, "index")[["alt"]], ref...

Dear all, I would like to use the values in vegaggr.BLMCMR02$colony str(vegaggr.BLMCMR02) `data.frame': 1678 obs. of 3 variables: $ vegtype : Factor w/ 27 levels "2010","2020",..: 3 4 5 19 4 5 19 5 $ colony : Factor w/ 406 levels "0","1","10","100",..: 1 1 1 1 2 2 2 $ Totvegproparea: num 0.00055 0.03956 0.95705 0.00284 0.05215 ... as the seq in a for loop (f...

Dear friends, when i use glm() to fit my data, i use glm(formula = snail ~ vegtype + mhveg + humidity + elevation + soiltem, *family = Gamma(link = inverse),* data =a,)) It shows: error in eval(expr, envir, enclos) : *gamma distribution don't allow negative value*. But i use result<-glm(formula = snail ~ vegtype + mhveg + humidity + elevation + soiltem, family = poisson...