[Robast-commits] r684 - in branches/robast-0.9/pkg/13.07.16 - Wrapper for RobAStBase, RobExtremes: . Bilder zum Auswahl
noreply at r-forge.r-project.org
noreply at r-forge.r-project.org
Mon Sep 2 17:26:28 CEST 2013
Author: pupashenko
Date: 2013-09-02 17:26:28 +0200 (Mon, 02 Sep 2013)
New Revision: 684
Added:
branches/robast-0.9/pkg/13.07.16 - Wrapper for RobAStBase, RobExtremes/Bilder zum Auswahl/
branches/robast-0.9/pkg/13.07.16 - Wrapper for RobAStBase, RobExtremes/Bilder zum Auswahl/.Rhistory
branches/robast-0.9/pkg/13.07.16 - Wrapper for RobAStBase, RobExtremes/Bilder zum Auswahl/AuswhalSkript.R
branches/robast-0.9/pkg/13.07.16 - Wrapper for RobAStBase, RobExtremes/Bilder zum Auswahl/ICAllPlot_GEV_optIC.png
branches/robast-0.9/pkg/13.07.16 - Wrapper for RobAStBase, RobExtremes/Bilder zum Auswahl/ICAllPlot_GEV_robIC.png
branches/robast-0.9/pkg/13.07.16 - Wrapper for RobAStBase, RobExtremes/Bilder zum Auswahl/ICAllPlot_GPD_optIC.png
branches/robast-0.9/pkg/13.07.16 - Wrapper for RobAStBase, RobExtremes/Bilder zum Auswahl/ICAllPlot_GPD_robIC.png
branches/robast-0.9/pkg/13.07.16 - Wrapper for RobAStBase, RobExtremes/Bilder zum Auswahl/ICAllPlot_Gamma_optIC.png
branches/robast-0.9/pkg/13.07.16 - Wrapper for RobAStBase, RobExtremes/Bilder zum Auswahl/ICAllPlot_Gamma_robIC.png
branches/robast-0.9/pkg/13.07.16 - Wrapper for RobAStBase, RobExtremes/Bilder zum Auswahl/ICAllPlot_Weibull_optIC.png
branches/robast-0.9/pkg/13.07.16 - Wrapper for RobAStBase, RobExtremes/Bilder zum Auswahl/ICAllPlot_Weibull_robIC.png
branches/robast-0.9/pkg/13.07.16 - Wrapper for RobAStBase, RobExtremes/Bilder zum Auswahl/cniperPointPlot_GEV.png
branches/robast-0.9/pkg/13.07.16 - Wrapper for RobAStBase, RobExtremes/Bilder zum Auswahl/cniperPointPlot_GPD.png
branches/robast-0.9/pkg/13.07.16 - Wrapper for RobAStBase, RobExtremes/Bilder zum Auswahl/cniperPointPlot_Gamma.png
branches/robast-0.9/pkg/13.07.16 - Wrapper for RobAStBase, RobExtremes/Bilder zum Auswahl/cniperPointPlot_Weibull.png
branches/robast-0.9/pkg/13.07.16 - Wrapper for RobAStBase, RobExtremes/Bilder zum Auswahl/infoPlot_GEV_optIC.png
branches/robast-0.9/pkg/13.07.16 - Wrapper for RobAStBase, RobExtremes/Bilder zum Auswahl/infoPlot_GEV_robIC.png
branches/robast-0.9/pkg/13.07.16 - Wrapper for RobAStBase, RobExtremes/Bilder zum Auswahl/infoPlot_GPD_optIC.png
branches/robast-0.9/pkg/13.07.16 - Wrapper for RobAStBase, RobExtremes/Bilder zum Auswahl/infoPlot_GPD_robIC.png
branches/robast-0.9/pkg/13.07.16 - Wrapper for RobAStBase, RobExtremes/Bilder zum Auswahl/infoPlot_Gamma_optIC.png
branches/robast-0.9/pkg/13.07.16 - Wrapper for RobAStBase, RobExtremes/Bilder zum Auswahl/infoPlot_Gamma_robIC.png
branches/robast-0.9/pkg/13.07.16 - Wrapper for RobAStBase, RobExtremes/Bilder zum Auswahl/infoPlot_Weibull_optIC.png
branches/robast-0.9/pkg/13.07.16 - Wrapper for RobAStBase, RobExtremes/Bilder zum Auswahl/infoPlot_Weibull_robIC.png
branches/robast-0.9/pkg/13.07.16 - Wrapper for RobAStBase, RobExtremes/Bilder zum Auswahl/outllyingPlot_GEV.png
branches/robast-0.9/pkg/13.07.16 - Wrapper for RobAStBase, RobExtremes/Bilder zum Auswahl/outllyingPlot_GPD.png
branches/robast-0.9/pkg/13.07.16 - Wrapper for RobAStBase, RobExtremes/Bilder zum Auswahl/outllyingPlot_Gamma.png
branches/robast-0.9/pkg/13.07.16 - Wrapper for RobAStBase, RobExtremes/Bilder zum Auswahl/outllyingPlot_Weibull.png
Log:
- Ticket # 149 - Vignette RobExtremes: Diagnostik
- Auswahl der Bilder
+ Outlyingnessplot
+ cniperPlot
+ ICplot mit Daten
+ infoPlot
Added: branches/robast-0.9/pkg/13.07.16 - Wrapper for RobAStBase, RobExtremes/Bilder zum Auswahl/.Rhistory
===================================================================
--- branches/robast-0.9/pkg/13.07.16 - Wrapper for RobAStBase, RobExtremes/Bilder zum Auswahl/.Rhistory (rev 0)
+++ branches/robast-0.9/pkg/13.07.16 - Wrapper for RobAStBase, RobExtremes/Bilder zum Auswahl/.Rhistory 2013-09-02 15:26:28 UTC (rev 684)
@@ -0,0 +1,512 @@
+#
+# # GEV
+# dev.new()
+# fam = GEVFamily(loc = 0, scale = 1, shape = 0.7)
+# X=distribution(fam)
+# x = r(X)(1000)
+# outlyingPlotWrapper(x,alpha=0.95,fam=fam, alpha.trsp = 100, withCall = FALSE)
+#
+# # Gamma
+# dev.new()
+# fam = GammaFamily()
+# X=distribution(fam)
+# x = r(X)(1000)
+# outlyingPlotWrapper(x,alpha=0.95,fam=fam, alpha.trsp = 100, withCall = FALSE)
+#
+# # Weibull
+# dev.new()
+# fam = WeibullFamily(scale = 0.5, shape = 1)
+# X=distribution(fam)
+# x = r(X)(1000)
+# outlyingPlotWrapper(x,alpha=0.95,fam=fam, alpha.trsp = 100, withCall = FALSE)
+#
+# ### infoPlotWrapper ########################################################################
+# GPD
+fam = GParetoFamily()
+IC <- optIC(model = fam, risk = asCov())
+Y=distribution(fam)
+data = r(Y)(1000)
+dev.new()
+infoPlotWrapper(IC, data, alpha.trsp=100, withCall = FALSE)
+rfam = InfRobModel(fam, ContNeighborhood(0.5))
+ICr <- optIC(model = rfam, risk = asBias())
+dev.new()
+infoPlotWrapper(ICr, data, alpha.trsp=100, withCall = FALSE)
+#
+#
+# # GEV
+# fam = GEVFamily()
+# IC <- optIC(model = fam, risk = asCov())
+# Y=distribution(fam)
+# data = r(Y)(1000)
+# # dev.new()
+# # infoPlotWrapper(IC, alpha.trsp=100, with.legend = TRUE, rescale = TRUE, withCall = TRUE)
+# dev.new()
+# infoPlotWrapper(IC, data, rescale = TRUE, withCall = FALSE)
+#
+# # Gamma
+# fam = GammaFamily()
+# IC <- optIC(model = fam, risk = asCov())
+# Y=distribution(fam)
+# data = r(Y)(1000)
+# # dev.new()
+# # infoPlotWrapper(IC)
+# dev.new()
+# infoPlotWrapper(IC, data, withCall = FALSE)
+#
+# # Weibull
+# fam = WeibullFamily()
+# IC <- optIC(model = fam, risk = asCov())
+# Y=distribution(fam)
+# data = r(Y)(1000)
+# # dev.new()
+# # infoPlotWrapper(IC, alpha.trsp=30, with.legend = TRUE, withCall = FALSE)
+# dev.new()
+# infoPlotWrapper(IC, data, withCall = FALSE)
+#
+# ### cniperPointPlotWrapper #################################################################
+# # GPD
+# dev.new()
+# fam = GParetoFamily()
+# cniperPointPlotWrapper(fam=fam, main = "GPD", lower = 0, upper = 10, withCall = FALSE)
+#
+# # GEV
+# dev.new()
+# fam = GEVFamily()
+# cniperPointPlotWrapper(fam=fam, main = "GEV", lower = 0, upper = 5, withCall = FALSE)
+#
+# # Gamma
+# dev.new()
+# fam = GammaFamily()
+# cniperPointPlotWrapper(fam=fam, main = "Gamma", lower = 0, upper = 5, withCall = FALSE)
+#
+# # Weibull
+# dev.new()
+# fam = WeibullFamily()
+# cniperPointPlotWrapper(fam=fam, main = "Weibull", withCall = FALSE)
+# GPD
+fam = GParetoFamily(loc = 0, scale = 0.5, shape = 0.1)
+IC <- optIC(model = fam, risk = asCov())
+Y=distribution(fam)
+data = r(Y)(1000)
+dev.new()
+infoPlotWrapper(IC, data, alpha.trsp=100, withCall = FALSE)
+rfam = InfRobModel(fam, ContNeighborhood(0.5))
+ICr <- optIC(model = rfam, risk = asBias())
+dev.new()
+infoPlotWrapper(ICr, data, alpha.trsp=100, withCall = FALSE)
+# GPD
+fam = GParetoFamily(loc = 0, scale = 0.1, shape = 0.1)
+IC <- optIC(model = fam, risk = asCov())
+Y=distribution(fam)
+data = r(Y)(1000)
+dev.new()
+infoPlotWrapper(IC, data, alpha.trsp=100, withCall = FALSE)
+rfam = InfRobModel(fam, ContNeighborhood(0.5))
+ICr <- optIC(model = rfam, risk = asBias())
+dev.new()
+infoPlotWrapper(ICr, data, alpha.trsp=100, withCall = FALSE)
+fam = GParetoFamily(loc = 0, scale = 0.5, shape = 0.1)
+IC <- optIC(model = fam, risk = asCov())
+Y=distribution(fam)
+data = r(Y)(1000)
+dev.new()
+infoPlotWrapper(IC, data, alpha.trsp=50, withCall = FALSE)
+rfam = InfRobModel(fam, ContNeighborhood(0.5))
+ICr <- optIC(model = rfam, risk = asBias())
+dev.new()
+infoPlotWrapper(ICr, data, alpha.trsp=50, withCall = FALSE)
+fam = GParetoFamily(loc = 0, scale = 0.5, shape = 0.1)
+IC <- optIC(model = fam, risk = asCov())
+Y=distribution(fam)
+data = r(Y)(1000)
+dev.new()
+infoPlotWrapper(IC, data, alpha.trsp=30, withCall = FALSE)
+rfam = InfRobModel(fam, ContNeighborhood(0.5))
+ICr <- optIC(model = rfam, risk = asBias())
+dev.new()
+infoPlotWrapper(ICr, data, alpha.trsp=30, withCall = FALSE)
+fam = GEVFamily()
+IC <- optIC(model = fam, risk = asCov())
+Y=distribution(fam)
+data = r(Y)(1000)
+dev.new()
+infoPlotWrapper(IC, data, alpha.trsp=100, rescale = TRUE, withCall = FALSE)
+rfam = InfRobModel(fam, ContNeighborhood(0.5))
+ICr <- optIC(model = rfam, risk = asBias())
+dev.new()
+infoPlotWrapper(ICr, data, alpha.trsp=100, rescale = TRUE, withCall = FALSE)
+GEVFamily
+fam = GEVFamily(loc = 0, scale = 0.5, shape = 1)
+IC <- optIC(model = fam, risk = asCov())
+Y=distribution(fam)
+data = r(Y)(1000)
+dev.new()
+infoPlotWrapper(IC, data, alpha.trsp=100, rescale = TRUE, withCall = FALSE)
+rfam = InfRobModel(fam, ContNeighborhood(0.5))
+ICr <- optIC(model = rfam, risk = asBias())
+dev.new()
+infoPlotWrapper(ICr, data, alpha.trsp=100, rescale = TRUE, withCall = FALSE)
+fam = GEVFamily(loc = 0, scale = 1, shape = 0.7)
+IC <- optIC(model = fam, risk = asCov())
+Y=distribution(fam)
+data = r(Y)(1000)
+dev.new()
+infoPlotWrapper(IC, data, alpha.trsp=100, rescale = TRUE, withCall = FALSE)
+rfam = InfRobModel(fam, ContNeighborhood(0.5))
+ICr <- optIC(model = rfam, risk = asBias())
+dev.new()
+infoPlotWrapper(ICr, data, alpha.trsp=100, rescale = TRUE, withCall = FALSE)
+fam = GammaFamily()
+IC <- optIC(model = fam, risk = asCov())
+Y=distribution(fam)
+data = r(Y)(1000)
+dev.new()
+infoPlotWrapper(IC, data, alpha.trsp=100, withCall = FALSE)
+rfam = InfRobModel(fam, ContNeighborhood(0.5))
+ICr <- optIC(model = rfam, risk = asBias())
+dev.new()
+infoPlotWrapper(ICr, data, alpha.trsp=100, withCall = FALSE)
+fam = GammaFamily(scale = 0.1, shape = 0.5)
+IC <- optIC(model = fam, risk = asCov())
+Y=distribution(fam)
+data = r(Y)(1000)
+dev.new()
+infoPlotWrapper(IC, data, alpha.trsp=100, withCall = FALSE)
+rfam = InfRobModel(fam, ContNeighborhood(0.5))
+ICr <- optIC(model = rfam, risk = asBias())
+dev.new()
+infoPlotWrapper(ICr, data, alpha.trsp=100, withCall = FALSE)
+fam = GammaFamily(scale = 0.3, shape = 1)
+IC <- optIC(model = fam, risk = asCov())
+Y=distribution(fam)
+data = r(Y)(1000)
+dev.new()
+infoPlotWrapper(IC, data, alpha.trsp=100, withCall = FALSE)
+rfam = InfRobModel(fam, ContNeighborhood(0.5))
+ICr <- optIC(model = rfam, risk = asBias())
+dev.new()
+infoPlotWrapper(ICr, data, alpha.trsp=100, withCall = FALSE)
+fam = GammaFamily(scale = 0.1, shape = 0.3)
+IC <- optIC(model = fam, risk = asCov())
+Y=distribution(fam)
+data = r(Y)(1000)
+dev.new()
+infoPlotWrapper(IC, data, alpha.trsp=100, withCall = FALSE)
+rfam = InfRobModel(fam, ContNeighborhood(0.5))
+ICr <- optIC(model = rfam, risk = asBias())
+dev.new()
+infoPlotWrapper(ICr, data, alpha.trsp=100, withCall = FALSE)
+fam = GammaFamily(scale = 0.1, shape = 0.7)
+IC <- optIC(model = fam, risk = asCov())
+Y=distribution(fam)
+data = r(Y)(1000)
+dev.new()
+infoPlotWrapper(IC, data, alpha.trsp=100, withCall = FALSE)
+rfam = InfRobModel(fam, ContNeighborhood(0.5))
+ICr <- optIC(model = rfam, risk = asBias())
+dev.new()
+infoPlotWrapper(ICr, data, alpha.trsp=100, withCall = FALSE)
+fam = GammaFamily(scale = 0.1, shape = 0.5)
+IC <- optIC(model = fam, risk = asCov())
+Y=distribution(fam)
+data = r(Y)(1000)
+dev.new()
+infoPlotWrapper(IC, data, alpha.trsp=100, withCall = FALSE)
+rfam = InfRobModel(fam, ContNeighborhood(0.5))
+ICr <- optIC(model = rfam, risk = asBias())
+dev.new()
+infoPlotWrapper(ICr, data, alpha.trsp=100, withCall = FALSE)
+fam = WeibullFamily()
+IC <- optIC(model = fam, risk = asCov())
+Y=distribution(fam)
+data = r(Y)(1000)
+dev.new()
+infoPlotWrapper(IC, data, alpha.trsp=100, withCall = FALSE)
+rfam = InfRobModel(fam, ContNeighborhood(0.5))
+ICr <- optIC(model = rfam, risk = asBias())
+dev.new()
+infoPlotWrapper(ICr, data, alpha.trsp=100, withCall = FALSE)
+IC <- optIC(model = fam, risk = asCov())
+fam = WeibullFamily(scale = 0.1, shape = 0.5)
+IC <- optIC(model = fam, risk = asCov())
+Y=distribution(fam)
+data = r(Y)(1000)
+dev.new()
+infoPlotWrapper(IC, data, alpha.trsp=100, withCall = FALSE)
+rfam = InfRobModel(fam, ContNeighborhood(0.5))
+ICr <- optIC(model = rfam, risk = asBias())
+dev.new()
+infoPlotWrapper(ICr, data, alpha.trsp=100, withCall = FALSE)
+fam = WeibullFamily(scale = 0.2, shape = 0.7)
+IC <- optIC(model = fam, risk = asCov())
+Y=distribution(fam)
+data = r(Y)(1000)
+dev.new()
+infoPlotWrapper(IC, data, alpha.trsp=100, withCall = FALSE)
+rfam = InfRobModel(fam, ContNeighborhood(0.5))
+ICr <- optIC(model = rfam, risk = asBias())
+dev.new()
+infoPlotWrapper(ICr, data, alpha.trsp=100, withCall = FALSE)
+dev.new()
+fam = GParetoFamily()
+cniperPointPlotWrapper(fam=fam, main = "GPD", lower = 0, upper = 10, withCall = FALSE)
+# ### cniperPointPlotWrapper #################################################################
+# GPD
+dev.new()
+fam = GParetoFamily()
+cniperPointPlotWrapper(fam=fam, main = "GPD", lower = 0, upper = 10, withCall = FALSE)
+# GEV
+dev.new()
+fam = GEVFamily()
+cniperPointPlotWrapper(fam=fam, main = "GEV", lower = 0, upper = 5, withCall = FALSE)
+# Gamma
+dev.new()
+fam = GammaFamily()
+cniperPointPlotWrapper(fam=fam, main = "Gamma", lower = 0, upper = 5, withCall = FALSE)
+# Weibull
+dev.new()
+fam = WeibullFamily()
+cniperPointPlotWrapper(fam=fam, main = "Weibull", withCall = FALSE)
+# Gamma
+dev.new()
+fam = GammaFamily(scale = 1, shape = 1)
+cniperPointPlotWrapper(fam=fam, main = "Gamma", lower = 0, upper = 5, withCall = FALSE)
+require(RobExtremes)
+require(distr)
+source("D:/Dropbox/My Mathematics/Researches Misha/Current Research/11.06 - KL PhD/PhD Thesis/Reports for Project/13.07.16 - Wrapper for RobAStBase, RobExtremes/cniperCont.R")
+source("D:/Dropbox/My Mathematics/Researches Misha/Current Research/11.06 - KL PhD/PhD Thesis/Reports for Project/13.07.16 - Wrapper for RobAStBase, RobExtremes/ICAllPlotWrapper.R")
+source("D:/Dropbox/My Mathematics/Researches Misha/Current Research/11.06 - KL PhD/PhD Thesis/Reports for Project/13.07.16 - Wrapper for RobAStBase, RobExtremes/outlyingPlotWrapper.R")
+source("D:/Dropbox/My Mathematics/Researches Misha/Current Research/11.06 - KL PhD/PhD Thesis/Reports for Project/13.07.16 - Wrapper for RobAStBase, RobExtremes/infoPlotWrapper.R")
+source("D:/Dropbox/My Mathematics/Researches Misha/Current Research/11.06 - KL PhD/PhD Thesis/Reports for Project/13.07.16 - Wrapper for RobAStBase, RobExtremes/cniperPointPlotWrapper.R")
+setwd("D:/Dropbox/My Mathematics/Researches Misha/Current Research/11.06 - KL PhD/PhD Thesis/Reports for Project/13.07.16 - Wrapper for RobAStBase, RobExtremes/Bilder zum Auswahl")
+windows.options(width=15, height=12)
+### ICAllPlotWrapper ########################################################################
+# GPD
+fam = GParetoFamily(loc = 0, scale = 0.5, shape = 0.1)
+IC <- optIC(model = fam, risk = asCov())
+Y=distribution(fam)
+y = r(Y)(1000)
+dev.new()
+png(filename = "ICAllPlot_GPD_optIC.png",
+width = 900, height = 700, units = "px", pointsize = 12,
+bg = "white")
+ICAllPlotWrapper(IC, y, alpha.trsp = 100, withCall = FALSE)
+dev.off()
+rfam = InfRobModel(fam, ContNeighborhood(0.5))
+ICr <- optIC(model = rfam, risk = asBias())
+dev.new()
+png(filename = "ICAllPlot_GPD_robIC.png",
+width = 900, height = 700, units = "px", pointsize = 12,
+bg = "white")
+ICAllPlotWrapper(ICr, y, alpha.trsp = 100, withCall = FALSE)
+dev.off()
+# GEV
+fam = GEVFamily(loc = 0, scale = 0.5, shape = 0.1)
+IC <- optIC(model = fam, risk = asCov())
+Y=distribution(fam)
+y = r(Y)(1000)
+dev.new()
+png(filename = "ICAllPlot_GEV_optIC.png",
+width = 900, height = 700, units = "px", pointsize = 12,
+bg = "white")
+ICAllPlotWrapper(IC, y, alpha.trsp = 100, rescale = TRUE, withCall = FALSE)
+dev.off()
+rfam = InfRobModel(fam, ContNeighborhood(0.5))
+ICr <- optIC(model = rfam, risk = asBias())
+dev.new()
+png(filename = "ICAllPlot_GEV_robIC.png",
+width = 900, height = 700, units = "px", pointsize = 12,
+bg = "white")
+ICAllPlotWrapper(ICr, y, alpha.trsp = 100, rescale = TRUE, withCall = FALSE)
+dev.off()
+# Gamma
+fam = GammaFamily(scale = 0.5, shape = 0.5)
+rfam = InfRobModel(fam, ContNeighborhood(0.5))
+IC <- optIC(model = fam, risk = asCov())
+Y=distribution(fam)
+y = r(Y)(1000)
+dev.new()
+png(filename = "ICAllPlot_Gamma_optIC.png",
+width = 900, height = 700, units = "px", pointsize = 12,
+bg = "white")
+ICAllPlotWrapper(IC, y, alpha.trsp = 100, withCall = FALSE)
+dev.off()
+rfam = InfRobModel(fam, ContNeighborhood(0.5))
+ICr <- optIC(model = rfam, risk = asBias())
+dev.new()
+png(filename = "ICAllPlot_Gamma_robIC.png",
+width = 900, height = 700, units = "px", pointsize = 12,
+bg = "white")
+ICAllPlotWrapper(ICr, y, alpha.trsp = 100, withCall = FALSE)
+dev.off()
+# Weibull
+fam = WeibullFamily(scale = 0.5, shape = 1)
+IC <- optIC(model = fam, risk = asCov())
+Y=distribution(fam)
+y = r(Y)(1000)
+dev.new()
+png(filename = "ICAllPlot_Weibull_optIC.png",
+width = 900, height = 700, units = "px", pointsize = 12,
+bg = "white")
+ICAllPlotWrapper(IC, y, alpha.trsp = 100, withCall = FALSE)
+dev.off()
+rfam = InfRobModel(fam, ContNeighborhood(0.5))
+ICr <- optIC(model = rfam, risk = asBias())
+dev.new()
+png(filename = "ICAllPlot_Weibull_robIC.png",
+width = 900, height = 700, units = "px", pointsize = 12,
+bg = "white")
+ICAllPlotWrapper(ICr, y, alpha.trsp = 100, withCall = FALSE)
+dev.off()
+### outllyingPlotWrapper ###################################################################
+# GPD
+dev.new()
+png(filename = "outllyingPlot_GPD.png",
+width = 900, height = 700, units = "px", pointsize = 12,
+bg = "white")
+fam = GParetoFamily(loc = 0, scale = 1, shape = 1)
+X=distribution(fam)
+x = r(X)(1000)
+outlyingPlotWrapper(x,alpha=0.95,fam=fam, alpha.trsp = 100, withCall = FALSE)
+dev.off()
+# GEV
+dev.new()
+png(filename = "outllyingPlot_GEV.png",
+width = 900, height = 700, units = "px", pointsize = 12,
+bg = "white")
+fam = GEVFamily(loc = 0, scale = 1, shape = 0.7)
+X=distribution(fam)
+x = r(X)(1000)
+outlyingPlotWrapper(x,alpha=0.95,fam=fam, alpha.trsp = 100, withCall = FALSE)
+dev.off()
+# Gamma
+dev.new()
+png(filename = "outllyingPlot_Gamma.png",
+width = 900, height = 700, units = "px", pointsize = 12,
+bg = "white")
+fam = GammaFamily()
+X=distribution(fam)
+x = r(X)(1000)
+outlyingPlotWrapper(x,alpha=0.95,fam=fam, alpha.trsp = 100, withCall = FALSE)
+dev.off()
+# Weibull
+dev.new()
+png(filename = "outllyingPlot_Weibull.png",
+width = 900, height = 700, units = "px", pointsize = 12,
+bg = "white")
+fam = WeibullFamily(scale = 0.5, shape = 1)
+X=distribution(fam)
+x = r(X)(1000)
+outlyingPlotWrapper(x,alpha=0.95,fam=fam, alpha.trsp = 100, withCall = FALSE)
+dev.off()
+### infoPlotWrapper ########################################################################
+# GPD
+fam = GParetoFamily(loc = 0, scale = 0.5, shape = 0.1)
+IC <- optIC(model = fam, risk = asCov())
+Y=distribution(fam)
+data = r(Y)(1000)
+dev.new()
+png(filename = "infoPlot_GPD_optIC.png",
+width = 900, height = 700, units = "px", pointsize = 12,
+bg = "white")
+infoPlotWrapper(IC, data, alpha.trsp=30, withCall = FALSE)
+dev.off()
+rfam = InfRobModel(fam, ContNeighborhood(0.5))
+ICr <- optIC(model = rfam, risk = asBias())
+dev.new()
+png(filename = "infoPlot_GPD_robIC.png",
+width = 900, height = 700, units = "px", pointsize = 12,
+bg = "white")
+infoPlotWrapper(ICr, data, alpha.trsp=30, withCall = FALSE)
+dev.off()
+# GEV
+fam = GEVFamily(loc = 0, scale = 1, shape = 0.7)
+IC <- optIC(model = fam, risk = asCov())
+Y=distribution(fam)
+data = r(Y)(1000)
+dev.new()
+png(filename = "infoPlot_GEV_optIC.png",
+width = 900, height = 700, units = "px", pointsize = 12,
+bg = "white")
+infoPlotWrapper(IC, data, alpha.trsp=100, rescale = TRUE, withCall = FALSE)
+dev.off()
+rfam = InfRobModel(fam, ContNeighborhood(0.5))
+ICr <- optIC(model = rfam, risk = asBias())
+dev.new()
+png(filename = "infoPlot_GEV_robIC.png",
+width = 900, height = 700, units = "px", pointsize = 12,
+bg = "white")
+infoPlotWrapper(ICr, data, alpha.trsp=100, rescale = TRUE, withCall = FALSE)
+dev.off()
+# Gamma
+fam = GammaFamily(scale = 0.1, shape = 0.5)
+IC <- optIC(model = fam, risk = asCov())
+Y=distribution(fam)
+data = r(Y)(1000)
+dev.new()
+png(filename = "infoPlot_Gamma_optIC.png",
+width = 900, height = 700, units = "px", pointsize = 12,
+bg = "white")
+infoPlotWrapper(IC, data, alpha.trsp=100, withCall = FALSE)
+dev.off()
+rfam = InfRobModel(fam, ContNeighborhood(0.5))
+ICr <- optIC(model = rfam, risk = asBias())
+dev.new()
+png(filename = "infoPlot_Gamma_robIC.png",
+width = 900, height = 700, units = "px", pointsize = 12,
+bg = "white")
+infoPlotWrapper(ICr, data, alpha.trsp=100, withCall = FALSE)
+dev.off()
+# Weibull
+fam = WeibullFamily(scale = 0.2, shape = 0.7)
+IC <- optIC(model = fam, risk = asCov())
+Y=distribution(fam)
+data = r(Y)(1000)
+dev.new()
+png(filename = "infoPlot_Weibull_optIC.png",
+width = 900, height = 700, units = "px", pointsize = 12,
+bg = "white")
+infoPlotWrapper(IC, data, alpha.trsp=100, withCall = FALSE)
+dev.off()
+rfam = InfRobModel(fam, ContNeighborhood(0.5))
+ICr <- optIC(model = rfam, risk = asBias())
+dev.new()
+png(filename = "infoPlot_Weibull_robIC.png",
+width = 900, height = 700, units = "px", pointsize = 12,
+bg = "white")
+infoPlotWrapper(ICr, data, alpha.trsp=100, withCall = FALSE)
+dev.off()
+# ### cniperPointPlotWrapper #################################################################
+# GPD
+dev.new()
+png(filename = "cniperPointPlot_GPD.png",
+width = 900, height = 700, units = "px", pointsize = 12,
+bg = "white")
+fam = GParetoFamily()
+cniperPointPlotWrapper(fam=fam, main = "GPD", lower = 0, upper = 10, withCall = FALSE)
+dev.off()
+# GEV
+dev.new()
+png(filename = "cniperPointPlot_GEV.png",
+width = 900, height = 700, units = "px", pointsize = 12,
+bg = "white")
+fam = GEVFamily()
+cniperPointPlotWrapper(fam=fam, main = "GEV", lower = 0, upper = 5, withCall = FALSE)
+dev.off()
+# Gamma
+dev.new()
+png(filename = "cniperPointPlot_Gamma.png",
+width = 900, height = 700, units = "px", pointsize = 12,
+bg = "white")
+fam = GammaFamily()
+cniperPointPlotWrapper(fam=fam, main = "Gamma", lower = 0, upper = 5, withCall = FALSE)
+dev.off()
+# Weibull
+dev.new()
+png(filename = "cniperPointPlot_Weibull.png",
+width = 900, height = 700, units = "px", pointsize = 12,
+bg = "white")
+fam = WeibullFamily()
+cniperPointPlotWrapper(fam=fam, main = "Weibull", withCall = FALSE)
+dev.off()
Added: branches/robast-0.9/pkg/13.07.16 - Wrapper for RobAStBase, RobExtremes/Bilder zum Auswahl/AuswhalSkript.R
===================================================================
--- branches/robast-0.9/pkg/13.07.16 - Wrapper for RobAStBase, RobExtremes/Bilder zum Auswahl/AuswhalSkript.R (rev 0)
+++ branches/robast-0.9/pkg/13.07.16 - Wrapper for RobAStBase, RobExtremes/Bilder zum Auswahl/AuswhalSkript.R 2013-09-02 15:26:28 UTC (rev 684)
@@ -0,0 +1,264 @@
+require(RobExtremes)
+require(distr)
+
+source("D:/Dropbox/My Mathematics/Researches Misha/Current Research/11.06 - KL PhD/PhD Thesis/Reports for Project/13.07.16 - Wrapper for RobAStBase, RobExtremes/cniperCont.R")
+
+source("D:/Dropbox/My Mathematics/Researches Misha/Current Research/11.06 - KL PhD/PhD Thesis/Reports for Project/13.07.16 - Wrapper for RobAStBase, RobExtremes/ICAllPlotWrapper.R")
+source("D:/Dropbox/My Mathematics/Researches Misha/Current Research/11.06 - KL PhD/PhD Thesis/Reports for Project/13.07.16 - Wrapper for RobAStBase, RobExtremes/outlyingPlotWrapper.R")
+source("D:/Dropbox/My Mathematics/Researches Misha/Current Research/11.06 - KL PhD/PhD Thesis/Reports for Project/13.07.16 - Wrapper for RobAStBase, RobExtremes/infoPlotWrapper.R")
+source("D:/Dropbox/My Mathematics/Researches Misha/Current Research/11.06 - KL PhD/PhD Thesis/Reports for Project/13.07.16 - Wrapper for RobAStBase, RobExtremes/cniperPointPlotWrapper.R")
+
+setwd("D:/Dropbox/My Mathematics/Researches Misha/Current Research/11.06 - KL PhD/PhD Thesis/Reports for Project/13.07.16 - Wrapper for RobAStBase, RobExtremes/Bilder zum Auswahl")
+windows.options(width=15, height=12)
+# ### ICAllPlotWrapper ########################################################################
+# # GPD
+# fam = GParetoFamily(loc = 0, scale = 0.5, shape = 0.1)
+# IC <- optIC(model = fam, risk = asCov())
+# Y=distribution(fam)
+# y = r(Y)(1000)
+# dev.new()
+# png(filename = "ICAllPlot_GPD_optIC.png",
+# width = 900, height = 700, units = "px", pointsize = 12,
+# bg = "white")
+# ICAllPlotWrapper(IC, y, alpha.trsp = 100, withCall = FALSE)
+# dev.off()
+#
+# rfam = InfRobModel(fam, ContNeighborhood(0.5))
+# ICr <- optIC(model = rfam, risk = asBias())
+# dev.new()
+# png(filename = "ICAllPlot_GPD_robIC.png",
+# width = 900, height = 700, units = "px", pointsize = 12,
+# bg = "white")
+# ICAllPlotWrapper(ICr, y, alpha.trsp = 100, withCall = FALSE)
+# dev.off()
+#
+# # GEV
+# fam = GEVFamily(loc = 0, scale = 0.5, shape = 0.1)
+# IC <- optIC(model = fam, risk = asCov())
+# Y=distribution(fam)
+# y = r(Y)(1000)
+# dev.new()
+# png(filename = "ICAllPlot_GEV_optIC.png",
+# width = 900, height = 700, units = "px", pointsize = 12,
+# bg = "white")
+# ICAllPlotWrapper(IC, y, alpha.trsp = 100, rescale = TRUE, withCall = FALSE)
+# dev.off()
+#
+# rfam = InfRobModel(fam, ContNeighborhood(0.5))
+# ICr <- optIC(model = rfam, risk = asBias())
+# dev.new()
+# png(filename = "ICAllPlot_GEV_robIC.png",
+# width = 900, height = 700, units = "px", pointsize = 12,
+# bg = "white")
+# ICAllPlotWrapper(ICr, y, alpha.trsp = 100, rescale = TRUE, withCall = FALSE)
+# dev.off()
+#
+# # Gamma
+# fam = GammaFamily(scale = 0.5, shape = 0.5)
+# rfam = InfRobModel(fam, ContNeighborhood(0.5))
+# IC <- optIC(model = fam, risk = asCov())
+# Y=distribution(fam)
+# y = r(Y)(1000)
+# dev.new()
+# png(filename = "ICAllPlot_Gamma_optIC.png",
+# width = 900, height = 700, units = "px", pointsize = 12,
+# bg = "white")
+# ICAllPlotWrapper(IC, y, alpha.trsp = 100, withCall = FALSE)
+# dev.off()
+#
+# rfam = InfRobModel(fam, ContNeighborhood(0.5))
+# ICr <- optIC(model = rfam, risk = asBias())
+# dev.new()
+# png(filename = "ICAllPlot_Gamma_robIC.png",
+# width = 900, height = 700, units = "px", pointsize = 12,
+# bg = "white")
+# ICAllPlotWrapper(ICr, y, alpha.trsp = 100, withCall = FALSE)
+# dev.off()
+#
+# # Weibull
+# fam = WeibullFamily(scale = 0.5, shape = 1)
+# IC <- optIC(model = fam, risk = asCov())
+# Y=distribution(fam)
+# y = r(Y)(1000)
+# dev.new()
+# png(filename = "ICAllPlot_Weibull_optIC.png",
+# width = 900, height = 700, units = "px", pointsize = 12,
+# bg = "white")
+# ICAllPlotWrapper(IC, y, alpha.trsp = 100, withCall = FALSE)
+# dev.off()
+#
+# rfam = InfRobModel(fam, ContNeighborhood(0.5))
+# ICr <- optIC(model = rfam, risk = asBias())
+# dev.new()
+# png(filename = "ICAllPlot_Weibull_robIC.png",
+# width = 900, height = 700, units = "px", pointsize = 12,
+# bg = "white")
+# ICAllPlotWrapper(ICr, y, alpha.trsp = 100, withCall = FALSE)
+# dev.off()
+#
+# ### outllyingPlotWrapper ###################################################################
+# # GPD
+# dev.new()
+# png(filename = "outllyingPlot_GPD.png",
+# width = 900, height = 700, units = "px", pointsize = 12,
+# bg = "white")
+# fam = GParetoFamily(loc = 0, scale = 1, shape = 1)
+# X=distribution(fam)
+# x = r(X)(1000)
+# outlyingPlotWrapper(x,alpha=0.95,fam=fam, alpha.trsp = 100, withCall = FALSE)
+# dev.off()
+#
+# # GEV
+# dev.new()
+# png(filename = "outllyingPlot_GEV.png",
+# width = 900, height = 700, units = "px", pointsize = 12,
+# bg = "white")
+# fam = GEVFamily(loc = 0, scale = 1, shape = 0.7)
+# X=distribution(fam)
+# x = r(X)(1000)
+# outlyingPlotWrapper(x,alpha=0.95,fam=fam, alpha.trsp = 100, withCall = FALSE)
+# dev.off()
+#
+# # Gamma
+# dev.new()
+# png(filename = "outllyingPlot_Gamma.png",
+# width = 900, height = 700, units = "px", pointsize = 12,
+# bg = "white")
+# fam = GammaFamily()
+# X=distribution(fam)
+# x = r(X)(1000)
+# outlyingPlotWrapper(x,alpha=0.95,fam=fam, alpha.trsp = 100, withCall = FALSE)
+# dev.off()
+#
+# # Weibull
+# dev.new()
+# png(filename = "outllyingPlot_Weibull.png",
+# width = 900, height = 700, units = "px", pointsize = 12,
+# bg = "white")
+# fam = WeibullFamily(scale = 0.5, shape = 1)
+# X=distribution(fam)
+# x = r(X)(1000)
+# outlyingPlotWrapper(x,alpha=0.95,fam=fam, alpha.trsp = 100, withCall = FALSE)
+# dev.off()
+
+# ### infoPlotWrapper ########################################################################
+# # GPD
+# fam = GParetoFamily(loc = 0, scale = 0.5, shape = 0.1)
+# IC <- optIC(model = fam, risk = asCov())
+# Y=distribution(fam)
+# data = r(Y)(1000)
+# dev.new()
+# png(filename = "infoPlot_GPD_optIC.png",
+# width = 900, height = 700, units = "px", pointsize = 12,
+# bg = "white")
+# infoPlotWrapper(IC, data, alpha.trsp=30, withCall = FALSE)
+# dev.off()
+#
+# rfam = InfRobModel(fam, ContNeighborhood(0.5))
+# ICr <- optIC(model = rfam, risk = asBias())
+# dev.new()
+# png(filename = "infoPlot_GPD_robIC.png",
+# width = 900, height = 700, units = "px", pointsize = 12,
+# bg = "white")
+# infoPlotWrapper(ICr, data, alpha.trsp=30, withCall = FALSE)
+# dev.off()
+#
+# # GEV
+# fam = GEVFamily(loc = 0, scale = 1, shape = 0.7)
+# IC <- optIC(model = fam, risk = asCov())
+# Y=distribution(fam)
+# data = r(Y)(1000)
+# dev.new()
+# png(filename = "infoPlot_GEV_optIC.png",
+# width = 900, height = 700, units = "px", pointsize = 12,
+# bg = "white")
+# infoPlotWrapper(IC, data, alpha.trsp=100, rescale = TRUE, withCall = FALSE)
+# dev.off()
+#
+# rfam = InfRobModel(fam, ContNeighborhood(0.5))
+# ICr <- optIC(model = rfam, risk = asBias())
+# dev.new()
+# png(filename = "infoPlot_GEV_robIC.png",
+# width = 900, height = 700, units = "px", pointsize = 12,
+# bg = "white")
+# infoPlotWrapper(ICr, data, alpha.trsp=100, rescale = TRUE, withCall = FALSE)
+# dev.off()
+#
+# # Gamma
+# fam = GammaFamily(scale = 0.1, shape = 0.5)
+# IC <- optIC(model = fam, risk = asCov())
+# Y=distribution(fam)
+# data = r(Y)(1000)
+# dev.new()
+# png(filename = "infoPlot_Gamma_optIC.png",
+# width = 900, height = 700, units = "px", pointsize = 12,
+# bg = "white")
+# infoPlotWrapper(IC, data, alpha.trsp=100, withCall = FALSE)
+# dev.off()
+#
+# rfam = InfRobModel(fam, ContNeighborhood(0.5))
+# ICr <- optIC(model = rfam, risk = asBias())
+# dev.new()
+# png(filename = "infoPlot_Gamma_robIC.png",
+# width = 900, height = 700, units = "px", pointsize = 12,
+# bg = "white")
+# infoPlotWrapper(ICr, data, alpha.trsp=100, withCall = FALSE)
+# dev.off()
+#
+# # Weibull
+# fam = WeibullFamily(scale = 0.2, shape = 0.7)
+# IC <- optIC(model = fam, risk = asCov())
+# Y=distribution(fam)
+# data = r(Y)(1000)
+# dev.new()
+# png(filename = "infoPlot_Weibull_optIC.png",
+# width = 900, height = 700, units = "px", pointsize = 12,
+# bg = "white")
+# infoPlotWrapper(IC, data, alpha.trsp=100, withCall = FALSE)
+# dev.off()
+#
+# rfam = InfRobModel(fam, ContNeighborhood(0.5))
+# ICr <- optIC(model = rfam, risk = asBias())
+# dev.new()
+# png(filename = "infoPlot_Weibull_robIC.png",
+# width = 900, height = 700, units = "px", pointsize = 12,
+# bg = "white")
+# infoPlotWrapper(ICr, data, alpha.trsp=100, withCall = FALSE)
+# dev.off()
+
+# # ### cniperPointPlotWrapper #################################################################
+# # GPD
+# dev.new()
+# png(filename = "cniperPointPlot_GPD.png",
+# width = 900, height = 700, units = "px", pointsize = 12,
+# bg = "white")
+# fam = GParetoFamily()
+# cniperPointPlotWrapper(fam=fam, main = "GPD", lower = 0, upper = 10, withCall = FALSE)
+# dev.off()
+#
+# # GEV
+# dev.new()
+# png(filename = "cniperPointPlot_GEV.png",
+# width = 900, height = 700, units = "px", pointsize = 12,
+# bg = "white")
+# fam = GEVFamily()
+# cniperPointPlotWrapper(fam=fam, main = "GEV", lower = 0, upper = 5, withCall = FALSE)
+# dev.off()
+#
+# # Gamma
+# dev.new()
+# png(filename = "cniperPointPlot_Gamma.png",
+# width = 900, height = 700, units = "px", pointsize = 12,
+# bg = "white")
+# fam = GammaFamily()
+# cniperPointPlotWrapper(fam=fam, main = "Gamma", lower = 0, upper = 5, withCall = FALSE)
+# dev.off()
+#
+# # Weibull
+# dev.new()
+# png(filename = "cniperPointPlot_Weibull.png",
+# width = 900, height = 700, units = "px", pointsize = 12,
+# bg = "white")
+# fam = WeibullFamily()
+# cniperPointPlotWrapper(fam=fam, main = "Weibull", withCall = FALSE)
+# dev.off()
Added: branches/robast-0.9/pkg/13.07.16 - Wrapper for RobAStBase, RobExtremes/Bilder zum Auswahl/ICAllPlot_GEV_optIC.png
===================================================================
(Binary files differ)
Property changes on: branches/robast-0.9/pkg/13.07.16 - Wrapper for RobAStBase, RobExtremes/Bilder zum Auswahl/ICAllPlot_GEV_optIC.png
___________________________________________________________________
Added: svn:mime-type
+ application/octet-stream
Added: branches/robast-0.9/pkg/13.07.16 - Wrapper for RobAStBase, RobExtremes/Bilder zum Auswahl/ICAllPlot_GEV_robIC.png
===================================================================
(Binary files differ)
Property changes on: branches/robast-0.9/pkg/13.07.16 - Wrapper for RobAStBase, RobExtremes/Bilder zum Auswahl/ICAllPlot_GEV_robIC.png
___________________________________________________________________
Added: svn:mime-type
+ application/octet-stream
Added: branches/robast-0.9/pkg/13.07.16 - Wrapper for RobAStBase, RobExtremes/Bilder zum Auswahl/ICAllPlot_GPD_optIC.png
===================================================================
(Binary files differ)
Property changes on: branches/robast-0.9/pkg/13.07.16 - Wrapper for RobAStBase, RobExtremes/Bilder zum Auswahl/ICAllPlot_GPD_optIC.png
___________________________________________________________________
Added: svn:mime-type
+ application/octet-stream
Added: branches/robast-0.9/pkg/13.07.16 - Wrapper for RobAStBase, RobExtremes/Bilder zum Auswahl/ICAllPlot_GPD_robIC.png
===================================================================
(Binary files differ)
Property changes on: branches/robast-0.9/pkg/13.07.16 - Wrapper for RobAStBase, RobExtremes/Bilder zum Auswahl/ICAllPlot_GPD_robIC.png
___________________________________________________________________
Added: svn:mime-type
+ application/octet-stream
Added: branches/robast-0.9/pkg/13.07.16 - Wrapper for RobAStBase, RobExtremes/Bilder zum Auswahl/ICAllPlot_Gamma_optIC.png
===================================================================
(Binary files differ)
Property changes on: branches/robast-0.9/pkg/13.07.16 - Wrapper for RobAStBase, RobExtremes/Bilder zum Auswahl/ICAllPlot_Gamma_optIC.png
___________________________________________________________________
Added: svn:mime-type
+ application/octet-stream
Added: branches/robast-0.9/pkg/13.07.16 - Wrapper for RobAStBase, RobExtremes/Bilder zum Auswahl/ICAllPlot_Gamma_robIC.png
===================================================================
(Binary files differ)
Property changes on: branches/robast-0.9/pkg/13.07.16 - Wrapper for RobAStBase, RobExtremes/Bilder zum Auswahl/ICAllPlot_Gamma_robIC.png
___________________________________________________________________
Added: svn:mime-type
+ application/octet-stream
Added: branches/robast-0.9/pkg/13.07.16 - Wrapper for RobAStBase, RobExtremes/Bilder zum Auswahl/ICAllPlot_Weibull_optIC.png
===================================================================
(Binary files differ)
Property changes on: branches/robast-0.9/pkg/13.07.16 - Wrapper for RobAStBase, RobExtremes/Bilder zum Auswahl/ICAllPlot_Weibull_optIC.png
___________________________________________________________________
Added: svn:mime-type
+ application/octet-stream
Added: branches/robast-0.9/pkg/13.07.16 - Wrapper for RobAStBase, RobExtremes/Bilder zum Auswahl/ICAllPlot_Weibull_robIC.png
===================================================================
(Binary files differ)
Property changes on: branches/robast-0.9/pkg/13.07.16 - Wrapper for RobAStBase, RobExtremes/Bilder zum Auswahl/ICAllPlot_Weibull_robIC.png
___________________________________________________________________
Added: svn:mime-type
+ application/octet-stream
Added: branches/robast-0.9/pkg/13.07.16 - Wrapper for RobAStBase, RobExtremes/Bilder zum Auswahl/cniperPointPlot_GEV.png
===================================================================
(Binary files differ)
Property changes on: branches/robast-0.9/pkg/13.07.16 - Wrapper for RobAStBase, RobExtremes/Bilder zum Auswahl/cniperPointPlot_GEV.png
___________________________________________________________________
Added: svn:mime-type
+ application/octet-stream
Added: branches/robast-0.9/pkg/13.07.16 - Wrapper for RobAStBase, RobExtremes/Bilder zum Auswahl/cniperPointPlot_GPD.png
===================================================================
(Binary files differ)
Property changes on: branches/robast-0.9/pkg/13.07.16 - Wrapper for RobAStBase, RobExtremes/Bilder zum Auswahl/cniperPointPlot_GPD.png
___________________________________________________________________
Added: svn:mime-type
+ application/octet-stream
[TRUNCATED]
To get the complete diff run:
svnlook diff /svnroot/robast -r 684
More information about the Robast-commits
mailing list