Eu tenho duas séries temporais de dados diários. Um é sign-upse o outro terminationsde assinaturas. Eu gostaria de prever o último usando as informações contidas nas duas variáveis. Observando o gráfico dessas séries, é óbvio que as terminações estão correlacionadas com os múltiplos de inscrições...
Versão curta: existe uma introdução aos escritos de Ronald Fisher (papéis e livros) sobre estatística que é voltada para aqueles com pouco ou nenhum conhecimento de estatística? Estou pensando em algo como um "leitor Fisher anotado" destinado a não estatísticos. Explico a motivação para esta...
Eu tenho lido sobre estimativa de máxima verossimilhança e estimativa máxima a posteriori e até agora encontrei exemplos concretos apenas com a estimativa de máxima verossimilhança. Encontrei alguns exemplos abstratos de estimativa máxima a posteriori, mas nada concreto ainda com números: S Pode...
Estou tentando usar o neuralnetpacote de R (documentação aqui ) para previsão. Aqui o que estou tentando fazer: library(neuralnet) x <- cbind(runif(50, min=1, max=500), runif(50, min=1, max=500)) y <- x[, 1] * x[, 2] train <- data.frame(x, y) n <- names(train) f <-...
Não sei como são chamadas essas tramas e, portanto, dei um título estúpido a essa pergunta. Digamos que eu tenha um conjunto de dados ordenado da seguinte maneira 4253 4262 4270 4383 4394 4476 4635 ... Cada número corresponde à quantidade de postagens que um determinado usuário contribuiu para...
Existe uma estratégia para escolher o número de árvores em um GBM? Especificamente, o ntreesargumento em R's gbmfunção. Não vejo por que você não deve definir ntreeso valor mais alto possível. Percebi que um número maior de árvores reduz claramente a variabilidade dos resultados de vários GBMs. Eu...
Usando Amelia em R, obtive vários conjuntos de dados imputados. Depois disso, realizei um teste de medidas repetidas no SPSS. Agora, quero reunir os resultados dos testes. Eu sei que posso usar as regras do Rubin (implementadas por meio de qualquer pacote de imputação múltipla em R) para agrupar...
Usei o seguinte código R para ajustar um modelo probit: p1 <- glm(natijeh ~ ., family=binomial(probit), data=data1) stepwise(p1, direction='backward/forward', criterion='BIC') Eu quero saber o que faz stepwisee backward/forwardfaz exatamente e como selecionar as
Eu tenho uma pergunta sobre a interpretação da chamada tsboot em R. Verifiquei a documentação do Kendall e do pacote de inicialização, mas não sou mais inteligente do que antes. Quando executo um bootstrap usando, por exemplo, o exemplo no pacote Kendall, onde a estatística de teste é a tau do...
Em um modelo hierárquico de dados , onde y ~ de Poisson ( λ ) λ ~ Gama ( α , β ) parece ser típico em prática a escolher os valores ( ct , β ) de tal modo que a média e variância de distribuição gama corresponder aproximadamente a média e variância dos dados y (por exemplo, Clayton e Kaldor,...
Estou executando uma simulação rápida para comparar diferentes métodos de cluster e atualmente enfrenta um problema ao tentar avaliar as soluções de cluster. Conheço várias métricas de validação (muitas encontradas em cluster.stats () em R), mas presumo que sejam melhor utilizadas se o número...
Eu apliquei alguns dados para encontrar a melhor solução de variáveis do modelo de regressão usando regressão de cumeeira em R. Eu usei lm.ridgee glmnet(quando alpha=0), mas os resultados são muito diferentes, especialmente quando lambda=0. Supõe-se que ambos os estimadores de parâmetros tenham...
Para não estatísticos como eu, é muito difícil capturar a ideia de VImétrica (variação de informação) mesmo depois de ler o artigo relevante de Marina Melia " Comparando agrupamentos - uma distância baseada em informação " (Journal of Multivariate Analysis, 2007). De fato, não estou familiarizado...
Gostaria de verificar Rse meus dados se encaixam nas distribuições log-normal ou Pareto. Como eu pude fazer isso? Talvez ks.testpudesse me ajudar a fazer isso, mas como eu poderia obter os parâmetros αα\alpha e kkk para a distribuição de Pareto para meus
Eu escrevi um programa para simular um embaralhamento de cartão overhand . Cada carta é numerada, com o naipe passando CLUBS, DIAMONDS, HEARTS, SPADESe a classificação de Dois a Dez, depois Valete, Rainha, Rei e Ás. Assim, o Two of Clubs tem um número 1, o Three of Clubs um 2 .... O Ace of Clubs...
Estou tentando calcular uma média e desvio padrão de 2 percentis para uma distribuição lognormal. Consegui realizar o cálculo de uma distribuição normal usando X = mean + sd * Ze resolvendo para média e sd. Acho que estou perdendo uma equação quando tento fazer o mesmo para uma distribuição...
Estou reproduzindo um exemplo dos modelos Generalizado, Linear e Misto . Meu MWE está abaixo: Dilution <- c(1/128, 1/64, 1/32, 1/16, 1/8, 1/4, 1/2, 1, 2, 4) NoofPlates <- rep(x=5, times=10) NoPositive <- c(0, 0, 2, 2, 3, 4, 5, 5, 5, 5) Data <- data.frame(Dilution, NoofPlates,...
Gostaria de corrigir manualmente um certo coeficiente, digamos , e depois ajustar os coeficientes a todos os outros preditores, mantendo β 1 = 1,0 no modelo.β1= 1,0β1=1.0\beta_1=1.0β1= 1,0β1=1.0\beta_1=1.0 Como posso conseguir isso usando R? Eu particularmente gostaria de trabalhar com o LASSO (...
Usando um biplot de valores obtidos através da análise de componentes principais, é possível explorar as variáveis explicativas que compõem cada componente principal. Isso também é possível com a Análise Discriminante Linear? Os exemplos fornecidos usam Os dados são "Dados de íris de Edgar...
Estou tentando avaliar o desempenho de um algoritmo de classificação de aprendizado de máquina supervisionado. As observações se enquadram em classes nominais (2 por enquanto, no entanto, eu gostaria de generalizar isso para problemas de várias classes), extraídos de uma população de 99...