similar to: Error unary operator

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

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

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

Dear all I would like to superpose some smoothing line through boxplot in ggplot > dput(ad) structure(list(konc.f = structure(c(1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 6L, 6L, 6L, 7L, 7L, 7L, 7L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 6L, 6L, 6L, 7L, 7L, 7L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 4L, 4L, 4L,

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

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 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, When I run a regression, I got the following response.  your help is appreciated.  Error in constants (lhs, cnames_symb) :   The hypothesis "Cost_PKL-QK.X_(Intercept)=0" is not well formed: contains bad coefficient/variable names.  In addition: Warning message: In constants(lhs, cnames_symb) : NAs introduced by coercion Dereje [[alternative HTML version deleted]]

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 all: I am having trouble dropping levels, got a few hints online?without success. Please consider the dataset below: ?I was under the inpression that subset(......drop=TRUE) would work but it doesn't library(ggplot2) ??? library(hmisc) x <- structure(list(first = c(38.2086, 43.1768, 43.146, 41.8044, 42.4232, 46.3646, 38.0813, 40.0745, 40.4889, 38.6246, 40.2826, 41.6056, 34.5353,

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 there, WinXP Prof-Client on a samba 3.0.2a-Domain The domain-login ist ok in general. But when a user logs in, he gets an editor-window with a desktop.ini. content: [.ShellClassInfo] LocalizedResourceName=@%SystemRoot%\System32\shell32.dll ,-21787 ) i think, this ist desktop.ini in //%L/profiles/%u/startmenu/autostart. Do you have an idea how to ban this editor window from desktop after

Dear list,   Using ggplot2 I could produce both boxplot and points in the same plot but instead of points I would like to label the different subjects with different colors or their idnumbers. Is there away to do it? Also how can I put three plots on the same graph with ggplot2? mfrow=c(3,1) did not do the job.   > dat    group time   id  freq 1      1   00 0018  5.21 2      1   00 3026  3.13

Hi When I run my script using ggplot and geom_smooth I get messages that I would like to suppress: p <- ggplot(dataSubset) p <- p + aes(x = as.Date(factor(key),format="%Y%m%d")) + geom_line() p <- p + geom_smooth(span=0.2,se=FALSE,size=0.7) The messages look like this: geom_smooth: method="auto" and size of largest group is <1000, so using loess. Use 'method