data413 <- read.csv("Table413.csv")
data413
attach(data413)

​

#散布図

plot(X,Y)
 

#回帰分析　関数 lm(変数１~変数２,データ)

# 変数１を変数２から予測する（順序に注意）

result <- lm(Y ~ X, data = data413)           
summary(result) 

​

#出力の中で以下が切片aと傾きb

#Coefficients:

#　              Estimate

#(Intercept)  -0.9435    

#   X               1.2536 

 

#lm(Y ~ X, data = data413)  だけを>から実行すると、

#Estimateのみが表示される

detach(data413)

#スピアマンの順位相関係数の無相関検定
data421 <- read.csv("Table421.csv")
data421
attach(data421)

cor.test(X, Y, method="spearman") 
detach(data421)

#ケンドールの順位相関係数の無相関検定
data422 <- read.csv("Table422.csv")
data422
attach(data422)

cor.test(X, Y, method="kendall") 
detach(data422)