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