Re: [R] Error in cor.default(x1, x2) : missing observations in cov/cor

From: Erik Iverson <iverson_at_biostat.wisc.edu>
Date: Thu, 28 Feb 2008 10:12:08 -0600

So you did have NAs in your data?

what do

 > any(is.na(x1))

and

 > any(is.na(x2))

give?

Ken Spriggs wrote:

> Thanks.  Your suggestion works. 
> 
> 
> Daniel Malter wrote:

>> Does the code below solve your problem? If you have NAs in the same rows,
>> you have to use "c" or "p" as use= parameters. Otherwise you get the error
>> you described.
>>
>> a=c(1,2,3,4,NA,6)
>> b=c(2,4,3,5,NA,7)
>>
>> which(is.na(a))==which(is.na(b))
>>
>> cor(a,b) ####Error
>> cor(a,b,use="all.obs") ####Error
>>
>> cor(a,b,use="complete.obs") #### Does it. AND
>> cor(a,b,use="pairwise.complete.obs") #### Does it too.
>>
>> -------------------------
>> 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 Ken Spriggs
>> Gesendet: Wednesday, February 27, 2008 4:34 PM
>> An: r-help_at_r-project.org
>> Betreff: [R] Error in cor.default(x1, x2) : missing observations in
>> cov/cor
>>
>>
>> Hello,
>>
>> I'm trying to do cor(x1,x2) and I get the following error:
>> Error in cor.default(x1, x2) : missing observations in cov/cor
>>
>> A few things:
>> 1. I've used cor() many times and have never encountered this error.
>> 2. length(x1) = length(x2)
>> 3. is.numeric(x1) = is.numeric(x2) = TRUE
>> 4. which(is.na(x1)) = which(is.na(x2)) = integer(0) {the same goes for
>> is.nan()}
>> 5. I also try cor(x1,x2, use = "all.obs") and get the same error.
>>
>> What can be going wrong?
>>
>>
>> --
>> View this message in context:
>> http://www.nabble.com/Error-in-cor.default%28x1%2C-x2%29-%3A-missing-observa
>> tions-in-cov-cor-tp15723848p15723848.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.
>>
>> ______________________________________________
>> 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.
>>
>>
>

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 Thu 28 Feb 2008 - 16:15:03 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 Thu 28 Feb 2008 - 16:30:17 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.

list of date sections of archive