Suponha que você ajuste uma regressão linear / logística g(y)=a0+a1⋅x1+a2⋅x2g(y)=a0+a1⋅x1+a2⋅x2g(y) = a_0 + a_1\cdot x_1 + a_2\cdot x_2 , com o objetivo de uma estimativa imparcial de a1a2a1a2\frac{a_1}{a_2} . Você está muito confiante de que tantoa1a1a_1quantoa2a2a_2 são muito positivos em relação...
Eu estava procurando maneiras de fazer um teste de razão de verossimilhança em R para comparar os ajustes do modelo. Eu primeiro o codifiquei, depois encontrei a anova()função padrão e também lrtest()no lmtestpacote. Quando verifiquei, no entanto, anova()sempre produz um valor p ligeiramente...
Para um dado problema de inferência, sabemos que uma abordagem bayesiana geralmente difere na forma e nos resultados de uma abordagem fequentista. Os freqüentistas (geralmente me incluem) frequentemente apontam que seus métodos não requerem um prévio e, portanto, são mais "orientados a dados" do...
Estou experimentando o algoritmo da máquina de aumento de gradiente através do caretpacote em R. Usando um pequeno conjunto de dados de admissões de faculdade, executei o seguinte código: library(caret) ### Load admissions dataset. ### mydata <-
Parece haver muita confusão na comparação entre usar glmnetdentro caretpara procurar uma lambda ideal e usar cv.glmnetpara fazer a mesma tarefa. Muitas perguntas foram feitas, por exemplo: Modelo de classificação train.glmnet vs. cv.glmnet? Qual é a maneira correta de usar glmnet com...
Recentemente, perguntas de navegação aleatória provocaram a memória de comentários de um de meus professores há alguns anos, alertando sobre o uso de proporções em modelos de regressão. Então comecei a ler sobre isso, levando eventualmente a Kronmal 1993. Quero ter certeza de que estou...
Atualmente, estou revisando alguns trabalhos e encontrei o seguinte, o que me parece errado. Dois modelos mistos são montados (em R) usando o lmer. Os modelos não são aninhados e são comparados por testes de razão de verossimilhança. Em suma, aqui está um exemplo reproduzível do que...
Eu tenho vários artigos que apresentam "OR" com IC95% (intervalos de confiança). Quero estimar a partir dos artigos o valor de P para o OR observado. Para isso, preciso de uma suposição sobre a distribuição OR. Que distribuição posso assumir / usar com
No exemplo a seguir > m = matrix(c(3, 6, 5, 6), nrow=2) > m [,1] [,2] [1,] 3 5 [2,] 6 6 > (OR = (3/6)/(5/6)) #1 [1] 0.6 > fisher.test(m) #2 Fisher's Exact Test for Count Data data: m p-value = 0.6699 alternative hypothesis: true odds ratio is not equal to 1 95 percent confidence...
O teste da razão de verossimilhança e a AIC são ferramentas para escolher entre dois modelos e ambos se baseiam na probabilidade de log. Mas, por que o teste da razão de verossimilhança não pode ser usado para escolher entre dois modelos não aninhados enquanto o AIC
Como um aplicativo de exemplo, considere seguir duas propriedades dos usuários do Stack Overflow: contagem de visualizações de reputação e perfil . Espera-se que, para a maioria dos usuários, esses dois valores sejam proporcionais: os usuários com alta reputação atraem mais atenção e, portanto,...
Y=μ+σGY=μ+σGY=\mu+\sigma Gμμ\muWWWGGGRnRn\mathbb{R}^nFFFH0:{μ∈U}H0:{μ∈U}H_0\colon\{\mu \in U\}U⊂WU⊂WU \subset W Como podemos saber que essa estatística fornece o teste mais poderoso paraH0(talvez depois de descartar casos particulares incomuns)? Isso não deriva do teorema de Neyman-Pearson, porque...
Alguém poderia me dizer quais são as condições de regularidade para a distribuição assintótica do teste da Razão de Verossimilhança? Em todo lugar que olho, está escrito 'Sob as condições de regularidade' ou 'sob as regularidades probabilísticas'. Quais são exatamente as condições? Que o primeiro...
Mais especificamente, por que os testes da razão de verossimilhança têm uma distribuição assintoticamente se os modelos estão aninhados, mas esse não é mais o caso dos modelos não aninhados? Entendo que isso se segue do teorema de Wilks, mas, infelizmente, não entendo sua prova...
Num artigo recente, Norton et al. (2018) [ 1 ] afirmam que[ 1 ][1]^{[1]} Odds ratio diferentes do mesmo estudo não podem ser comparados quando os modelos estatísticos que resultam em estimativas de odds ratio possuem variáveis explicativas diferentes, pois cada modelo possui um fator de escala...
Estou estudando como construir um intervalo de confiança de 95% para odds ratio a partir dos coeficientes obtidos na regressão logística. Portanto, considerando o modelo de regressão logística, log(p1−p)=α+βxlog(p1−p)=α+βx \log\left(\frac{p}{1 - p}\right) = \alpha + \beta x \newcommand{\var}{\rm...
Deixe- fff , ggg e hhh ser densidades e suponha que você tem xi∼hxi∼hx_i \sim h , i∈Ni∈Ni \in \mathbb{N} . O que acontece com a razão de verossimilhança ∏i=1nf(xi)g(xi)∏i=1nf(xi)g(xi) \prod_{i=1}^n \frac{f(x_i)}{g(x_i)} comon→∞n→∞n \rightarrow \infty? (Ele converge? Para quê?) Por exemplo,...
Suponha que onde X i ∼ N ( 0 , σ 2 ) são independentes.X=X1+X2+⋯+XnX=X1+X2+⋯+Xn X = X_1 + X_2+\cdots+ X_n Xi∼N(0,σ2)Xi∼N(0,σ2)X_i \sim N(0,\sigma^2) Minha pergunta é: qual distribuição faz Z=X2X21+X22+⋯+X2nZ=X2X12+X22+⋯+Xn2 Z = \frac{X^2}{X_1^2 + X_2^2 + \cdots + X_n^2} Segue? Eu sei daqui que a...
Seja uma distribuição uniforme e uma distribuição normal. O que pode ser dito sobre ? Existe uma distribuição para isso?Y XXXXYYYXYXY\frac X Y Descobri que a proporção de dois normais com zero médio é
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,...