surface plots using wireframe, color at high res, grid at low res?

Hello,

I'm working on some simple copula plots for a poster presentation, 
and by simple I mean simple, gaussian, etc., etc., etc.

I know I might be taking this a bit too far but I'm trying to see if 
I can overlay one surface plot over another.  I want to use color 
but I like the addition of a black wireframe to hold together a 
better sense of shape.

Right now, I have it set up to create a 300x300 colored surface 
using col="transparent", and I'd like to add a n x n (10, 20, 
whatever) on top of this.  Right now I'm using wireframe and I'm 
using my employer's official colors for gradient coloring.  My code 
is below.  I apologize for the code length but it does execute the 
way I want from a new session.

I'm also curious as to why I'm some grid effects or lighting effects 
(unsure of which) that I'm using on my personal 3.0.1 copy on 
Mountain Lion OSX, that I wasn't getting on 2.14.1 on Windows XP... 
something to figure out for another time.  This is a much more 
secondary concern for which I have a solution (compile the final 
graphic at work).  For now I'm just trying to see if its possible to 
lay a smaller resolution grid over the colored high resolution grid.

I will say that other help documents have been very useful to get 
the code along to this point.

Sincerely,

Patrick Joyce

Work e-mail:  patrick.m.joyce at census.gov


library(RColorBrewer)
library(colorspace)


library(lattice)

normal.copula=function(u1,u2,rho){

x1=qnorm(u1)
x2=qnorm(u2)

1/(2*pi)/sqrt(1-rho^2)*exp(-1/2*(x1^2+x2^2-rho*x1*x2)/(1-rho^2))/dnorm(x1)/dnorm(x2)

}

new.census.red="#AB0537"
new.census.blue="#0A6FB7"

ncred=hex2RGB(new.census.red)
ncblue=hex2RGB(new.census.blue)
ncwhite=hex2RGB("#FFFFFF")

census.color.function=function(alpha.cen){
out.hex=ifelse(alpha.cen<.5,hex(mixcolor(2*alpha.cen,ncblue,ncwhite),fixup=TRUE),ifelse(alpha.cen==.5,"#FFFFFF",hex(mixcolor(2*alpha.cen-1,ncwhite,ncred),fixup=TRUE)))

out.hex

}

rho.boat=.7

n.cuts=300
u1.corrd=(1:n.cuts-.5)/n.cuts
u2.corrd=(1:n.cuts-.5)/n.cuts
shareframe=data.frame(expand.grid(u1.corrd,u2.corrd))
shareframe$z.dir=normal.copula(shareframe$Var1,shareframe$Var2,rho.boat)
names(shareframe)=c("u1.dir","u2.dir","z.dir")

plot.new()
wireframe(shareframe$z.dir~ shareframe$u1.dir*shareframe$u2.dir,
  xlab = expression(U[1]), ylab = expression(U[2]), zlab="Density",main=" ",

scales = list(arrows=FALSE,col=1),
  drape = c(TRUE),

par.settings=list(axis.line=list(col="transparent")),col="transparent",colorkey 
= FALSE,
  screen = list(z = 195,y=180, x = 105),col.regions = 
census.color.function((1:10000-.5)/10000))


        [[alternative HTML version deleted]]