Perguntas com a marcação «model»

Nota: esta pergunta é um repost, pois minha pergunta anterior teve que ser excluída por razões legais. Ao comparar o PROC MIXED do SAS com a função lmedo nlmepacote no R, deparei-me com algumas diferenças bastante confusas. Mais especificamente, os graus de liberdade nos diferentes testes...

Recentemente, li sobre Bayes empírico (Casella, 1985, Uma introdução à análise de dados empíricos de Bayes) e parecia muito com o modelo de efeitos aleatórios; em que ambos têm estimativas reduzidas à média global. Mas eu não li completamente ... Alguém tem alguma idéia sobre a semelhança e as...

Estou tendo problemas de desempenho usando o MCMCglmmpacote no R para executar um modelo de efeitos mistos. O código fica assim: MC1<-MCMCglmm(bull~1,random=~school,data=dt,family="categorical" , prior=list(R=list(V=1,fix=1), G=list(G1=list(V=1, nu=0))) , slice=T, nitt=iter, ,burnin=burn,...

Não consegui descobrir como executar a regressão linear em R in para um design de medida repetida. Em uma pergunta anterior (ainda sem resposta), me foi sugerido não usar, lmmas usar modelos mistos. Eu usei lmda seguinte maneira: lm.velocity_vs_Velocity_response <-

Eu quero calcular a função de sobrevivência prevista para um modelo de riscos proporcionais de Cox com termos de fragilidade [usando o pacote de sobrevivência]. Parece que quando termos de fragilidade estão no modelo, a função de sobrevivência prevista não pode ser calculada. ## Example...

Dado o mesmo conjunto de covariáveis ​​e família de distribuição, como posso comparar modelos com diferentes funções de link? Acho que a resposta correta aqui é "AIC / BIC", mas não tenho 100% de certeza. É possível ter modelos aninhados se eles tiverem um link

Eu tenho um conjunto de dados que representa decaimento exponencial. Eu gostaria de ajustar uma função exponencial a esses dados. Eu tentei log transformando a variável de resposta e, em seguida, usando menos quadrados para ajustar uma linha; usando um modelo linear generalizado com uma função de...

Recebi um conjunto de dados que contém o número de prêmios ganhos por alunos de uma escola secundária, onde os preditores do número de prêmios ganhos incluem o tipo de programa em que o aluno estava matriculado e a pontuação no exame final em matemática. Fiquei me perguntando se alguém poderia me...

Estou familiarizado com notações como: ondeβ0j=β0+uj, eyij=β0+βixij+uj+eij=β0j+βixij+eijyij=β0+βixij+uj+eij=β0j+βixij+eij\begin{align} y_{ij} &= \beta_0 + \beta_i x_{ij} + u_j + e_{ij}\\ &= \beta_{0j} + \beta_i x_{ij} + e_{ij} \end{align}β0j=β0+ujβ0j=β0+uj\beta_{0j}=\beta_{0}+u_j...

Eu tenho um conjunto de dados com 8000 clusters e 4 milhões de observações. Infelizmente, meu software estatístico, Stata, roda lentamente ao usar sua função de dados em painel para regressão logística: xtlogitmesmo com uma subamostra de 10%. No entanto, ao usar a logitfunção não painel , os...

O que explica as diferenças nos valores de p no seguinte aove nas lmchamadas? A diferença é apenas devido a diferentes tipos de cálculos de somas de

Estou executando um modelo linear generalizado, no qual tenho que especificar uma família diferente da normal. Qual é a distribuição esperada de resíduos? Por exemplo, os resíduos devem ser distribuídos

Gostaria de entender o que o código a seguir está fazendo. A pessoa que escreveu o código não trabalha mais aqui e é quase completamente indocumentada. Me pediram para investigá-lo por alguém que pensasse " é um modelo de regressão logística bayesiana " bglm <- function(Y,X) { # Y is a vector...

Quais são os parâmetros start, etastart, mustartna função GLM () ? Eu estive pesquisando os documentos e a Internet, mas não encontrei uma explicação clara do que isso significa. Assemelha-se a "valores iniciais" bayesianos para as cadeias, mas duvido que isso esteja relacionado, pois a função glm...

Considere um modelo de efeitos não observados linear do tipo: onde é uma característica não observada, mas invariante no tempo e é um erro, e indexe observações individuais e tempo, respectivamente. A abordagem típica em uma regressão de efeitos fixos (EF) seria remover por meio de manequins...

Eu tenho um conjunto de dados de séries temporais no qual estou tentando ajustar um Modelo de Markov oculto (HMM) para estimar o número de estados latentes nos dados. Meu pseudo-código para fazer isso é o seguinte: for( i in 2 : max_number_of_states ){ ... calculate HMM with i states ......

Estou ajustando um modelo de efeitos aleatórios glmera alguns dados comerciais. O objetivo é analisar o desempenho de vendas por distribuidor, levando em consideração a variação regional. Eu tenho as seguintes variáveis: distcode: ID do distribuidor, com cerca de 800 níveis region: ID geográfico...

Fiquei com a impressão de que a função lmer()no lme4pacote não produzia valores de p (veja lmervalores de p e tudo isso ). Eu tenho usado valores MCMC gerados p vez como por esta pergunta: efeito significativo no lme4modelo misto e esta pergunta: Não é possível encontrar p-valores na saída do...

Correlacionei dados e estou usando um modelo de efeitos mistos de regressão logística para estimar o efeito no nível individual (condicional) de um preditor de interesse. Eu sei que, para modelos marginais padrão, a inferência nos parâmetros do modelo usando o teste de Wald é consistente para os...

Ao usar splines cúbicos naturais (isto é, restritos), as funções básicas criadas são altamente colineares e, quando usadas em uma regressão, parecem produzir estatísticas VIF (fator de inflação de variação) muito altas, sinalizando multicolinearidade. Quando se está considerando o caso de um modelo...