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...
Estou confuso com declarações em uma página da UCLA sobre regressão logística de efeitos mistos. Eles mostram uma tabela de coeficientes de efeitos fixos ao se ajustar a esse modelo e o primeiro parágrafo abaixo parece interpretar os coeficientes exatamente como uma regressão logística normal. Mas...
Eu tenho uma amostra de 1.449 pontos de dados que não estão correlacionados (r ao quadrado 0,006). Ao analisar os dados, descobri que, dividindo os valores das variáveis independentes em grupos positivo e negativo, parece haver uma diferença significativa na média da variável dependente para...
Eu vi essa notação para mínimos quadrados comuns aqui . minw∥Xw−y∥22minw‖Xw−y‖22 \min_w \left\| Xw - y \right\|^2_2 Eu nunca vi as barras duplas e as 2 na parte inferior. O que significam esses símbolos? Eles têm terminologia específica para
Tanto bayesglm()(no pacote arm R) quanto várias funções no pacote MCMCpack têm como objetivo fazer uma estimativa bayesiana de modelos lineares generalizados, mas não tenho certeza de que eles estejam computando a mesma coisa. As funções do MCMCpack usam a cadeia de Markov Monte Carlo para obter...
Vamos supor que eu tenha um vetor - de variáveis dependentes e um vetor - de variável independente. Quando é plotado contra , vejo que há uma relação linear (tendência ascendente) entre os dois. Agora, isso também significa que há uma tendência de queda linear entre e...
Quais são as diferenças no uso de Modelos Lineares Generalizados, como Determinação Automática de Relevância (ARD) e regressão de Ridge, versus modelos de séries temporais como Box-Jenkins (ARIMA) ou suavização exponencial para previsão? Existem regras práticas sobre quando usar o GLM e quando usar...
Estou executando glms em R (modelos lineares generalizados). Eu pensei que conhecia pvalues - até que vi que chamar um resumo para um glm não lhe desse um pvalor predominante representativo do modelo como um todo - pelo menos não no lugar onde os modelos lineares o fazem. Gostaria de saber se...
Existem poucas explicações que descrevem como interpretar os coeficientes de regressão linear após diferenciar uma série temporal (para eliminar uma raiz unitária). É tão simples que não há necessidade de declarar formalmente? (Estou ciente dessa pergunta , mas não tinha certeza de quão geral foi...
Tenho uma pergunta sobre a validade do uso do RMSE (Root Mean Squared Error) para comparar diferentes modelos logísticos. A resposta é ou 0ou 1e as previsões são probabilidades entre 0- 1? A maneira aplicada abaixo também é válida com as respostas binárias? # Using
Em R, por que as configurações padrão qqplot(linear model)usam os resíduos padronizados no eixo y? Por que R não usa os resíduos
Usando R, ajustei um modelo linear para uma variável de resposta única a partir de uma mistura de preditores contínuos e discretos. Isso é super básico, mas estou tendo problemas para entender como funciona um coeficiente para um fator discreto. Conceito: Obviamente, o coeficiente da variável...
Sei que essa é provavelmente uma pergunta básica ... Mas não acho a resposta. Estou ajustando um GLM com uma família Poisson e tentei dar uma olhada nas previsões, no entanto, o deslocamento parece ser levado em
Estou tentando ajudar um aluno de um colega. O aluno observou e contou o comportamento das aves (número de chamadas) em uma configuração experimental. O número de chamadas atribuíveis a uma ave observada específica durante cada experimento não pôde ser determinado, mas foi possível contar o número...
Estou simulando ensaios de Bernoulli com um aleatório entre grupos e, em seguida, encaixo o modelo correspondente em o pacote:logitθ∼N(logitθ0,12)logitθ∼N(logitθ0,12)\text{logit}\, \theta \sim {\cal N}(\text{logit}\, \theta_0, 1^2)lme4 library(lme4) library(data.table) I <- 30 # number of...
O seguinte modelo logístico multinível com uma variável explicativa no nível 1 (nível individual) e uma variável explicativa no nível 2 (nível do grupo): π 0 j = γ 00 + γ 01 z j + u 0 j … ( 2 ) π 1 j = γ 10 + γ 11 z j + u 1 j … ( 3
Questões: Modelos lineares impróprios são utilizados na prática ou algum tipo de curiosidade é descrita de tempos em tempos em revistas científicas? Em caso afirmativo, em que áreas eles são usados? Existem outros exemplos de tais modelos? Finalmente, os erros padrão, valores- , etc. retirados do...
No meu campo, a maneira usual de plotar dados emparelhados é como uma série de segmentos finos de linhas inclinadas, sobrepondo-os à mediana e ao IC da mediana para os dois grupos: No entanto, esse tipo de gráfico se torna muito mais difícil de ler à medida que o número de pontos de dados se...