[R-br] Calculo de correlcao entre coluas que apresentam NAS de varias matrizes

Quarta Março 18 14:38:34 BRT 2015

Boa tarde,

Eu tenho tres matrizes de tres empresas diferentes que estao dispostas da
seguinte maneira:

Empresa AA
data            9:31
            data
09:31 09:32 09:33 09:34 09:35 09:36 09:37 09:38 09:39 09:40 09:41 09:42
09:43 09:44  09:45  09:46  09:47

2004-07-01 700 200 7000 9400 1300 1500 6300 1800 4300 1300 4800 3900 2500
4300 4300

 2004-07-02 NA NA 68900 2800 4900 1200 12300 NA NA 6500 NA 2900 8900 11600
26900

 2004-07-06 800 4400 13700 10700 4700 4000 18700 9000 2100 2300 700 8900
10500 3900 8400

 2004-07-07 NA 5600 1000 7300 5600 11200 3900 5800 5700 4500 34400 4100 500
15300 500

 2004-07-08 NA NA NA NA NA NA NA NA NA NA 154100 42400 16700 4700 5900

 2004-07-09 NA 96900 1000 20400 18300 18000 22500 4100 6000 1000 1900 3700
600 2900 14600

 2004-07-12 1400 100 4500 3000 1300 4200 6500 4500 1600 1600 5600 7600 1900
4600 8300

 2004-07-13 NA 5600 3400 16900 9000 5700 16500 1800 21000 6300 5900 12500
51000 7500 13100

 2004-07-14 NA NA NA NA NA NA 164400 35000 41600 56600 27800 36500 55300
29900 13100

 2004-07-15 NA 8100 2000 6100 4900 4500 2500 1400 5200 4100 5400 4000 13400
3200 11300

 2004-07-16 282300 18100 30400 NA 15900 40800 9000 5300 5200 27600 7600
12100 8400 1000 8400

 2004-07-19 28700 1700 2900 9500 400 2800 1100 6200 11500 17900 11100 16000
6800 5200 7100

Onde as datas vao de 2004-07-01 ate 2004-12-31(ao todo 128)  e os horarios
vao de 09:31 ate 16:04( ao todo 397)

Eu quero calcular a correlacao entre essas tres empresas para cada hora
fixa,ou seja, para a empresa AA eu calculo com ela mesma e com a AIG e a
AXP para cada hora fixa, e por isso montei este programa:

A<-read.table(file="/home/usuario/Documentos/AAparaCalculodaMedia.txt",
header=TRUE)
AA.<-as.matrix(AA)
##write.table(AA.,"AA..txt")
AXP<-read.table("/home/usuario/Documentos/AXPparaCalculodaMedia.txt",header=TRUE)
AXP.<-as.matrix(AXP)
##write.table(AXP.,"AXP..txt")
AIG<-read.table("/home/usuario/Documentos/AIGparaCalculodaMedia.txt",header=TRUE)
AIG.<-as.matrix(AIG)

Symbols<-c("AA.","AIG.","AXP.")
matrizAA<-matrix()
j<-1
for (j in 1:length(Symbols)){
  empresa<-paste0(Symbols)#,"[",j,"]")

  for (i in 2:397){
    empresa1<-as.numeric(paste0(empresa,"[,",i,"]"))
    head(empresa1)
    matrizAA[,i-1]<-cor(as.numeric(AA.[,i]), as.numeric(empresa1),
use="pairwise.complete.obs")

  }

}

Mas da o erro: Erro em cor(as.numeric(AA.[, i]), as.numeric(empresa1), use
= "pairwise.complete.obs") :
  incompatible dimensions
Além disso: Mensagens de aviso perdidas:
NAs introduzidos por coerção

Eu achava que usando o use = "pairwise.complete.obs" as correlacoes seriam
calculadas apenas entre pares de numeros, mas acho que isso nao esta
acontecendo.
Os tres arquivos vou deixar em dropbox.

Muito obrigada,

Michelle
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20150318/b2c5038c/attachment.html>