Estou reproduzindo um exemplo dos modelos Generalizado, Linear e Misto . Meu MWE está abaixo: Dilution <- c(1/128, 1/64, 1/32, 1/16, 1/8, 1/4, 1/2, 1, 2, 4) NoofPlates <- rep(x=5, times=10) NoPositive <- c(0, 0, 2, 2, 3, 4, 5, 5, 5, 5) Data <- data.frame(Dilution, NoofPlates,...
Recentemente, ajustei 4 modelos de regressão múltiplos para os mesmos dados de previsão / resposta. Dois dos modelos eu me encaixo com a regressão de Poisson. model.pois <- glm(Response ~ P1 + P2 +...+ P5, family=poisson(), ...) model.pois.inter <- glm(Response ~ (P1 + P2 +...+ P5)^2,...
Estou trabalhando em um conjunto de dados. Depois de usar algumas técnicas de identificação de modelos, criei um modelo ARIMA (0,2,1). Usei a detectIOfunção no pacote TSAem R para detectar um outlier inovador (IO) na 48ª observação do meu conjunto de dados original. Como faço para incorporar esse...
Entendo que, se eu tenho dois modelos A e B e A está aninhado em B, dados alguns dados, posso ajustar os parâmetros de A e B usando o MLE e aplicar o teste de razão de verossimilhança generalizada de log. Em particular, a distribuição do teste deve ser com n graus de liberdade, onde n é a diferença...
Pelo que tenho lido, entre outros no site do grupo de consultoria de estatística da UCLA, os testes de razão de verossimilhança e testes de wald são bastante semelhantes ao testar se dois modelos de glm mostram uma diferença significativa na adequação de um conjunto de dados (desculpe-me se minha...
Eu tenho um conjunto de dados em que testo diferenças significativas entre três populações com relação a cerca de 50 variáveis diferentes. Faço isso usando os testes de Kruskal-Wallis, por um lado, e por testes de razão de verossimilhança do modelo GLM aninhado (por e sem população como uma...
Estou usando a glmfitfunção no MATLAB. A função retorna apenas o desvio e não a probabilidade do log. Entendo que o desvio é basicamente o dobro da diferença entre as probabilidades de log dos modelos, mas o que não entendo é que estou apenas usando glmfitpara criar um modelo, mas de alguma forma...
Se houver 0 na tabela de contingência e estivermos ajustando modelos poisson / loglinear aninhados (usando a glmfunção de R ) para um teste de razão de verossimilhança, precisamos ajustar os dados antes de ajustar os modelos glm (por exemplo, adicione 1/2 a todos as contagens)? Obviamente, alguns...
Atualmente, estou tentando implementar um método usado em um artigo popular intitulado "Acabei de executar dois milhões de regressões". A idéia básica por trás disso é que existem certos casos em que não é óbvio quais controles devem ser incluídos no modelo. Uma coisa que você pode fazer nesse caso...
Eu tenho dados com um pico duplo que estou tentando modelar e há sobreposição suficiente entre os picos que não posso tratá-los independentemente. Um histograma dos dados pode ser algo como isto: Eu criei dois modelos para isso: um usa duas distribuições de Poisson e o outro usa duas...
Sejam e duas variáveis aleatórias independentes com os respectivos PDFs:XXXYYY f(x;θi)={1θie−x/θi0<x<∞,0<θi<∞0elsewheref(x;θi)={1θie−x/θi0<x<∞,0<θi<∞0elsewheref \left(x;\theta_i \right) =\begin{cases} \frac{1}{\theta_i} e^{-x/ {\theta_i}} \quad...
Estou analisando um determinado conjunto de dados e preciso entender como escolher o melhor modelo que se ajusta aos meus dados. Estou usando R. Um exemplo de dados que tenho é o seguinte: corr <- c(0, 0, 10, 50, 70, 100, 100, 100, 90, 100, 100) Esses números correspondem à porcentagem de...
Estou curioso para saber como o pacote lmerTest em R, especificamente a função "rand", lida com testes de efeitos aleatórios. Considere o exemplo do pdf lmerTest no CRAN que usa o conjunto de dados incorporado "cenouras": #import lme4 package and lmerTest package library(lmerTest) #lmer model...
Eu tenho duas variáveis aleatórias independentes de poisson, e , com e . Desejo testar versus a alternativa .X1X1X_1X2X2X_2X1∼Pois(λ1)X1∼Pois(λ1)X_1 \sim \text{Pois}(\lambda_1)X2∼Pois(λ2)X2∼Pois(λ2)X_2 \sim \text{Pois}(\lambda_2)H0:λ1=λ2H0:λ1=λ2H_0:\, \lambda_1 = \lambda_2H1:λ1≠λ2H1:λ1≠λ2H_1:\,...
Estou trabalhando em um projeto em que observamos o comportamento de uma tarefa (por exemplo, tempo de resposta) e modelamos esse comportamento em função de várias variáveis manipuladas experimentalmente, bem como de várias variáveis observadas (sexo do participante, QI do participante,...