search for: jeanwaijang

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 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 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 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 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, In plyr I am using coding I have used before but an error seems to be apparent. Has the function "join" been replaced by another word? I tried this before and it didn't work and then worked for me. But now it is not working again. Please advise, Jean > data.to.analyze <- join(responses, envir, by="QuadratEvent") Error: could not find function

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

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

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

Hi everyone, I just did a GLM binomial regression. I am wondering how do I know when my responses are fixed? I have found that an interaction between two factors is the most significant predictor amongst other factors by comparing the AIC and add1 with each factor. Does this mean that the interaction between these two factors are my fixed responses? I would like to do a quasibinomial

I've originally made 48 GLM binomial models and compare the AIC values. But dispersion was very large: Example: Residual deviance: 8811.6 on 118 degrees of freedom I was suggested to do a quasibinomial afterwards but found that it did not help the dispersion factor of models and received a warning: Residual deviance: 3005.7 on 67 degrees of freedom AIC: NA Number of Fisher Scoring

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

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