Examinei o Google e este site e ainda estou confuso sobre a função lmer na biblioteca lme4. Tenho alguns dados coletados de diferentes enfermarias psiquiátricas, que possuem uma estrutura multinível. Para simplificar, selecionarei duas variáveis de nível 2 e duas de nível 1, embora eu realmente...
Uma ocorrência não tão incomum ao lidar com modelos mistos máximos complexos (estimar todos os efeitos aleatórios possíveis para dados e modelo) é perfeita (+1 ou -1) ou correlação quase perfeita entre alguns efeitos aleatórios. Para o propósito da discussão, vamos observar o seguinte modelo e...
Considere um experimento com vários participantes humanos, cada um medido várias vezes em duas condições. Um modelo de efeitos mistos pode ser formulado (usando a sintaxe lme4 ) como: fit = lmer( formula = measure ~ (1|participant) + condition ) Agora, digamos que eu queira gerar intervalos de...
Suponha que, em um estudo de 15 sujeitos, a variável resposta (res) seja modelada com duas variáveis explicativas, uma (nível) é categórica com 5 níveis e a outra (tempo de resposta: TR) é contínua. Com o lmer no pacote lme4 do R, tenho: fm1 <- lmer(res ~ level * RT + (level-1 | subject),...
Essas perguntas podem parecer estúpidas, mas ... é correto que efeitos aleatórios possam ser aplicados apenas a variáveis categóricas (como identificação individual, identificação de população, ...), por exemplo, digamos que é uma variável categórica:xixix_i ~ β x iyiyiy_iβxiβxi\beta_{x_i} ~...
Eu tenho o seguinte tipo de dados. Avaliei 10 indivíduos, cada um repetido 10 vezes. Tenho matriz de relação 10x10 (relação entre todas as combinações dos indivíduos). set.seed(1234) mydata <- data.frame (gen = factor(rep(1:10, each = 10)), repl = factor(rep(1:10, 10)), yld = rnorm(10, 5,...
Meu colega e eu estamos ajustando uma variedade de modelos de efeitos mistos lineares e não-lineares em R. Somos solicitados a realizar validação cruzada nos modelos ajustados, para que se possa verificar se os efeitos observados são relativamente generalizáveis. Normalmente, essa é uma tarefa...
Em R, como especifico o modelo lmer sem efeito fixo global? Por exemplo, se eu disser algo como lmer(y ~ (1 | group) + (0 + x | group), data = my_df) o modelo ajustado será yij=a+αi+βixijyij=a+αi+βixijy_{ij} = a + \alpha_{i} + \beta_{i} x_{ij} Como encaixo o
Estou tentando lidar com uma análise de tempo para evento usando resultados binários repetidos. Suponha que o tempo do evento seja medido em dias, mas, no momento, discretizamos o tempo em semanas. Quero aproximar um estimador de Kaplan-Meier (mas permitir covariáveis) usando resultados binários...
Os seguintes enxertos são retirados deste artigo . Eu sou novato no bootstrap e estou tentando implementar o bootstrap paramétrico, semiparamétrico e não paramétrico para o modelo misto linear com o R bootpacote. Código R Aqui está o meu
Estou tentando analisar o efeito de Year na variável logInd para um grupo específico de indivíduos (eu tenho 3 grupos). O modelo mais simples: > fix1 = lm(logInd ~ 0 + Group + Year:Group, data = mydata) > summary(fix1) Call: lm(formula = logInd ~ 0 + Group + Year:Group, data =...
Estou tentando produzir um modelo misto linear, o código R é o seguinte. lme (Average.payoff ~ Jogo + Tipo + Outros.Tipo + Jogo: Tipo + Jogo: Outros.Tipo + Tipo: Outros.Tipo aleatório = ~ 1 | Assuntos, método = "REML", dados = Assuntosm1) -> lme1 O termo de resposta Average.payoff é...
Em um estudo longitudinal, os resultados de unidades são medidos repetidamente nos pontos de tempo com um total de ocasiões de medição fixas (fixo = medições nas unidades são feitas ao mesmo tempo). i t mYeu tYitY_{it}Euiitttmmm As unidades são atribuídas aleatoriamente a um tratamento, , ou a um...
Eu tenho uma pergunta básica sobre abordagens para modelar a média usando critérios de TI para ponderar modelos em um conjunto de candidatos. A maioria das fontes que li sobre a média do modelo defende a média das estimativas do coeficiente de parâmetros com base nos pesos do modelo (usando um...
Estou ajustando um glmermodelo no lme4pacote R. Estou procurando uma tabela anova com o valor p mostrado nela, mas não consigo encontrar nenhum pacote que se encaixe nela. É possível fazer isso em R? O modelo que estou ajustando é da forma: model1<-glmer(dmn~period*teethTreated+(1|fullName),...
Devo incluir efeitos aleatórios em um modelo, mesmo que não sejam estatisticamente significativos? Eu tenho um desenho experimental de medidas repetidas, no qual cada indivíduo experimenta três tratamentos diferentes em ordem aleatória. Eu gostaria de controlar quaisquer efeitos individuais e de...
Exemplos: Eu tenho uma frase na descrição do trabalho: "Java senior engineer in UK". Eu quero usar um modelo de aprendizado profundo para prever em duas categorias: English e IT jobs. Se eu usar o modelo de classificação tradicional, ele poderá prever apenas 1 rótulo com softmaxfunção na última...
Estou com problemas para obter resultados equivalentes entre um aovmodelo de medidas repetidas entre dentro e um lmermodelo misto. Meus dados e script têm a seguinte aparência data=read.csv("https://www.dropbox.com/s/zgle45tpyv5t781/fitness.csv?dl=1") data$id=factor(data$id) data id FITNESS TEST...
A distinção entre variáveis estritamente aleatórias (que devem ser modeladas como tais) e variáveis não aleatórias que alguns argumentam que poderiam ser modeladas como aleatórias se for um modelo hierárquico / multinível, é pouco clara para mim. Bates e Bolker exemplificam efeitos aleatórios...
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...