[R-br] Problema nas dimensões. Cálculo Matricial

[Luciano Ramos Gonçalves]

Mais uma dúvida simples:
O problema é no final do código ( em vermelho) :

datas <- read.table("ddd.txt", header=F)
datas <- as.Date(datas[,1], format="%m/%d/%y")
i <- read.table("iii.txt", header=F)
i <- i[,1]
cupons <- read.table("cc1.txt", header=F)
cupons <- as.Date(cupons[,1], format="%m/%d/%Y")
fluxo <- read.table("ff1.txt", header=F)
fluxo <- fluxo[,1]
falta <- outer (cupons, datas, "-") *.69
falta <- matrix(falta,nrow(falta),ncol(falta))
falta <- ifelse(falta < 0, NA, falta)
VP <- cbind((fluxo)/((1+i)^(falta)))
plot(colSums(VP, na.rm=T), type="l")

i tem 297 elementos, datas tem 297 elementos
cupons tem 14 elementos, e fluxo tem 14 elementos;
Logo, falta e VP tem 297 colunas e 14 linhas
quanto ploto VP, aparecem 14 vales no gráfico. Acho que algo de errado 
em (1+i)^(falta).
O que quero fazer é,(1+i[1])^(falta[,1]), (1+i[2])^(falta[,2]), ... , 
(1+i[n])^(falta[,n]), ou seja, elevar (1+ i), para cada i["n"] em todos 
os elementos da coluna "n" em "falta". Assim, acabaria também com uma 
matriz 14 por 297. E quero dividir os 14 elementos de "fluxo" pelos 14 
elementos de cada coluna na matriz gerada. Deu pra entender?

Obrigado,

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