Dear Mr. Ligges,
my apologies! As I can see in your E-Mail the quote of mine is really just
unreadable! Anyhow, this was not produced by me writing messy messages, this
is due to the Java application behind the client not being able to cope with
external tab-characters on my browser when parsing the message while sending(I
tipped the message on an external editor because of the code lines and copy
pasted it). Of course I also never would add any html files!
Anyhow, I send you my original message again (new = True did not help me.
Assuming I used it the way you suggested it produced the same error). I send
it to myself first, so I am sure, this mess won't happen again!
My previous Message:
#########
Dear Madame or Sir,
I am having a problem in combining density-smoothed scatterplot matrices with
a plot of kernel destiny estimations of each dimension plotted on the
respective diagonals.
I have tried following approach using the package "sm" for the kernel density
estimation, as well as "MASS" respectively:
pairs
(myTable[, 1:4],
panel=function(x,y, ...)
{
xy=kde2d(x,y)
image(xy, add=TRUE)
}
,diag.panel = function(x,...)
{
plot(density(x,...), add = TRUE)
#boxplot(x, add = TRUE)
}
)
where myTable is an ordinary table I am using the first four columns of.
Unfortunately this produces an error stating that plot has created a new plot.
I suppose this is the same situation when forgetting "add = TRUE" with
plotting functions such as boxplot (as in my code used with "boxplot(x, add =
TRUE)" ) or any other (e.g. smoothScatter). Analyzing the partial plot, which
can be seen before R rejects my code, I come to the conclusion that add might
not be a valid argument for "plot".
My questions therefore is, whether there is an option to plot the result of
the 1d kernel density estimation to the diagonal of the scatterplot produced
by pairs or which plot function is valid in this case. So basically I would
like to generate the same output I am able to produce with boxplots on the
diagonal with a density estimation.
Additional question: I bet there is an opportunity to include both the
boxplot and the kernel density estimation on each field of the diagonal (as I
saw before in a published paper), but I just could not figure out how to
"overdraw". Maybe you could give me a slight idea (reference would be fine,
too) on how to achieve this. I would be highly thankful!
Kind regards and many thanks in advance,
-Patrick Grossman
##############
Please excuse the inconvenience and note my thanks.
Warmest regards,
-Patrick Grossman
----- Ursprüngliche Nachricht -----
Von: Uwe Ligges <lig...@statistik.tu-dortmund.de>
Datum: Samstag, 2. Juli 2011, 16:36
Betreff: Re: [R] Error when using plot in diag.panel argument of pairs
An: Patrick Großmann <patrick.grossm...@uni-bielefeld.de>
Cc: r-help@r-project.org
> If you send a better readable message with really reproducible
> code, we
> may provide a solution, this way it is really hard to read your
> code.
> And yes, please do not send html mail. My guess is that
> par(new=TRUE)
> will help, but I do not want to construct the examples at first.
>
>
>
>
> On 02.07.2011 14:37, "Patrick Großmann" wrote:
> > Dear Madame or Sir,I am having a problem in
> combining density-smoothed scatterplot matrices with a plot of
> kernel destiny estimations of each dimension plotted on the
> respective field of the diagonal.I have tried following approach
> using the package "sm" for the kernel density estimation, as
> well as "MASS" respectively:pairs(myTable[,
> 1:4],panel=function(x,y, ...){
> xy=kde2d(x,y) image(xy,
> add=TRUE)},diag.panel = function(x,...){plot(density(x), add =
> TRUE)#boxplot(x, add = TRUE)})where myTable is an ordinary table
> I am using the first four columns of. Unfortunately this
> produces an error stating that plot has created a new plot. I
> suppose this is the same situation when forgetting "add =
> TRUE" with plotting functions such as boxplot (as in my code
> used with "boxplot(x, add = TRUE)" ) or any other (e.g.
> smoothScatter). Analyzing the partial plot, which can be
> seen before R rejects my code, I come to the conclusion that add
> might not be a valid argument for "plot
> ". My questions therefore is, whether there is an option to plot
> the result of the 1d kernel density estimation to the diagonal
> of the scatterplot produced by pairs or which plot function is
> valid in this case. So basically I would like to generate the
> same output I am able to produce with boxplots on the diagonal
> with a density estimation.Additional question: I bet there is an
> opportunity to include both the boxplot and the kernel density
> estimation on each field of the diagonal (as I saw before in a
> published paper), but I just could not figure out how to
> "overdraw". Maybe you could give me a slight idea (reference
> would be fine, too) on how to achieve this. I would be highly
> thankful!Kind regards and many thanks in advance,-Patrick
> Grossman --
> > Patrick Großmann - HiWi zur Netzwerkdokumentation
> > Network Operation Center NOC
> > Hochschulrechenzentrum HRZ
> Universität Bielefeld
> >
> > Büro:
> > UHG V0-251
> > Telefon:
> > +49 521 106-12618
> > Fax:
> > +49 521 106-2969
> > Telefon Sekretariat:
> > +49 521 106-4951
> >
> > http://www.uni-bielefeld.de/hrz/
> >
> >
> >
> > [[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.
--
Patrick Großmann - HiWi zur Netzwerkdokumentation
Network Operation Center NOC
Hochschulrechenzentrum HRZ Universität Bielefeld
Büro:
UHG V0-251
Telefon:
+49 521 106-12618
Fax:
+49 521 106-2969
Telefon Sekretariat:
+49 521 106-4951
http://www.uni-bielefeld.de/hrz/
[[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.
