[Distr-commits] r190 - branches/distr-2.0/pkg/distrMod/R

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

Author: stamats
Date: 2008-07-22 20:09:51 +0200 (Tue, 22 Jul 2008)
New Revision: 190

Modified:
   branches/distr-2.0/pkg/distrMod/R/SimpleL2ParamFamilies.R
Log:
corrected some minor bugs

Modified: branches/distr-2.0/pkg/distrMod/R/SimpleL2ParamFamilies.R
===================================================================
--- branches/distr-2.0/pkg/distrMod/R/SimpleL2ParamFamilies.R	2008-07-22 13:18:07 UTC (rev 189)
+++ branches/distr-2.0/pkg/distrMod/R/SimpleL2ParamFamilies.R	2008-07-22 18:09:51 UTC (rev 190)
@@ -372,11 +372,11 @@
 ExpScaleFamily <- function(rate = 1, trafo){ 
     L2ScaleFamily(loc = 0, scale = 1/rate, name = "Exponential scale family", 
                   centraldistribution = Exp(rate = 1),
-                  LogDeriv = function(x)  x-1,  
+                  LogDeriv = function(x)  1,
                   L2derivDistr.0 = (Exp(rate = 1)-1)*rate,
                   FisherInfo.0 = 1, 
                   distrSymm = NoSymmetry(), 
-                  L2derivSymm = FunSymmList(NonSymmetric()), 
+                  L2derivSymm = FunSymmList(EvenSymmetric(SymmCenter = 1/rate)), 
                   L2derivDistrSymm = DistrSymmList(NoSymmetry()),
                   trafo = trafo)
 }
@@ -389,8 +389,8 @@
     L2ScaleFamily(loc = 0, scale = exp(meanlog),  
                   name = "lognormal scale family", 
                   centraldistribution = Lnorm(meanlog = 0, sdlog = sdlog),
-                  LogDeriv = function(x) log(x)/sdlog^2/x,  
-                  L2derivDistr.0 = Norm(mean=0, sd=1/sdlog),
+                  LogDeriv = function(x) log(x)/x/sdlog^2 + 1/x,
+                  L2derivDistr.0 = Norm(mean=0, sd=1/sdlog^2),
                   FisherInfo.0 = 1/sdlog^2, 
                   distrSymm = NoSymmetry(), 
                   L2derivSymm = FunSymmList(NonSymmetric()), 
@@ -405,10 +405,10 @@
 GumbelLocationFamily <- function(loc = 0, scale = 1, trafo){ 
     L2LocationFamily(loc = loc, scale = scale, 
                      name = "Gumbel location family", 
-                     centraldistribution = Gumbel(loc = 0, scale = scale),
-                     LogDeriv = function(x)  (1-exp(-(x-loc)/scale))/scale,  
+                     centraldistribution = Gumbel(loc = 0),
+                     LogDeriv = function(x) 1 - exp(-x),
                      L2derivDistr.0 = (1 - Exp(rate = 1))/scale,
-                     FisherInfo.0 = 1/scale^2, 
+                     FisherInfo.0 = 1, 
                      distrSymm = NoSymmetry(), 
                      L2derivSymm = FunSymmList(NonSymmetric()), 
                      L2derivDistrSymm = DistrSymmList(NoSymmetry()),