similar to: quantile regression using copulas

Hi r-users, I try to fit the t copula using the gamma marginals.  But I got error message which I don't really understand. Thank you for any help given. myCop.t <- ellipCopula(family = "t", dim = 2, dispstr = "toep", param = 0.5, df = 8) myCop.t myMvd <- mvdc(copula = myCop.t, margins = c("gamma", "gamma"), paramMargins = list(list(mean = 0, sd

I am using R 2.5.0 on windows XP and trying to fit copula. I see the following code works for some users, however my code crashes on the chol. Any suggestions? > mycop <- tCopula(param=0.5, dim=8, dispstr="ex", df=5) > x <- rcopula(mycop, 1000) > myfit <- fitCopula(x, mycop, c(0.6, 10), optim.control=list(trace=1), method="Nelder-Mead")

Hi r-users,   I really need help in fitting the t-copula.  I try to reproduce the example given by Jun Yan in “Enjoy the joy of copula” but I’m not sure how to correct the error based on the error message.  I tried so many ways but still could not get it working.   loglik.marg <- function(b, x) sum(dgamma(x, shape = b[1], scale = b[2], log = TRUE))   ctrl <- list(fnscale = -1)   #dat <-

Hi r-users, I really need help in fitting the t-copula. I try to reproduce the example given by Jun Yan in "Enjoy the joy of copula" but I'm not sure how to correct the error based on the error message. I tried so many ways but still could not get it working. loglik.marg <- function(b, x) sum(dgamma(x, shape = b[1], scale = b[2], log = TRUE)) ctrl <- list(fnscale

Hi, I'd like to know why using the program "R" I can't add a number of margins> 3, I have a problem with the graphics. Post here my script: > myCop.norm <- ellipCopula(family = "normal", dim = 3, param = 0.4) > myMvd <- mvdc(copula = myCop.norm, margins = c("norm", "norm","norm"), > paramMargins = list(list(mean = 0, sd

Hi, I want to run a quantile regression (Y=a+bX+e) using normal and t copula for my dissertation. I 've read the documentation of "copula" and "copBasic". However, I still have difficulty to deal with my data. Details are as following: I've already loaded xls data into r using "XLConnect" package. excel.file<-file.path("Q:/dailyvstoxx.xls")

Is anybody using the Copula package in R? The particular problem I'm facing is that R is not acknowledging the fitCopula command/function when I load the package and (try to) run something very simple: fit1 <- fitCopula(x1 = list(u11,u12,u13,u14,u15,u16,u17,u18), tCopula, optim.control = list(NULL), method = "BFGS") Anybody also using it, successfully or unsuccessfully?

I am porting some FORTRAN to R in which an Inf triggers an if(). The trigger is infinite on exp(lgamma(OVER)). What is the canonical R style of determining OVER when exp(OVER)== Inf? The code structure that I am porting is best left intact--so I need to query R somehow to the value of OVER that causes exp(lgamma(OVER)) to equal Inf. On my system, exp(lgamma(171)) is about first to equal Inf.

I upgraded to R-2.2.1 on two PPC G5 computers today. Further I want to work with the tseries package for the first time. As root with R CMD INSTALL tseries_0.10-0.tar.gz I get the following gcc-3.3 -bundle -flat_namespace -undefined suppress -L/usr/local/lib - o tseries.so arma.o bdstest.o boot.o dsumsl.o garch.o ppsum.o tsutils.o -framework vecLib -L/usr/local/lib/gcc/powerpc-apple-

Dear list, I have a question concerning the above mentioned methods in the pls package with respect to the loadings matrix produced by the call. In some work I am doing I have found that the values produced are nearly of the same magnitude but of opposite sign. When I use the example data (sensory) I find this result reproduced. I am prepared to work this through but I have a feeling that

CRAN et al., I would like to add an extented introduction or other arbitrary sections to my package lmomco. I have been shipping inst/doc/Introduction.Rd. I would like to have this content inserted to the front of the PDF build for the CRAN. The R-exts.pdf seems to be a little silent on this subject? For my purposes, I have been doing this R CMD Rd2dvi --pdf

Hi, Could someone please suggest where I might find some instructions / tutorials / FAQs that describe how to create a frequency distribution and cumulative frequency distribution in R using different class withs. I have about a 2-million observations (distances between points ranging from sub-millimetre to about 400km, and I want to get a feel for how they are distributed). I'd like the

Hi all! I would like to import a matlab or gauss code to R. Could you help me? Bye, Sebasti?n. 2007/6/23, r-help-request en stat.math.ethz.ch <r-help-request en stat.math.ethz.ch>: > Send R-help mailing list submissions to > r-help en stat.math.ethz.ch > > To subscribe or unsubscribe via the World Wide Web, visit >

Dear all, I wanna to do multiple comparison test for proportions (multiple chi squre ?), could someone tell me how in R, thank you!

I have question (curiosity) regarding returned values of R's qcauchy () function, for nonexceedance probability (F). It seems the ideal returned range of cauchy distribution should be [-Inf,Inf]. For F=0 > qcauchy(0) [1] -Inf but for F=1 > qcauchy(1) [1] 8.16562e+15 It seems to me that the proper return value should be Inf??? For default (location=0,scale=1) quantile function of

I have a function XOLD() from a nearly verbatim port of legacy FORTRAN in a package. I have remplemented this function as XNEW() using much cleaner native R and built-in functions of R. I have switched the package to the XNEW(), but for historical reasons would like to retain the XOLD() somewhere in the package directory structure. An assertion through a README or other will point to

I am looking at using R to analyze time series data containing a tidal component. I need to remove the tidal signal to extract the time series of the phenomena I seek to study. A browse of R-project search engines has not been too fruitful? I've found 'hoa' and 'Rwave', but need further help getting started. THANKS. -wa

Hi, I found this document, but it concerns S+. If it could interest you'll see: http://faculty.washington.edu/ezivot/book/QuanCopula.pdf Cordially Vito You wrote: Dear R: Is there a function or a reference to simulate Gumbel copulas, please? Thanks in advance! Sincerely, Erin Hodgess mailto: hodgess at gator.uhd.edu R version 2.0.1 windows ===== Diventare costruttori di soluzioni

I am trying to make a loglikelihood function using copulas. I am trying to use mvdc to find the density function. When I run this I got the error that the pdf and cdf of my function tobit doesn't exist. Can somebody guide me where my mistake is? dtobit <- function(beta,sigma, x, y) {ifelse(y>0, dnorm(y,x%*%beta, sigma),(1-pnorm((x%*%beta)/sigma)))} ptobit <- function(beta,sigma, x,

Dear All, I'm now looking for R-code on how to find the Goodness of Fits for Archimedean Copulas. If anyone have a guide for this problems please lets me know. Your prompt action is much appreciated. Regards, Ummul [[alternative HTML version deleted]]