Duncan,
First off, I admit it is not clear to me what you are trying to
achieve and more importantly, why? by "why" I mean 1) I don't see
the
advantage of writing one general panel function for completely
different situations (one/multiple smoothers, grouping levels
etc.) 2)
your intended result as I understand it seems rather cluttered,
google
<chartjunk>. 3) I am unfamiliar with locfit package, but are we
reinventing the wheel here ? i.e. will modifying settings in
xyplot(y
~x, xx, groups = Farm, type=c('p','smooth')) achieve the same ?
With your initial reproducible example (thank you) it was easy to
eliminate the errors, but clearly the resulting plots are not
what you
intended (continue inline):
On Thu, Apr 19, 2012 at 4:23 PM, Duncan Mackay <mac...@northnet.com.au >
wrote:
<snip>
3. What I want to be able to add in the above is extra lines with
different
values of nn.
I think I will have to modify panel.Locfit so that it goes
through
different values of nn in each of the panels and groups if I want
different
colours for extra lines with different nn values
Yes you could. There are several options:
add group.number to the arguments of panel.locfit and use it to
make
nn a vector, along the lines of
panel.foo <- function(x,y,group.number,theta,...){
smpar <- theta[group.number]
panel.loess(x,y,smpar,...)
panel.xyplot(x,y,...)
}
xyplot (y ~ x ,xx ,group
=Farm,theta=c(4,1,.4),panel=panel.superpose,panel.groups=panel.foo)
# or
xyplot(y~x|Farm,xx,group=Padd,theta=c(.6,1),
panel=panel.superpose,panel.groups=panel.foo)
Here you will need to modify the Farm group to 6 levels - 3*two
smoothers.
You could make nn a list and loop over it inside the panel
function.
Looks like you tried something like that with specifying 2
panel.Locfit, one suggestion to your code:
panel.Locfit(x,y,...) # default 0.7
panel.Locfit(x,y,nn=0.9) # i.e. remove the
... to avoid clashes
Finally, use ?trellis.focus to plot the second smoother "post-hoc".
also the latticeExtra package has many useful tools to create
layers
of the same (or different) plot with different settings.
4 Produce an extra line for a fit for all the groups in 1/2+
panels.
As for 3 but I do not know how to group all the x and y's for
each
of the
panes using panel.groups
Why does it matter ? seems you have failed to learn the lesson from
the first post - the same functionality applies to 1 as to multiple
panels. Does each panel have a different grouping structure ? use
packet.number() for panels similar to group.number idea.
I need to do this and then scale up for a panel function to
include
confidence bands
than expand the xlim,ylim or scales in ?xyplot
For the record making Farm and Padd factors. With 1 panel and
groups =
Farm
works with the extra line the same colour for its group
a similar situation for the three panels when conditioned by
Farm and
groups
= Pad
????
Like I said I am a little lost on this problem but I hope this
helps
giving some direction.
Cheers
xyplot(y ~x, xx,
groups = Farm,
par.settings = list(strip.background = list(col =
"transparent"),
superpose.line = list(col =
c("black","grey"),
lwd =
c(1,2,3),
lty =
c(2,1,3)),
superpose.symbol = list(cex = c(0.8,
0.7,0.7),
col =
c("red","black","blue"),
pch =
c(20,4,16))
),
auto.key=list(lines=T,points = T,rectangles=F),
panel = panel.superpose,
panel.groups=function(x,y, ...){
panel.xyplot(x,y,...)
panel.Locfit(x,y,...) # default 0.7
panel.Locfit(x,y,nn=0.9,...)
}
) ## xyplot
Regards
Duncan
At 02:12 20/04/2012, you wrote:
On Thu, Apr 19, 2012 at 2:30 AM, Duncan Mackay
<mac...@northnet.com.au>
wrote:
Hi
xyplot(y ~x|Farm,xx,
groups = Padd,
panel = panel.superpose,
panel.groups=function(x,y, ...){
panel.Locfit(x,y,...)
panel.xyplot(x,y,...)
}
) ## xyplot
The above works nicely and also without par.setting giving
lattice
defaults.
The par.setting is handy for a lot of graphs that I do.
But when I tried a 1 panel plot I get the error message.
xyplot(y ~x,xx,
groups = Farm,
auto.key=TRUE,
panel = function(x,y, ...){
panel.Locfit(x,y,...)
panel.xyplot(x,y,...)
}
)
These two plots are NOT THE SAME. Did you want the same as the
first
but with groups being Farm and Padd ignored ? in that case you
(again)
need a panel.groups:
xyplot(y ~x,xx,
groups = Farm,
auto.key=TRUE,
panel = panel.superpose,panel.groups=function(x,y,...){
panel.Locfit(x,y,...)
panel.xyplot(x,y,...)
}
)
If I want to plot another curve with different smoothing
but gives an error message without par.settings if i want to add
panel.Locfit(x,y,nn= 0.9,lwd =
c(1,2,3), ...)
Error using packet 1
formal argument "Iwd" matched by multiple actual arguments
It is all in the way you initially specified how to pass the
arguments
for panel.Locfit. This works without error:
xyplot(y ~x,xx,
groups = Farm,
auto.key=TRUE,lwd=1:3,
panel = panel.superpose,panel.groups=function(x,y,nn,...){
panel.Locfit(x,y,nn=.9,...)
panel.xyplot(x,y,...)
}
)
HTH
I also need to plot a smoothed line for all groups trying
groups,
subscripts
and panel.groups as arguments without success
Any solutions to solve the above will be gratefully received and
faithfully
applied.
Duncan
sessionInfo()
R version 2.15.0 (2012-03-30)
Platform: i386-pc-mingw32/i386 (32-bit)
locale:
[1] LC_COLLATE=English_Australia.1252
LC_CTYPE=English_Australia.1252
LC_MONETARY=English_Australia.1252 LC_NUMERIC=C
LC_TIME=English_Australia.1252
attached base packages:
[1] datasets utils stats graphics grDevices grid
methods
base
other attached packages:
[1] locfit_1.5-7 R.oo_1.9.3 R.methodsS3_1.2.2
foreign_0.8-49
chron_2.3-42 MASS_7.3-17 latticeExtra_0.6-19
RColorBrewer_1.0-5
[9] lattice_0.20-6
loaded via a namespace (and not attached):
[1] tools_2.15.0
