R 복수 데이터의 총당 단일 회귀 분석
4895 단어 R
0 소개
꽤 초보적인 내용입니다. 비망록도 겸해.
R에 있어서 복수 요인간의 단회귀 분석을 정리해 실시할 때의 흐름을 정리합니다.
이번은 심리학 분야에서 이용되기 위해서 정리된 psych의 corr.test 함수를 이용했습니다.
상관 계수만 알면 좋은 경우에는 cor 함수를 사용하면 됩니다.
1 이전 준비
사용하는 패키지와 데이터 세트를 R로 읽습니다.
준비하는 데이터 세트는 아래 이미지와 같은 형식입니다.
각 열에는 외부 요인이 있으며 각 행에는 샘플 데이터가 들어 있습니다.
(이번 사용한 것은 199 샘플×22 요인)
setwd("/Directory1/Directory2") #ディレクトリの指定
install.packages("psych") #psychパッケージのインストール
library("psych")#パッケージの呼び出し
dat1 <- read.csv("for_cor.csv", header=TRUE,row.names=1,fileEncoding="CP932")
# データセットの読み込み
# RはEXCELで通常保存される[.xlsx]の拡張子ファイルを扱うことは難しいので[.csv]形式でのデータセットを用いる
#列名、行名が不要の場合にはheader=FALSE, row.names <- NULL
2 상관 분석 수행
디폴트에서는 "pearson"의 상관 계수가 출력됩니다만 "spearman"또는 "kendall"에 의한 상관 계수의 산출도 가능합니다.
기타 adjust 소개 여러 파라미터를 변경할 수 있습니다.
자세한 내용은 [4 인용]에서.
result1 <- corr.test(dat1) #dat1についての相関分析結果をresult1に格納
result1 #結果表示(必要なら)
3 데이터 출력
각 항목에 대해 별도의 파일로 출력하고 싶었으므로 다음을 수행했습니다.
names(result1) #相関分析により得られたデータ項目の確認(必要なら)
cor1 <- result1$r #result1のr(相関係数)に関するデータのみを取り出しcor1に格納
write.csv(cor1, file = "cor1.csv") #相関係数のデータをcsvファイルに出力
corr.test(cbind(dat4))$p
pv1 <- result1$p #result1のp(p-value)に関するデータのみを取り出しpv1に格納
write.csv(pv1, file = "pvalue.csv") #p-valueに関するデータをcsvファイルに出力
p-value에 대해서는 출력되는 대칭 행렬의 상측이 보정 후, 하측이 보정 전의 값이 출력되는 점에 주의.
4 인용
다음 사이트를 참고로 하고 있습니다.
htps //w w. r 도쿠 멘들 온. rg/파c게이 s/psych/ゔぇr 시온 s/2.0.9/와 피cs/코르 st
Reference
이 문제에 관하여(R 복수 데이터의 총당 단일 회귀 분석), 우리는 이곳에서 더 많은 자료를 발견하고 링크를 클릭하여 보았다
https://qiita.com/t_tac/items/3a3470430a7b8f9cdded
텍스트를 자유롭게 공유하거나 복사할 수 있습니다.하지만 이 문서의 URL은 참조 URL로 남겨 두십시오.
우수한 개발자 콘텐츠 발견에 전념
(Collection and Share based on the CC Protocol.)
사용하는 패키지와 데이터 세트를 R로 읽습니다.
준비하는 데이터 세트는 아래 이미지와 같은 형식입니다.
각 열에는 외부 요인이 있으며 각 행에는 샘플 데이터가 들어 있습니다.
(이번 사용한 것은 199 샘플×22 요인)
setwd("/Directory1/Directory2") #ディレクトリの指定
install.packages("psych") #psychパッケージのインストール
library("psych")#パッケージの呼び出し
dat1 <- read.csv("for_cor.csv", header=TRUE,row.names=1,fileEncoding="CP932")
# データセットの読み込み
# RはEXCELで通常保存される[.xlsx]の拡張子ファイルを扱うことは難しいので[.csv]形式でのデータセットを用いる
#列名、行名が不要の場合にはheader=FALSE, row.names <- NULL
2 상관 분석 수행
디폴트에서는 "pearson"의 상관 계수가 출력됩니다만 "spearman"또는 "kendall"에 의한 상관 계수의 산출도 가능합니다.
기타 adjust 소개 여러 파라미터를 변경할 수 있습니다.
자세한 내용은 [4 인용]에서.
result1 <- corr.test(dat1) #dat1についての相関分析結果をresult1に格納
result1 #結果表示(必要なら)
3 데이터 출력
각 항목에 대해 별도의 파일로 출력하고 싶었으므로 다음을 수행했습니다.
names(result1) #相関分析により得られたデータ項目の確認(必要なら)
cor1 <- result1$r #result1のr(相関係数)に関するデータのみを取り出しcor1に格納
write.csv(cor1, file = "cor1.csv") #相関係数のデータをcsvファイルに出力
corr.test(cbind(dat4))$p
pv1 <- result1$p #result1のp(p-value)に関するデータのみを取り出しpv1に格納
write.csv(pv1, file = "pvalue.csv") #p-valueに関するデータをcsvファイルに出力
p-value에 대해서는 출력되는 대칭 행렬의 상측이 보정 후, 하측이 보정 전의 값이 출력되는 점에 주의.
4 인용
다음 사이트를 참고로 하고 있습니다.
htps //w w. r 도쿠 멘들 온. rg/파c게이 s/psych/ゔぇr 시온 s/2.0.9/와 피cs/코르 st
Reference
이 문제에 관하여(R 복수 데이터의 총당 단일 회귀 분석), 우리는 이곳에서 더 많은 자료를 발견하고 링크를 클릭하여 보았다
https://qiita.com/t_tac/items/3a3470430a7b8f9cdded
텍스트를 자유롭게 공유하거나 복사할 수 있습니다.하지만 이 문서의 URL은 참조 URL로 남겨 두십시오.
우수한 개발자 콘텐츠 발견에 전념
(Collection and Share based on the CC Protocol.)
result1 <- corr.test(dat1) #dat1についての相関分析結果をresult1に格納
result1 #結果表示(必要なら)
각 항목에 대해 별도의 파일로 출력하고 싶었으므로 다음을 수행했습니다.
names(result1) #相関分析により得られたデータ項目の確認(必要なら)
cor1 <- result1$r #result1のr(相関係数)に関するデータのみを取り出しcor1に格納
write.csv(cor1, file = "cor1.csv") #相関係数のデータをcsvファイルに出力
corr.test(cbind(dat4))$p
pv1 <- result1$p #result1のp(p-value)に関するデータのみを取り出しpv1に格納
write.csv(pv1, file = "pvalue.csv") #p-valueに関するデータをcsvファイルに出力
p-value에 대해서는 출력되는 대칭 행렬의 상측이 보정 후, 하측이 보정 전의 값이 출력되는 점에 주의.
4 인용
다음 사이트를 참고로 하고 있습니다.
htps //w w. r 도쿠 멘들 온. rg/파c게이 s/psych/ゔぇr 시온 s/2.0.9/와 피cs/코르 st
Reference
이 문제에 관하여(R 복수 데이터의 총당 단일 회귀 분석), 우리는 이곳에서 더 많은 자료를 발견하고 링크를 클릭하여 보았다 https://qiita.com/t_tac/items/3a3470430a7b8f9cdded텍스트를 자유롭게 공유하거나 복사할 수 있습니다.하지만 이 문서의 URL은 참조 URL로 남겨 두십시오.
우수한 개발자 콘텐츠 발견에 전념 (Collection and Share based on the CC Protocol.)