Foro de discusión y soporte para usuarios de R

Hola a todos,
Tengo un script para construir una matriz de "scatter plots", en donde se indican también los valores del coeficiente de correlación de spearman y respectivos valores de p para cada par de variables (además de otros detales como la función de densidad de cada variable y smoothers) .
Quisiera añadir a este script un comando que me permitiera incluir en la grafica final, además de los valores de rho y de p, el numero de observaciones usadas en cada correlación. ¿Alguién tiene alguna idea de como lo podria hacer?

Muchas gracias!

Dulce

Aqui tenéis el script:

#scatter plot matrix with density function in diagonal and smoother

panel.cor <- function(x, y,method="spearman",digits=2,...) {
points(x,y,type="n");
usr <- par("usr"); on.exit(par(usr))
par(usr = c(0, 1, 0, 1));
correl <- cor.test(x, y,method=method);
r=correl$estimate;
pval=correl$p.value;
color="black";
if (pval<0.05) color="red";

txt <- format(r,digits=2)
pval <- format(pval,digits=2)
txt <- paste("rho=", txt, "\np=",pval,sep="")
text(0.5, 0.5, txt,col=color)
}

scatterplot.matrix (~ variable1 + variable2 + variable3 ,
data=mydata, main="Mydata", smooth=TRUE, lower.panel=panel.cor, pch=20, cex=0.5, col=c("red","black"),cex.labels=1, font.labels=2, lwd=0.5)