Hello, I'd like to automate this script a bit more and cycle several parameters(both the species and the metric). For example where AnnualDepth occurs, I need to process about 12 metrics so instead of writing this entire script 12 times once for each metric I'd like to be able to automatically get another metric.
Any suggestion will be greatly appreciated. Currently running Windows XP, R 2.11.1 ############################################################################### Marsh <- cbind(SoilVegHydro, vegcode) AnnualDepth <- Marsh[,'meanAnnualDepthAve'] cattail_0 <- Marsh[,'cattail'] == '0' # no need to run for 8 species, automate if possible cattail_1 <- Marsh[,'cattail'] == '1' # yes need to run for 8 species spbase.rate.d1 <- sum(cattail_1)/(sum(cattail_1) + sum(cattail_0) ) annualDepth.density <- density(AnnualDepth) # this line needs either interactively defined or automatically cycle thru a number of metrics cattail.d0.density <- density(AnnualDepth[cattail_0]) cattail.d1.density <- density(AnnualDepth[cattail_1]) approxfun (cattail.d0.density$x, cattail.d0.density$y) -> cattail.d0.f approxfun (cattail.d1.density$x, cattail.d1.density$y) -> cattail.d1.f p.d.given.AnnualDepth <- function(AnnualDepth, spbase.rate.d1) { p1 <- cattail.d1.f(AnnualDepth) * spbase.rate.d1 p0 <- cattail.d0.f(AnnualDepth) * (1 - spbase.rate.d1) p1/(cattail_0+cattail_1) } x <- 1:1292 y <- p.d.given.AnnualDepth(x, spbase.rate.d1) plot (x, y, type='l', col='red', xlab='Mean Annual Depth', main=c("Cattail"), ylab='estimated\nProbability(cattail|AnnualDepth)') plot (cattail.d0.density, col ='red', lty= 1, main = "") lines(cattail.d1.density, col = 'blue', lty=1) lines(annualDepth.density , col='black', lty = 1) legend(2000, 0.0023, c("No Cattail", "Cattail", "Mean Annual Depth"), col=c("red", "blue", "black"),lty=c(1)) ##################################################################################################### Steve Friedman Ph. D. Spatial Statistical Analyst Everglades and Dry Tortugas National Park 950 N Krome Ave (3rd Floor) Homestead, Florida 33034 steve_fried...@nps.gov Office (305) 224 - 4282 Fax (305) 224 - 4147 ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.