search for: soilveghydro

Hello, vegMeans <- by(SoilVegHydro[3:37] , SoilVegHydro$Physiogomy, mean) vegSD <- by(SoilVegHydro[3:37] , SoilVegHydro$Physiogomy, sd) write.table(vegMeans, file="A:\\Work_Area\\Steve\\Hydrology_Data\\data\\vegMeans.txt") Error in as.data.frame.default(x[[i]], optional = TRUE, stringsAsFactors = stringsAsFactors) :...

hello, I am trying to develop a triangle plot but am having difficultly assigning the row.names to the 3 columns in the data.frame Here is what I've done, attach(SoilVegHydro) dim(SoilVegHydro) 1292 39 # now take 3 variables from main data.frame for plotting dat <- cbind.data.frame(TP, meanAnnualDepthAve, BulkDensity) # These are variables held in the data frame SoilVegHydro row.names(dat) <- paste(row.names(SoilVegHydro$Physiogomy), rep(c(1,2,3), rep(12...

...ing 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.d...

...I'd like this comparison to be based on probabilities. An example of the graph I'm seeking can be seen here: http://www.statmethods.net/graphs/density.html (scroll to the bottom and examine the Comparing Groups via Kernel Density). A graph (attached) produced using my data: > attach(SoilVegHydro) > vegtype.f <- factor(Physiogomy) > sm.density.compare(meanAnnualDepthAve, vegtype.f, xlab="Mean Annual Depth") # create color fill patterns >colfill <- c(2:(2+length(levels(vegtype.f)))) > legend(locator(1), levels(vegtype.f), fill=colfill) (See attached fil...