Estou tentando resolver um problema que lida com a imputação de dados ausentes de um estudo de dados em painel (não tenho certeza se estou usando o 'estudo de dados em painel' corretamente - como o aprendi hoje.) Tenho dados de contagem total de mortos nos anos 2003 até 2009, todos os meses,...
Estou lendo um texto, "Estatística Matemática e Análise de Dados", de John Rice. Estamos preocupados com a aproximação do valor esperado e variância da variável aleatória . Podemos calcular o valor e a variação esperados da variável aleatória e conhecemos a relação . Portanto, é possível aproximar...
Eu tenho usado rpart.controlpara minsplit=2, e obteve os seguintes resultados de rpart()função. Para evitar o ajuste excessivo dos dados, preciso usar as divisões 3 ou 7? Não devo usar divisões 7? Por favor deixe-me saber. Variáveis realmente usadas na construção de árvores: [1] ct_a ct_b...
Uso R há vários anos, mas sempre usei uma GUI como Tinn-R, JGR e, mais recentemente, R-Studio. Embora eu goste da interface do R-Studio, sinto que, para programas mais longos com poucos / nenhum gráfico, eu poderia codificar com mais eficiência, aproveitando os atalhos oferecidos por um editor de...
Preciso desenhar 20 distribuições em um único gráfico em R, e isso não me parece bom (confuso) com boxplot regular (20 caixas), mesmo com boxwex = 0,3. Você poderia me sugerir como plotar um tipo de boxplot em R para as 20 distribuições, com pontos para mediana e apenas uma linha em vez de caixa,...
Eu tenho as seguintes séries temporais obtido usando os dados postados abaixo. Para um tamanho de janela deslizante de 10, estou tentando calcular a divergência KL entre o PMF de valores na janela deslizante atual e o PMF da história com o objetivo final de traçar o valor da divergência KL ao...
Eu preciso fazer inferência sobre um parâmetro positivo . Para acomodar a positividade, reparametrizei . Usando a rotina do MLE, calculei a estimativa de pontos e se para . A propriedade invariância do MLE me fornece diretamente uma estimativa de pontos para , mas não sei como calcular se para ....
Fiz uma ANOVA descobrindo, por exemplo, uma interação entre gênero e série do que gostaria de saber em que séries meninos e meninas diferem, mas em muitos casos encontro valores-p (ajustados) de 0 e 1. Como / por que isso é possível? Não parece certo ... as.factor(gender) 1 16 16.2 2.6377 0.104396...
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. As funções que eu estou familiarizado incluem a escala da base R, a nova...
Entendo que não devemos permitir que o mesmo conjunto de dados que estamos analisando direcione / defina como são as distribuições anteriores em uma análise bayesiana. Especificamente, não é apropriado definir distribuições anteriores para uma análise bayesiana com base em estatísticas resumidas do...
Enquanto verificava se duas variáveis estavam correlacionadas, observei que a aplicação da correlação de Pearson produziu números tão baixos quanto 0,1, indicando que não há correlação. Existe algo que eu possa fazer para fortalecer essa afirmação? O conjunto de dados (subconjunto por causa dos...
Estou revisando o livro 'Introductory Time Series with R' de Cowpertwait e Metcalfe. Na página 36, diz que as linhas estão em: . Eu li aqui no fórum R que as linhas estão em . ±1,96/ √−1/n±2/n−−√−1/n±2/n-1/n \pm 2/\sqrt{n}±1.96/n−−√±1.96/n\pm 1.96/\sqrt{n} Eu executei o seguinte código: b =...
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 vetor x <- c(1,2,3,4,5,5,5,6,6,6,6, 7,7,7,7,7,7,7,7,7,7,7,7,7,7,7,7,7,7,7, 7,7,7,7,7,7,7,7,8,8,8,8,9,9,9,10) (meu vetor real tem um comprimento> 10.000) e gostaria de encontrar os intervalos em que estão os 90% da densidade. É quantile(x, probs=c(0.05,0.95), type=5)o mais...
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,...
Existe alguma diferença entre os seguintes termos ou são iguais? Viés Viés sistemático Erros sistemáticos Se houver algumas diferenças, por favor, explique-as. Esses erros podem ser reduzidos quando se aumenta o tamanho da amostra? UPDATE: Meu campo de interesse é inferência estatística....
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 visualizar alguns dados do consumidor, que tem quatro categorias. Os usuários podem alternar entre diferentes categorias. Gostaria de visualizar os últimos três ou quatro comutadores para cada indivíduo. Então começamos com um gráfico com uma coluna com 4 proporções empilhadas....
Tenho gráficos de caixa de 13 grupos que mostro em um gráfico. Os grupos têm populações desequilibradas e não são normalmente distribuídas. Quero mostrar quais pares são estatisticamente semelhantes (ou seja, têm valor de p kruskal.test <0,05) colocando a, b, c, etc. em cima das caixas...
Se eu tenho uma variável com 4 níveis, em teoria eu preciso usar 3 variáveis fictícias. Na prática, como isso é realmente realizado? Uso 0-3, uso 1-3 e deixo o 4 em branco? Alguma sugestão? NOTA: Eu vou trabalhar na R. ATUALIZAÇÃO: O que aconteceria se eu apenas usasse uma coluna que use 1-4...