similar to: Q - scatterplot, plot function & trellis linear regressions???

Hi R-listers, I have developed 47 GLM models with different combinations of interactions from 1 variable to 5 variables. I have manually made each model separately and put them into individual tables (organized by the number of variables) showing the AIC score. I want to compare all of these models. 1) What is the best way to compare various models with unique combinations and different number

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

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

Hi R-listers, I am trying to compare sections of the beach separated from the HTL to the Veg (east to west), separated into indices (-5 to 30m), HTLIndex. Cross parallel (north to south) are major beach sections (Rayos 1, 2, 3, 4 and MNB). I am thinking to do an ANOVA for each independent rectangle of beach (not exactly but will be treated as). The HTL from (0-5m) in Rayo 1 is one rectangle of

Hi R-listers, I am trying to do an ANOVA for the following scatterplot and received the following error: library(car) scatterplot(HSuccess ~ Veg, data = data.to.analyze, xlab = "Vegetation border (m)", ylab = "Hatching success (%)") anova(HSuccess ~ Veg, data=data.to.analyze) Error in UseMethod("anova") : no applicable

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

After looking at 48 glm binomial models I decided to try the quasibinomial with the top model 25 (lowest AIC). To try to account for overdispersion (residual deviance 2679.7/68 d.f.) After doing so the dispersion factor is the same for the quasibinomial and less sectors of the beach were significant by p-value. While the p-values in the binomial were more significant for each section of the

Hi R listers, I am trying to group distances of nests to the vegetation into classes that are define by (0-5m, 6-10m, 11-15m, 16-20m, 21-25m, 26-30m, 31-35m, 36-40m, 41-45m, 46-50m, 51-55m, 56-60m). Each row is a nest and all the distances to the vegetation is in a column. In plyr, I have tried - below script but I think I am going about this the wrong way and am not successful. #Veg index

Hello all, turns out i'm having a bad R week. I am at my wits end with a function that I am trying to write. When I run the lines of code outside of a function, I get the desired output. When I wrap the lines of code into a function it doesn't work as expected. Not sure what is going on here. I suspected that the syntax of the if statement with the for loop was the culprit, but when I

Hello all, I am going slightly mad trying to create a table for running co-correspondence analysis. What I have is seed bank and vegetation data, and my aim is to see if the vegetation found in a site (containing several seed bank samples) can predict the composition of a seed bank sample within that site. So for this I need two tables with matching rows. I have created an empty matrix,

Hi R-users: [R 2.2 on OSX 10.4.3] I have a (sparse) vegetation data frame with 500 rows (sampling units) and 177 columns (plant species) where the data represent % cover. I need to summarize the cover data by returning the names of the most dominant and the second most dominant species per plot. I reduced the data frame to omit cover below 5%; this is what it looks like stacked. I have

Hi R-listers, I am in plyr and using the script below and when I try to pull up the .csv file with head() it does not show up. I am not sure where I am going wrong. I have also checked the file box for turtlehatch.csv and it still does not show up when I head(turtlehatch.csv). Am I missing a command? Please advise, Jean > devel.index <- function(values, weights=c(1, 2, 3, 4, 5, 6)) { +

Hi everybody, I know this might be very off topic and it took me quite a while to up my courage to post this?. But I remember a thread some time ago about how we can find the packages we need to do specific tasks in R if we don?t know before hand which ones actually do it. Now all the packages are listed alphabetically on the web site. Since I am not very advanced in writing my own functions I

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 + Aeventexhumed + HTLIndex:Aeventexhumed,

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,

Hello, We (two mac users) have been attempting to install rgdal from "http://www.kyngchaos.com/software:frameworks", given that it is not available as a binary on the CRAN (binaries) of the "Package Installer". I have also tried to solve this problem by looking on the net for an old question, and though I have found it, the answers do not help very much because our porgraming

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

Hi, I would like to construct a transition matrix from a data frame with annual transitions of marked plants. plant<-c(1:6) class<-c("seed","seed", "seed", "veg", "rep", "rep") fate<-c("dead", "veg","veg","rep", "rep", "veg") trans<-data.frame(plant, class, fate)

Hello I've created a 3d scatterplot, and had no problems labeling the points. However, I've been really struggling to change the color of the points based on a factor (see 'group' below). Is such a thing possible? My data look like this: food group x y z apple fruit 0.216 -0.112 -0.893 orange fruit 0.814 0.097 0.460 broccoli veg -0.239 0.240 -0.425 banana

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