[R-br] Cálculo da correlação entre vetores que apresentam Nas

25 Mai 2015

      Boa tarde,

Eu estou calculando a correlação entre vetores que apresentam NAs. Eu usei
a função use="pairwise.complete.obs" mas da o erro:

Mensagens de aviso perdidas:
In cor(AA[, i], get(Symbols[j])[, i], use = "pairwise.complete.obs") :
  the standard deviation is zero

Eu preciso que ele calcule a correlação entre os pares de vectores que
contém valores sem propagar os NAs. Alguém saberia me dizer se preciso usar
outra função ou se estou fazendo algo errado?

Abaixo segue o código e um exemplo de um dos banco de dados usados.
Muito obrigada!

AA<-read.table("/Users/bau/MatrizCorrelação25%/AA_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",
header=TRUE)
AXP<-read.table("/Users/bau/MatrizCorrelação25%/AIG_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE)
AIG<-read.table("/Users/bau/MatrizCorrelação25%/AXP_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE)
BA<-read.table("/Users/bau/MatrizCorrelação25%/BA_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE)
CAT<-read.table("/Users/bau/MatrizCorrelação25%/C_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE)
C<-read.table("/Users/bau/MatrizCorrelação25%/CAT_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",
header=TRUE)
DD<-read.table("/Users/bau/MatrizCorrelação25%/DD_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",
header=TRUE)
DIS<-read.table("/Users/bau/MatrizCorrelação25%/DIS_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE)
GE<-read.table("/Users/bau/MatrizCorrelação25%/GE_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",
header=TRUE)
GM<-read.table("/Users/bau/MatrizCorrelação25%/GM_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",
header=TRUE)
HD<-read.table("/Users/bau/MatrizCorrelação25%/HD_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE)
HON<-read.table("/Users/bau/MatrizCorrelação25%/HON_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE)
HPQ<-read.table("/Users/bau/MatrizCorrelação25%/HPQ_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE)
IBM<-read.table("/Users/bau/MatrizCorrelação25%/IBM_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE)
INTCc<-read.table("/Users/bau/MatrizCorrelação25%/INTC_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",
header=TRUE)
INTC<-INTCc[,91:485]
JNJ<-read.table("/Users/bau/MatrizCorrelação25%/JNJ_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE)
JPM<-read.table("/Users/bau/MatrizCorrelação25%/JPM_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE)
KO<-read.table("/Users/bau/MatrizCorrelação25%/KO_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE)
MCD<-read.table("/Users/bau/MatrizCorrelação25%/MCD_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE)
MMM<-read.table("/Users/bau/MatrizCorrelação25%/MMM_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE)
MO<-read.table("/Users/bau/MatrizCorrelação25%/MO_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",
header=TRUE)
MRK<-read.table("/Users/bau/MatrizCorrelação25%/MRK_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE)
MSFTt<-read.table("/Users/bau/MatrizCorrelação25%/MSFT_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE)
MSFT<-MSFTt[,91:485]
PFE<-read.table("/Users/bau/MatrizCorrelação25%/PFE_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",
header=TRUE)
PG<-read.table("/Users/bau/MatrizCorrelação25%/PG_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE)
SBC<-read.table("/Users/bau/MatrizCorrelação25%/SBC_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE)
UTX<-read.table("/Users/bau/MatrizCorrelação25%/UTX_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE)
VZ<-read.table("/Users/bau/MatrizCorrelação25%/VZ_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE)
WMT<-read.table("/Users/bau/MatrizCorrelação25%/WMT_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE)
XOM<-read.table("/Users/bau/MatrizCorrelação25%/XOM_Volume25medioParaCalculoCorrelaçãoREORGANIZADO.txt",header=TRUE)

matrizAA<-matrix(nrow=395, ncol=31)
Symbols<-c("AA","AIG","AXP","BA","C","CAT","DD","DIS","GE","GM","HD","HON","HPQ","IBM","INTC","JNJ","JPM","KO","MCD","MMM","MO","MRK","MSFT","PFE","PG","SBC","UTX","VZ","WMT","XOM")

for(j in 1:30)
{
  empresa<-Symbols[j]

  for (i in 2:395){

    matrizAA[i-1,j+1]<-cor(AA[,i],get(Symbols[j])[,i],
use="pairwise.complete.obs")
    ##show(get(Symbols[j])[,i])

  }
}
matrizAA[,1]<-c(1:395)
write.table(matrizAA,"/Users/bau/MatrizCorrelação25%/AAcorrelaçãoParaMatriz25.txt")
...
show(INTC)
      X09.31  X09.32  X09.33  X09.34  X09.35 X09.36 X09.37 X09.38  X09.39
X09.40
1      95991      NA      NA      NA      NA     NA     NA     NA      NA
  NA
2         NA      NA      NA      NA      NA     NA     NA     NA      NA
  NA
3     165142  288968      NA      NA  108902     NA     NA     NA      NA
  NA
4     118675      NA      NA  106849      NA     NA     NA     NA      NA
  NA
5     143731      NA      NA      NA      NA     NA     NA     NA      NA
  NA
6     182431      NA  320616  121284  102200 110047 177755     NA   89746
  NA
7     385596  233253  102160   98434      NA 107502 125117     NA      NA
  NA
8     164819      NA      NA      NA   82725     NA     NA     NA      NA
  NA
9    1037347 1391429 1168113 2144977  996927 555185 399196 665563  911069
594260
10    375979  229744  123207  142073  389064 164749     NA     NA      NA
208184
11    734943  258833  133831  228149  173972 148795 108966 182753  109337
211194
12    302175   95693      NA  131569      NA     NA 148767 185680  103596
  NA
13    301830      NA      NA      NA      NA     NA     NA     NA   96945
111397
14    380442  175182      NA      NA  154135 107285     NA  94698   82688
  NA
15    257121  122878      NA      NA  113273     NA 170090     NA      NA
  NA
16    135805      NA      NA      NA      NA  96105 230759 178512   98562
152584
17    143954  105240      NA  135897      NA 162295     NA  89958  175599
171491
18    110088      NA      NA  232589   91542     NA  90570 118529      NA
176981
19    133450      NA      NA      NA  103523     NA     NA     NA      NA
  NA
20    151424  535320  124938  202910  150950 228580     NA     NA      NA
  NA
21    174497  146076  379007  230754      NA 259161 262723 141953      NA
209095
22    202101  125722  178210  241326  222106     NA     NA     NA      NA
234913
23    193474   81358  126753  157524      NA 140525     NA 132483  122650
 94888
24    226467      NA  106886      NA      NA     NA     NA     NA      NA
  NA
25     94991      NA      NA      NA      NA 100226     NA 138565      NA
  NA
    X09.41  X09.42 X09.43 X09.44 X09.45  X09.46 X09.47 X09.48  X09.49
X09.50 X09.51
1       NA      NA     NA     NA     NA      NA     NA 125230      NA
NA     NA
2       NA      NA     NA     NA     NA      NA     NA     NA      NA
NA     NA
3    95210      NA  83615     NA     NA      NA     NA     NA   81950
NA     NA
4       NA      NA     NA     NA 525228      NA     NA     NA      NA
122300     NA
5       NA      NA     NA     NA     NA      NA     NA     NA      NA
NA     NA
6       NA      NA     NA 118085     NA      NA     NA     NA      NA
NA     NA
7   101348      NA     NA  85359     NA      NA     NA  84274   95450
 90829     NA
8       NA      NA     NA     NA     NA      NA     NA     NA      NA
 92065     NA
9   475210 1588546 828340 591364 849922  604366 592396 500498  231407
622671 633772
10      NA      NA  82904 115703 133982  600720     NA     NA  139664
142184  92034
11   91429   84681 288650     NA 182264  134730     NA 234712   82140
NA 138098
12      NA      NA     NA 103916 130810  164900  94492 282866  117227
 94789     NA
13      NA      NA     NA     NA     NA      NA     NA     NA      NA
NA     NA
14      NA  290261  82617     NA     NA   91485     NA     NA      NA
NA     NA
15      NA      NA     NA     NA     NA      NA 114630 120390      NA
150594     NA
16  191759      NA     NA     NA     NA  111460 212266 135440      NA
NA 143878
17  163178      NA     NA 181344  94675      NA 124844     NA      NA
NA  89650
18      NA      NA     NA     NA     NA      NA     NA     NA      NA
NA     NA
19  168241      NA 204800     NA     NA      NA     NA     NA      NA
NA     NA
20  382726  197013 261650     NA 193270  191328 162645     NA      NA
187996 134218
21  257324  248162 391465 212684 125445  101610 228834     NA      NA
226410 128451
22      NA      NA 209705     NA     NA      NA     NA     NA      NA
NA  88528
23      NA  145010     NA     NA     NA   86450 106323     NA  127235
131685     NA
24      NA      NA     NA     NA     NA      NA     NA     NA      NA
NA     NA
25      NA      NA     NA     NA     NA      NA     NA     NA      NA
NA     NA
    X09.52 X09.53 X09.54 X09.55 X09.56 X09.57 X09.58 X09.59 X10.00 X10.01
X10.02
1       NA     NA     NA     NA     NA     NA 189095     NA     NA     NA
  NA
2       NA     NA 212714 124657     NA 125427     NA     NA     NA     NA
  NA
3       NA     NA     NA     NA     NA     NA 108377     NA     NA     NA
  NA
4       NA     NA     NA     NA     NA     NA     NA     NA     NA     NA
  NA
5       NA     NA     NA     NA     NA     NA     NA     NA 111310     NA
  NA
6   135987     NA     NA     NA     NA     NA     NA     NA     NA     NA
  NA
7       NA 143816  86290  91866     NA     NA 186791     NA     NA     NA
  NA
8       NA     NA     NA     NA     NA     NA     NA     NA     NA     NA
  NA
9   398724 348269 209850 192490 427251 419138 879114 231852 289164 164084
560886
10  152704     NA 115500     NA  91505     NA 211354 107150 135024     NA
110234
11   88621 212824 194126 147968     NA     NA 113810     NA     NA     NA
  NA
12      NA  85738     NA  84984     NA     NA     NA     NA     NA     NA
  NA
13   94627     NA     NA     NA     NA     NA     NA     NA     NA     NA
  NA
14  107388     NA     NA     NA 160841  97684 291855     NA     NA     NA
  NA
15      NA 148828     NA 370458     NA     NA     NA     NA 135246     NA
  NA
16      NA 187592     NA     NA  87650     NA     NA     NA     NA     NA
  NA
17      NA 138967 117446     NA     NA     NA     NA     NA     NA     NA
  NA
18  230329 261154     NA     NA 153873     NA     NA     NA     NA     NA
  NA
19      NA     NA     NA     NA     NA     NA     NA     NA     NA     NA
  NA
20  295116 163870 145028 231132 154587 166770 286393 143550 190173     NA
146736
21      NA     NA     NA     NA     NA     NA  91722 325090     NA     NA
  NA
22      NA     NA     NA     NA     NA     NA 112804     NA     NA     NA
  NA
23      NA     NA     NA     NA     NA     NA     NA     NA     NA     NA
 81763
24      NA 125616     NA     NA     NA     NA     NA     NA     NA     NA
  NA
25      NA     NA     NA     NA     NA     NA  84216     NA     NA     NA
  NA

[R-br] Cálculo da correlação entre vetores que apresentam Nas

Michelle Bau Graczyk