similar to: masked by GlobalEnv ???

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

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

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,

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

Hi All, I am a beginner to R and a prof helped me with some script. I am having trouble understanding the below line. Is it finding the file turtlehatch.csv? I do not have my working directory set to this file. If so I tried ("file.choose") and it did not work. R did nothing in response, no error message with this entry and the below entry. Can you please help me interpret? :

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 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 R-listers, Q1) What is the difference between the scatterplot and plot function? Q2) I am able to make a graph with the scatterplot function: scatterplot(DevelopIndex ~ Veg, + data = Turtle, + xlab = "Vegetation border (m)", + ylab = "Embryonic development index") And have been successful. But I do not know if the lines are for:

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

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

I was able to make a scatterplot but ... 1) what does the "86" mean? The "86" shows up on the graph as well. > scatterplot (Shells/TotalEggs ~ Sector, data = data.to.analyze) [1] "86" 2) Also how do you change the Y axis title? I don't want it to read Shells/TotalEggs, instead I would like it to read Average Hatching Rate (%). 3) What does this error mean?

Hi All! I'm trying to convert serial numbers in Excel to dates in R. For each single "thedate" entry, I get a correct answer. But if I try using the for loop, I get bizarre numbers in "mynewdata". thedate<-as.matrix(40548:40759,ncol=1) exdate<-function(){ mynewdate<-NULL for(i in 1:nrow(thedate)){

Hi all, I create some objects in local(), and want to pass them to GlobalEnv. How can I do this? Thanks! xiagao1982 2011-05-02 [[alternative HTML version deleted]]

Full_Name: Mark Bravington Version: 1.5.0 OS: Windows 2000 Submission from: (NULL) (140.79.2.3) In R1.3.1, I used the following code inside a function to set a "path" attribute to .GlobalEnv: env_ pos.to.env( 1) attr( env, 'path')_ .Path # .Path is a named char vector of length 1 .Internal( lib.fixup( env, .GlobalEnv)) # adds the attribute And this works fine. But in

Hi First, sorry if I get the terminology wrong, I am still quite new to the concept of using environments and workspaces. Say I have a statement in a package SIM like sim <- TYPE where the variable TYPE is initialized in the package to e.g. "exponential" (SIM::TYPE == "exponential"). Now, I want to give the user the option of specifying the variable TYPE, but to the

Hello, Suppose I have an expression, E, which accesses some variables present in an environment V. I do this via eval(E,envir=V) however all assignments end up in V. I would like the results of assignments in E to end up in the .GlobalEnv ? Or at least the calling environment. Is there a quick way to this instead of iterating over all objects E and assigning into .GlobalEnv? Thank you Saptarshi