data911a <- read.csv("Table911a.csv")
data911a
attach(data911a)

cor(hour,score) # 時間と成績の相関

cor(data911a[2:4])# ３変数間の相関行列

​

​

#スキルの影響を除いた場合の時間と成績の偏相関

#パッケージを使わずに計算する場合（公式に代入）

rxy=-0.511 
rxz=-0.716
ryz=0.755
rxy.z=(rxy-rxz*ryz)/(sqrt(1-rxz^2)*sqrt(1-ryz^2))

rxy.z
 

​

#psych パッケージを使用する場合
install.packages("psych")
library("psych")

partial.r(data911a[1:4],c(2,3),4)

​

​