[Distr-commits] r937 - branches/distr-2.6/pkg/distrMod/R

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

Author: ruckdeschel
Date: 2014-06-27 00:52:18 +0200 (Fri, 27 Jun 2014)
New Revision: 937

Modified:
   branches/distr-2.6/pkg/distrMod/R/returnlevelplot.R
Log:
noch eine Kleinigkeit rausgenommen

Modified: branches/distr-2.6/pkg/distrMod/R/returnlevelplot.R
===================================================================
--- branches/distr-2.6/pkg/distrMod/R/returnlevelplot.R	2014-06-26 22:18:28 UTC (rev 936)
+++ branches/distr-2.6/pkg/distrMod/R/returnlevelplot.R	2014-06-26 22:52:18 UTC (rev 937)
@@ -97,19 +97,25 @@
     ycl <- p2rl(yc.o)
 
     ### extend range somewhat
-    xyall <- sort(unique(c(yc.o,x,
+#    pyn <- p(y)(10^(seq(-1, 3.75 + log10(npy), by = 0.1)))
+    xyall <- force(sort(unique(c(yc.o,x,
                     q(y)(c(seq(0.01, 0.09, by = 0.01),(1:9)/10,
-                         0.95, 0.99, 0.995, 0.999)),
-                         10^(seq(-1, 3.75 + log10(npy), by = 0.1))
-                         )))
-    rxyall <- (max(xyall)-min(xyall))*0.6
+                         0.95, 0.99, 0.995, 0.999))
+                         ))))
+    rxyall  <- (max(xyall)-min(xyall))*0.6
     rxymean <- (max(xyall)+min(xyall))/2
 
-    xyallc <- seq(rxymean-rxyall,rxymean+rxyall, length.out=300)
+    xyallc  <- seq(from=rxymean-rxyall,to=rxymean+rxyall, length.out=300)
+    print(xyallc)
+    pxyall  <- p(y)(xyallc)
+    print(pxyall)
+
     pxyallc <- p2rl(xyallc)
-    xyallc <- xyallc[pxyallc>0.00001 & pxyallc<0.99999]
-    pxyallc <- pxyallc[pxyallc>0.00001 & pxyallc<0.99999]
+     xyallc <-  xyallc[pxyall>0.00001 & pxyall<0.99999]
+    pxyallc <- pxyallc[pxyall>0.00001 & pxyall<0.99999]
 
+    print(cbind(pxyallc,xyallc))
+
     if("support" %in% names(getSlots(class(y))))
        ycl <- sort(jitter(ycl, factor=jit.fac))