Dear ashz,
Unfortunately, much of you want is not possible with the current
implementation of pairs.panels.
Since pairs.panels is adapted from the help file of pairs, you might
try pairs and then add in the panel functions that do what you want.
That the lm option does it what it does met a need of mine for a demo
of the difference of regression slopes of X on Y versus Y on X. Your
request is very reasonable and I will implement it in the next
revision.
To get the pairwise complete rather than pairwise, you can preprocess
your data file with na.omit
e.g., cfc <- na.omit(cfcap[8:11])
pairs.panesl(cfc)
At 1:37 PM -0700 10/3/10, ashz wrote:
Hi,
I used the pairs.panels() in pkg:psych and it is helpful. It saves time.
but if I use this line:
pairs.panels(cfcap[8:11], scale = FALSE, lm=TRUE,ellipses=TRUE, digits = 2
)
The results are:
- The upper.panel does not show the pearson r but the lm data. Furthermore,
can I use the pairwise.complete.obs method for the upper.panel. Can it be
fixed?
- Can I remove the histograms?
Not as a call, but you change pairs.panels to draw just the densities
by substituting
"panel.hist.density" <-
function(x,...) {
usr <- par("usr"); on.exit(par(usr))
par(usr = c(usr[1:2], 0, 1.5) )
h <- hist(x, plot = FALSE)
breaks <- h$breaks; nB <- length(breaks)
y <- h$counts; y <- y/max(y)
#rect(breaks[-nB], 0, breaks[-1], y,col=hist.col) # <---
comment this line out
tryd <- try( d <- density(x,na.rm=TRUE,bw="nrd",adjust=1.2),silent=TRUE)
if(class(tryd) != "try-error") {
d$y <- d$y/max(d$y)
lines(d)}
}
in place of the current panel.hist.density function.
- Can I control the eliipse alpha?
Not yet. Good idea.
Your requests are all very reasonable and will be added to my wish
list of additions to pairs.panels. This will not happen for several
weeks, however.
Bill
Thanks a lot.
--
View this message in context:
http://r.789695.n4.nabble.com/Scatterplot-matrix-Pearson-linear-correlation-and-Density-Ellipse-tp2763552p2953521.html
Sent from the R help mailing list archive at Nabble.com.
______________________________________________
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.
