เป็นไปได้หรือไม่ที่จะหาค่า p-pearson correlation ใน R?
เพื่อค้นหาความสัมพันธ์ลูกแพร์สันฉันมักจะทำเช่นนี้
col1 = c(1,2,3,4)
col2 = c(1,4,3,5)
cor(col1,col2)
# [1] 0.8315218
แต่ฉันจะหาค่า p ของสิ่งนี้ได้อย่างไร
เป็นไปได้หรือไม่ที่จะหาค่า p-pearson correlation ใน R?
เพื่อค้นหาความสัมพันธ์ลูกแพร์สันฉันมักจะทำเช่นนี้
col1 = c(1,2,3,4)
col2 = c(1,4,3,5)
cor(col1,col2)
# [1] 0.8315218
แต่ฉันจะหาค่า p ของสิ่งนี้ได้อย่างไร
คำตอบ:
คุณสามารถใช้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
หากคุณต้องการเฉพาะค่า P:
> cor.test(col1,col2)$p.value
[1] 0.1684782
สิ่งต่อไปนี้จะทำตามที่คุณถาม:
library(Hmisc) # You need to download it first.
rcorr(x, type="pearson") # type can be pearson or spearman
นี่คือกรอบข้อมูล x และ rcorr คืนค่าสหสัมพันธ์ซึ่งเป็นไปได้ที่จะสร้างจากกรอบข้อมูล "x"
หรือคุณสามารถคำนวณสถิติด้วยตัวคุณเอง:
โดยที่มีความสัมพันธ์แบบลูกแพร์ที่ประเมินจากข้อมูลและ n คือขนาดตัวอย่าง
cor
(?cor
) ระบุอย่างชัดเจนcor.test
(ภายใต้ "ดูเพิ่มเติม")