similar to: Fit Student Copula

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

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?

Hi R-users, I would like to use the copula package.? I? the package plus the mvtnorm and try to run the example given, but I got the following message: install.packages(repos=NULL,pkgs="c:\\Tinn-R\\copula_0.8-3.zip") norm.cop <- normalCopula(c(0.5, 0.6, 0.7), dim = 3, dispstr = "un") t.cop <- tCopula(c(0.5, 0.3), dim = 3, dispstr = "toep", df = 2, df.fixed =

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

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

Hi all, Many thank in advance for helping me.? I tried to fit Expectation Maximization algorithm for mixture data. I must used one of numerical method to maximize my function. I built my code but I do not know how to make the optim function run over a different value of the parameters.? That is, For E-step I need to get the value of mixture weights based on the current (initial) values of

Hi, How can I generate a random signal that's correlated with a given signal at a given correlation (say 0.7)? I've been looking at rmvnorm etc but don't seem to figure it out. Thanks ----- Yasir H. Kaheil Columbia University -- View this message in context: http://www.nabble.com/generating-a-random-signal-with-a-known-correlation-tp18932541p18932541.html Sent from the R help

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. #################################

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

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

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"

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 >

Hello everyone, Here is the setup. z is a 119 x 15 matrix, m_index is a 119 x 5 matrix What I am trying to do is return the results from fitCopula by sequentially binding all 15 columns of z to the first column of m_index, (cbind(z[,1],m_index[,1]),(cbind(z[,2],m_index[,1]), etc. Unfortunately, my code below only binds z[,1] and m_index[,1] and return the same result 14 times. Any ideas on

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,

Hi, I'm a new user of R. I'm using package copula implemented in R. I want to know how to test goodness of fit of student's t copula for 3 dimensional cases using real world data. In manual gof test has been performed on the copula families generated from rcopula function that i understood as random data generated from copula? How to perform the test with real/observed data matrix.

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