Hi Rolf, Hi Mark, Hi List,
I have not digested Rolf's response yet. It may well answer my problems.
In the meantime I have some reproducible code which actually shows my
problem:
patches <-
structure(list(URBAN_AREA = structure(c(2L, 19L, 23L, 2L, 19L,
23L, 2L, 19L, 23L, 2L, 19L, 23L), .Label = c("CENTRAL AUCKLAND ZONE",
"CHRISTCHURCH", "DUNEDIN", "HAMILTON ZONE", "HASTINGS ZONE",
"INVERCARGILL", "LOWER HUTT ZONE", "mean", "NAPIER ZONE", "NELSON",
"NEW PLYMOUTH", "NORTHERN AUCKLAND ZONE", "PALMERSTON NORTH",
"PORIRUA ZONE", "ROTORUA", "SD", "SE", "SOUTHERN AUCKLAND ZONE",
"TAURANGA", "WANGANUI", "WELLINGTON ZONE", "WESTERN AUCKLAND ZONE",
"WHANGAREI"), class = "factor"), NO_PATCHES = c(11L, 16L, 21L,
87L, 192L, 324L, 164L, 417L, 773L, 679L, 757L, 3083L), MEAN_AREA =
c(9.623631225,
15.29089619, 149.2063532, 14.1676, 247.5262, 28.611, 11.5698,
221.0022, 37.3725, 11.918, 133.5804, 25.6759), AREA.ZONE = c(13683,
3666, 1558, 64830, 41103, 22581, 123819, 90107, 57627, 264735,
223963, 174456), Buffer.zone = c(0L, 0L, 0L, 5L, 5L, 5L, 10L,
10L, 10L, 20L, 20L, 20L)), .Names = c("URBAN_AREA", "NO_PATCHES",
"MEAN_AREA", "AREA.ZONE", "Buffer.zone"), class = "data.frame",
row.names = c(2L,
15L, 19L, 22L, 36L, 40L, 42L, 56L, 60L, 62L, 76L, 80L))
library(lattice)
Region = factor(patches$URBAN_AREA)
lpatches = patches
for (i in 2:4) lpatches[,i] = log10(patches[,i])
# zone not transformed
lpatches.pca = princomp(lpatches[,-c(1,17)],cor=FALSE)
x = lpatches.pca$scores[,1]
y = lpatches.pca$scores[,2]
zz = as.character(patches$Buffer.zone/5)
table(zz)
plsy <- trellis.par.get("plot.symbol")
# only 0 or 1 used as plotting symbol
plsy$pch = as.character(rep(1:6,2))
trellis.par.set("plot.symbol",plsy)
xyplot(y ~ x |Region)
# only 1,2,3,4 used as plotting symbol
I actually wish 0,1,2, or 4 to be used - to indicate the zone coded in zz.
Cheers, Murray
PS The xyplots produced on R 2.7.0 for Mac OS X, the first one also on
an older Windows version.
Rolf Turner wrote:
Murray:
I'm not at all sure that I understand what you're driving at --- but
does this do something like what you want?
require(lattice)
set.seed(260808)
n = 50
x = rnorm(n)
y = rnorm(n)
z = ceiling(runif(n,0,4))
g = runif(n,0,6)
G = factor(ceiling(g))
print(xyplot(y ~ x | G,pch=G,
panel=function(x,y,...,subscripts,pch) {
panel.xyplot(x,y,pch=pch[subscripts])
}
))
cheers,
Rolf
######################################################################
Attention: This e-mail message is privileged and confidential. If you
are not the intended recipient please delete the message and notify
the sender. Any views or opinions presented are solely those of the
author.
This e-mail has been scanned and cleared by MailMarshal
www.marshalsoftware.com
######################################################################
Dr Murray Jorgensen http://www.stats.waikato.ac.nz/Staff/maj.html Department of Statistics, University of Waikato, Hamilton, New Zealand Email: maj at waikato.ac.nz majorgensen at ihug.co.nz Fax 7 838 4155 Phone +64 7 838 4773 wk Home +64 7 825 0441 Mobile 021 139 5862