R help - Aug 2007 - Using sunflowerplot to add points in a xyplot panel

Hi,

I use panel.points to add points to a xyplot graphic. But I like to use the 
sunflowerplot to plot my points because this is very superimposed. It is 
possible to use this? I try but it dont work directly. It may be need to put 
this function inside a panel.??? 

Thanks
Ronaldo
-- 
Where there's a will, there's a relative.
--
> Prof. Ronaldo Reis J?nior
|  .''`. UNIMONTES/Depto. Biologia Geral/Lab. de Ecologia
| : :'  : Campus Universit?rio Prof. Darcy Ribeiro, Vila Mauric?ia
| `. `'` CP: 126, CEP: 39401-089, Montes Claros - MG - Brasil
|   `- Fone: (38) 3229-8187 | ronaldo.reis em unimontes.br | chrysopa em
gmail.com
| http://www.ppgcb.unimontes.br/ | ICQ#: 5692561 | LinuxUser#: 205366

Ronaldo Reis Junior said the following on 8/14/2007 7:08
AM:
> Hi,
> 
> I use panel.points to add points to a xyplot graphic. But I like to use the
> sunflowerplot to plot my points because this is very superimposed. It is 
> possible to use this? I try but it dont work directly. It may be need to
put
> this function inside a panel.??? 
> 
> Thanks
> Ronaldo
You'll need to write your own panel function. Here's one shot at it. 
Most of the code is from ?sunflowerplot with added touches for lattice 
capability.

HTH,

--sundar

panel.sunflowerplot <- function(x, y, number, log = "", digits = 6,
rotate = FALSE,
                                 cex.fact = 1.5, size = 1/8, seg.col = 
2, seg.lwd = 1.5, ...) {
   n <- length(x)
   if(missing(number)) {
     x <- signif(x, digits = digits)
     y <- signif(y, digits = digits)
     orderxy <- order(x, y)
     x <- x[orderxy]
     y <- y[orderxy]
     first <- c(TRUE, (x[-1] != x[-n]) | (y[-1] != y[-n]))
     x <- x[first]
     y <- y[first]
     number <- diff(c((1:n)[first], n + 1))
   } else {
     if(length(number) != n)
       stop("'number' must have same length as 'x' and
'y'")
     np <- number > 0
     x <- x[np]
     y <- y[np]
     number <- number[np]
   }
   n <- length(x)
   n.is1 <- number == 1
   cex <- trellis.par.get("plot.symbol")$cex
   if(any(n.is1))
     lpoints(x[n.is1], y[n.is1], cex = cex, ...)
   if(any(!n.is1)) {
     lpoints(x[!n.is1], y[!n.is1], cex = cex/cex.fact, ...)
     i.multi <- (1:n)[number > 1]
     ppin <- par("pin")
     pusr <- unlist(current.panel.limits())
     xr <- size * abs(pusr[2] - pusr[1])/ppin[1]
     yr <- size * abs(pusr[4] - pusr[3])/ppin[2]
     i.rep <- rep.int(i.multi, number[number > 1])
     z <- numeric()
     for (i in i.multi) z <- c(z, 1:number[i] + if (rotate) 
stats::runif(1) else 0)
     deg <- (2 * pi * z)/number[i.rep]
     lsegments(x[i.rep], y[i.rep], x[i.rep] + xr * sin(deg),
               y[i.rep] + yr * cos(deg), col = seg.col, lwd = seg.lwd)
   }
}

library(lattice)
xyplot(Petal.Width ~ Petal.Length, iris, panel = panel.sunflowerplot)

Ronaldo Reis Junior <chrysopa <at> gmail.com> writes:
> 
> Hi,
> 
> I use panel.points to add points to a xyplot graphic. But I like to use the
> sunflowerplot to plot my points because this is very superimposed. It is 
> possible to use this? I try but it dont work directly. It may be need to
put
> this function inside a panel.??? 
sunflower is "old style R graphics", and xyplot is lattice/trellis.
these two
normally don't get along together too well. Paul Murrell's grid which is
the
base of lattice, can do some magic to bring them together.

Dieter


## Modfied after Cris Bergstresser and Paul Murrell
##http://finzi.psych.upenn.edu/R/Rhelp02a/archive/73780.html
library(grid);
library(lattice);

sunpanel <- function(x, y, subscripts, ...) {
  pushViewport(viewport(x = 0.5, y = 0.5, just = "center"));
  sunflowerplot(x, y, axes = FALSE, xlab = "", ylab = "");
  popViewport(1);
} 

x = round(runif(100, 1, 5));
y = round(runif(100, 1, 5));
print(xyplot(y ~ x, aspect = "iso", panel = sunpanel));