similar to: Problem w/ axis and distortion in a plotting function

Dear R users, has anybody made an attempt to create the soil texture triangle graph in R? For an example see here: http://www.teachingkate.org/images/soiltria.gif I would like to get the lines in black and texture labels in gray to allow for plotting my texture results on top. Any examples or suggestions are very welcome! Thanks in advance, Sander. --

Dear list: I encountered some problems using the function point.in.polygon() of the sp package, when trying to determine whether some points lye inside, outside, on the border or on a vertice of a polygon. I have a list of point I know should lye right on the border of a polygon, but some of them are not classified as such by point.in.polygon() (see the example code below). To make a long story

Dear R users, If I manage to sort out this graph, it is certainly a candidate for the new R graph gallery (http://addictedtor.free.fr/graphiques/displayGallery.php)! I created the following lattice graph: library(lattice) tmp <- expand.grid(geology = c("Sand","Clay","Silt","Rock"), species =

Hi, I am trying to create a triangular plot to show the 'composition' of a set of items with three variables (historically the percent sand, silt and clay in soil). So far I have tried the 'soil texture triangle plot' in the package plotrix and the 'ternary or triangular plots' in the package cwhtool (cwhmisc). Both have strengths and weaknesses, but neither has a

Hello, Using the dataset below I produce a soil.texture plot (R code for this included at the end). One will notice I have plotted the points based on whether or not they are in a region called 'upstream' or 'downstream'. I'm curious if there is a way to somehow extract counts of the number of points (red points and blue points) falling within each classification (e.g. silty

Hi, I am interested in ggplot2 and I found this lattice code very interesting (http://addictedtor.free.fr/graphiques/graphcode.php?graph=48). Code: library(lattice) lattice.options(default.theme = canonical.theme(color = FALSE)) tmp <- expand.grid(geology = c("Sand","Clay","Silt","Rock"), species = c("ArisDiff",

I am trying to reproduce an analysis done in SAS from which Type III Sums of squares are used. How can you get Type III sums of squares in an aov similar to what SAS provides? With best wishes and kind regards I am Sincerely, Corey A. Moffet Instructor Department of Range, Wildlife, and Fisheries Management Mail Stop 2125 Texas Tech University Lubbock, Texas 79409-2125 (806) 742-2843 Voice

Dear R graphics gurus, Another question about lattice graphics. This time I would like to plot means and confidence intervals by group factor in a lattice graph. I can not find any working lattice examples. Maybe a custom panel function is the answer, but that is a bit beyond me for now. The individual plots within the lattice graph could look like this: # Example with confidence intervals

I have some data that was created for import into a Lotus 1-2-3 spreadsheet and on of the columns is time. The time is akin to Julian were the value 1 is mapped "01-Jan-00 12:00:00 AM" in Lotus 1-2-3. Is there a function in an R package that can convert this numeric vector to a POSIXct vector? With best wishes and kind regards I am Sincerely, Corey A. Moffet Instructor Department

Hello, I've been looking through ?phantom and ?expression and this forum for examples of how I might be able to manipulate some of the names that appear on the y-axis of the barplot below. For example, the "gw" in "ECgw" would appear as a subscript...or "qr" would be the theta symbol followed by subscript "r". My attempts haven't even come close

Is there a function/package that will allow you to test the hypothesis beta1 = x in a simple linear regression, where x is a constant? The AOV tests if beta1 in different from 0 but what if I want to know if beta1 one is different from 1? With best wishes and kind regards I am Sincerely, Corey A. Moffet Instructor Department of Range, Wildlife, and Fisheries Management Mail Stop 2125 Texas

Dear R-help: Does anyone know of a package in R that will do Welch's modified F-test for heterogeneous error variances? Are there other statistical techniques available in R that test the equality of means when homoscedastisity is violated? 't.test' does this in the pairwise sense when var.equal = TRUE. With best wishes and kind regards I am Sincerely, Corey A. Moffet Support

Dear R-Help, Does anyone know of a way to use CGIwithR in Windows? David Firth the author of CGIwithR states "It ought to be possible to adapt it for use also in conjunction with web servers running under Windows NT ... but the author has no plans for that." I have developed some applications in R that I would like to make available through a web interface using CGIwithR but at this

I am having difficulty with plot.formula. For example: ------------------- dat <- data.frame(y = c(1,1,1,2,2,2,3,3,3,4,4,4), x1 = as.ordered(c(1,1,1,1,1,1,2,2,2,2,2,2)), x2 = as.ordered(c(1,1,1,2,2,2,1,1,1,2,2,2))) plot(y ~ x1 + x2, data = dat) # Works fine "interactively" asking: # "Hit <Return> to see next plot:" #

Dear R users I'd like to know if it is possible to create conditional plots in the soiltexture package. I have data from the particle size distribution of soils derived from two parent materials and I want to use the parent material as conditioner. Thanks a lot. Alessandro Samuel-Rosa Post-Graduation Program in Soil Science Federal University of Santa Maria Av. Roraima, n? 1000, CEP

o <- (structure(list(sand.silt = structure(list(statistic = structure(185, .Names = "W"), parameter = NULL, p.value = 0.0478835773838087, null.value = structure(0, .Names = "location shift"), alternative = "two.sided", method = "Wilcoxon rank sum test with continuity correction", data.name = ".column by site"), .Names =

hello, I am very new to R. My current data set is a mix of values and categories. It is a geoscience data set, with values per rock sample. Case in point, each sample belongs to a lithology class, and each sample has several physical property measurements (density, porosity...). I want to be able to plot these physical properties for all samples in each lithology class. this is how i'm doing

I would like to preform a t.test to each of the measured variables (sand.silt etc.) with a mean and sd for each of the treatments (up or down), and out put this as a table.... I am having a hard time starting- maybe it is to close to lunch. Any suggestions would be greatly appreciated. Stephen Sefick x <- (structure(list(sample. = structure(c(1L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 2L, 3L,

Dear List, I'm using gam in a multiple imputation framework -- specifying the knot locations, and saving the results of multiple models, each of which is fit with slightly different data (because some of it is predicted when missing). In MI, coefficients from multiple models are averaged, as are variance-covariance matrices. VCV's get an additional correction to account for how

Hi everyone! I'm trying to lear how to call external C code in R but even the R help web-page example is drive me crazy. I copy-paste the example at: http://stat.ethz.ch/R-manual/R-patched/doc/manual/R-exts.html#Interface%20functions%20.C%20and%20.Fortran here is the example: void convolve(double *a, int *na, double *b, int *nb, double *ab) { int i, j, nab = *na + *nb - 1; for(i = 0;