similar to: Help needed with lattice graph!

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

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",

Dear all, I have a data frame containing the results of an experiment. Like this: a<-seq(1,4,by=1) b<-seq(1,2,by=1) test<-expand.grid(b,a,a) colnames(test)<-c("replicates","bins", "groups") test$abc <- rnorm(32) test$def <- rnorm(32) test$ghi <- rnorm(32) test The following code snippet aggregates the data for one variable and then draws a

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 All, I wonderif there is a simple way to draw a regression line in the xyplot: more specifically, let: age <- c(20:30, 31:40 ) age.cut <- cut(age, breaks = 2 ) y<- rnorm(20) x <- rnorm(20,4,1) xyplot(y ~ x| age.cut, xlab="x", ylab="y") How to draw (in the plot given by xyplot) the two regression lines (y ~ x) corresponding to the two category

I have a function for plotting soil texture that I am having a problem with. The function is: plot.psd <- function(sand, clay, ...) { conv.ter <- function(x, y) { x.con <- 100 - x - y*cos(1/3*pi) y.con <- y*sin(1/3*pi) data.frame(x=x.con, y=y.con) } plot(conv.ter(sand, clay), xlim = c(0,100), ylim = c(0,100*sin(1/3*pi)), axes = FALSE, xlab="",

On Mon, Oct 24, 2011 at 9:17 AM, Devang Patel <dpatel at apple.com> wrote: > > On Oct 23, 2011, at 12:03 AM, Talin wrote: > > Just a follow up on this - I am still having problems, I never did figure > out a solution. (I've been running with debug off for the last month so that > I could get work done.) > > Here's what I am seeing: I am definitely calling

Hi, If I have an instance of DISubprogram, can I get the debug info of local variables of the function, including parameters? I tried to use the getVariables() function defined in DISubprogram, but it seemed to return an empty DIArray node when I ran my pass alone using opt. Do I need to enable any other analysis passes in order to populate the data? My related snippet of code is like the

Dear all, I am not fluent in R and am struggling to 1) apply a lm to a weight-size dataset, thus the model has to run separately for each species, each year; 2) extract coefs, r-squared, n, etc. The data look like this: year sps cm w 2009 50 16 22 2009 50 17 42 2009 50 18 45 2009 51 15 45 2009 51 16 53 2009 51 17 73 2010 50 15 22 2010 50 16 41 2010 50 16 21 2010

Hi everyone, I'm trying to import .sps (SPSS portable file) file. the read.spss function (library foreign) doesn't allow to import such files. should I import in spss and then save as sav file? there is not other solutions available? what I mostly like from spss file is that they have variable labels. want is really wish to keep are the variable.labels from the spss file; so, if there is a

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

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 =

Dear all, I get stuck when i try to export the data into SPSS format/file using "write.foreign()" Do you know how to do it exactly? What i have done is 1) First i type the following code in R: df<-data.frame(id,year,res1) names(df)<-c("idcode","year","resarrvl")

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

On Oct 27, 2011, at 10:39 PM, Talin wrote: > On Mon, Oct 24, 2011 at 9:17 AM, Devang Patel <dpatel at apple.com> wrote: > > On Oct 23, 2011, at 12:03 AM, Talin wrote: > >> Just a follow up on this - I am still having problems, I never did figure out a solution. (I've been running with debug off for the last month so that I could get work done.) >> >>

Hello, I have a Salicru SPS One 700VA device, I want to manage it with Debian Whezzy. I read some information in the net about how to configure but I can not get it working. Could you help with this? This is the log when I connect the USB: [ 3194.688204] usb 4-4: new low-speed USB device number 2 using ohci_hcd [ 3194.857239] usb 4-4: New USB device found, idVendor=0665, idProduct=5161 [

As previously noted, I need to replace a UPS and want to do my product research in public so others can benefit. The information I elicit on this thread will be merged into the UPS-HOWTO. I went to my local computer-equiuipment big-box store, a place called MicroCenter in Ardmore, Pennsylvania. I looked at the range of "UPS" devices they sell and was *completely* confused. Bear in

Dear R Users, I have an xyplot() where different plotting symbols are used for subgroups (originally used within S-Plus, but hopefully it?s also applicable to R users). How can I fit separate regression lines for every subgroup? So far, I can only plot the overall fitted line. The code looks like this: trellis.device() sps<-trellis.par.get("superpose.symbol") sps$pch<-1:7

Dear R Lists, Can anyone help me add site IDs (site: 1~50) directly to my xyplot. I have 50 sites and collected observations from the sites at 13 different time points. I want to look at the change of my observations in each site. I was able to make a graph using xyplot, however, I can't find how to add site IDs within the graph not adding them in the legend. Overall, I would like to be able