Considere um problema estatístico em que você tem uma responsevariável que deseja descrever condicional em uma explanatoryvariável e uma nestedvariável, em que a variável aninhada surge apenas como uma variável significativa para valores específicos da variável explicativa . Nos casos em que a...
Devo confessar que nunca ouvi falar desse termo em nenhuma das minhas aulas, graduação ou pós-graduação. O que significa uma regressão logística ser bayesiana? Estou procurando uma explicação com uma transição da logística regular para a logística bayesiana semelhante à seguinte: Esta é a equação...
Estou tentando modelar as intensidades médias de parasitas que afetam um hospedeiro em R usando um modelo binomial negativo. Eu recebo 50 ou mais avisos que dizem: In dpois(y, mu, log = TRUE) : non-integer x = 251.529000 Como posso lidar com isso? Meu código fica assim: mst.nb =...
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,...
Eu construí um modelo de glm em R e o testei usando um grupo de teste e treinamento, por isso estou confiante de que funciona bem. Os resultados de R são: Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) -2.781e+00 1.677e-02 -165.789 < 2e-16 *** Coeff_A 1.663e-05 5.438e-06...
No GLM, assumindo um escalar e para a distribuição subjacente com pdf Pode ser mostrado que . Se a função de link g (\ cdot) satisfizer o seguinte, g (\ mu) = \ theta = X '\ beta em que X' \ beta é o preditor linear, então g (\ cdot) será chamada de função de link canônico para este...
Eu tenho 100.000 observações (9 variáveis indicadoras fictícias) com 1000 positivos. A regressão logística deve funcionar bem neste caso, mas a probabilidade de corte me intriga. Na literatura comum, escolhemos um ponto de corte de 50% para prever 1s e 0s. Não posso fazer isso, pois meu modelo...
Qual é a diferença estatística entre fazer uma regressão linear em R com o formulaconjunto para em y ~ x + 0vez de y ~ x? Como interpreto esses dois resultados
Estou tentando executar um modelo para estimar quão bem as doenças catastróficas, como tuberculose, AIDS etc afetam os gastos com hospitalização. Eu tenho "por custo de hospitalização" como variável dependente e vários marcadores individuais como variáveis independentes, quase todos fictícios,...
Acabei de ser apresentado à distribuição Tweedie (veja isto ou isto ), mas estou tendo dificuldade em encontrar qual é a função de link para um modelo linear generalizado
Estou executando um glmm com uma variável de resposta binomial e um preditor categórico. O efeito aleatório é dado pelo design aninhado usado para a coleta de dados. Os dados são assim: m.gen1$treatment [1] sucrose control protein control no_injection ..... Levels: no_injection control sucrose...
Minha pergunta revela meu fraco entendimento da regressão de Poisson e GLMs em geral. Aqui estão alguns dados falsos para ilustrar minha pergunta: ### some fake data x=c(1:14) y=c(0, 1, 2, 3, 1, 4, 9, 18, 23, 31, 20, 25, 37, 45) Algumas funções personalizadas para retornar psuedo-R2: ###...
Eu tenho uma pergunta sobre se deve ou não usar um deslocamento. Suponha um modelo muito fácil, onde você deseja descrever o número (geral) de gols no hóquei. Então você tem gols, número de jogos jogados e uma variável dummy "atacante" que é igual a 1 se o jogador for atacante e 0 caso contrário....
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...
Eu tenho uma pergunta sobre modelos lineares generalizados (GLM). Minha variável dependente (DV) é contínua e não é normal. Então eu log transformá-lo (ainda não é normal, mas melhorou). Quero relacionar o DV com duas variáveis categóricas e uma covariável contínua. Para isso, quero conduzir um...
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...
No Rque pudermos "peso antes" uma glmregressão via pesos parâmetro. Por exemplo: glm.D93 <- glm(counts ~ outcome + treatment, family = poisson(), weights=w) Como isso pode ser feito em um JAGSou BUGSmodelo? Encontrei algum artigo discutindo isso, mas nenhum deles fornece um exemplo. Estou...
No código a seguir, realizo uma regressão logística em dados agrupados usando glm e "manualmente" usando mle2. Por que a função logLik em R me fornece uma probabilidade de log logLik (fit.glm) = - 2.336 que é diferente do logLik (fit.ml) = - 5.514 que recebo manualmente? library(bbmle) #successes...
Estou simplesmente tentando recalcular com dnorm () a probabilidade de log fornecida pela função logLik de um modelo lm (em R). Funciona (quase perfeitamente) para um grande número de dados (por exemplo, n = 1000): > n <- 1000 > x <- 1:n > set.seed(1) > y <- 10 + 2*x +...
Sou relativamente novo em R e estou tentando ajustar um modelo aos dados que consistem em uma coluna categórica e uma coluna numérica (número inteiro). A variável dependente é um número contínuo. Os dados têm o seguinte formato: predCateg, predIntNum, ResponseVar Os dados são mais ou menos...