Estou executando um modelo de regressão de Poisson com 1 variável de resposta e 6 covariáveis. A seleção de modelos usando o AIC resulta em um modelo com todas as covariáveis e 6 termos de interação. O BIC, no entanto, resulta em um modelo com apenas 2 covariáveis e sem termos de interação. É...
Estive usando o ajuste de modelo caret, mas depois executei novamente o modelo usando o gbmpacote. Entendo que o caretpacote usa gbme a saída deve ser a mesma. No entanto, apenas um teste rápido usando data(iris)mostra uma discrepância no modelo de cerca de 5% usando RMSE e R ^ 2 como métrica de...
Estou realizando uma regressão de Poisson com o objetivo final de comparar (e tendo a diferença de) as contagens médias previstas entre dois níveis de fator no meu μ^1- μ^2μ^1-μ^2\hat{\mu}_1-\hat{\mu}_2 , mantendo outras variáveis do modelo (que são todos binária ) constante. Eu queria saber se...
Na modelagem de dados de contagem de solicitações em um ambiente de seguro, comecei com Poisson, mas depois notei superdispersão. Um Quasi-Poisson modelou melhor a maior relação média-variância que o Poisson básico, mas notei que os coeficientes eram idênticos nos modelos de Poisson e...
Eu me deparei com várias questões práticas ao modelar dados de contagem de pesquisas experimentais usando um experimento dentro do assunto. Descrevo brevemente o experimento, os dados e o que fiz até agora, seguidos pelas minhas perguntas. Quatro filmes diferentes foram exibidos para uma amostra...
Gostaria de testar em que regressão se encaixa melhor nos meus dados. Minha variável dependente é uma contagem e possui muitos zeros. E eu precisaria de ajuda para determinar qual modelo e família usar (poisson ou quasipoisson ou regressão de poisson inflada a zero) e como testar as...
O mgcvpacote para Rpossui duas funções para ajustar as interações do produto tensorial: te()e ti(). Entendo a divisão básica do trabalho entre os dois (ajustando uma interação não linear versus decompondo essa interação em efeitos principais e uma interação). O que não entendo é o porquê te(x1,...
Corri modelos binomiais negativos e quase-Poisson com base em uma abordagem de teste de hipóteses. Meus modelos finais, usando os dois métodos, têm covariáveis e interações diferentes. Parece que não há padrões quando planto meus resíduos nos dois casos. Assim, eu queria saber qual teste eu...
Estou analisando os dados de rastreamento ocular de um experimento projetado. Uma versão simplificada dos meus dados tem esta aparência (você pode obter os dados de dput () aqui ), head(lookDATA) participant fixationImage fixationCount 1 9 Automobile 81 2 9 Bird 63 3 9 Chair 82 4 9 Dog 64 5 9...
Eu tenho um GLMM do formulário: lmer(present? ~ factor1 + factor2 + continuous + factor1*continuous + (1 | factor3), family=binomial) Quando uso drop1(model, test="Chi"), obtenho resultados diferentes dos que utilizo Anova(model, type="III")na embalagem do carro ou summary(model). Estes dois...
Estou tentando entrar um pouco nas estatísticas, mas estou preso a alguma coisa. Meus dados são os seguintes: Year Number_of_genes 1990 1 1991 1 1993 3 1995 4 Agora, quero construir um modelo de regressão para poder prever o número de genes para um determinado ano com base nos dados. Fiz isso...
Quero calcular as taxas de incidência para apresentar as taxas de risco, a fim de apresentar medidas de risco relativas e absolutas. Vi em outros estudos que essas taxas de incidência podem ser calculadas usando modelos de Poisson com tempo de acompanhamento no modelo como um deslocamento. Então,...
Eu tenho esses dados: set.seed(1) predictor <- rnorm(20) set.seed(1) counts <- c(sample(1:1000, 20)) df <- data.frame(counts, predictor) Fiz uma regressão de Poisson poisson_counts <- glm(counts ~ predictor, data = df, family = "poisson") E uma regressão binomial...
Eu já vi muito que discute se uma regressão básica de Poisson é uma versão aninhada de uma regressão de Poisson inflada a zero. Por exemplo, este site argumenta que sim, uma vez que o último inclui parâmetros extras para modelar zeros adicionais, mas inclui os mesmos parâmetros de regressão de...
Estou interessado em estimar uma taxa de risco ajustada, análoga a como se estima uma taxa de chances ajustada usando regressão logística. Alguma literatura (por exemplo, isso ) indica que o uso da regressão de Poisson com erros padrão de Huber-White é uma maneira baseada em modelo para fazer...
O Guia para Iniciantes do Zuur 2013 do GLM e GLMM sugere a validação de uma regressão de Poisson, plotando os resíduos de Pearsons em relação aos valores ajustados. Zuur afirma que não devemos ver os resíduos se espalhando à medida que os valores ajustados aumentam, como o gráfico anexado...
Estou examinando tendências (entre 1998 e 2011) nas taxas de mortalidade entre pacientes com doença de Crohn. Cada paciente (caso) foi incluído entre 1998 e 2011. Na inclusão, cada paciente foi comparado a um controle saudável com a mesma idade e sexo. Estou analisando tendências nas taxas de...
Eu estava tentando modelar uma variável dependente de contagem com exposição desigual. Glms clássicos usariam log (exposição) como deslocamento, também o gbm, mas o xgboost não permite o deslocamento até agora ... Tentando encontrar uma desvantagem neste exemplo em validação cruzada ( onde o...
Exposição e deslocamento são duas técnicas frequentemente usadas em regressões de Poisson por atuários para prever a frequência da reivindicação. Pelo meu entendimento, deslocamento e exposição são as mesmas coisas, então não entendo por que existem dois termos para descrever a mesma coisa. É...
No "Guia do praticante para modelos lineares generalizados" no parágrafo 1.83, afirma-se que: "No caso particular de um GLM multiplicativo de Poisson, pode ser demonstrado que a reivindicação de modelagem conta com um termo de compensação igual ao logaritmo da exposição produziu resultados...