search for: rmvdc

Displaying 7 results from an estimated 7 matches for "rmvdc".

Did you mean: mvdc
2013 May 03
0
Empirica Copula
...paramMargins=list(list(rate=0.5),list(rate=0.5))) # Use own data for bivariate CDF construction myCDF2<- mvdc(gmb, margins=c("SAR","per"), paramMargins=list(list(rate=.5),list(rate=.5))) # Generate (bivariate) random numbers from that, and visualize x <- rMvdc(1000, myCDF2) And I get error message everytime as: > x <- rMvdc(1000, myCDF2) Error in qSAR(x, rate = 0.5) : unused argument(s) (rate = 0.5) It works fine with myCDF and generate bivariate data: x <- rMvdc(1000, myCDF2) But my problem is simulated data (using myCDF) does not show the...
2008 Jul 30
2
Sampling two exponentials
Hi all, I am going to sample two variables from two exponential distributions, but I want to specify a covariance structure between these two variables. Is there any way to do it in R? Or is there a "Multivariate Exponential" thing corresponding to the multivariate normal? Thanks in advance. Sincerely, Yanwei Zhang Department of Actuarial Research and Modeling Munich Re America Tel:
2010 Jun 10
0
error message fitting tcopula
...ellipCopula(family = "t", param = 0.5,dim = 2, dispstr = "un", df = 8) myCop.t   myMvd <- mvdc(copula = myCop.t, margins = c("gamma", "gamma"), paramMargins = list(list(shape = 1.5, scale = 38), list(shape = 1.7, scale = 50))) myMvd n <- 200 dat <- rmvdc(myMvd, n)   mm <- apply(dat, 2, mean) vv <- apply(dat, 2, var) rho <- rcorr(dat,type="spearman")[[1]]; round(rho,2) rbind(mm,vv)   b1.0 <- c(mm[1]^2/vv[1], vv[1]/mm[1]) b2.0 <- c(mm[2]^2/vv[2], vv[2]/mm[2]) a.0  <- sin(cor(dat[, 1], dat[, 2], method = "kendall"...
2010 Jun 10
0
error message in fitting tcopula
...opula(family = "t", param = 0.5,dim = 2, dispstr = "un", df = 8) myCop.t myMvd <- mvdc(copula = myCop.t, margins = c("gamma", "gamma"), paramMargins = list(list(shape = 1.5, scale = 38), list(shape = 1.7, scale = 50))) myMvd n <- 200 dat <- rmvdc(myMvd, n) mm <- apply(dat, 2, mean) vv <- apply(dat, 2, var) rho <- rcorr(dat,type="spearman")[[1]]; round(rho,2) rbind(mm,vv) b1.0 <- c(mm[1]^2/vv[1], vv[1]/mm[1]) b2.0 <- c(mm[2]^2/vv[2], vv[2]/mm[2]) a.0 <- sin(cor(dat[, 1], dat[, 2], method = "kend...
2006 Oct 08
2
Generating bivariate or multivariate data with known parameter values
Greetings, I'm interested in generating data from various bivariate or mulitivariate distributions (e.g. gamma, t, etc), where I can specify the parameter values, including the correlations among the variables. I haven't been able to dig anything up on the faq, but I probably missed something. A nudge in the right direction would be appreciated. David --
2007 Jul 16
3
R and Copula
hi, first I want to say that I'm new here, and new with copula and R. That is the reason why I'm writing, if somebody can help me. I have to make an example of Copula. On internet I've found this forum and that copula can calculate with R. Can somebody help me with the thing how can I start and where can read about these stuffs. Thank to all who can help! -- View this message
2006 May 12
3
Maximum likelihood estimate of bivariate vonmises-weibulldistribution
Thanks Dimitris!!! That's much clearer now. Still have a lot of work to do this weekend to understand every bit but your code will prove very useful. Cheers, Aziz -----Original Message----- From: Dimitrios Rizopoulos [mailto:Dimitris.Rizopoulos at med.kuleuven.be] Sent: May 12, 2006 4:35 PM To: Chaouch, Aziz Subject: RE: [R] Maximum likelihood estimate of bivariate