# Re: [R] ellipse in pairs plot

From: Kaveh Vakili <kaveh.vakili_at_ulb.ac.be>
Date: Sat, 15 Nov 2008 00:28:42 +0100 (CET)

Dear list moderator,

I'm sending you a solution i found to the question i asked the list yesterday (i have not been able to figure out how to post it to the list directly).

elly<-function(w,col=array(1,dim(w)),pch=array(1,dim(w))){

```	p<-dim(w)
op<-par(mfrow=c(3,3),mar=c(0,0,0,0)+1,oma=c(1,1,1,1))
c1<-rep(1:p,p)
c2<-sort(c1)
dan<-covMcd(w)
for (i in 1:p^2) {
if(c1[i]!=c2[i]){
plot(w[,c(c1[i],c2[i])],xlab="",ylab="",xaxt="n",yaxt="n",main="",col=col,pch=pch)
grid()
ava<-c(dan\$cov[c1[i],c1[i]],dan\$cov[c1[i],c2[i]],dan\$cov[c1[i],c2[i]],dan\$cov[c2[i],c2[i]])
ellipse(mu=dan\$center[c(c1[i],c2[i])],sigma=matrix(ava,2,2),col=max(col)+1)
} else {
plot(density(w[,c2[i]],from=min(w[,c2[i]]),to=max(w[,c2[i]])),xlab="",ylab="",xaxt="n",yaxt="n",main="",xlim=range(w[,c2[i]]),type="h",col="light grey")
}
if(c1[i]!=c2[i]){
if(c1[i]==1){
axis(2)
}
if(c2[i]==1){
axis(3)
}
if(c1[i]==p){
axis(4)
}
if(c2[i]==p){
axis(1)
}
}
if(i==1){
axis(3)
}
if(i==p^2){
axis(1)
}
}
par(op)
```

}

>Hi dear list,
>
>I'm scrambling with this:
>
>i have a p-dimensional data frame for which i want to have a pairs() plot, with ellipse added in each lower.panel.
>
>I want the ellipse's parameter (mu, sigma) to be extracted from the appropriate entries in a pre-computed matrix ("dab" in the code below).
>
>Can it be done ?
>
>
>
>melly<-function (x, y, col = par("col"), bg = NA, pch = par("pch"),cex = 1){
> points(x, y, pch = pch, col = col, bg = bg, cex = cex)
> ok <- is.finite(x) & is.finite(y)
> if (any(ok))
> ellipse(mu=dab\$center[?],sigma=dab\$cov[?,?])
>}
>
>pairs(wf[,1:5],diag.panel=bagwag,col=wc,pch=wp,cex=.5,lower.panel=melly)
>
>
> no attachments have been sent
>
>

R-help_at_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. Received on Sat 15 Nov 2008 - 00:26:23 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Sat 15 Nov 2008 - 01:30:25 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.