Rutils/maybe-not-useful/pretty.sqrt.r

In manfredo89/ED2io: Manages Input/Output for Ecosystem Demography 2 Model

#==========================================================================================#
#==========================================================================================#
#     Function that creates a pretty scale in square root.                                 #
#------------------------------------------------------------------------------------------#
pretty.sqrt    = function(x,n=10){
   vlevels     = pretty(x=sign(x)*sqrt(abs(x)),n=n)
   nact        = length(vlevels)
   vlevels     = abs(vlevels)*vlevels

   power.base  = 10^(floor(log10(x=abs(vlevels))))
   npow        = length(unique(power.base))
   ndigits     = ceiling(log10(nact/npow))
   vlevels     = ifelse( vlevels==0
                       , 0.
                       , round(vlevels / power.base,digits=ndigits) * power.base
                       )#end ifelse
   return(vlevels)
}#end function
#==========================================================================================#
#==========================================================================================#






#==========================================================================================#
#==========================================================================================#
#     Function that creates a pretty scale in cubic root.                                  #
#------------------------------------------------------------------------------------------#
pretty.cbrt    = function(x,n=10){
   vlevels     = pretty(x=cbrt(x),n=n)
   nact        = length(vlevels)
   vlevels     = vlevels^3

   power.base  = 10^(floor(log10(x=abs(vlevels))))
   npow        = length(unique(power.base))
   ndigits     = ceiling(log10(nact/npow))
   vlevels     = ifelse( vlevels==0
                       , 0.
                       , round(vlevels / power.base,digits=ndigits) * power.base
                       )#end ifelse
   return(vlevels)
}#end function
#==========================================================================================#
#==========================================================================================#