Dear R Users, I have searched through the archives but I am still struggling to find a way to process the below dataset. I have a dataset that has stratum and plot identifier. Within each plot there is variable (Top) stating the number of measurments that should be used to to calculate the mean to the largest "top" elements within one of the vectors (X). I would like to process this summary statistic by groups. At this stage I have been trying to use the "summarize" function within the Hmisc library but I am getting the following error "Error in eval(expr, envir, enclos) : numeric 'envir' arg not of length one In addition: Warning message: no finite arguments to max; returning -Inf".
Any suggetsions on how I can fix this would be greatly appreciated. Kind regards Andrew test <- read.table("test.csv", header=TRUE, sep=",") #function to calculate mean of "top" elements within a plot > g<-function(y) { + top_no <-max(y$top) + weight <- with(y,as.numeric(x>=x[order(x,decreasing=TRUE)[top_no]])) + wtd.mean(y$x,weight) + } > g(test) [1] 172.6667 #call to summarize function - use function g and summarise by stratum plot > test.2 <- with(test,summarize(test$x,llist(test$Stratum,test$plot),g)) Error in eval(expr, envir, enclos) : numeric 'envir' arg not of length one In addition: Warning message: no finite arguments to max; returning -Inf >traceback() 9: eval(substitute(expr), data, enclos = parent.frame()) 8: with.default(y, as.numeric(x >= x[order(x, decreasing = TRUE)[top_no]])) 7: with(y, as.numeric(x >= x[order(x, decreasing = TRUE)[top_no]])) 6: FUN(X, ...) 5: summarize(test$x, llist(test$Stratum, test$plot), g) 4: eval(expr, envir, enclos) 3: eval(substitute(expr), data, enclos = parent.frame()) 2: with.default(test, summarize(test$x, llist(test$Stratum, test$plot), g)) 1: with(test, summarize(test$x, llist(test$Stratum, test$plot), g)) The version im running on is $platform [1] "i386-pc-mingw32" $arch [1] "i386" $os [1] "mingw32" $system [1] "i386, mingw32" $status [1] "" $major [1] "2" $minor [1] "2.0" $year [1] "2005" $month [1] "10" $day [1] "06" $"svn rev" [1] "35749" $language [1] "R" > Stratum plot id top x 1 1 1 2 12 1 1 2 2 41 1 1 3 2 12 1 1 4 2 43 1 1 5 2 12 1 1 6 2 14 1 1 7 2 43 1 1 8 2 12 1 2 1 4 42 1 2 2 4 12 1 2 3 4 432 1 2 4 4 12 1 2 5 4 12 1 2 6 4 14 1 2 7 4 41 1 2 8 4 1 2 1 1 2 12 2 1 2 2 41 2 1 3 2 12 2 1 4 2 43 2 1 5 2 12 2 1 6 2 14 2 1 7 2 43 2 1 8 2 12 2 2 1 3 42 2 2 2 3 12 2 2 3 3 432 2 2 4 3 12 2 2 5 3 12 2 2 6 3 14 2 2 7 3 41 2 2 8 3 1 [[alternative HTML version deleted]] ______________________________________________ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html