2011-09-26 40 views

Trả lời

13

Hãy thử cor(A,B, use = "pairwise.complete.obs"). Điều đó sẽ bỏ qua các NA trong các quan sát của bạn.

Để thống kê nghiêm ngặt, bạn cũng nên xem xét số lượng các mục bị thiếu trong dữ liệu của mình và xem liệu thiếu sót ở mức giả định ngẫu nhiên.

Chỉnh sửa 1: Hãy xem ?cor để xem các tùy chọn khác cho thông số use.

4

Bạn có thể xem xét sử dụng hàm rcorr trong gói Hmisc.

Tốc độ rất nhanh và chỉ bao gồm các quan sát hoàn chỉnh theo cặp. Các đối tượng quay trở lại chứa một ma trận

  1. điểm số tương quan
  2. với số lượng quan sát sử dụng cho mỗi giá trị tương quan
  3. của một p-giá trị cho mỗi tương quan

Một số mã ví dụ có sẵn here:

Các vấn đề liên quan