[R-SIG-Finance] fit.gausscopula function

Sylvain BARTHELEMY barth at tac-financial.com
Tue Aug 21 12:39:19 CEST 2007


Dear Xiaochen,

I think that there are conflicting versions of packages on your machine, as
the R code you provide works well on two of my machines, with R 2.5 and
QRMlib v2.5.1. 

You should update all your R packages by using the "update packages" menu.

Regards.

---
Sylvain Barthélémy
Research Director, TAC
www.tac-financial.com | www.sylbarth.com

-----Message d'origine-----
De : r-sig-finance-bounces at stat.math.ethz.ch
[mailto:r-sig-finance-bounces at stat.math.ethz.ch] De la part de Xiaochen Sun
Envoyé : mardi 21 août 2007 11:15
À : r-sig-finance at stat.math.ethz.ch
Objet : [R-SIG-Finance] fit.gausscopula function

I do not know why, fit.gausscopula( ) function in QRMlib doesn't work.....
 
If you run the example code in help file, Here is the message I've got:
 
> data(ftse100);
> data(smi);
> TS1 <- cut(ftse100, "1990-11-08", "2004-03-25");
> TS1Augment <- alignDailySeries(TS1, method="before");
> TS2Augment <- alignDailySeries(smi, method="before");
> INDEXES.RAW <- merge(TS1Augment,TS2Augment);
> #Cleanup:
> rm(TS1, TS1Augment, TS2Augment);
> INDEXES <- mk.returns(INDEXES.RAW);
> PARTIALINDEXES <- cut(INDEXES, "1993-12-31", "2003-12-31");
> #Now create a data matrix from the just-created timeSeries 
> data <- seriesData(PARTIALINDEXES);
> #Keep only the data items which are non-zero for both smi and ftse100
> data <- data[data[,1]!=0 & data[,2] !=0,];
> # Construct pseudo copula data. The 2nd parameter is MARGIN=2 
> #when applying to columns and 1 applied to rows. Hence this says to
> #apply the 'edf()' empirical distribtion function() to the columns
> #of the data. 
> Udata <- apply(data,2,edf,adjust=1);
> mod.gauss <- fit.gausscopula(Udata); 
Error in dmnorm(Qdata, rep(0, d), P, logvalue = TRUE) : 
        unused argument(s) (logvalue = TRUE)

However it did work once:
X3X2 <- cbind(x3,x2)
U3U2 <- apply(X3X2,2,edf,adjust=1)
mod.gaussU3U2 <- fit.gausscopula(U3U2)
mod.gaussU3U2

> mod.gaussU3U2[1]
$P
          [,1]      [,2]
[1,] 1.0000000 0.2770152
[2,] 0.2770152 1.0000000
I can not figure out:(((( 
 
Could anyone tell me why?
 
Cheers
Mc

_______________________________________________
R-SIG-Finance at stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-sig-finance
-- Subscriber-posting only. 
-- If you want to post, subscribe first.



More information about the R-SIG-Finance mailing list