rcoder
2008-Aug-15 13:44 UTC
[R] cor() btwn columns in two matrices - no complete element pairs
Hi everyone, I'm trying to calculate correlation coefficients between corresponding columns in two matrices with identical dimensions but different data. The problem is that the matrices contain NAs in different locations. I am using the following code to try to calculate correlations between complete sets of data: #Code start maxcol<-ncol(mat1) for (i in 1:maxcol) { corr_results[1,i]<-cor(mat1[,i],mat2[,i], use="complete.obs") } #Code end ...but I get the following error message: Error in cor(mat1[,i], mat2[,i], use="complete.obs") : no complete element pairs Is there something I'm not including in the 'cor' parentheses? I apologise for not including the true original data frames. Thanks, rcoder -- View this message in context: http://www.nabble.com/cor%28%29-btwn-columns-in-two-matrices---no-complete-element-pairs-tp18998875p18998875.html Sent from the R help mailing list archive at Nabble.com.
rcoder
2008-Aug-15 14:13 UTC
[R] cor() btwn columns in two matrices - no complete element pairs
Is there any way to increment the loop when the error occurs, returning an NA in the results matrix. I have included an if starement to only calculate cor() when columns are !all(is.na()), but I still exit out with this error, albeit after some loops through the matrix - i.e. when the condition is satisfied, but the data rows do not overlap at all. rcoder wrote:> > Hi everyone, > > I'm trying to calculate correlation coefficients between corresponding > columns in two matrices with identical dimensions but different data. The > problem is that the matrices contain NAs in different locations. I am > using the following code to try to calculate correlations between complete > sets of data: > > #Code start > maxcol<-ncol(mat1) > for (i in 1:maxcol) > { > corr_results[1,i]<-cor(mat1[,i],mat2[,i], use="complete.obs") > } > #Code end > > ...but I get the following error message: > > Error in cor(mat1[,i], mat2[,i], use="complete.obs") : > no complete element pairs > > Is there something I'm not including in the 'cor' parentheses? I apologise > for not including the true original data frames. > > Thanks, > > rcoder >-- View this message in context: http://www.nabble.com/cor%28%29-btwn-columns-in-two-matrices---no-complete-element-pairs-tp18998875p18999356.html Sent from the R help mailing list archive at Nabble.com.
Daniel Malter
2008-Aug-16 18:02 UTC
[R] cor() btwn columns in two matrices - no complete element pairs
Hi, the argument use should be: use="p" (for pairwise complete obs.) Daniel ------------------------- cuncta stricte discussurus ------------------------- -----Urspr?ngliche Nachricht----- Von: r-help-bounces at r-project.org [mailto:r-help-bounces at r-project.org] Im Auftrag von rcoder Gesendet: Friday, August 15, 2008 9:45 AM An: r-help at r-project.org Betreff: [R] cor() btwn columns in two matrices - no complete element pairs Hi everyone, I'm trying to calculate correlation coefficients between corresponding columns in two matrices with identical dimensions but different data. The problem is that the matrices contain NAs in different locations. I am using the following code to try to calculate correlations between complete sets of data: #Code start maxcol<-ncol(mat1) for (i in 1:maxcol) { corr_results[1,i]<-cor(mat1[,i],mat2[,i], use="complete.obs") } #Code end ...but I get the following error message: Error in cor(mat1[,i], mat2[,i], use="complete.obs") : no complete element pairs Is there something I'm not including in the 'cor' parentheses? I apologise for not including the true original data frames. Thanks, rcoder -- View this message in context: http://www.nabble.com/cor%28%29-btwn-columns-in-two-matrices---no-complete-e lement-pairs-tp18998875p18998875.html Sent from the R help mailing list archive at Nabble.com. ______________________________________________ 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.