similar to: xlab at top in xyplot()

Dear r-helpers, In an xyplot I have xyplot(y ~ x | c, groups = g, scales = list( x = list( at = v1, labels = c('A', 'B', 'C', 'D'), rot = 45 ) ) ) g consists of two groups. How do I define a different set of x ticks and labels for the second panel? I'm looking for something like xyplot(y ~ x | c, groups = g, scales = list( x =

I am using xYplot to plot data and CIs. How do I add several lines to the figure? _____________________________ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400 Charlottesville, VA 22904-4400 Parcels: Room 102 Gilmer Hall McCormick Road Charlottesville, VA 22903 Office: B011 +1-434-982-4729 Lab: B019

Dear R-helpers, I have an xyplot of the following kind: xYplot(y ~ x | p, groups = factor(gg, levels = c('1', '5', '2', '4', '3')), strip = strip.custom(var.name = 'g', factor.levels = c(1', '5', '2', '4', '3'), strip.levels = T, strip.names = T, sep = ' = ', shingle.intervals = NULL), data =

Dear R-helpers, How do I tell panel.polygon what greoup and panel it applies to whithin xyplot(y ~ x | c, groups = g panel = function(x, y, groups, ...){ panel.polygon(x = xpol[c, g], y = ypol[c, g], default.units = 'native') panel.xYplot(x, y, groups, ...) llines(x = c(1, 6), y = c(-24.283333, 35.941667), lwd = 2, lty = 3, col = 4) } x[c, g] and y[c, g] describe

Dear R-helpers, I would like the bands in an xYplot() to be gray without dark edges. I've tried playing around with col.bands, but that seems not to affect the edges, although changin col.fill does change the fill as I would expect. I would appreciate pointers. _____________________________ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box

Dear R-helpers, x <- rep(1:2, 4) y <- c(2, 4, 3, 5, 1, 3, 2, 4) w <- factor(rep(1:2, each = 4)) v <- rep(1:2, each = 2, 2) xyplot(y ~ x | w, groups = v, type = 'b') How do I tell xyplot to use four colors and four plot characters and four line types? And how do I set up an appropriate key in each? _____________________________ Professor Michael Kubovy University of Virginia

Dear r-helpers, Does anyone have a straightforward example of putting together three unrelated (expect for a common y-axis) xyplot() figures in what would be in base graphics a par(mfrow = c(1, 3)) arrangement? _____________________________ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400 Charlottesville, VA 22904-4400 Parcels: Room 102

Hi, Does anyone know how to do a coord_trans() in which the y-axis is tranformed into (for example) -1000/y? Thanks, _____________________________ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400 Charlottesville, VA 22904-4400 Parcels: Room 102 Gilmer Hall McCormick Road Charlottesville, VA 22903 Office: B011

Dear r-helpers, In basic graphics, I have a figure with x and y axes suppresed. I would like to move the xlab and the ylab closer to the axes. Do I have to use mtext()? _____________________________ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400 Charlottesville, VA 22904-4400 Parcels: Room 102 Gilmer Hall McCormick Road

Dear r-helpers, tst <- data.frame(as.numeric(x <- 1:20), f <- rep(c('hi','lo'), times = 10)) with(tst, bwplot(f ~ x, panel = function(x, y){panel.bwplot(x, y, pch = '|', stats = boxplot.stats, fill = 8, varwidth = T)})) I can't figure out from the documentation how to tell stats that I would like to see notches or (even bands). Here is what I've

splom() doesn't complain here, but writes no names in the diagonal boxes. What am I missing? I believe that I need to add something like grid.text(x, ...) to the diagonal panel, but I don't know how to get it cycle through the column labels. And should varname.col = 'blue', varname.cex = 1 be inside the diag.panel() function? splom(szw[, n], pscales = 0, diag.panel

The following plot is a first approximation to what I need: *********************************** mu1 <- 0 mu2 <- 5 s <- 1 x <- seq(-2.5, 7.5, length = 41) y <- seq(-2.5, 2.5, length = 41) f <- function(x,y){ term1 <- 1/(2*pi*sqrt(s*s)) term2 <- -1/2 term3 <- (x - mu1)^2/s term4 <- (y - mu1)^2/s term5 <- (x - mu2)^2/s term1*(.5 * exp(term2*(term3 + term4)) + .5 *

Dear r-helpers, When I do an xYplot and display the result in a JavaGD() window, the font is sans-serif (presumably Helvetica). When I send the figure to a pdf, I get a serif font (presumably times). How do I insure that the font in the pdf is indeed the default sans serif? > sessionInfo() R version 2.4.1 (2006-12-18) i386-apple-darwin8.8.1 locale: C attached base packages: [1]

What am I doing wrong? (please cc me when replying) yy <- 1:10 xx <- yy*2 xYplot(yy~xx,ylab="") #this plots as it should text(.4,5,expression(paste(log, frac(p(b),p(a)) ))) Error in text.default(0.4, 13, expression(paste(log, frac(p(b), p (a))))) : plot.new has not been called yet This is the case even when I write: xxyy <- xYplot(yy~xx,ylab="") xxyy <-

I would like to suppress the text 'Scatter Plot Matrix' that appears under the plot. Could someone please suggest how? _____________________________ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400 Charlottesville, VA 22904-4400 Parcels: Room 102 Gilmer Hall McCormick Road Charlottesville, VA 22903 Office: B011

Dear R-helpers, How do I transform v <- c('insd-otsd', 'sppr-unsp') into c('insd--otsd', 'sppr--unsp') ? _____________________________ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400 Charlottesville, VA 22904-4400 Parcels: Room 102 Gilmer Hall McCormick Road Charlottesville, VA 22903

df1 <- data.frame(a = LETTERS[1:2], b = LETTERS[3:4], c = 1:2) I am looking for an idiom that swaps the elements of df$a and df$b when (e.g.) df$c == 2, resulting in df2 <- data.frame(a = LETTERS[c(1, 4)], b = LETTERS[c(3, 2)], c = 1:2) _____________________________ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400 Charlottesville, VA

Dear R-helpers, It is not clear to me how you get measurement.error SD when you have a single dataset, and it is not clear to me how sensitive SIMEX is to errors in the estimates of measurement error. Could someone please point me to the relevant literature? _____________________________ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400

library(glmmADMB) #Example for glmm.admb data(epil2) glmm.admb(y~Base*trt+Age +Visit,random=~Visit,group="subject",data=epil2,family="nbinom") Gives: Error in glmm.admb(y ~ Base * trt + Age + Visit, random = ~Visit, group = "subject", : The function maximizer failed ****************** R version 2.4.1 RC (2006-12-14 r40181) powerpc-apple-darwin8.8.0 locale: C

How does one define PDFs as yet undefined in R, such as the ex- gaussian, the sum of two RVs, one exponential, one Gaussian? The PDF would then be the convolution of an exponential PDF, dexp(), and a normal, dnorm(). Kindly cc me in your reply to r-help. Thanks, _____________________________ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400