การหาค่า p-pearson ในสหสัมพันธ์ p


26

เป็นไปได้หรือไม่ที่จะหาค่า p-pearson correlation ใน R?

เพื่อค้นหาความสัมพันธ์ลูกแพร์สันฉันมักจะทำเช่นนี้

col1 = c(1,2,3,4)
col2 = c(1,4,3,5)
cor(col1,col2)
# [1] 0.8315218

แต่ฉันจะหาค่า p ของสิ่งนี้ได้อย่างไร


1
ความช่วยเหลือในcor( ?cor) ระบุอย่างชัดเจนcor.test (ภายใต้ "ดูเพิ่มเติม")
Glen_b

คำตอบ:


32

คุณสามารถใช้cor.test:

col1 = c(1,2,3,4) 
col2 = c(1,4,3,5)
cor.test(col1,col2) 

ซึ่งจะช่วยให้ :

# Pearson's product-moment correlation   
# data:  col1 and col2   
# t = 2.117, df = 2, p-value = 0.1685   
# alternative hypothesis: true correlation is not equal to 0   
# 95 percent confidence interval:   
#  -0.6451325  0.9963561   
# sample estimates:   
#       cor    
# 0.8315218    

ข้อมูลเพิ่มเติมเกี่ยวกับสถิติและพารามิเตอร์เพิ่มเติมได้ที่หน้าอย่างเป็นทางการ:

https://stat.ethz.ch/R-manual/R-patched/library/stats/html/cor.test.html


ขอบคุณ, คะแนนไม่เพียงพอที่จะ upvote, ใช้มันแม้ว่า
tubby


7

สิ่งต่อไปนี้จะทำตามที่คุณถาม:

 library(Hmisc) # You need to download it first.
 rcorr(x, type="pearson") # type can be pearson or spearman

นี่คือกรอบข้อมูล x และ rcorr คืนค่าสหสัมพันธ์ซึ่งเป็นไปได้ที่จะสร้างจากกรอบข้อมูล "x"

หรือคุณสามารถคำนวณสถิติด้วยตัวคุณเอง:

t=ρ^1ρ^2n2

โดยที่มีความสัมพันธ์แบบลูกแพร์ที่ประเมินจากข้อมูลและ n คือขนาดตัวอย่างρ^


ขอบคุณ แต่ x คืออะไร ฉันคิดว่ามันคือการต่อกันของ col1 และ col2 เพราะเราต้องการเวกเตอร์สองตัวเพื่อคำนวณความสัมพันธ์ของเพียร์สัน แต่คุณสามารถบอกฉันว่า x คืออะไร?
tubby

มันเป็นกรอบข้อมูลให้ดูการปรับปรุงของฉัน
Repmat
โดยการใช้ไซต์ของเรา หมายความว่าคุณได้อ่านและทำความเข้าใจนโยบายคุกกี้และนโยบายความเป็นส่วนตัวของเราแล้ว
Licensed under cc by-sa 3.0 with attribution required.