INDELFIT(1)

✨ INDELFIT β€” the girlies are modeling indels ✨

okay so bestie, welcome to indelFit (that's section 1 of the manual, iykyk) πŸ’… she lives in the User Commands girlie ecosystem and she has ONE job: figuring out the parameters of a simple indel model when you already know the indel history. no cap.


🌸 NAME 🌸

indelFit β€” estimate the parameters of a simple indel model, given an indel history. literally that's the whole thing. she's focused. she's booked. she's serving parameter estimation. πŸ’–


πŸ¦‹ DESCRIPTION πŸ¦‹

So here's the vibe: you hand indelFit an indel history and a tree, and she estimates the parameters of a simple indel model for you. no drama, no fuss. πŸŽ€

By default (and this is the important part, don't skip it) she writes estimates of three parameters to stdout:

  • an insertion rate alpha 🌱
  • a deletion rate beta πŸ‚
  • a parameter tau ✨ which is roughly equal to the inverse of the expected indel length β€” modulo the adjustments required to make the probabilities sum to one, bestie, we don't make the rules
Not me explaining where you even get these inputs: an indel history can be obtained with indelHistory, and a tree with phyloFit. 🫢

πŸ’… OPTIONS πŸ’…

  • --alpha, -a <val> β€” starting value of alpha, the insertion rate. default is 0.02. 🌱
  • --beta, -b <val> β€” starting value of beta, the deletion rate. default is 0.04. πŸ‚
  • --tau, -t <val> β€” starting value of tau, the length parameter. default is 0.05. ✨
  • --lnl, -L β€” compute the log likelihood of the model only, without optimizing any parameters. the options above can be used to set the parameters as desired. it's giving "just tell me the score, don't touch anything" 😌
  • --columns, -c β€” output the log likelihood of each column. every single one. she's thorough like that. πŸ“Š
  • --features, -f <file.gff> β€” estimate separate parameters for each type of feature in the given file, as well as for background sites. slay, segmentation queen πŸ‘‘
  • --reference, -r <name> β€” (for use with --features) name of species (node of tree) defining the coordinate frame of features. by default, the coordinate frame of the entire alignment is assumed. 🧭
  • --log, -l <file> β€” write the log of the optimization to the specified file. receipts, bestie, always keep the receipts 🧾
  • --help, -h β€” display this help message and exit. the classic. the icon. the moment. 🫢

and that's the whole page, angel πŸ’œ go forth and estimate those rates, the model believes in you 🌈⭐🧚

this was indelFit version 1.4, dated May 2016 β€” yeah she's a little vintage, but vintage is a vibe, and she still eats πŸ’–