Esta pergunta descreve a diferença básica entre um histograma uniforme e um não uniforme. E essa pergunta discute a regra básica para escolher o número de posições de um histograma uniforme que otimiza (em certo sentido) o grau em que o histograma representa a distribuição a partir da qual as...
Estou tentando ajustar um modelo de tempo discreto no R, mas não sei como fazê-lo. Eu li que você pode organizar a variável dependente em linhas diferentes, uma para cada observação no tempo, e usar a glmfunção com um link logit ou cloglog. Neste sentido, tem três colunas: ID, Event(1 ou 0, em...
Estou trabalhando em um conjunto de dados. Depois de usar algumas técnicas de identificação de modelos, criei um modelo ARIMA (0,2,1). Usei a detectIOfunção no pacote TSAem R para detectar um outlier inovador (IO) na 48ª observação do meu conjunto de dados original. Como faço para incorporar esse...
Aqui está um exemplo de caso: Eu tenho uma população de 10.000 itens. Cada item tem um ID exclusivo. Escolho aleatoriamente 100 itens e registro os IDs Coloquei os 100 itens de volta na população Escolho aleatoriamente 100 itens novamente, registro os IDs e os substituo. No total, repito essa...
Eu tenho alguns dados e estava tentando ajustar uma curva suave para ele. No entanto, não quero impor muitas crenças anteriores ou pré-concepções muito fortes (exceto as implícitas pelo restante da minha pergunta) sobre ela, ou quaisquer distribuições específicas. Eu só queria ajustá-lo com alguma...
Apesar de várias tentativas de ler sobre o bootstrap, pareço sempre bater em uma parede de tijolos. Gostaria de saber se alguém pode dar uma definição razoavelmente não técnica de bootstrapping? Sei que não é possível neste fórum fornecer detalhes suficientes para que eu possa entendê-lo...
No meu conjunto de dados, tenho cinco grupos (ordinais) com uma quantidade x de medida. Como a homoscedasticidade é violada, realizei o teste do qui-quadrado de Friedman para verificar se existem diferenças estatísticas entre os grupos: fried = stats.friedmanchisquare(*[grp for idx, grp in...
Eu tenho que analisar um planejamento fatorial com cinco fatores (um deles aninhado em outro) e respostas numéricas. Eu gostaria de realizar uma ANOVA não paramétrica, mas é claro que não posso usar nem o teste de Kruskall Wallis nem Friedman (eu repliquei medidas). Existe um comando ou um código...
Eu escrevi um script que testa os dados usando o wilcox.test, mas quando obtive os resultados, todos os valores de p eram iguais a 1. Li em alguns sites que você poderia usar o jitter antes de testar os dados (para evitar empates, como eles disseram), Eu fiz isso e agora tenho um resultado...
Estou analisando dados extremamente não lineares para os quais os modelos ARMA / ARIMA não funcionam bem. No entanto, vejo alguma autocorrelação e suspeito ter melhores resultados para a autocorrelação não linear. 1 / existe um equivalente do PACF para correlação de classificação? (em R?) 2 /...
Digamos que definimos uma distância, que não é uma métrica , entre N itens. Com base nessa distância, usamos um cluster hierárquico aglomerativo . Podemos usar cada um dos algoritmos conhecidos (ligação única / máxima / média etc.) para obter resultados significativos? Em outras palavras, qual é...
Estou usando o PROC GLM no SAS para ajustar uma equação de regressão da seguinte forma Y= b0 0+ b1 1X1 1+ b2X2+ b3X3+ b4tY=b0 0+b1 1X1 1+b2X2+b3X3+b4t Y = b_0 + b_1X_1 + b_2X_2 + b_3X_3 + b_4t O gráfico QQ dos vermelhos resultantes indica desvio da normalidade. Qualquer transformação de não é...
O que seria um bom livro para estatísticas não paramétricas. Não apenas a introdução, mas o nível avançado. Também estou procurando algo que eu possa usar para aprender e não para referência. Em particular, estou procurando um livro que possa conter conceitos básicos por trás de métodos não...
Para uma tabela 2 por 2 contingência, alguns disseram que o teste exato de Fisher usa a contagem no (1,1) célula na tabela como a estatística de teste, e sob hipótese nula, vontade tem uma distribuição hipergeométrica. X 1 , 1X1 , 1X1,1X_{1,1}X1 , 1X1,1X_{1,1} Alguns disseram que sua estatística...
Gostaria de entender esse tópico, mas aprender com os white papers e tutoriais é difícil, porque existem muitas lacunas que geralmente são preenchidas em livros didáticos. Se for importante, tenho um histórico matemático relativamente forte, como fiz meu Ph.D. em matemática aplicada (CFD para ser...
Digamos que temos os seguintes dados: set.seed(123) data <- data.frame(x = c(rnorm(50, 1, 1), rnorm(50, 5, 2)), y = c(rep('A', 50), rep('B', 50))) Qual produz o seguinte boxplot ( boxplot(data$x ~ data$y)): Agora, digamos que eu queira testar se as duas amostras têm os mesmos parâmetros de...
Eu já vi na literatura publicada (e postada aqui) que a eficiência relativa assintótica do teste de posto assinado de Wilcoxon é de pelo menos 0,864 quando comparada ao teste t. Também ouvi dizer que isso se aplica apenas a grandes amostras, embora alguns livros não mencionem isso (o que é...
Gostaria de codificar um kmeans agrupado em python usando pandas e scikit learn. Para selecionar o bom k, gostaria de codificar a estatística Gap de Tibshirani e al 2001 ( pdf ). Gostaria de saber se poderia usar o resultado de inércia_ do scikit e adaptar a fórmula estatística de gap sem precisar...
As bolsas geralmente exigem análise de energia para suportar um tamanho de amostra proposto. Na proteômica (e na maioria das ômicas), existem centenas de milhares de características / variáveis medidas em 10 das amostras (talvez 100, mas improvável). Além disso, sabe-se que algumas dessas...
Sabemos que, formalmente, as suposições de um teste não podem ser testadas, porque se escolhermos qual teste usar com base nos resultados do teste, o teste composto resultante terá propriedades desconhecidas (taxas de erro do tipo I e II). Eu acho que essa é uma das razões pelas quais o tipo "Six...