[Robast-commits] r655 - branches/robast-0.9/pkg/RobExtremesBuffer

[noreply at r-forge.r-project.org]

Author: ruckdeschel
Date: 2013-04-18 14:41:15 +0200 (Thu, 18 Apr 2013)
New Revision: 655

Modified:
   branches/robast-0.9/pkg/RobExtremesBuffer/investigateL2derivGEV.R
Log:
enhanced investigateL2derivGEV.R -> now also plots L2deriv[scale] x Lderiv[shape] as function in observations (relevant for integration) -> highly dubious ... 

Modified: branches/robast-0.9/pkg/RobExtremesBuffer/investigateL2derivGEV.R
===================================================================
--- branches/robast-0.9/pkg/RobExtremesBuffer/investigateL2derivGEV.R	2013-04-18 12:39:30 UTC (rev 654)
+++ branches/robast-0.9/pkg/RobExtremesBuffer/investigateL2derivGEV.R	2013-04-18 12:41:15 UTC (rev 655)
@@ -51,6 +51,10 @@
         plot(z, pf0(L2deriv(G)[[1]]@Map[[i]](q(G)(z))),
               type="l",xlab="x",ylab=if(i==1) "shape" else "scale",axes=F)
         tus()}
+  plot2 <- function(){
+        plot(z, pf0(L2deriv(G)[[1]]@Map[[1]](q(G)(z)))*pf0(L2deriv(G)[[1]]@Map[[2]](q(G)(z))),
+              type="l",xlab="x",ylab="shape x scale",axes=F)
+        tus()}
   windows()
   plot(G)
   if(p(G)(0)<=0){
@@ -58,9 +62,10 @@
     plot(G,log="x")
   }
   windows()
-  if(ncol(trafo(G))>1) {par(mfrow=c(2,1))
+  if(ncol(trafo(G))>1) {par(mfrow=c(3,1))
   ploti(1);
   ploti(2)
+  plot2()
   par(mfrow=c(1,1))}else ploti(1)
 }
 plotG(G=NormLocationScaleFamily())