Estou tentando implementar o algoritmo MCMC hit and run, mas estou com um pouco de dificuldade para entender como fazê-lo. A ideia geral é a seguinte: Para gerar um salto de proposta no MH, nós: Gere uma direção partir de uma distribuição na superfície da esfera unitáriadddOO\mathcal{O}...
Desenvolvi o pacote ez para R como um meio de ajudar as pessoas a fazer a transição de pacotes de estatísticas como SPSS para R. Isso é (espero) alcançado, simplificando a especificação de vários tipos de ANOVA e fornecendo resultados semelhantes ao SPSS (incluindo tamanhos e suposições de efeitos...
Como devo definir uma fórmula de modelo em R, quando uma (ou mais) restrições lineares exatas que vinculam os coeficientes está disponível. Como exemplo, digamos que você saiba que b1 = 2 * b0 em um modelo de regressão linear
Como a modelagem de regressão geralmente é mais "arte" do que ciência, muitas vezes me pego testando muitas iterações de uma estrutura de regressão. Quais são algumas maneiras eficientes de resumir as informações dessas várias execuções de modelo na tentativa de encontrar o "melhor" modelo? Uma...
Estou tentando implementar uma análise de "ponto de mudança" ou uma regressão multifásica usando nls()em R. Aqui estão alguns dados falsos que eu criei . A fórmula que eu quero usar para ajustar os dados é: y= β0 0+ β1x + β2max ( 0 , x - δ)y=β0 0+β1x+β2max(0 0,x-δ)y = \beta_0 + \beta_1x +...
Em geral, padronizo minhas variáveis independentes em regressões, a fim de comparar adequadamente os coeficientes (dessa forma, eles têm as mesmas unidades: desvios-padrão). No entanto, com dados longitudinais / em painel, não tenho certeza de como devo padronizar meus dados, especialmente se...
Estou explorando diferentes métodos de classificação para um projeto em que estou trabalhando e estou interessado em experimentar as Florestas Aleatórias. Estou tentando me educar ao longo do curso e gostaria de receber qualquer ajuda fornecida pela comunidade do CV. Dividi meus dados em conjuntos...
Fechadas. Esta questão está fora de tópico . No momento, não está aceitando respostas. Deseja melhorar esta pergunta? Atualize a pergunta para que ela esteja no tópico de Validação cruzada. Fechado há 4 anos . Eu tenho um quadro de dados que contém alguns...
Bloqueado . Esta pergunta e suas respostas estão bloqueadas porque a questão está fora do tópico, mas tem um significado histórico. No momento, não está aceitando novas respostas ou interações. Graças à pergunta do Tormod (postada aqui ), me deparei com o enredo Parallel...
Estou tentando gerar amostras aleatórias de um pdf personalizado usando R. Meu pdf é: fX( x ) = 32( 1 - x2) , 0 ≤ x ≤ 1fX(x)=32(1-x2),0 0≤x≤1f_{X}(x) = \frac{3}{2} (1-x^2), 0 \le x \le 1 Gerei amostras uniformes e tentei transformá-lo em minha distribuição personalizada. Eu fiz isso encontrando o...
Me pedem para escrever uma introdução à estatística e estou lutando para mostrar graficamente a maneira como o valor-p e o poder se relacionam. Eu vim com este gráfico: Minha pergunta: existe uma maneira melhor de exibir isso? Aqui está o meu código R x <- seq(-4, 4, length=1000) hx <-...
Eu tentei várias vezes "ir sozinho" - mas com sucesso limitado. Sou um usuário ocasional do SPSS e tenho alguma experiência com SAS. Gostaria de receber um ponteiro ou dois de alguém que tenha antecedentes semelhantes e agora use
Eu tenho um modelo linear clássico, com 5 possíveis regressores. Eles não estão correlacionados entre si e têm uma correlação bastante baixa com a resposta. Cheguei a um modelo em que três dos regressores têm coeficientes significativos para sua estatística t (p <0,05). A adição de uma ou das...
Olá gurus estatísticos e assistentes de programação R, Estou interessado em modelar capturas de animais em função das condições ambientais e do dia do ano. Como parte de outro estudo, tenho contagens de capturas em ~ 160 dias em três anos. Em cada um desses dias, tenho temperatura, precipitação,...
Fui encarregado de mudar um dos nossos grandes modelos estocásticos atuais do SAS para um novo idioma. Pessoalmente, prefiro uma linguagem compilada tradicional, mas o PI quer que eu verifique o R, que nunca usei. Nossa motivação para tirar o modelo do SAS é (1) muitas pessoas não têm acesso a ele...
Aqui está o meu experimento: Estou usando a findPeaksfunção no pacote quantmod : Quero detectar picos "locais" dentro de uma tolerância 5, ou seja, os primeiros locais após a série cronológica caírem dos picos locais em 5: aa=100:1 bb=sin(aa/3) cc=aa*bb plot(cc, type="l") p=findPeaks(cc,...
Estou usando plm()para estimar modelos de efeitos fixos do formulário y ~ x + time + time:fixed_trait onde fixed_traité uma variável que varia entre indivíduos, mas é constante dentro dos indivíduos. O ponto de interagir timecom fixed_traité permitir o efeito defixed_trait variar ao longo do...
Preciso prever as seguintes 4 variáveis para a 29ª unidade de tempo. Eu tenho aproximadamente 2 anos de dados históricos, onde 1, 14 e 27 são todos do mesmo período (ou época do ano). No final, estou fazendo uma decomposição no estilo Oaxaca-Blinder em WWW , , e .w c pw dWdwdw cWcwcppp time W wd...
Suponha que w1,w2,…,wnw1,w2,…,wnw_1,w_2,\ldots,w_n e x1,x2,...,xnx1,x2,...,xnx_1,x_2,...,x_n são cada iid extraído de algumas distribuições, com wiwiw_i independente de xixix_i . O wiwiw_i é estritamente positivo. Você observa todo o wiwiw_i , mas não o xixix_i ; em vez disso você observa...
Costumo lidar com dados confusos da pesquisa, que exigem muita limpeza antes que qualquer estatística possa ser feita. Eu costumava fazer isso "manualmente" no Excel, às vezes usando fórmulas do Excel e, às vezes, verificando as entradas uma a uma. Comecei a executar cada vez mais essas tarefas...