similar to: Interpreting script - write.csv

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

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

Hello, Is there a way to put multiple data.frames or objects into a single .csv ? I know I can use sink() to do this for a text file. Is there way to do it for a .csv? e.g. for sink() sink("c5.k.s.p.92.LN.sl.vs.mean.txt") symnum(c5.k.92.LN.sl.vs.mean) symnum(c5.s.92.LN.sl.vs.mean) symnum(c5.p.92.LN.sl.vs.mean) sink() Could I put these all into a single .csv, in a similar format to

Please see the following example. I can not write '"' to a csv file successfully. Could somebody let me if it is possible to write '"' to a csv file with the default setting of write.csv? my_home$ Rscript main_quote.R > x=rbind( + "\"A\"" + , "\"B\"" + ) > x [,1] [1,] "\"A\"" [2,]

Hello all, Have a little annoying issue with wreite.csv. It always writes out observation number. Is there a way to avoid this. Thanks, DL -- View this message in context: http://r.789695.n4.nabble.com/Have-write-csv-write-csv-without-observation-number-tp4030672p4030672.html Sent from the R help mailing list archive at Nabble.com.

hi i want to ask a simple question: exists a samba VFS module for interpreting a file on access? (mostly the same as a dynamic file on an webserver) if not, i want to suggest it ;-) a simple example: a share with an single MS-DOS .bat file on it; this .bat file is (for example) a bash script file; if the file is accessed via SMB than smbd executes this script and handles the stdout output 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 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, +

Hello, I am pretty new to R. I am working on neural network classifiers and I am feeding the nnet input from different regions of interest (fMRI data). The script that I am using is this: library (MASS) heap_lda <- data.frame(as.matrix(t(read.table(file="R_10_5runs_matrix9.txt")))*100000,syll = c(rep("heap",3),rep("hoop",3),rep("hop",3))) library(nnet)

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