search for: meanannualdepthav

Hello, I'm trying to calculate non-parametric probabilities using the np package and having some difficulties. OS is Windows, R version 2.11.1 Here is what I've done so far. library(np) veg <- data.frame(factor(Physiogomy), meanAnnualDepthAve, TP) attach(veg) : for clarification dim(veg) returns 1292 3 fy.x <- npcdens(veg$factor.Physiogomy ~ veg$meanAnnualDepthAve, nmulti=1) # this works, but I haven't found any information explaining what the nmulti=1 term is doing? Does this set the number of levels in the factor? My data...

...ems that it did not appear on the list. If you've rec'd 2 m I'm trying to calculate non-parametric probabilities using the np package and having some difficulties. OS is Windows, R version 2.11.1 Here is what I've done so far. library(np) veg <- data.frame(factor(Physiogomy), meanAnnualDepthAve, TP) attach(veg) : for clarification dim(veg) returns 1292 3 fy.x <- npcdens(veg$factor.Physiogomy ~ veg$meanAnnualDepthAve, nmulti=1) # this works, but I haven't found any information explaining what the nmulti=1 term is doing? Does this set the number of levels in the factor? My data...

...nodes, my analysis produces 55 terminal nodes, so unless someone is really interested I'm only using the iris data as an example. This is the model I am sending to ctree. The graph produced is very informative, but I need the information from the plot(Marshct) Marshct <- ctree(Physiogomy ~ meanAnnualDepthAve + meanAnnualDepthWetAve + medAnnualDepthAve + medianAnnualDepthWetAve + medianAnnualDepthDryAve + continHydroWetAve + DCHperiodAverage + DCHydroDryAve + threeDayWaterDepthMinAve + threeDayWaterDepthMaxAve + sevenDayDryFreqAve + sevenDayWaterDepthMaxAve + sevenDayWaterDepthMinAve + seventeenDayWate...

...o, 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(1292, 3)), sep =" ") # following the syntax from the help triangle.plot page this is returned when the last line is submitted. ro...

...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) )...

...aph 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 file: MeanAnnualDepth_Density.png) Thanks Steve Steve Friedman Ph. D. Spatial Statis...