similar to: Value at Risk using a volatility model?

Ok, I try it again with plain text, with a simple R code example and just sending it to the r list and you move it to sig finance if it is necessary. I try to be as detailed as possible. I want to fit a distribution to my financial data using a volatility model to estimate the VaR. So in case of a normal distribution, this would be very easy, I assume the returns to follow a normal distribution

Dear everyone, i`m a german economics student, writing my master´s thesis about "Multivariate Volatility Models". After having read about theoretical aspects of Multivariate GARCH ans Stochastic Volatility Models, I would like to compare DCC-GARCH and DC-SV with help of an empirical application. I figuered out that one has to use MCMC-simulation-methods for that. Some days ago I

Dear All, i want to use garch model in return of stock. and the data should presence volatility cluster (Heteroscedasticity). Do you know how to test volatility cluster (the presence of heteroscedasticity) in series data of stock return in R? Is it using Langrange Multiplier (LM) ARCH test? what package i should use? I really need the help. Thanks for the attention. Eko A P

Hi, as i suppose to know in a stationary GARCH(1,1) model the sum of alpha and beta has to be smaller than 1. But if i use the garchfit() function from the package fGarch for my timeseries the sum is bigger than 1. The adf.test tells me a p-value smaller than 0.01 instead. What does this mean for me? Can i trust in the coefficients in this case? mfg user84 -- View this message in context:

Dear All, I have loaded package(tseries), but when I run predict.garch(...) R tells me could not find function "predict.garch", however ?predict.garch shows me something. I am confused about this. How can I forecast garch volatility? I have tried: predict(...,n.ahead=...),give me fitted value predict(...,n),give me NA,NA

Hey, I'm trying to implement a GARCH model with Johnson-Su innovations in order to simulate returns of financial asset. The model should look like this: r_t = alpha + lambda*sqrt(h_t) + sqrt(h_t)*epsilon_t h_t = alpha0 + alpha1*epsilon_(t-1)^2 + beta1 * h_(t-1). Alpha refers to a risk-free return, lambda to the risk-premium. I've implemented it like this: #specification of the model

This script below pulls yahoo data via a function in quantmod, then massages the data around to forumalate a 3D graph with RGL library, attached is a ggplot to show the data i'm trying to create a surface with in separate line geoms . the issue is that the 3D graph looks very ugly and cut up because of the limited quantities of points on the front month expirations.. can anyone tell me whats

I think the rule is that you can do anything as long as you don't complain. If you want to complain, you must follow the instructions. -- Jari Oksanen in Re: [Rd] Keeping up to date with R-devel -- Patrick Burns pburns at pburns.seanet.com twitter: @burnsstat @portfolioprobe http://www.portfolioprobe.com/blog http://www.burns-stat.com (home of: 'Impatient R' 'The R

Dear R People: There is a function from the fBasics library to get the probability and quantiles for the hyperbolic probability function. Is there one that will estimate parms of the hyperbolic probability function from a data set, please? Thanks in advance! Sincerely, Erin Hodgess mailto: hodgess at gator.uhd.edu R Version 2.0.1 windows

R-help is all about solving R problems. So here ya go: http://www.portfolioprobe.com/2011/09/12/solve-your-r-problems/ -- Patrick Burns pburns at pburns.seanet.com twitter: @portfolioprobe http://www.portfolioprobe.com/blog http://www.burns-stat.com (home of 'Some hints for the R beginner' and 'The R Inferno')

I created a EM algorithm for Generalized hyperbolic distribution. I want to estimate mutheldaplus, sigmatheldaplus, betasigmaplus in my code. After getting use these value , then my iteration have to be begin of this code. But I can not to do iteration part. Can you help me use my code and get iteration ? Do know any useful code for EM algorithm for Generalized Hyperbolic library(QRMlib)

Dear all, I am using GARCH (1,1) model to simulate volatility. But seems that I am missing something about how it works in R. The following code produces negative results, though vola cannot be. What is wrong here? library("fSeries") library("fGarch") spec = garchSpec(model = list(omega = 0.01, alpha = 0.13, beta = 0.86)) gat <- garchSim(spec, n = 10) Thanks a lot!

Hi team, I am facing issues with PearsonIV distribution fitting in R. I am applying Hyperbolic and PearsonIV distributions on the equity returns in UK over a period of 30 years. For the same data set i am getting strikingly different results under which Hyperbolic distribution does produce negative percentiles of the return after fitting but PearsonIV distribution does not. I think there is an

How to use the package generalized hyperbolic distribution in order to estimate the four parameters in the NIG-distribution? I have a data material with stock returns that I want to fit the parameters to. -- View this message in context: http://r.789695.n4.nabble.com/Generalized-Hyperbolic-distribution-tp3504369p3504369.html Sent from the R help mailing list archive at Nabble.com.

Dear all, We seem to have found a "strange" behaviour in the hyperbolic tangent function tanh on Windows. When running tanh(356 + 0i) the Windows result is NaN + 0.i while on Mac the result is 1 + 0i. It doesn't seem to be a floating point error because on Mac it is possible to run arbitrarily large numbers (say tanh(

Hello! I am optimizing my code in R and for this I need to know a bit more about the internals. It would help tremendously if someone could link me to a page with O()-complexities of all the operations. In this particular case, I need something like a linked list with O(1) insertLast/First ability. I can't preallocate a vector since I do not know the final size of the list ahead of time. The

Dear R forum helpers,I have following datatrans <- data.frame(currency_transacted = c("EURO", "USD", "USD", "GBP", "USD", "AUD"), position_amt = c(10000, 25000, 20000, 15000, 22000, 30000))date <- c("12/31/2010", "12/30/2010", "12/29/2010", "12/28/2010", "12/27/2010",

Hi, I have limited experience on R and recently started using REXcel. Although I have been able to run both simple functions (like mean etc) and some complex ones (like Principal Component analysis, PCA) using RExcel, I am facing some problems while running EGARCH model. For this I have downloaded the 'betategarch' package for R to run EGARCH with student t dist. Although the package has

Dear all Currently, I'm implementing the generalized hyperbolic distribution into Splus. Unfortunately the Bessel function is not implemented in Splus. In R the Bessel function does exist but it is an internal function and I'm not able to look at the code. Is there any possibility to see the code of the Bessel function in R or does anybody has an implementation of the Bessel function in

Hello, I'm very interested in using financial time series data, but I'm a beginner of R programming. I'd like to fully understand internal logics on several time-series related packages such as quantmod, xts, zoo, chron, etc. So, I read some books, 'R Cookbook' and 'Art of R Programming' and another simple tutorials. But I still can't understand grammars of the