Estou lidando com um modelo linear hierárquico bayesiano , aqui a rede que o descreve. YYY representa as vendas diárias de um produto em um supermercado (observado). XXX é uma matriz conhecida de regressores, incluindo preços, promoções, dia da semana, clima e feriados. SSS é o nível de...
Eu tenho alguns dados em [0,1] que gostaria de analisar com uma regressão beta. Obviamente, algo precisa ser feito para acomodar os valores 0,1. Não gosto de modificar dados para ajustar-se a um modelo. Também não acredito que a inflação zero e 1 seja uma boa ideia, porque acredito que, neste caso,...
Eu estive analisando a modelagem de efeitos mistos usando o pacote lme4 em R. Estou usando principalmente o lmercomando, por isso vou fazer minha pergunta por meio de código que usa essa sintaxe. Suponho que possa ser uma pergunta fácil geral, é correto comparar dois modelos construídos no lmeruso...
Se f1,…,fkf1,…,fkf_1,\ldots,f_k são densidades conhecidas das quais posso simular, ou seja, para as quais um algoritmo está disponível. e se o produto ∏i=1kfi(x)αiα1,…,αk>0∏i=1kfi(x)αiα1,…,αk>0\prod_{i=1}^k f_i(x)^{\alpha_i}\qquad \alpha_1,\ldots,\alpha_k>0 é integrável, existe uma abordagem...
Sabemos que um teste t emparelhado é apenas um caso especial de ANOVA unidirecional de medidas repetidas (ou dentro do sujeito), bem como modelo linear de efeitos mistos, que pode ser demonstrado com a função lme (), o pacote nlme em R como mostrado abaixo. #response data from 10 subjects under...
Digamos que eu tenha um conjunto de dados de "cateter renal". Estou tentando modelar uma curva de sobrevivência usando um modelo de Cox. Se eu considerar um modelo de Cox: preciso da estimativa do risco da linha de base. Usando a função interna do pacote R , posso fazer isso da seguinte maneira:h (...
Eu tenho uma tabela de dados T1, que contém quase mil variáveis (V1) e cerca de 200 milhões de pontos de dados. Os dados são escassos e a maioria das entradas é NA. Cada ponto de dados possui um par único de identificação e data para distinguir um do outro. Eu tenho outra tabela T2, que contém...
Isso parece tão elementar, mas eu sempre fico preso neste momento ... A maioria dos dados com os quais eu lido não é normal e a maioria das análises é baseada em uma estrutura GLM. Para minha análise atual, tenho uma variável de resposta que é "velocidade de caminhada" (metros / minuto). É fácil...
Estou tentando usar lmedo nlmepacote para replicar resultados de aovANOVAs de medidas repetidas. Eu fiz isso para um experimento de medidas repetidas de fator único e para um experimento de dois fatores com um fator entre sujeitos e um fator dentro dos sujeitos, mas estou tendo problemas para fazer...
Considere os seguintes dados de uma maneira bidirecional no design de assuntos: df <- "http://personality-project.org/r/datasets/R.appendix4.data" df <- read.table(df,header=T) head(df) Observation Subject Task Valence Recall 1 1 Jim Free Neg 8 2 2 Jim Free Neu 9 3 3 Jim Free Pos 5 4 4 Jim...
Tenho uma planilha com algumas pontuações da equipe. Primeira equipe com 10 pontos ganha. Existem 2 jogadores em cada equipe. Os jogadores jogam com diferentes companheiros de equipe o tempo todo, embora não sejam escolhidos perfeitamente aleatoriamente. Nenhuma pontuação individual é...
Se eu fixar os valores dos nós observados de um MRF, ele se tornará um
Eu sou interessante em aprender sobre como classificar e classificar indivíduos em um grupo que apenas interagem / competem de maneira pareada (ou seja, sistemas como o sistema de classificação ELO para xadrez). Existem métodos de acesso ou métodos mais precisos e avançados por aí? Existem...
Digamos que estamos trabalhando em um modelo de efeitos aleatórios de alguns dados de contagem ao longo do tempo e queremos controlar algumas tendências. Normalmente, você faria algo como: lmer(counts ~ dependent_variable + (1+t+I(t^2)|ID), family="poisson") para incluir uma forma quadrática para...
Basicamente, o que eu quero saber é como as diferentes estruturas de covariância são aplicadas e como os valores dentro dessas matrizes são calculados. Funções como lme () permitem escolher qual estrutura gostaríamos, mas eu adoraria saber como elas são estimadas. Considere o modelo linear de...
Em resposta a essa pergunta, sobre se meu projeto, onde eu aleatoriamente apresentei aos participantes fotos de diferentes categorias, era um exemplo em que eu deveria usar uma ANOVA de medidas repetidas, obtive a resposta de que deveria usar um modelo misto, com uma das opções razões são que eu...
Por favor, considere estes dados: dt.m <- structure(list(id = c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12), occasion = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L), .Label = c("g1", "g2"), class =...
Digamos que precisamos dos GLMMs mod1 <- glmer(y ~ x + A + (1|g), data = dat) mod2 <- glmer(y ~ x + B + (1|g), data = dat) Esses modelos não são aninhados no sentido usual de: a <- glmer(y ~ x + A + (1|g), data = dat) b <- glmer(y ~ x + A + B + (1|g), data = dat) então não podemos...
Eu tenho um gráfico de valores residuais de um modelo linear em função dos valores ajustados, onde a heterocedasticidade é muito clara. No entanto, não tenho certeza de como devo proceder agora, porque, tanto quanto eu entendo essa heterocedasticidade, torna meu modelo linear inválido. (Isso está...
Preciso de alguns conselhos sobre dois dilemas principais da minha pesquisa, que é um estudo de caso de três grandes produtos farmacêuticos e inovação. O número de patentes por ano é a variável dependente. Minhas perguntas são Quais são os critérios mais importantes para um bom modelo? O que é...