Changes in version 0.3.8 - Relatório Quarto: distance (cluster::daisy só na Gower), discriminant (MASS::lda ou MASS::qda), logistic (stats::glm ou nnet::multinom) e os jackknife registram no resultado a ferramenta que rodou (atributo trama_ferramentas). - Relatório Quarto exportado: multi/pca (resumo, autovalores e cargas), multi/fa (resumo e cargas de padrão), multi/dist (resumo e, até 12 indivíduos, a matriz) e multi/cluster (resumo e os membros de cada grupo) ganham report — tabelas em Markdown no lugar do print cru da lista. Pede trama.data (>= 0.4.4). Changes in version 0.3.7 - Templates regenerados com o formato atual do documento (saídas nomeadas e grupos do trama 0.5.10). Pede trama (>= 0.5.10). Changes in version 0.3.6 - Cards que são gráfico seguem o tema padrão do projeto (antes ficavam no tema embutido, escuro). Requer trama 0.5.9. Changes in version 0.3.5 - Testes passam no R CMD check do pacote instalado (antes liam arquivos da árvore do repositório). Changes in version 0.3.1 - Parâmetros que só valem para certa escolha de outro parâmetro agora declaram trama::tr_when() e somem do card quando não se aplicam. Exige trama (>= 0.2.0). Changes in version 0.3.0 Versão sobe de 0.2.0 para 0.3.0: a integração com a coesão das coleções muda o que os blocos devolvem (a discriminante e a logística saem como models/fit, e classificar, confundir, ROC, PR e coeficientes passam aos blocos da trama.models). Integração 9.1b (coesão × main) - multi/logistic ganha metodo = "firth" (main); os coeficientes (com IC perfilado e p da razão de verossimilhanças, ou Wald) saem pelo models/coefficients, e fluxos com multi/logistic_coefficients (inclusive com nivel, v2, ou intervalo, v4) abrem migrados para ele. - multi/roc e multi/pr_curve migram para models/roc e models/pr_curve, com DeLong, Hand & Till e as curvas por grupo portadas; multi/confusion migra para models/confusion com tabela. - Jackknife por grupo (grupo) e confianca (v2) ficam aqui; as duas migrações (da coleção e por versão do nó) convivem. 0.2.0 Versão sobe de 0.1.0: o resultado padrão mudou (intervalo perfilado como padrão na multi/logistic_coefficients). Blocos novos - multi/pr_curve: curva precisão-revocação dos classificadores (LDA/QDA e logística), por deixa-um-fora como a multi/roc, com a precisão média (AP = Σ ΔR·P) e a área de Davis & Goadrich (2006, doi:10.1145/1143844.1143874) integrada em forma fechada (Keilwagen, Grosse & Grau 2014, doi:10.1371/journal.pone.0092209), e a prevalência como linha do acaso (Saito & Rehmsmeier 2015, doi:10.1371/journal.pone.0118432). Três ou mais grupos: cada um contra os outros. A mesma conta da ml/pr_curve, copiada (a coleção não depende da trama.ml). Validação: exemplo à mão da ml; yardstick::average_precision 1.4.0 a 1e-10 e PRROC::pr.curve 1.4 (auc.integral) a 1e-8 no pima por deixa-um-fora (AP 0,721, área 0,719, acaso 0,333) e nos três cultivares dos vinhos com empates. - multi/mardia: teste de normalidade multivariada de Mardia (1970, doi:10.1093/biomet/57.3.519) — assimetria b1,p (χ², e com a correção de amostra pequena de Mardia 1974) e curtose b2,p (z), na tabela toda ou dentro de cada grupo. Covariância de divisor n, como no artigo e em MVN::mardia(use_population = TRUE). Validação: psych::mardia reescalado (divisor n − 1) a 1e-10 em quatro tabelas; estatísticas e p iguais ao código do MVN::mardia 6.3 a 1e-13. Os pressupostos de normalidade da discriminante, do M de Box, da fatorial e do KMO/Bartlett passam a apontá-lo. Logística de Firth - multi/logistic ganha metodo = c("ml", "firth") (padrão ml, sem mudança de resultado). firth (só binária) maximiza a verossimilhança penalizada de Firth (1993, doi:10.1093/biomet/80.1.27): coeficientes finitos mesmo com separação. multi/logistic_coefficients (versão 2) devolve a coluna nova intervalo (Wald ou perfilado); no Firth, IC da verossimilhança penalizada perfilada e p da razão de verossimilhanças penalizadas (Heinze & Schemper 2002, doi:10.1002/sim.1047). EP pela inversa da informação de Fisher em β̂ (o logistf usa (X'W(1 + h)X)⁻¹ e dá EP menores). Validação: logistf 1.26.1 no sex2 — coeficientes a 1e-6, limites e p a 1e-4 (também com confiança 0,9 e num exemplo com separação completa); nos vinhos B × C quase separados, onde o logistf não converge em quatro dos oito limites, os limites conferem pela definição (χ²₁ do perfil a 1e-6) e batem com os quatro em que ele converge. Intervalo da AUC - multi/roc (versão 2): a AUC sai com o intervalo de DeLong, DeLong & Clarke-Pearson (1988, doi:10.2307/2531595), no subtítulo (dois grupos) ou na legenda (cada grupo contra os outros), com o param novo confianca (0,95). Validação: pROC::ci.auc(method = "delong") 1.18 a 1e-8 no aSAH (três escores, com empates; também a 90%) e nas probabilidades de deixa-um-fora da logística do pima (AUC 0,849, IC 0,816–0,882). Métricas da matriz de confusão - multi/confusion ganha tabela = c("matriz", "métricas") (padrão matriz, sem mudança). métricas devolve acurácia, acurácia balanceada (Brodersen et al. 2010, doi:10.1109/ICPR.2010.764), kappa de Cohen (1960, doi:10.1177/001316446002000104) e precisão, revocação e F1 por grupo (Sokolova & Lapalme 2009); precisão de grupo nunca previsto é NA. Validação: contas à mão; exemplo 20/5/10/15 (κ = 0,4, o exemplo da página "Cohen's kappa" da Wikipédia — não é fonte primária); irr::kappa2 e psych::cohen.kappa a 1e-12 nas previsões cruzadas da LDA do iris. Aviso de n pequeno no Mardia - multi/mardia: com menos de 20 linhas na tabela ou num grupo, avisa (tr_multi_warning_small_n) que o teste é assintótico e pouco confiável em n pequeno. 20 é a regra do MVN::mardia para usar a assimetria de amostra pequena (conferida no código-fonte do MVN), não um limiar da literatura. Resultado igual. z de Wald e p perfilado na multi/logistic_coefficients - A ajuda explica que o z é sempre de Wald e, com intervalo = "perfilado", o p_valor é da razão de verossimilhanças (penalizadas no Firth): os dois podem discordar. O EP do Firth (inversa da informação de Fisher em β̂) não foi conferido no texto de Heinze & Schemper (2002), sem acesso aberto; o resumo só diz que Wald existe e que os testes da razão de verossimilhanças penalizadas são muitas vezes preferíveis. IC da AUC indisponível em vez de recusa - multi/roc (versão 4): com AUC 0 ou 1 (variância de DeLong zero) ou menos de dois casos numa classe, a AUC sai com "IC indisponível" e a nota vai na legenda do gráfico, como a ml/roc (IC NA com nota). Antes, AUC 1 dava o intervalo degenerado 1–1 (a setosa na LDA da iris) e uma classe pequena derrubava o gráfico inteiro com tr_multi_error_small_group. Com três ou mais grupos, o grupo sem IC não afeta as curvas dos outros. Validação: bordas à mão (AUC 0, 1, classe com um caso); pROC inalterado nos casos regulares. Jackknife por grupo com grupos desiguais - Os quatro jackknifes (versão 3), com grupo de tamanhos diferentes: viés, corrigida e pseudovalor saem NA, e o intervalo centra na estimativa da amostra toda (θ̂ ± t(G − 1) · EP), com a coluna nota. Antes centrava na corrigida de G réplicas, que supõe grupos iguais (Shao & Tu 1995, sec. 2.3); a correção ponderada para grupos desiguais não foi implementada por falta de oráculo conferido. O EP JK1 não muda (vale com grupos desiguais; survey). Menos de 5 grupos avisa (tr_multi_warning_few_groups). Muda o resultado só com grupos desiguais; grupos iguais e deixa-um-fora ficam iguais. Validação: contas à mão (5 grupos iguais e desiguais; escala log). confianca nos jackknifes - multi/jackknife_pca, _fa, _discriminant e _logistic (versão 2): o param do nível do intervalo passa de nivel a confianca, como na multi/logistic_coefficients e na multi/roc. Fluxo salvo com nivel abre migrado (tr_node(migracoes = )). Na chamada R, nivel = deixa de ser aceito. Resultado igual. confianca na multi/logistic_coefficients - multi/logistic_coefficients (versão 3): o param do nível do intervalo passa de nivel a confianca, a convenção da coleção (multi/roc) e das irmãs. Fluxo salvo com nivel abre migrado: o nó declara a migração da v3 (tr_node(migracoes = )), o núcleo troca nivel por confianca ao abrir o flow ou colar o template e avisa "migrado de v2 para v4". Sem alias: o nó exige que o fn seja a função exportada e que cada argumento dela seja um param declarado, então nivel = também deixa de ser aceito na chamada R. Resultado igual. Intervalo perfilado na logística ML - multi/logistic_coefficients (versão 4) ganha intervalo = c("perfilado", "Wald"), com perfilado como padrão: na logística ML binária, o IC da verossimilhança perfilada (Venables & Ripley 2002, sec. 7.2, doi:10.1007/978-0-387-21706-2) e o p da razão de verossimilhanças. O de Wald supõe a log-verossimilhança quadrática; Hosmer, Lemeshow & Sturdivant (2013, sec. 1.4, doi:10.1002/9781118548387) preferem o perfilado em amostra pequena, e com n grande os dois coincidem. Muda o resultado padrão (os limites e o p da ML binária; coeficientes e EP iguais); intervalo = "Wald" reproduz o anterior. Na multinomial fica Wald (o nnet::multinom não tem perfil e não há implementação de referência para conferir um próprio), dito na coluna intervalo. No Firth, Wald passa a ser opção também. Validação: confint do glm (o perfil do MASS, no stats desde o R 4.4) com grade fina a 1e-4 (2e-6 observado) no pima inteiro e numa subamostra de 64, a 95% e 90%; a grade padrão erra ~6e-4 na subamostra pela spline, e os limites do trama conferem pela definição (desvio perfilado = χ²₁ a 1e-6); p contra drop1(test = "LRT") a 1e-10. No pima (glicose, imc, pedigree), a razão de chances do pedigree 3,71 tem IC 1,88–7,45 perfilado e 1,86–7,39 de Wald. AUC multiclasse de Hand & Till - multi/roc (versão 3): com três ou mais grupos, o subtítulo traz a AUC multiclasse M de Hand & Till (2001, Machine Learning 45(2):171–186, doi:10.1023/A:1010920819831) — média, sobre os pares de grupos, de [A(i|j) + A(j|i)]/2, cada A só com os casos do par. Não depende das proporções dos grupos. Sem intervalo. Validação: pROC::multiclass.roc 1.19.1 a 1e-10 (LDA da iris por deixa-um-fora, M = 0,998133, diferença 1e-16; logística dos vinhos; 4 grupos com empates) e conta à mão. Jackknife por grupo - Os quatro blocos de jackknife (multi/jackknife_pca, _fa, _discriminant, _logistic) ganham o param grupo (padrão em branco, sem mudança): com uma coluna de conglomerado, cada réplica tira o grupo inteiro (jackknife apagar-um-grupo; Shao & Tu 1995, The Jackknife and Bootstrap, doi:10.1007/978-1-4612-0795-5; Kott 2001, Journal of Official Statistics 17(4):521–526, sem DOI no Crossref). EP = √((G − 1)/G · Σ(θ₍g₎ − θ̄)²), intervalo com t(G − 1); viés e pseudovalores com G (exatos com grupos iguais). Para dados em conglomerados, onde deixar uma linha fora subestima a variância. Recusa faltante no grupo, menos de 2 grupos e, nos classificadores, a própria coluna do grupo previsto. Validação: EP igual ao do survey 4.5 com os grupos como UPAs e réplicas JK1 (as.svrepdesign(type = "JK1")) a 1e-10 na média (0,1244884, diferença 6e-17) e na razão (0,006755280, diferença 1e-17), 9 grupos desiguais; exemplo de 3 grupos à mão; réplica da PCA igual ao ajuste sem o grupo.