R help - Nov 2005 - obtaining p values & CI for a non-parametric correlation matrix

Hello,

I tried to obtain a nonparametric correlation matrix with the following syntax:

by(mydat, mydat$GROUP, function(subset) { matrix <- cbind(subset$WK1FREQ, 
subset$WK2FREQ, subset$WK3FREQ, subset$WK4FREQ, subset$BESTDAY) cor(matrix, 
method="kendall",use="pairwise.complete.obs") })
Error: syntax error

In addition to the matrix, I was hoping to obtain p-value and CI (rounded 
to 3 decimal points).

Any suggestions are much appreciated,

Bob Green