研究を加速するツール Exploratory
Rベースで基本機能は学生は無料
Exploratory勝手に普及委員会 委員長の
夏一番です。
Exploratoryはサンフランシスコ発
アメリカ留学中に作成したスクリプトを公開します
今回は相関係数を返すスクリプト
f.cor<-function(df, x,y,digits=3,method="pearson"){
## MakeP function
makep<-function(x){
##P >= 0.1
pvalue<-ifelse(x>=0.1,as.character(round(x,digits=2)),
## P >= 0.01
ifelse(x<0.1 & x>=0.01, as.character(round(x,digits=3)),
## 0.001 - P - 0.01
ifelse(x<0.01 & x>=0.001,as.character(round(x,digits=3)),
## P <0.001
ifelse(x<0.001,"< 0.001",NA))))
return(pvalue)}
factor1<-as.character(substitute(x))
factor2<-as.character(substitute(y))
library(broom)
library(dplyr)
res<-cor.test(x,y,method=method,data=df)
res<-res %>% tidy %>%
mutate_at(vars(estimate,conf.low,conf.high),funs(round(.,digits)))%>%
mutate(P.value=(makep(p.value)))%>%
select(estimate,P.value,method,alternative)%>%
cbind(factor1,factor2,.)
res
}
これで相関係数とP valueが返されます