Find agent skills by outcome

- `SANN` (simulated annealing): convergence code is always 0 — it never fails. `maxit` = total function evals (default 10000)

no other stopping criterion.

- Default method is Nelder-Mead (no gradients

robust but slow). Poor for 1D — use `Brent` or `optimize()`.

- Returns object with `predict`

`print`

- `predict.glm` with `type = response` gives predictions on the response scale (e.g.

probabilities for logistic); `type = link` (default) gives on the link scale.

- `lowess` parameter is `f` (span

default 2/3); `loess` parameter is `span` (default 0.75).

- `lowess` is the older function; returns `list(x

y)` — cannot predict at new points.

- Direction: `direction = both` (default)

`forward`

Creative

- Self-starting models (`SSlogis`

`SSasymp`

Creative

- `aov` is a wrapper around `lm` that stores extra info for balanced ANOVA. For unbalanced designs

Type I SS (sequential) are computed — order of terms matters.

- For Type III SS

use `car::Anova()` or set contrasts to `contr.sum`/`contr.helmert`.

- Quasi-families (`quasibinomial`

`quasipoisson`) allow overdispersion — no AIC is computed.

- `binomial` accepts response as: 0/1 vector

logical

- `weights` in `glm` means prior weights (not frequency weights) — for frequency weights

use the cbind trick or offset.

- `anova(glm_obj

test = Chisq)` for deviance-based tests; `F` is invalid for non-Gaussian families.

- Common families: `binomial(link = logit)`

`poisson(link = log)`

- Default `family = gaussian(link = identity)` — `glm()` with no `family` silently fits OLS (same as `lm`

but slower and with deviance-based output).

- Column names from `model.matrix` can be surprising: e.g.

`factorLevelName` concatenation.

- `terms` object attributes: `order` (interaction order per term)

`intercept`

Creative

- `model.matrix` creates the design matrix including dummy coding. Default contrasts: `contr.treatment` for unordered factors

`contr.poly` for ordered.

- `I()` is required to use arithmetic operators literally: `y ~ x + I(x^2)`. Without `I()`

`^` means interaction crossing.

Coding & Debugging

- `URLencode(url

reserved = FALSE)` by default does NOT encode reserved chars (`/`

- Captures output from `cat`

`print`

> Non-obvious behaviors

gotchas

- `as.is = TRUE` (recommended): keeps characters as character

not factor.