similar to: accessing column and row numbers inside splom in lattice

I thought one nice addition to a splom figure would be to have the scatterplots in the upper triangle and a color-coordinated correlation matrix on the bottom. So I tried my hand at customizing panel.pairs(), and was rebuffed. Many times. Four hours of fruitless debugging later, I turn to you for help: panel.pairs(z=teststatfull[,6:12], pscales=0, panel.subscripts=FALSE, subscripts=,

Hello, I have a customized pairs () fonction as follows that displays correctely my data. ------------------------------------------------------------------------ panel.cor1 <- function (x, y, digits=2, prefix="") { usr <- par("usr"); on.exit(par(usr)) par(usr = c(0, 1, 0, 1)) r <- cor(x, y,use="pairwise.complete.obs",

Dear ExpeRts, I have the scatter plot matrix as given below. I would like the different "sub-plots" in the scatter plot matrix to be colored differently. How do I get all points shown in the upper-left plot (on position (1,1) in the scatter plot matrix) to be plotted in blue, and the points shown in the plot to the right (on position (1,2) in the scatter plot matrix) to be plotted in

Dear R-experts, I have found a splom-modification online which is given below. This works perfectly, but I would like to have a matrix of given correlation values to be used in the lower triangular part (lower.panel) of the splom-plot instead of calculated correlation values. Here is the matrix I would like to use (it can be any other convenient data structure):

Dear expeRts, I would like to create a scatter plot matrix with splom(). The lower panel should contain some additional information about the samples shown in the upper panel plot, see the splom() call below. Now two questions came up: (1) The lower panels show "tau" and "alpha" on top of each other. How can I plot *three* expressions on top of each other? I tried several

Hi, I am trying to produce a 30x30 lattice grid. The 30 variables naturally fall into three groups of ten and I would like to add thicker axis lines around these to highlight this. However, I can only do this separately, so far, and I don't know how to superimpose the grid with the thicker lines (separating the 10x10 groups) onto the overall grid of all 30x30 plots. I enclose a

Hello list, I am looking to create a figure for my dataset using splom, where there is a splom subplot for each level of a factor within, for example, a 2x2 layout. For each subplot, I wish to put the r-value between each variable pair in the lower panel. The code I have thus far is below, using the iris dataset so that it is reproducible. This almost does what I want, but for some reason

Sorry Folks, I'm sure I could suss out the answer myself but I need it soon ... ! 1. Given a set of 4 variables X,Y,Z,W in a dataframe DF, I make a scatter-plot matrix using splom(DF). 2. I do all regressions of U on V using lm(U~V), where U and V are all 12 different ordered pairs from X,Y,Z,W. 3. Now I would like to superpose the regression lines from (2) onto the corresponding

I'm interested in doing multinomial logistic regression with equality constraints on some of the parameter values. For example, with categorical outcomes Y_1 (baseline), Y_2, and Y_3, and covariates X_1 and X_2, I might want to impose the equality constraint that \beta_{2,1} = \beta_{3,2} that is, that the effect of X_1 on the logit of Y_2 is the same as the effect of X_2 on the

Hi all, I use the excellent packages rjags and cacheSweave, and unfortunately seem to have found an incompatibility between them. Below are a minimal .Rnw file and corresponding JAGS model file which illustrate the problem: *** JAGS model file; name=j.bug *** model { mu ~ dnorm(0,1.0E-5) sigma ~ dunif(0,100) for(i in 1:length(y)) { y[i] ~ dnorm(mu,sigma) } } *** *** .Rnw

Dear R Users, When plotting with "splom" I tried to use the pscales=list(...) feature Unfortunately it didn't work at all. Instead the scales always were suppressed. So I looked at the source of panel.pairs and found: "draw <- is.numeric(pscales) && pscales != 0" which rather has to be: "draw <- is.list(pscales) || (is.numeric(pscales) &&

I'm interested in writing a function that constructs a new plot on the current graphics device if no plot exists there yet, but adds lines to the existing plot if a plot is already there. How can I do this? It seems to me that the exists() function might be co-opted to do this, but it's not obvious how. Many thanks, Roger -- Roger Levy Email: rlevy at ucsd.edu

I'm writing a book using Sweave which includes Unicode characters which R's postscript and pdf graphics devices don't handle properly. The cairo graphics device does handle these characters properly. Thus I'd like to be able to instruct Sweave to use cairo (either sometimes or always) to produce graphics when it processes R code. Is this possible? Many thanks, Roger Levy --

Dear list, I want to plot several variables with splom and in the main diagonal, instead of the variable names, I'd like to plot an histogram of corresponding variables. Searching I did not find the correct syntax, only some tips in an old post in the list, but this comments help to plot only density lines instead of histograms. I had some code, but it fails to plot (I've commented the

Hi, I'd like to be able to generate a splom plot in R and then use my mouse to click on one of the sub-panels (panel.pairs, specifically) and have R return either the coordinates of that sub-panel, or even better, the names of the corresponding variables plotted in that sub-panel. Here's an example to work with: > library(lattice) > splom(~iris[1:4], groups = Species, data = iris,

Hi, Kindly let me know if I posted a wrong question in the forum. I want to draw a splom plot with different symbols in plot. My command is as follows: splom(~ log10(splomData[2:3]), groups = programs, data = splomData, panel = panel.superpose, key = list(title = paste(splomLoop,"Programs of Hog Analysis (Sorted by LR(GB))"), panel = panel.superpose, columns =

I'm not seeing how to plot the quantities associated with two values of a factor by reading the ?xyplot help page and Deepayan's book. Perhaps I need to generate many different subsets from the data frame, but that would require _many_ new data frames. The structure of the data frame (for a single stream) has sites (factor), a date (as.Date), parameters (factor), and quantities

Dear expeRts, how can I decrease the space between the tick marks and the corresponding labels in an splom? See here: library(lattice) U <- matrix(runif(4000), ncol = 8) splom(U, axis.text.cex = 0.2) # => space between the [small] tick labels and tick marks is/seems to be too large I checked ?panel.pairs but could not find an option for that. Cheers, Marius

All, I have long variable names that are being fed through splom (R 1.8.1). I'd like to resize the text printed on the diagonals to better display the names (unless anyone can suggest another approach - creative use of varnames). I've looked at the code, R-Help, ?splom and the Trellis User's Guide to no avail. Any suggestions? Thanks in advance, Max Kuhn, Ph.D. Becton Dickinson

Using the mtcars dataset, how to define the grouping variable to be valid only for the upper or lower panel? The following doesn't work: # Code start splom(~data.frame(mpg, disp, hp, drat, wt, qsec), data=mtcars, pscales=0, auto.key=list(columns=3), upper.panel = function(...){ panel.grid(...) panel.xyplot(groups=cyl,...) } ) # Code end -- View