similar to: Error - function join()

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

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

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

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 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, 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 just tried updating my R and now I can't even open it and it is prompting me to relaunch then relaunch just reappears. And it will not open R. I am afraid I may have lost my scripts. What should I do? I am running a MacBook OS X Version 10.5.8 1) Restore the entire system to an earlier date? 2) Drag R from applications and reload R program (updated version)? Will my

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

Dear R-listers, I am a newbie with R and I am struggling with something I consider very basic. I wish to produce a table (to import in a latex file) of summary statistics, but for as much as I've been looking around and trying various alternatives (plyr, reporttools, pastecs and Hmisc) I haven't found what I am looking for. Probably I am doing something wrong, but I can't figure out

Hi R-listers, I am having trouble plotting the values for the R2 line, lowess smooth, smothered conditional spread, outlier identification. I have tried the function text and legend but was unsuccessful. Please see below. Your help would be appreciated. Thanks in advance. Jean library(car) > SPHSHTL <- scatterplot(HSuccess ~ HTL, + data = data.to.analyze, +

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

R-listers, I'm working on a project where I need to get the lowercase of the county variable in a dataset alike to the example below (only difference is that the full dataset has all the states and counties in the southeast United States). I keep getting this strange error with the lowerize function (which didn't occur the first few times I use the code below), and oddly enough the error

I retrieve for a few hundred times a group of time series (10-15 ts with 10000 values each), on every group I do some calculation, graphs etc. I wonder if there is a faster method than what presented below to get an appropriate timeseries object. Making a query with RODBC for every group I get a data frame like this: > X ID DATE VALUE 14 3 2000-01-01 00:00:03 0.5726334

(moved from R-help) I tried this on Ubuntu with R-2.14.1 built from source, and I do not get the segfault problem. (I don't at the moment have a debian binary R, or I would confirm whether I get the segfault problem.) My sessioninfo() is reporting additional information about namespace imports: > library(ggplot2) Loading required package: reshape Loading required package: plyr

I want to learn how to use the reshape package. The reshape package is not included in the Ubuntu repositories, so I attempted to install reshape with: install.packages("reshape") This is what I got for output: Warning in install.packages("reshape") : argument 'lib' is missing: using '/home/andy/R/i486-pc-linux-gnu-library/2.7'