similar to: problem with table.CAPM in PerformanceAnalytics

Hey Guys, i want to do a CAPM-GARCH model. I didn?t find anything posted online. (If there is something - shame on me - i didn?t find it.) My Problem: What is the difference if I let the residuals ?e? follow a garch process ? How do I do my regression analysis now? I began reading about regression analyis with heteroscedasticity, but didn?t get it. So i started programming. First

Hi This is my function. It serves an HTML page after the calculations. I'm connecting to a MSSQL DB using pyodbc. def CAPM(self,client): r=self.r cds="1590" bm="20559" d1 = [] v1 = [] v2 = [] print"Parsing GET Params" params=client.g[1].split("&") for items in

Hi all, I am also looking for interesting statistical experiments about testing and estimating CAPM, APT, Fama models, etc. using R using financial series data... please give me some pointers... I have been searching the R archives for the past a few hours and I vaguely got to know that there are programs do these interesting statistical things, but I just could not find where are they... I have

I'm trying to setup a snow grid using sockets (Windows 7).  On the test grid (my computer and another) I have an SSHD server up and running, can connect OK via public key authentication.  Running "makeSOCKcluster" on just my local machine works OK.  Running it to the other computer fails.  My SSH logs show the connection accepted, authenticated and then dropped.   I

Dear all, I want to do piecewise CAPM linear regression in R: RRiskArb−Rf = (1−δ)[αMktLow+βMktLow(RMkt−Rf)] + δ[αMkt High +βMkt High(RMkt −Rf )] where δ is a dummy variable if the excess return on the value-weighted CRSP index is above a threshold level and zero otherwise. and at the same time add the restriction: αMkt Low + βMkt Low · Threshold = αMkt High + βMkt High · Threshold to ensure

Inthe code below I was trying to to obtain the GMM estimates for CAPM (REGRESSION) for 36 stocks each have 180 observations,however it only gives me one output rather than 36. In SAS i would just put in a *By statement*. I have a variable TICKER that categorize them into 36 groups. *How can I obtain all 36 output instead of just one.* **

Hi, I am trying to calculate the std dev of returns of YHOO so far i got: getSymbols("YHOO") retYHOO <- Return.calculate(Cl(YHOO)) > sd(retYHOO) YHOO.Close NA but i received an NA....can any assist? tks! -- View this message in context: http://www.nabble.com/Newbie-question-re-stddev%2C-quantmod-and-performanceanalytics-tp25001293p25001293.html Sent from the R help

Hi All, I am trying to plot 4 graphs on to 1 page using layout(...), or par(mfcol = c(...)); with the function QQplot from the package PerformanceAnalytics. The problem is that, no matter what order I use, it only plots 3 graphs on to 1 page and the last QQplot is shunted to the next page. Also, this only happens to the QQplot, i.e. there is no problem with 4 Histograms. set.seed(1033) data

We are pleased to announce the availability on CRAN of PerformanceAnalytics version 0.9.6. This is a feature and bugfix release. http://cran.r-project.org/src/contrib/Descriptions/PerformanceAnalytics.html PerformanceAnalytics is a library of econometric functions for performance and risk analysis. This library aims to aid practitioners and researchers in utilizing the latest research in

We are pleased to announce the availability on CRAN of PerformanceAnalytics version 0.9.6. This is a feature and bugfix release. http://cran.r-project.org/src/contrib/Descriptions/PerformanceAnalytics.html PerformanceAnalytics is a library of econometric functions for performance and risk analysis. This library aims to aid practitioners and researchers in utilizing the latest research in

Dear R users, I am trying to carry out MLE of the time-varying CAPM using the FKF package. My approach so far has been to try and adapt the example given in the help file found using ?fkf which demonstrates the MLE of an ARMA(2,1) model. When I attempt to run my R code (given below) I get the following error: Error in fkf(a0 = sp$a0, P0 = sp$P0, dt = sp$dt, ct = sp$ct, Tt = sp$Tt, : Some of

Hi all I am looking for an add-in. I am currently working on something and I use daily data of closing stock prices. As not all companies are traded daily (e.g. on monday, then on thursday etc) at the stock exchange, there is satistically a problem. There are some papers which explain the approach to handle infrequent trading of a stock or non synchronous data and beta estimation (Dimson, 1979;

# Using Moving/Rolling Windows, here we do an OLS Regression with corrections for #Heteroscedasticity and Autocorrelations (HAC) using Newey West Method. This code is a #extension of Ajay Shah?s code for moving windows simple OLS regression. # The easiest way to adjust for Autocorrelations and Heteroscedasticity in the OLS residuals is to #use the coeftest function that is included in the

Hi, I have a problem with a regression I try to run. I did an estimation of the market model with daily data. You can see to output below: /> summary(regression_resn) Time series regression with "ts" data: Start = -150, End = -26 Call: dynlm(formula = ror_resn ~ ror_spi_resn) Residuals: Min 1Q Median 3Q Max -0.0255690 -0.0030378 0.0002787

Hi all ARCH modelling I have a problem now on how to proceed with further steps in my analysis. I did a linear OLS regression with my daily data of stock and index returns. There is now the problem of arch in my error terms. Thus I used the following r command: garch(resid_desn, order=c(0,2)) ## This ARCH(2) process seems to fit the best after trial and error. Consequently, I get there three

Dear Experts, I am trying to do a time-stratified case-crossover analysis on air pollution data and number of myocardial infarctions. In order to avoid model selection bias, I started with a simple simulation. I'm still not sure if my simulation is right. But the results I get from the "ts-case-crossover" are much more variable than those from a glm. Is this: a. Due to

Dear Experts, I am trying to do a time-stratified case-crossover analysis on air pollution data and number of myocardial infarctions. In order to avoid model selection bias, I started with a simple simulation. I'm still not sure if my simulation is right. But the results I get from the "ts-case-crossover" are much more variable than those from a glm. Is this: a. Due to the simple

Hi R Helpers, I recently tried to take advantage of the ability to download all the tickers in the S&P 500 using the functionality of tidyquant, but it threw an error. For summary, the set of commands that I ran was library(tidyquant) tq_index_options() tq_index("SP500") sessionInfo() R feedback including error message and sessionInfo are provided below. Guidance would be

I have hundreds of megabytes of price data time series, and perl scripts that extract it to tab delimited files (I have C++ programs that must analyse this data too, so I get Perl to extract it rather than have multiple connections to the DB). I can read the data into an R object without any problems. thedata = read.csv("K:\\Work\\SignalTest\\BP.csv", sep = "\t", header =

Hi, I am trying to install a package (PortfolioAnalytics) of R 2.15.2 on Win 7. I have tried the following instructions on : http://stackoverflow.com/questions/11105131/cannot-install-r-forge-package-using-install-packages I used svn checkout svn://svn.r-forge.r-project.org/svnroot/returnanalytics/ to get the package and then copy it to my Win 7. Then, on Win 7 I run: >