similar to: Copula package

is() does not catch parent S3 classes: > library(splines) > temp <- bs(1:99, df=5) > class(temp) [1] "bs" "basis" > is(temp, "basis") [1] FALSE In contrast, is() does catch parent S4 classes: > library(copula) > norm.cop <- ellipCopula("normal", param = c(0.5, 0.6, 0.7), + dim = 3, dispstr = "un")

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

Hello everybody, I have a big problem that I do not manage to solve ! I will be very grateful if you can solve this ! I want to fit a t Copula with the copula package : > student.cop <- ellipCopula("t", param = c(0.5, 0.6, 0.7), dim = 3, dispstr = "un",df=5) > x<-rcopula(student.cop,1000) > fit <- fitCopula(x, student.cop, c(0.5,0.5,0.5,5)) And there is an

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")

Hola, Quiero saber si es posible programar una cópula donde las funciones marginales son multivariantes, siguiendo el esquema del package ''Copula''. Es decir, Copula(F(x,y), G(w,z))) En el caso de funciones marginales univariantes, un ejemplo de la normal multivariante quedaria de la siguiente forma,

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 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,

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

Hey all, I'm trying to construct a 7-dimensional normal copula using the copula package. I'd like to supply as parameter a randomly generated correlation matrix (that I'll convert to a vector so I can feed it to the normalCopula function). What order do the pairwise correlations inside that vector have to be in? In other words: What does the normalCopula function "expect"?

I'm a new user of R and a novice user in copula R package. I want to fit 3-dimensional t copula for my trivariate data. So I used the command t.cop <- tCopula(c(0.785,0.283,0.613),dim=3,dispstr="un",df=6,df.fixed = TRUE) where c(0.785,0.283,0.613) is the correlation pattern of my data with 0.785 pearson correlation between variable 1-2, 0.283 correlation between 1-3 and 0.613

Hi all, I need to test gof of 3-dimensional t copula for my trivariate observed data set. So I used the command t.cop <- tCopula(c(0.785,0.283,0.613),dim=3,dispstr="un",df=6,df.fixed = TRUE) where c(0.785,0.283,0.613) is the correlation pattern of my data with 0.785 pearson correlation between variable 1-2, 0.283 correlation between 1-3 and 0.613 is the correlation between variable

Hello, I am trying to do a fit a loglikelihood function with Multivariate distribution via copulas with fitMdvc. The problem is that it doesn't recognize that my beta is a vector of km parameter and when I try to run it it say that the length of my initial values is not the same as the parameter. Can somebody guide me where my mistake is. Thanks, Elisa. #################################

>>>>> Soumen Banerjee <soumen08 at gmail.com> >>>>> on Sat, 21 Apr 2018 17:22:56 +0800 writes: > Hello! I am trying to fit a copula to some data in R and > I get the error mentioned above. This is the code for a > reproducible example - (not really reproducible: You did not set the random seed, so the data is different every time;

Dear R helper, does anyone have an idea on why R.2.4.0 draws the surface for the two command lines below and the next time it renders the error message below for exactly the same command lines: > norm.cop <- normalCopula(0.5) > persp(norm.cop, dcopula) Error in ceiling(length.out) : Non-numeric argument to mathematical function. I will appreciate any help from anyone thanks, Dominique

-- Hello, I am attempting to fit monthly stock returns to possible copula functions using the copula package in R. Below is my code (mat2 is a 2x119 matrix of the two stock returns): my.cop <- normalCopula(param=.3, dim = 2, dispstr = "un") myfit <- fitCopula(mat2,my.cop, start=.65, optim.control= list(NULL), method = "BFGS") myfit Unfortunately, I continue to receive

Hello R users, A beginners question which I could not find the answer to in earler posts. My thought process: Here "z" is a 119 x 15 data matrix Step 1: start at column one, bind every column with column 1 Step2: use the new matrix, "test", in the fitCopula package Step3: store each result in myfit, bind each result to "answer" Step4: return "answer"

I am using curve3d in the emdbook package to graph a gaussian copula density function generated via the copula package. Unfortunately, it appears that emdbook masks dmvnorm from the package mvtnorm in a way that prohibits copula from generating the gaussian copula. (Sounds very confusing!) For example, > library(copula) > f<-function(x,y) dcopula(normalCopula(0),c(x,y)) >

Hello! I am trying to fit a copula to some data in R and I get the error mentioned above. This is the code for a reproducible example - library(copula) data = matrix(data=runif(600),nrow=200,ncol=3) data[,2] = 2*data[,1] data[,3] = 3*data[,1] fr_cop = frankCopula(dim=3) fit_fr_cop = fitCopula(fr_cop,pobs(data),method = "mpl") #Error Here The error says : Error in fitCopula.ml(copula, u

Hi, I am trying to fit t copula to some data, and I am using the following function in the library(QRMlib). Udatac <- apply(datac, 2, edf,adjust=1) tcopulac <- fit.tcopula.rank(Udatac) But the error message come out "Error in fit.tcopula.rank(Udatac) : Non p.s.d. covariance matrix" Could anyone give me some advice? In fact, I am not sure what the "adjust=1" is used for.