similar to: Data format question for triangle.plot package ade4

Hello, This is a re-submittal of question I submitted last week, but haven't rec'd any responses. I need to extract the probabilities used to construct the barplots displayed as part of the graph produced by plot("ctree"). For example, library(party) iris.ct <- ctree(Species ~ . , data = iris) plot(iris.ct) Instead of a simple example with only 4 terminal nodes, my

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) : cannot coerce class '"by"' into a

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

Hello, I would like to produce a series of graphs comparing the probability distributions for 8 factors against a continuous metric. The kind of graph I'm hoping to produce would look like the density comparison graphs (library sm) using the function sm.density.compare. However, instead of calculating the density distributions for comparisons, I'd like this comparison to be based on

Hello Trying to develop variograms and kriged surfaces from a point file. Here is what I've done so far. library(gstat) # also loads library(sp) library(lattice) soilpts$x <- soilpts$UTM_X soilpts$y <- soilpts$UTM_Y soil.dat <- subset(soilpts, select=c(x, y, Area, BulkDensity, LOI, TP, TN, TC, Total_Mg)) dim(soil.dat) [1] 1292 7 coordinates(soil.dat) <- ~ x+y

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

Hello, I tried to post this earlier, but it seems 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

Hello, I need to chenge axes orirentation in triangle plot. (function triangle.plot in ade4 package) I want to plot elasticities of some species in demographic triangle, where axes values commnly increace "clockwise". If some better imangination is needed, see http://www.open.ac.uk/science/biosci/personalpages/j.silvertown/pdfs/Silvertown%20et%20al.%201993.pdf I am sorry if I just

Dear all, mydata<-data.frame(x1=c(1,4,6),x2=c(3,1,2),x3=c(2,1,3)) cor(mydata) x1 x2 x3 x1 1.0000000 -0.5960396 0.3973597 x2 -0.5960396 1.0000000 0.5000000 x3 0.3973597 0.5000000 1.0000000 I wonder if it is possible to fill only lower triangle of this correlation matrix? Using 'dist' doesn't seem to be useful as it doesnt allow to convert this table

Hello, I am using triangle.class in ade4 package and a factor. I would like to use this graphic parameter (pch) to set a different kind of point for each group in my factor (lakes). I have tried pch=1:4, I guess it does what it is supposed to do but not what I want. I can't use colors. Thank you in advance, -- View this message in context:

Hello, I am trying to run some AMOVA analyses with the amova function in the package ade4. When running the example dataset provided in ade4, I noticed a difference between the published results from the same data (Excoffier et al. 1992) and what ade4 calculates. Below are the data for "within sample/population" from ade4 and from the haplotypic distance matrix in the paper:

Hello, The ade4 package is now available on CRAN. ade4 is a multivariate data analysis package oriented towards statistical ecology. The main functions include one-table (PCA, COA, MCA, PCO, etc.), two-tables (PCAIV, Co-inertia, CCA, etc.) and k-tables multivariate analysis methods, with several Monte-Carlo tests and many graphical display methods. A complete documentation with many example data

Hello, The ade4 package is now available on CRAN. ade4 is a multivariate data analysis package oriented towards statistical ecology. The main functions include one-table (PCA, COA, MCA, PCO, etc.), two-tables (PCAIV, Co-inertia, CCA, etc.) and k-tables multivariate analysis methods, with several Monte-Carlo tests and many graphical display methods. A complete documentation with many example data

Dear List, I tried to conduct cca using csv data but failed. The message said that data.frame is required. Please kindly share how to convert a csv-imported file to a data.frame. Thank you. Elaine code rm(list=ls()) spec <-read.csv("c:/migration/M_R_20100718_winterM_spec_vegan.csv",header=T, row.names=1) dim(spec) spec[1,] envi

Dear R-Users, ade4TkGUI is a new package available on CRAN. It implements a Tcl/Tk graphical user interface (GUI) for the ade4 package. Only the most basic functions of ade4 have a GUI in this first version : classical one-table data analysis methods (PCA, COA, MCA, PCO, etc.), one table with groups of rows (BGA, WGA, DA), and two-tables analysis methods (Coinertia analysis, CCA, PCAIV).

Dear R-Users, ade4TkGUI is a new package available on CRAN. It implements a Tcl/Tk graphical user interface (GUI) for the ade4 package. Only the most basic functions of ade4 have a GUI in this first version : classical one-table data analysis methods (PCA, COA, MCA, PCO, etc.), one table with groups of rows (BGA, WGA, DA), and two-tables analysis methods (Coinertia analysis, CCA, PCAIV).

Dear List, I used spec and envi for cca in ade4. (both are data.frame) However, there is a message telling that error in if (nf > rank) nf <- rank R missing value in TRUE/FALSE Please kindly help how to modify the code below. Thank you. Elaine code rm(list=ls()) spec <-read.csv("c:/migration/M_R_20100718_winterM_spec_vegan.csv",header=T, row.names=1) dim(spec) spec[1,]

Hello I'm using the package ade4 to obtain classification from a .txt file. I use the following commands: cronquist <- read.table("cronquist.txt", h = T, row.names = 8) cronquist <- as.taxo(cronquist[7:1]) cro.phy <- taxo2phylog(cronquist) in which cronquist.txt is a file with a classification of 218 genus of tree species from french Guiana. We I try to use the

Un saludo a todos. Estoy haciendo un ACC y utilizo el paquete Ade4. Cuando quiero realizar un plot del resultado utilizo plot(). el resultado son 6 gráficos: loadings, correlation, inertia axes, scores and predictions, eigenvalues y species. Lo que quería saber es cómo representar en una sola ventana (los 6 gráficos a parecen en la misma) el gráfico correspondiente a Species. Gracias por vuestra

Was there ever an answer to this question about how to do amova with factorial data? t hundley [[alternative HTML version deleted]]