similar to: xyplot() and controlling panel.polygon()

I'm looking for a way to plot lines on top of a levelplot(), where the lines are borders between cells of different values. The clines() function provides contours suitable for continuous data. I am dealing with discrete values, spatial clusters of nodes where each cluster has an integer value, and I want to plot the borderlines between these areas. So, in a levelplot having nodes with

Dear r-help list members, I'm experiencing problems getting type="b" (or "o" or "c") to work in llines(). Try, e.g., the following scaled-down example: x <- factor(c("a", "b", "a", "b")) y <- c(1,2,1,2) z <- factor(c("A", "A", "B", "B")) symbols <- 1:2 lines <- 1:2

Another question about lattice: Is there a way to use plotmath in the labels and strips? E.g. xyplot(y~x|group,data=test.df, panel=function(x,y,subscripts,df) { llines(x,y,col="black",lwd=2) llines(x,df[subscripts,"lowerCI"],lty=2,col="#aaaaaa") llines(x,df[subscripts,"upperCI"],lty=2,col="#aaaaaa") },

Dear R-Users, how to include error bars within lattice? How should the panel = function(x,y,...){ looks like? Does panel.arrows works here as well? I appreciate any help on this. Regards, Mario AT [[alternative HTML version deleted]]

If I have a contourplot (in the lattice package) and I want to add straight lines to it, how do I do this? I see that there are llines() and lsegement() functions for lattice plots, but they don't seem to do anything in this case: library(lattice) library(KernSmooth) x=rnorm(10000) y=x+rnorm(x,0,.5) a=bkde2D(cbind(x,y),.7) z=as.vector(a$fhat) grid=expand.grid(x=a$x1,y=a$x2) grid$z=z

Hi, I have a data frame like this: ScoreDoseSex 2.81Dose1M 1.81Dose1M 1.22Dose1M 0.81Dose1M 0.49Dose1M 0.22Dose1M 0.00Dose1M -0.19Dose1M -0.17Dose1F -0.32Dose1F -0.46Dose1F -0.58Dose1F -0.70Dose1F -0.81Dose1F -0.91Dose1F -1.00Dose1F -1.77Dose2M -1.85Dose2M -1.93Dose2M -2.00Dose2M -2.07Dose2M -2.14Dose2M -2.20Dose2M -2.26Dose2M -2.32Dose2F -2.38Dose2F -2.17Dose2F -2.49Dose2F -2.54Dose2F

Hi, I was taken off guard by the following behavior in a lattice plot. I frequently want to add a predicted curve defined at more points than in the formula expression of xyplot. There have been numerous examples of how to do this on r-help, but I still often struggle to make this work. I just realized that specifying one of the axes on a log scale does not guarantee that the added data for a

Hi all, I have been previously been using scatterplot3d package to create some graphs but unfortunately it does not allow me to rotate the plot on all three axis. The cloud() function in the lattice package does allow me to do so. When I was using scatterplot3d I was using a script (Shown Below) to calculate the mean, quartiles and range limits for all three axis and I was representing that on the

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 lattice wizards, I am trying to figure out how to plot predicted values in xyplot, where the intercept, but not the slope, varies among conditioning factor levels. I am sure it involves the groups, but I have been unsuccessful in my search in Pinhiero and Bate, in the help files, or in the archive, or in my attempts on my own. My example follows: FACT is a factor with levels a,b,c

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, I would like to use the effect() function (actually a slightly modified version of it) on the output of the lmer() function in the lme4 package. But the effects package requires the nlme pacvkage, which is incompatible with lme4. Workaround? ______________________________________________ Professor Michael Kubovy University of Virginia Department of Psychology for mail add: for FedEx or

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

Dear Friends, Here is an algorithm for finding concave hulls: http://get.dsi.uminho.pt/local/ Has anyone implemented such an algorithm in R? RSiteSearch('concave hull') didn't reveal one (I think). _____________________________ Professor Michael Kubovy University of Virginia Department of Psychology Postal Address: P.O.Box 400400, Charlottesville, VA 22904-4400 Express Parcels

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

R folks, I can't seem to find the instructions in the help files for the lattice package that explain how to add lines, such as with lines() or ?, to a levelplot. I'd be grateful if someone could point me in the proper direction. Cheers, Jeff ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ Jeff Jorgensen Center for Limnology jcjorgensen at wisc.edu

Dear r-helpers, xx <- c(0.000, 0.210, 0.714, 0.514, 1.000, 0.190, 0.590, 0.152) yy <- c(0.000, 0.265, 0.256, 0.521, 0.538, 0.761, 0.821, 1.000) aa <- c(19, 19, 19, 21, 19, 21, 21, 21) x0 <- xx[c(1, 1, 2, 2, 2, 3, 3, 4, 4, 4, 5, 6, 6, 7)] y0 <- yy[c(1, 1, 2, 2, 2, 3, 3, 4, 4, 4, 5, 6, 6, 7)] x1 <- xx[c(2, 3, 3, 4, 6, 4, 5, 5, 6, 7, 7, 7, 8, 8)] y1 <- yy[c(2, 3, 3, 4, 6, 4, 5,

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