[R] quantile curves
Aishe Brei
aishe.brei at hotmail.de
Thu Sep 23 11:21:00 CEST 2010
Dear List,
I fitted a bivariate extreme value distribution to my matrix (data). Now
I would like to add ‘quantile curves’ to the plot, similar to
those calculated with
qcbvnonpar()
But I need the
probability of the upper right corner instead of the lower left corner.
(So that the curves in the plot below would look like they were turned
by 180 degrees).
Is there a package that could help me? Or is
there an option within the ‘evd’ package that I didn’t see?
My
data:
a<-c(6.8044, 6.4422, 6.0900, 6.1778, 6.2778, 6.3978,
6.2156, 6.1512, 6.2712, 6.3112, 6.1590, 6.2368, 6.1968, 6.1746, 6.2202,
6.5836, 6.1970, 6.2870, 6.1126, 6.3460, 6.0616, 6.2016, 6.3294, 6.1494,
6.0350, 6.2006, 6.1106, 6.3506, 6.3706, 6.0484, 6.3140, 6.0618, 6.0418,
6.4618, 6.1996, 6.3196, 6.0996, 6.2174, 6.1574, 6.0774, 6.1874, 6.0230,
6.0608, 6.3386, 6.0686, 6.3042, 6.2042, 6.6242, 6.0720, 6.1276, 6.1410,
6.1010, 6.2688, 6.5788, 6.3566, 6.1766, 6.1444, 6.6200, 6.0178, 6.1856,
6.1634)
b<-c(20, 16, 12, 16, 10, 26, 11, 5, 12, 16, 6, 12, 9, 15, 13, 22, 10,
15, 10, 9, 5, 8, 30, 8, 5, 8, 7, 23, 18, 5, 8, 6, 5, 28, 6, 9, 7, 8, 6,
5, 6, 6, 6, 15, 6, 8, 10, 26, 7, 5, 8, 8, 6, 11, 8, 8, 11, 16, 6, 12,
12)
data<-cbind(a,b)
bifit<-fbvevd(data,model="log")
plot(bifit,which=5)
Thanks a lot,
Aishe
More information about the R-help
mailing list