Por exemplo, no manual BUGS ou no próximo livro de Lee e Wagenmakers ( pdf ) e em muitos outros lugares, é usado um tipo de notação que me parece muito flexível, pois pode ser usada para descrever sucintamente a maioria dos modelos estatísticos. Um exemplo dessa notação é o seguinte: yEu∼ Binomial...
No modelo misto, assumimos que os efeitos aleatórios (parâmetros) são variáveis aleatórias que seguem distribuições normais. Parece muito semelhante ao método bayesiano, no qual todos os parâmetros são assumidos como aleatórios. Então, o modelo de efeito aleatório é um caso especial do método...
O capítulo 2.4 da ESL parece classificar a regressão linear como "baseada em modelo", porque assume , enquanto nenhuma aproximação semelhante é declarada para os vizinhos k-mais próximos. Mas os dois métodos não fazem suposições sobre ?f( x ) ≈ x ⋅ βf(x)≈x⋅βf(x) \approx x\cdot\betaf( X...
Eu sou novo na análise de sobrevivência e aprendi recentemente que existem diferentes maneiras de fazer isso, com um determinado objetivo. Estou interessado na implementação real e na adequação desses métodos. Fui apresentado aos tradicionais Riscos Proporcionais de Cox , modelos de tempos de...
Estou analisando os resultados de um experimento de tempo de reação em R. Executei uma ANOVA de medidas repetidas (1 fator dentro do sujeito com 2 níveis e 1 fator entre sujeitos com 2 níveis). Eu executei um modelo misto linear semelhante e queria resumir os resultados anteriores na forma de...
Qual é a estrutura de variância-covariância padrão para efeitos aleatórios no glmerou lmerno lme4pacote? Como se especifica outra estrutura de variância-covariância para efeitos aleatórios no código? Não encontrei nenhuma informação sobre isso na
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...
Na regressão logística, um odds ratio de 2 significa que o evento é 2 vezes mais provável, dado um aumento de uma unidade no preditor. Na regressão de Cox, uma taxa de risco de 2 significa que o evento ocorrerá duas vezes mais a cada momento, dado um aumento de uma unidade no preditor. Não são...
Eu estou analisando alguns dados comportamentais usando lme4em R, principalmente após excelentes tutoriais do Bodo de Inverno , mas eu não entendo se eu estou lidando com interações corretamente. Pior, ninguém mais envolvido nesta pesquisa usa modelos mistos, então estou um pouco à deriva quando se...
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 tenho o seguinte modelo m_plotequipado com lme4::lmerefeitos aleatórios cruzados para participantes ( lfdn) e itens ( content): Random effects: Groups Name Variance Std.Dev. Corr lfdn (Intercept) 172.173 13.121 role1 62.351 7.896 0.03 inference1 24.640 4.964 0.08 -0.30 inference2...
Antes de fazer essa pergunta, pesquisei em nosso site e encontrei muitas perguntas semelhantes (como aqui , aqui e aqui ). Mas acho que essas perguntas relacionadas não foram bem respondidas ou discutidas, portanto, gostaria de levantar essa questão novamente. Eu acho que deve haver uma grande...
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...
Espero sinceramente que tenha formulado esta pergunta de forma que ela possa ser respondida definitivamente - caso contrário, informe-me e tentarei novamente! Eu também acho que devo usar R para essas análises. Eu tenho várias medidas plant performance (Ys)que eu suspeito que foram influenciadas...
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...
Eu tenho que criar gráficos (semelhantes aos gráficos de crescimento) para crianças de 5 a 15 anos (apenas 5,6,7 etc; não existem valores fracionários como 2,6 anos) para uma variável de saúde que não seja negativa, contínua e o intervalo de 50 a 150 (com apenas alguns valores fora desse...
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...
Estou interessado em modelar a captura total de peixes usando gam em mgcv para modelar efeitos aleatórios simples para navios individuais (que fazem viagens repetidas ao longo do tempo na pesca). Eu tenho 98 sujeitos, então pensei em usar gam em vez de gamm para modelar os efeitos aleatórios. Meu...
Dadas três variáveis, ye x, que são positivas contínuas e z, que são categóricas, tenho dois modelos candidatos dados por: fit.me <- lmer( y ~ 1 + x + ( 1 + x | factor(z) ) ) e fit.fe <- lm( y ~ 1 + x ) Espero comparar esses modelos para determinar qual modelo é mais apropriado. Parece-me...