similar to: zoo and cex

I have insect data from twelve sites and like most environmental data it is non-normal mostly. I would like to preform an anova and a means seperation like tukey's HSD in a nonparametric sense (on some sort of central tendency measure - median?). I am searching around at this time on the internet. Any suggestions, books, etc. would be greatly appreciated. -- Stephen Sefick Research

R version 2.8.0 (2008-10-20) i386-pc-mingw32 locale: LC_COLLATE=English_United States.1252;LC_CTYPE=English_United States.1252;LC_MONETARY=English_United States.1252;LC_NUMERIC=C;LC_TIME=English_United States.1252 attached base packages: [1] stats graphics grDevices utils datasets methods base other attached packages: [1] StreamMetabolism_0.01 chron_2.3-24 zoo_1.5-4 loaded

library(zoo) d<-(structure(c(1.39981554315924, 0.89196314359498, 0.407816250252697, 0.823496839063978, 1.14429021220358, 1.23971035967413, 0.960868900583432, 0.927685306209829, 1.22072345292821, 0.249842897450642, 1.00879641624694, 0.925372139878243, 0.317259909172362, 0.382677149697482), index = structure(c(11808, 11869, 11961, 11992, 12084, 12173, 12265, 12418, 12600, 12631, 12753, 12996,

I have a time series from 1933-2005 of precipitation at Fayetteville NC. I get the following error messages when I plot the zoo series. Any help would be appreciated. If you need the data I can dput it or send the csv. I didn't include it here because I didn't want to clog up anybodies email account. I know that this is not reproducible, and I will send along the file if needed.

I am was going to look at the as.yearmon function in the zoo package and write a as.day function to aggregate a time series of 96 observations per day into the mean for each day, but I don't know how to look at the code so that I can convert it into something I can use. On top of that I believe that it is probably an S3 method and I haven't quite gotten that far in my programming

I would like to merge zoo objects that are stored in a list into one big zoo object with one index for all of the observations. I have created the list (74 dataframes) with the code below, and have tried the do.call(merge, foo) in the call and the output is not what I expected. Any help would be greatly appreciated. Stephen Sefick ###################################################level logger

library(zoo) library(chron) t1 <- chron("1/1/2006", "00:00:00") t2 <- chron("1/31/2006", "23:45:00") deltat <- times("00:15:00") tt <- seq(t1, t2, by = times("00:15:00")) d <- sample(33:700, 2976, replace=TRUE) sin.zoo <- zoo(d,tt) #there are ninety six reading in a day d.max <- rollapply(sin.zoo, width=96, FUN=max)

# chron library(chron) fmt.chron <- function(x) { chron(sub(" .*", "", x), gsub(".* (.*)", "\\1:00", x)) } z1 <- read.zoo(SC3.csv, sep = ",", header = TRUE, FUN = fmt.chron) z2 <- read.zoo(SC2.csv, sep = ",", header = TRUE, FUN = fmt.chron) z3<-c(z1, z2) write.table(z3, sep="," , "SC.csv") How do you

z1 = read.zoo("chemmgL.csv", sep=",", header=T, format="%m/%d/%y") I would like to the entire date field in the plot - 1/1/07 all that I get now is 2007 I do not include data because it is to large of a data set for an email -- Let's not spend our time and resources thinking about things that are so little or so large that all they really do for us is puff us

#this is a subset of a larger data frame and I am okay with subsetting it as there are redundant time stamps, but I would like to create a zoo object out of this and I am having a hard #time figuring out how to do this the date structure is year and then month x <- structure(list(Yearmonth = structure(c(12L, 24L, 1L, 13L, 14L, 3L, 15L, 4L, 16L, 5L, 17L, 6L, 18L, 7L, 19L, 8L, 20L, 9L, 21L,

#Is there a way to break the below zoo object into non-NA data frames algorithmically #this is a small example of a much larger problem. #It is really no even necessary to have the continuous chunks #end up as zoo objects but it is important to have them end #up with the index column. #thanks for all of your help in advance, and #if you need anything else please let me know library(zoo) ind.

I have a matrix with data that runs from 1/1/06 00:01:00-1/31/08 23:46:00. I have read in the data with this fmt.chron <- function(x) { chron(sub(" .*", "", x), gsub(".* (.*)", "\\1:00", x)) } x <- read.zoo(file.choose(), sep=",", header=T, FUN=fmt.chron) plotted with this plot(x[,(seq(3, by=9, length.out=12))],

the comma seperated file is 37Mb, and I get the below message: it is zoo object read in this way: # chron > library(chron) > fmt.chron <- function(x) { + chron(sub(" .*", "", x), gsub(".* (.*)", "\\1:00", x)) + } > z1 <- read.zoo("all.csv", sep = ",", header = TRUE, FUN = fmt.chron) and then the plot is done with:

These are date and times in the format YYYYMMDDhhmmss. I would like to take this column and make a chron object form them. I have tried a couple of the split family of functions but they need character input here is the data: date.time <- c(19851001001500, 19851001003000, 19851001004500, 19851001010000, 19851001011500, 19851001013000, 19851001014500, 19851001020000, 19851001021500,

#I am putting a test together for an introductory biology class and I would like to put different cross hatching inside of each bar for the bar plot below color <- c("Brightly Colored", "Dull", "Neither") lizards <- c(277, 70, 3) liz.col <- data.frame(color, lizards) qplot(color, lizards, data=liz.col, geom="bar", ylab="Observed Matings",

NOTE: I will provide data if necessary, but I didn't want clutter everyones mailbox All: I have a time series with level and temperature data for 11 sites for each of three bases. I will have to do this more than once is what I am saying here. OK, The time series are zoo objects with index values in chron format. The problem is that the date and times should be at even 15 min intervals,

plot(1,2, ylab= paste("insects", expression(m^2), sep=" ")) I get insects m^2 I would like m to the 2 what is the problem? -- Let's not spend our time and resources thinking about things that are so little or so large that all they really do for us is puff us up and make us feel like gods. We are mammals, and have not exhausted the annoying little problems of being

Is there any easy way to pull out the row indexes for a logical matching statment? #################example code######################################### foo <- data.frame(name=c(rep("A", 25), rep("B", 25), rep("C", 25), rep("A", 25)), stuff=rnorm(100), and=rnorm(100), things=rnorm(100)) #this is what I want but I would like the row indexes

I would like to rotate the axis labels 45 deg. -- Stephen Sefick Let's not spend our time and resources thinking about things that are so little or so large that all they really do for us is puff us up and make us feel like gods. We are mammals, and have not exhausted the annoying little problems of being mammals. -K. Mullis

I have tried to compile this from source. I don't know what Endianess is, but it is probably not debian power pc. Am I would of luck with this package? Stephen Sefick * Installing *source* package ?ifultools? ... ** libs gcc -std=gnu99 -I/usr/local/lib/R/include -I"../inst/include/" -D"MUTIL_STATIC" -D"DEF_TF" -D"INTERRUPT_ENABLE"