Estou ajustando vários modelos do formulário .. glm(DV ~ I(1/IV), family = Gamma(link = "log") .. e estou procurando maneiras de comparar os modelos obtidos para diferentes variáveis. Gostaria de saber se o valor alfa é de alguma utilidade prática? Para as três parcelas abaixo, os valores alfa...
A regressão logística binomial possui assíntotas superiores e inferiores de 1 e 0, respectivamente. No entanto, os dados de precisão (apenas como exemplo) podem ter assíntotas superiores e inferiores muito diferentes de 1 e / ou 0. Posso ver três soluções possíveis para isso: Não se preocupe se...
Eu tenho dois grupos de 10 participantes que foram avaliados três vezes durante um experimento. Para testar as diferenças entre os grupos e nas três avaliações, executei um ANOVA de desenho misto 2x3 com group(controle, experimental), time(primeiro, segundo, três) e group x time. Ambos timee...
Da Wikipedia O modelo linear geral (GLM) é um modelo linear estatístico. Pode ser escrito como 1 Y=XB+U,Y=XB+U, \mathbf{Y} = \mathbf{X}\mathbf{B} + \mathbf{U}, onde YYY é uma matriz com uma série de medidas multivariadas, XXX é uma matriz que pode ser uma matriz de design, BBB é uma matriz...
Minha pergunta geral é: por que usar em bayesglmvez de outros métodos de classificação? Nota: Estou interessado apenas em previsão. Eu tenho uma quantidade decente de dados (~ 100.000 obs). Sinto que o tamanho da amostra é grande o suficiente para que os parâmetros de uma regressão logística...
Estou tentando criar um modelo usando regressão binomial negativa (GLM binomial negativo). Eu tenho um tamanho de amostra relativamente pequeno (maior que 300) e os dados não são redimensionados. Percebi que existem duas maneiras de medir a qualidade do ajuste - uma é desvio e a outra é a...
Existe alguma maneira de executar o LASSO com regressão binomial negativa em R? Estou executando uma regressão binomial negativa no meu conjunto de dados porque os dados são muito dispersos para impor a regressão de poisson. Enquanto isso, também estou enfrentando algum problema de...
Eu tenho o seguinte modelo linear generalizado. O objeto glmDVé modelado como uma proporção de sucessos em relação ao total de tentativas. Os objetos x_isão variáveis contínuas. Como isso se parece na notação matemática? winp.glm = glm(glmDV ~ x1 + x2 + x3 + x4 + x5 + x6 + x7, data=myData,...
Então, aqui estou estudando modelos lineares generalizados. Sei que essa pergunta é bastante ingênua e simples, mas não sei exatamente por que a função canônica do link é tão útil. Alguém poderia me fornecer uma intuição sobre esse
Eu tenho dados que descrevem com que frequência um evento ocorre durante uma hora ("número por hora", nph) e quanto tempo os eventos duram ("duração em segundos por hora", dph). Estes são os dados originais: nph <- c(2.50000000003638, 3.78947368414551, 1.51456310682008, 5.84686774940732,...
Eu quero ajustar um modelo linear por R com family=binomial(link="identity"), no entanto, a família binomial não tem vínculo de identidade. O que devo
Por que um GLM prevê a média e não o modo de um sinal? Isso não contradiz a própria base do GLM, ou seja, a máxima probabilidade? As equações a serem resolvidas para os parâmetros do modelo em um GLM são baseadas na maximização da probabilidade, conforme descrito pela distribuição de probabilidade...
Eu tenho uma regressão logística binária com apenas um preditor de fator fixo binário. A razão de eu não fazer isso como um quadrado de Chi ou o teste exato de Fisher é que também tenho vários fatores aleatórios (existem vários pontos de dados por indivíduo e os indivíduos estão em grupos, embora...
Fechadas. Esta questão está fora de tópico . No momento, não está aceitando respostas. Deseja melhorar esta pergunta? Atualize a pergunta para que ela esteja no tópico de Validação cruzada. Fechado há 12 meses . Estou tentando ajustar o modelo linear...
Primeiro, deixe-me dar uma base; Resumirei minhas perguntas no final. A distribuição Beta, parametrizada por sua média e ϕ , possui Var ( Y ) = V ( μ ) / ( ϕ + 1 ) , onde V ( μ ) = μ ( 1 - μ ) é a função de variância.μμ\muϕϕ\phiVar(Y)=V(μ)/(ϕ+1)Var(Y)=V(μ)/(ϕ+1)\operatorname{Var}(Y) =...
Considere a seguinte regressão binomial: # Create some data set.seed(10) n <- 500 x <- runif(n,0,100) y <- x + rnorm(n,sd=100) < 0 # Fit a binomial regression model model <- glm(y ~ x, family="binomial") summary(model) A summaryfunção retorna um valor-p de 1.03e-05. Ao usar...
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...
Eu estava lendo esta página em Princeton.edu . Eles estão executando uma regressão logística (com R). Em algum momento, eles calculam a probabilidade de obter um desvio residual maior do que aquele que obtiveram em uma com graus de liberdade iguais aos graus de liberdade do modelo. Copiando e...
Agora, estou passando para os GLMs depois dos modelos padrão. No modelo padrão, y = Xb + epsilon e epsilon é assumido como sendo normalmente distribuído. Isso significa que podemos escrever y - Xb = epsilon e então podemos minimizar o lhs usando alguma norma adequada, dada a suposição de...
Considere o seguinte exemplo simples: library( rms ) library( lme4 ) params <- structure(list(Ns = c(181L, 191L, 147L, 190L, 243L, 164L, 83L, 383L, 134L, 238L, 528L, 288L, 214L, 502L, 307L, 302L, 199L, 156L, 183L), means = c(0.09, 0.05, 0.03, 0.06, 0.07, 0.07, 0.1, 0.1, 0.06, 0.11, 0.1,...