получить все возможные корреляции между двумя наборами данных - PullRequest
Новая
       10

получить все возможные корреляции между двумя наборами данных

0 голосов
/ 30 января 2019

Я пытаюсь найти несколько данных между каждой корреляцией

corr.test

У меня есть два набора данных df1 и df2 

df1<- structure(list(col1A = c(1.64, 0.03, 0, 4.202, 2.981, 0.055, 
0, 0.002, 0.005, 0, 0.002, 0.649, 2.55, 2.762, 6.402), col2A = c(2.635, 
0.019, 0, 5.638, 3.542, 0.793, 0.259, 0, 0.046, 0.004, 0.017, 
0.971, 3.81, 3.104, 5.849), col3A = c(0.91, 0.037, 0, 5.757, 
3.916, 0.022, 0, 0, 0.003, 0, 0.262, 0.136, 2.874, 3.466, 5.003
), col4A = c(1.027, 0.021, 0, 4.697, 2.832, 0.038, 0.032, 0.001, 
0.003, 0, 0, 0.317, 2.743, 3.187, 6.455)), class = "data.frame", row.names = c(NA, 
-15L))

вторые данные похожиниже 

 df2<-structure(list(col1 = c(2.172, 0, 0, 4.353, 4.581, 0.001, 0.027, 
0, 0.002, 0, 0, 0.087, 2.129, 4.317, 5.849), col2 = c(2.093, 
0, 0, 4.235, 3.166, 0, 0, 0.006, 0.01, 0, 0, 0.475, 0, 2.62, 
5.364), col3 = c(3.322, 0, 0, 4.332, 4.018, 0.049, 0.169, 0.004, 
0.02, 0, 0.032, 1.354, 2.944, 4.323, 5.44), col4 = c(0.928, 0.018, 
0, 3.943, 3.723, 0.02, 0, 0, 0, 0, 0.075, 0.136, 3.982, 3.875, 
5.83)), row.names = c("A", "AA", "AAA", "Aab", "buy", "yuyn", 
"gff", "fgd", "kil", "lilk", "hhk", "lolo", "fdd", "vgfh", "nghg"
), class = "data.frame")

Я хочу получить все возможные корреляции между двумя и извлечь все значения p и скорректированные значения p

Я использую 

library(psych)
corr.test(df1,df2, use = "pairwise",method="pearson",adjust="holm",alpha=.05,ci=TRUE,minlength=5)

, это не даетмне любое значение р.также я не могу контролировать любые перестановки для вычисления скорректированного значения p.

Я думал использовать следующее 

x <-df1[,1]
y <-df2[,2] 
corr_init <- cor(x,y) # original correlation
N <- 1000 # initialize number of permutations
count <- 0 # counts correlation greater than corr_init
for (i in 1:N) {
y_perm <- permute(y)
  if (cor(y_perm,x) > corr_init) count <- count+1
  }
p <- count/N #final p

, но затем я должен сделать это один за другим, и все же мне нужно извлечь каждый столбец и проверить ...

Мне интересно, есть ли лучший способ рассчитать всю корреляцию между двумя данными, получить значения R, значения p и P, скорректированные с определенным числом рандомизации?

1 Ответ

0 голосов
/ 30 января 2019

Это можно сделать с помощью пакета Hmisc: 

library(Hmisc)

df1_cor_matrix <- rcorr(as.matrix(df1), type = "pearson")
df2_cor_matrix <- rcorr(as.matrix(df2), type = "pearson")

Затем вы можете извлечь коэффициенты с помощью следующего: 

df1_coef <- df1_cor_matrix$r
df2_coef <- df2_cor_matrix$r

Вы можете извлечь значения pиспользуя следующее: 

df1_p_values <- df1_cor_matrix$P
df2_p_values <- df2_cor_matrix$P

Вы можете получить скорректированные значения p, используя функцию rcorr.adjust: 

rcorr.adjust(df1_cor_matrix, type = "pearson")
rcorr.adjust(df2_cor_matrix, type = "pearson")
...