Na regressão linear (perda ao quadrado), usando matriz, temos uma notação muito concisa para o objetivo minimizar ∥ A x - b ∥ 2minimizar __UMAx-b__2\text{minimize}~~ \|Ax-b\|^2 Onde é a matriz de dados, x são os coeficientes eb é a resposta.UMAUMAAxxxbbb Existe uma notação matricial...
Gostaria de combinar as saídas do lmer (realmente glmer) com um exemplo binomial de brinquedo. Eu li as vinhetas e acredito que entendo o que está acontecendo. Mas aparentemente eu não. Depois de ficar preso, consertei a "verdade" em termos dos efeitos aleatórios e fui atrás da estimativa dos...
Estou fazendo uma regressão de riscos proporcionais de Cox em R usando coxph, que inclui muitas variáveis. Os resíduos de Martingale parecem ótimos, e os resíduos de Schoenfeld são ótimos para QUASE todas as variáveis. Existem três variáveis cujos resíduos de Schoenfeld não são planos e a...
O fenômeno de 'super dispersão' em um GLM surge sempre que usamos um modelo que restringe a variação da variável de resposta, e os dados exibem uma variação maior do que a restrição do modelo permite. Isso ocorre geralmente ao modelar dados de contagem usando um Poisson GLM, e pode ser...
Li várias vezes que efeitos aleatórios (BLUPs / modos condicionais para, digamos, sujeitos) não são parâmetros de um modelo linear de efeitos mistos, mas podem ser derivados dos parâmetros estimados de variância / covariância. Por exemplo, Reinhold Kliegl et al. (2011) afirmam: Efeitos...
Estou interessado em alterar as hipóteses nulas usando glm()em R. Por exemplo: x = rbinom(100, 1, .7) summary(glm(x ~ 1, family = "binomial")) testa a hipótese de que p=0.5p=0.5p = 0.5 . E se eu quiser alterar o valor nulo para ppp = algum valor arbitrário glm()? Sei que isso também pode ser...
Então, eu tenho brincado com SVMs e me pergunto se isso é uma boa coisa a fazer: Eu tenho um conjunto de recursos contínuos (0 a 1) e um conjunto de recursos categóricos que converti em variáveis fictícias. Nesse caso em particular, codifico a data da medição em uma variável dummy: Existem 3...
Um dos problemas que sempre tive com modelos mistos é descobrir visualizações de dados - do tipo que pode acabar em um papel ou pôster - depois de obter os resultados. No momento, estou trabalhando em um modelo de efeitos mistos de Poisson com uma fórmula semelhante à seguinte: a <-...
Tenho dados sobre os funcionários de uma grande empresa italiana há dez anos e gostaria de ver como a diferença de gênero nos ganhos entre homens e mulheres mudou ao longo do tempo. Para esse propósito, eu OLS agrupado: onde é o log de ganhos por ano, inclui covariáveis que diferem por...
Quais são os benefícios de especificar uma estrutura de covariância em um GLM (em vez de tratar todas as entradas fora da diagonal na matriz de covariância como zero)? Além de refletir o que se sabe dos dados, ele faz melhorar a qualidade do ajuste? melhorar a precisão preditiva em dados...
Em um artigo recente, eu instalei modelos de efeitos fixos de três vias. Como um dos fatores não era significativo (p> 0,1), eu o removi e reinstalei o modelo com dois efeitos fixos e uma interação. Acabei de receber comentários dos árbitros, para citar: Esse tempo não foi um fator...
Quero incluir o tempo gasto fazendo alguma coisa (semanas de amamentação, por exemplo) como uma variável independente em um modelo linear. No entanto, algumas observações não se envolvem no comportamento. Codificá-los como 0 não é realmente certo, porque 0 é qualitativamente diferente de qualquer...
Quero desenvolver um modelo de previsão (Cox PH) para mortalidade por todas as causas em um conjunto de dados de participantes dos quais (quase) todos morreram no final do acompanhamento (por exemplo, 1 ano). Em vez de prever o risco absoluto de morrer em um determinado ponto no tempo, eu gostaria...
Na análise de dados em painel, usei modelos de vários níveis com efeitos aleatórios / mistos para lidar com problemas de auto-correlação (ou seja, as observações são agrupadas em indivíduos ao longo do tempo) com outros parâmetros adicionados para ajustar algumas especificações de tempo e choques...
Alguém pode me dizer como R estimar o ponto de interrupção em um modelo linear por partes (como um parâmetro fixo ou aleatório), quando eu também precisar estimar outros efeitos aleatórios? Incluí um exemplo de brinquedo abaixo que se encaixa em uma regressão de taco de hóquei / taco quebrado com...
Gostaria de saber a diferença entre análise de dados em painel e análise de modelo misto. Que eu saiba, tanto os dados do painel quanto os modelos mistos usam efeitos fixos e aleatórios. Se sim, por que eles têm nomes diferentes? Ou eles são sinônimos? Eu li o seguinte post, que descreve a...
Eu só queria saber se a regressão de Poisson tem um termo de erro? Uma regressão de Poisson pode ter efeitos aleatórios e um termo de erro? Estou confuso sobre esse ponto. Na regressão logística, não há termo de erro porque sua variável de resultado é binária. Esse é o único modelo de glm que não...
Eu apenas comecei a aprender a usar Stan e rstan. A menos que eu sempre tenha ficado confuso sobre como o JAGS / BUGS funcionou, pensei que você sempre tivesse que definir uma distribuição anterior de algum tipo para cada parâmetro no modelo a ser extraído. Parece que você não precisa fazer isso em...
Eu gostaria de empregar os dados de contagem como covariáveis ao ajustar um modelo de regressão logística. Minha pergunta é: Eu viole alguma suposição dos modelos logísticos (e, mais geralmente, do linear generalizado), empregando variáveis inteiras não negativas de contagem como variáveis...
Ao especificar um modelo de efeitos mistos cruzados, estou tentando incluir interações. No entanto, recebo a seguinte mensagem de erro: Error in lme.formula(rate ~ nozzle, random = ~nozzle | operator, data = Flow) : nlminb problem, convergence error code = 1 message = iteration limit reached...