[R-br] Cálculo da correlação entre vetores que apresentam Nas
Michelle Bau Graczyk
mbgraczyk em gmail.com
Segunda Maio 25 15:47:32 BRT 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
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20150525/0da93ff0/attachment.html>
Mais detalhes sobre a lista de discussão R-br