search for: aeventexhumed

Hi R-listers, I am receiving an error - see below. Aeventexhumed is the event in which nesting occured, so it is defined by A, B, C. I thought as a factor was ok, tried to change it to as.character but it still gave me the same error. Is there something I should do about this error or just ignore it? Please advise, Jean >summary(aov(EDI ~ HTLIndex + Aevent...

Hi R-listers, I am trying to make a trellis boxplot with the HSuccess (y-axis) in each Rayos (beach sections) (x-axis), for each Aeventexhumed (A, B, C) - nesting event. I am not able to do so and keep receiving: Error in eval(expr, envir, enclos) : object 'Rayos' not found Please advise, Jean require(plyr) resp <- read.csv("ABC Arribada R File Dec 12 Jean Jang.csv") envir <- read.csv("Responses En...

I am wondering if anyone can tell me what the error I'm receiving means below. I thought it said that Aeventexhumed should be converted to a factor, so I tried to do so and received the following error. Please advise. J --------------------------------------------------------- > data.to.analyze.glm <- glm(cbind(MaxHatch, TotalEggs-MaxHatch) ~ > Aeventexhumed, family=binomial, data=data.to.analyze) Wa...

Hi R-listers, I am wondering if the function aov() in plyr is appropriate for two different types of tests: 1) > summary(aov(EDI ~ VegIndex, data=data.to.analyze)) AND 2) > summary(aov(HSuccess ~ VegIndex + Aeventexhumed + VegIndex:Aeventexhumed, data=data.to.analyze)) the later inclusive of an interaction of the two explanatory variables. This is for obtaining P value and F statistics. Please advise, Jean -- View this message in context: http://r.789695.n4.nabble.com/aov-usage-tp4645380.html Sent from the R...

Hi All, I am wondering if anyone can tell me what the warning message below the model means? J add1(DTA.glm,~ Aeventexhumed + Veg + Berm + HTL + Estuary + Rayos) Single term additions Model: cbind(MaxHatch, TotalEggs - MaxHatch) ~ Aeventexhumed + Veg + Berm + HTL Df Deviance AIC <none> 488.86 4232.9 Estuary 1 454.96 4201.0 Rayos 3 258.80 4008.9 Warning messages: 1: In model.matrix.de...

...Each Rayo is a little bit different in size but for management purposes, it is identified by these sections (leave alone or use), so maybe ANOVA is appropriate? 3) Script suggestions greatly appreciated. Please advise, Jean head() included below ---------------- summary(aov(HSuccess ~ HTLIndex*Aeventexhumed, data=data.to.analyze)) Df Sum Sq Mean Sq F value Pr(>F) HTLIndex 6 0.691 0.11524 3.023 0.00937 ** Aeventexhumed 2 0.106 0.05308 1.393 0.25329 HTLIndex:Aeventexhumed 11 0.376 0.03418 0.897 0.54645 Residuals 98 3.73...

Hi R - listers, I am receiving an error. Does anyone know what this means? J ggplot(subset(foo, Rayos != "Rayos.NA"), aes(x=HTL, y=DevelopIndex, colour=TotalEggs)) +geom_point() +geom_jitter() + facet_grid(Aeventexhumed ~ Rayos) + geom_smooth(method="lm", fill=NA) + ylim(c(0, 7)) Error in +geom_smooth(method = "lm", fill = NA) : invalid argument to unary operator [[alternative HTML version deleted]]

...p the most simplest model ideally. Even though adding another variable would lower the AIC, how do I interpret it is worth it to include another variable in the model? How do I know when to stop? Definitions of Variables: HTL - distance to high tide line (continuous) Veg - distance to vegetation Aeventexhumed - Event of exhumation Sector - number measurements along the beach Rayos - major sections of beach (grouped sectors) TotalEggs - nest egg density Example of how all models were created: Model2.glm <- glm(cbind(Shells, TotalEggs-Shells) ~ Aeventexhumed, data=data.to.analyze, family=binomial) Mo...

...ot;) > abline(EDIHTL) Q-3 I am trying to make a trellis of linear regressions. I want 3 windows displaying 3 events of linear regression (two-way interactions). I have been trying this: require(lattice) > library(lattice) > trellis.par.set(col.whitebg()) > scatterplot(DevelopIndex~Veg|Aeventexhumed, + xlab = "Vegetation border", + ylab = "Embryonic development index") Error in eval(expr, envir, enclos) : object 'DevelopIndex' not found Please advise, many thanks. Jean -- View this message in context: http://r.789695.n4.nabble.com/Q-scatte...

Hi R Listers, I am trying to upload a data file and I received this message. It seems that I am still able to make graphs and Aeventexhumed still works in the analysis. Can I ignore this message or do I need to do something about this? Jean > require(plyr) Loading required package: plyr > turtlehatch <- read.csv(file.choose()) > attach(turtlehatch) The following object(s) are masked _by_ '.GlobalEnv': Aevent...

Hi R-listers, I am using the package plyr. I am just trying to get the hatching success mean of each nesting event and have typed in the following and received the below results: > tapply(HSuccess, Aeventexhumed, mean) A B C 0.2156265 0.1288559 NaN What can I do about NaN? I should be able to get a result for event C because I was able to make a boxplot graph with this information showing all three nesting events (A, B, and C). Please advise, Jean -...

...significant for each section of the beach. -- telling me more about the beach. Is this ok? Can I just look at the binomial glm model 25 and look at its p-values for beach sections and forget about the quasibinomial model 25? J Call: glm(formula = cbind(Shells, TotalEggs - Shells) ~ Sector:Veg:Aeventexhumed, family = quasibinomial, data = data.to.analyze) -- View this message in context: http://r.789695.n4.nabble.com/binomial-vs-quasibinomial-tp4364371p4364371.html Sent from the R help mailing list archive at Nabble.com.

I am trying to create a frequency distribution and I am a bit confused. Here are the commands I have entered: > data <- read.csv(file="40609_sortedfinal.csv",head=TRUE,sep=",") > NumberOfActionsByStatus = data$STATUS > NumberOfActionsByUser = data$ETS_LOGIN > NumberOfBidOffer = data$BID_OFFER > NumberOfActionsByUser.freq = table(NumberOfActionsByUser) >