GBEPFIT(1)

✨ GBEPFIT ✨

okay so bestie, welcome to GBEPFIT(1) β€” that's section 1 of the User Commands manual, i.e. the "you type this in a terminal" zone. Sit down, hydrate, we're fitting a density πŸ’…


🌸 NAME 🌸

gbepfit β€” fit a symmetric power exponential density. that's it, that's the whole name, and she's iconic


πŸ¦‹ SYNOPSIS πŸ¦‹

gbepfit [options]

πŸ’– DESCRIPTION πŸ’–

So what she does: gbepfit fits a symmetric power exponential density, and you get to pick your whole method β€” either maximum-likelihood or the method of moments. Pick your fighter, both are valid 🫢


πŸ‘‘ OPTIONS πŸ‘‘

  • -O β€” output type (default 0). The vibe check on what you actually get back:
    • 0 β€” parameter b a m and log-likelihood
    • 1 β€” the estimated distribution function computed on the provided points
    • 2 β€” the estimated density function computed on the provided points
    • 3 β€” parameters b a m and their standard errors
  • -x β€” set initial conditions b,a,m (default 2,1,0). this is literally where she starts her little journey 🧚
  • -m β€” the mode is not estimated but is set to the value provided. you're locking her in, no negotiation βœ‹
  • -s β€” number of intervals to explore at each iteration (default 10)
  • -V β€” verbosity level (default 0). how much she talks to you:
    • 0 β€” just the final result
    • 1 β€” headings and summary table
    • 2 β€” intermediate steps results
    • 3 β€” intermediate steps internals
    • 4+ β€” details of optim. routine
  • -M β€” active estimation steps. The value is the sum of (default 7). yes bestie, it's additive, you add them up:
    • 1 β€” initial estimation based on method of moments
    • 2 β€” global optimization not considering lack of smoothness in m
    • 4 β€” local optimization taking non-smoothness in m into consideration
  • -G β€” set global optimization options. Fields are step,tol,iter,eps,msize,algo. Empty field implies default (default .1,1e-2,100,1e-3,1e-5,3)
  • -I β€” set local optimization options. Fields are step,tol,iter,eps,msize,algo. Empty field implies default (default .01,1e-4,200,1e-4,1e-5,5)
  • -F β€” input fields separators (default " \t")
  • -h β€” this help. hi πŸ‘‹

🌈 the optimization parameters, explained 🌈

  • step β€” initial step size of the searching algorithm
  • tol β€” line search tolerance iter: maximum number of iterations
  • eps β€” gradient tolerance : stopping criteria ||gradient||<eps
  • msize β€” simplex max size : stopping criteria ||max edge||<msize
  • algo β€” optimization methods: 0 Fletcher-Reeves, 1 Polak-Ribiere, 2 Broyden-Fletcher-Goldfarb-Shanno, 3 Steepest descent, 4 Nelder-Mead simplex, 5 Broyden-Fletcher-Goldfarb-Shanno v.2, 6 Nelder-Mead simplex ver. 2, 7 Nelder-Mead simplex rnd init.

πŸ“ EXAMPLES πŸ“

gbepfit -m 1 -M 6 <file

what this does: estimate a and b with m=1 and skipping initial method of moments estimation. she ate, no cap πŸ”₯


πŸ’œ AUTHOR πŸ’œ

Written by Giulio Bottazzi. king behavior πŸ‘‘


πŸ™ˆ REPORTING BUGS πŸ™ˆ

Found a bug? not the vibe. Report bugs to <gbutils@googlegroups.com>

Package home page <http://cafim.sssup.it/~giulio/software/gbutils/index.html>


Copyright Β© 2001-2018 Giulio Bottazzi This program is free software; you can redistribute it and/or modify it under the terms of the GNU General Public License (version 2) as published by the Free Software Foundation;

This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details.

okay real talk for one sec πŸ’€ β€” that "WITHOUT ANY WARRANTY" bit is not a joke and not decoration, it means exactly what it says, so don't come for her if it doesn't do the thing.

and that's the whole page, bestie ✨ this one's gbepfit version 6.0.beta2, dated February 2019 β€” so if you're reading this in the future, yes, she's an elder but she still serves. mwah πŸ’‹