XTS data de diferentes fontes. Usando R para calcular beta

Question

Aug 03, 2012, 09:27 PM

XTS data de diferentes fontes. Usando R para calcular beta

Sou um pouco novo para R. Imagino que meu erro seja trivial para os experientes.

Eu estou tentando escrever um programa R que irá calcular beta para um número de ações. Os símbolos das ações são lidos deInput.csve os dados são baixados do yahoo. O código, em seguida, percorre um cálculo beta para cada ação e gera um csv resumindo as regressões.

Consegui que o código funcionasse quando uma única taxa livre de risco era assumida em todos os períodos, mas acredito que possa precisar usar a taxa livre de risco real em cada mês ao normalizar os retornos excedentes. Estou tendo problemas com este passo. O xts é baixado com êxito do FRED (GS20), mas quando o retorno é subtraído do retorno de segurança e de mercado, ele produz um xts de comprimento zero. Isso mata o programa.

Eu acredito que isso pode ser porque as datas estão em diferentes formatos entre FRED e Yahoo. Notei que o comando getSymbols ignorou as datas de saída. Qualquer ajuda seria apreciada.

require(PerformanceAnalytics)
require(quantmod)
require(car)

setwd("R Projects/Beta Test")

proxy <- read.csv("Input.csv",header=FALSE)[,1]
summary <- as.data.frame(matrix(0, ncol = 5, nrow = 0))

mar <- getSymbols("^GSPC", src = "yahoo", from = as.Date("2006-01-01"),
                  to = as.Date("2011-12-31"),auto.assign=FALSE)
riskFree <- getSymbols("GS20", src = "FRED", from = as.Date("2006-12-01"),
                       to = as.Date("2011-12-31"),auto.assign=FALSE)

for (n in proxy){

    sec <- getSymbols(n, src = "yahoo", from = as.Date("2006-01-01"),
                      to = as.Date("2011-12-31"),auto.assign=FALSE)

    #Monthly Returns
    #ERROR PRODUCED HERE
    sec.xsmonthly <- monthlyReturn(to.monthly(sec),type="log") - riskFree
    mar.xsmonthly <- monthlyReturn(to.monthly(mar),type="log") - riskFree

    sec.reg <- lm(sec.xsweekly ~ mar.xsmonthly + lag(mar.xsmonthly,-1))

    summary[n,1] <- coef(sec.reg)[1]
    summary[n,2] <- coef(sec.reg)[2]
    summary[n,3] <- coef(sec.reg)[3]
    summary[n,5]<-summary(sec.reg)$r.squared
}

summary[,4] <- summary[,2]+summary[,3]

colnames(summary) <- c("Alpha","Beta","Beta-Lag","Sum Beta","R-Squared")
write.csv(summary,file="output.csv")

leaveComments