mycor <- function (x,...)                
#
# Ausgabe einer Korrelatiosnmatrix sowie einer Matrix der dazugehörigen p-Werte
#
{nvar <- dim(x)[2]     # Anzahl Variablen
 corr<-cor(x,...)      # Korrelationsmatrix
 pval<-corr            # provisorische Matrix der p-Werte
 for (i in 1:nvar)     # Zeilen
  { for (j in 1:nvar)  # Spalten
    { pval[i,j]<-cor.test(x[,i],x[,j],...)$p.value }
   }
  list(c=corr,p=pval)
}