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...
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...
Eu tenho dois modelos de regressão logística em R feitos com glm(). Ambos usam as mesmas variáveis, mas foram criados usando diferentes subconjuntos de uma matriz. Existe uma maneira fácil de obter um modelo médio que forneça os meios dos coeficientes e depois usá-lo com a função predict...
Estou investigando a interação entre duas variáveis ( e ). Existe uma grande correlação linear entre essas variáveis com . Pela natureza do problema, não posso dizer nada sobre a causa (se causa ou o contrário). Eu gostaria de estudar os desvios da linha de regressão, a fim de detectar valores...
Se houver 0 na tabela de contingência e estivermos ajustando modelos poisson / loglinear aninhados (usando a glmfunção de R ) para um teste de razão de verossimilhança, precisamos ajustar os dados antes de ajustar os modelos glm (por exemplo, adicione 1/2 a todos as contagens)? Obviamente, alguns...
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
Foi-me pedido para analisar alguns dados de um ensaio clínico, procurando dois métodos para medir a pressão arterial. Eu tenho dados de 50 indivíduos, cada um com entre 2 e 57 medidas usando cada método. Eu estou querendo saber a melhor forma de proceder. Obviamente, preciso de uma solução que...
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),...
Dada a seguinte configuração experimental: Várias amostras são coletadas de um sujeito e cada amostra é tratada de várias maneiras (incluindo um tratamento de controle). O que é principalmente interessante é a diferença entre o controle e cada tratamento. Eu posso pensar em dois modelos simples...
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} ~...
A Wikipedia sugere que uma maneira de analisar a confiabilidade entre avaliadores é usar um modelo de efeitos aleatórios para calcular a correlação intraclasse . O exemplo da correlação intraclasse fala sobre olhar
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 analisando se a abundância está relacionada ao tamanho. O tamanho é (é claro) contínuo, no entanto, a abundância é registrada em uma escala tal que A = 0-10 B = 11-25 C = 26-50 D = 51-100 E = 101-250 F = 251-500 G = 501-1000 H = 1001-2500 I = 2501-5000 J = 5001-10,000 etc... A a Q ... 17...
Digamos que eu gere a probabilidade de um resultado com base em um determinado fator e trace a curva desse resultado. Existe alguma maneira de extrair a equação para essa curva de R? > mod = glm(winner~our_bid, data=mydat, family=binomial(link="logit")) > summary(mod) Call: glm(formula =...
Qual é a diferença entre um modelo GLM (regressão logística) com uma variável de resposta binária que inclui sujeito e tempo como covariáveis e o modelo GEE análogo que leva em consideração a correlação entre medições em vários momentos? Meu GLM se parece com: Y(binary) ~ A + B1X1(subject id)...
Criei minha própria versão ligeiramente aprimorada do termplot que utilizo neste exemplo, você pode encontrá-la aqui . Eu publiquei anteriormente no SO, mas quanto mais penso nisso, acredito que isso provavelmente esteja mais relacionado à interpretação do modelo de riscos proporcionais de Cox do...
Como eu usaria a saída de uma análise de componentes principais (PCA) em um modelo linear generalizado (GLM), assumindo que o PCA é usado para seleção de variáveis para o GLM? Esclarecimento: Eu quero usar o PCA para evitar o uso de variáveis correlacionadas no GLM. No entanto, o PCA me...
Existem métodos para corrigir o viés no modelo de risco proporcional de Cox causado por amostra selecionada não aleatoriamente (algo como a correção de Heckman)? Antecedentes : Digamos que a situação seja a seguinte: - Nos primeiros dois anos, todos os clientes são aceitos. - Após esses dois...