Hi Patrick
You should do better by using a specific graphic device and then send
the file to the printer
...
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")
trellis.device(device = postscript,
file = "copula1.ps",
title = "",
colormodel = "srgb")
wireframe(shareframe$z.dir~ shareframe$u1.dir*shareframe$u2.dir,
scales = list(arrows=FALSE,col=1),
drape = c(TRUE),
par.settings=list(axis.line=list(col="transparent")),
col="transparent",
colorkey = FALSE,
paper = "special",
horizontal = "FALSE",
height = 18,
width =18,
screen = list(z = 195,y=180, x = 105),
col.regions = census.color.function((1:10000-.5)/10000),
xlab = expression(U[1]), ylab = expression(U[2]), zlab="Density",main=" "
)
dev.off()
see ?trellis.device and ?Devices
You will have to modify the arguments to suit the graph and your device.
If you use a device like png use a high resolution value
HTH
Duncan
Duncan Mackay
Department of Agronomy and Soil Science
University of New England
Armidale NSW 2351
Email: home: mac...@northnet.com.au
At 07:55 6/07/2013, you wrote:
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.jo...@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]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
