Eu tenho um conjunto de dados mensais de temperatura da superfície do mar (SST) e quero aplicar alguma metodologia de cluster para detectar regiões com padrões SST semelhantes. Eu tenho um conjunto de arquivos de dados mensais em execução de 1985 a 2009 e quero aplicar o cluster a cada mês como...
Gostaria de saber se existe uma fórmula de tamanho de amostra como a fórmula de Lehr que se aplica a um teste F? A fórmula de Lehr para testes t é n=16/Δ2n=16/Δ2n = 16 / \Delta^2 , onde ΔΔ\Delta é o tamanho do efeito ( por exemplo, Δ=(μ1−μ2)/σΔ=(μ1−μ2)/σ\Delta = (\mu_1 - \mu_2) / \sigma ). Isso...
Estou muito interessado no potencial da análise estatística para simulação / previsão / estimativa de funções, etc. No entanto, eu não sei muito sobre isso e meu conhecimento matemático ainda é bastante limitado - eu sou um estudante de graduação em engenharia de software. Estou procurando um...
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 . Gostaria de plotar quatro gráficos de barras...
Estou tentando executar uma Anova de medidas repetidas em R, seguida de alguns contrastes específicos nesse conjunto de dados. Eu acho que a abordagem correta seria usar Anova()o pacote do carro. Vamos ilustrar minha pergunta com o exemplo retirado do ?Anovauso dos OBrienKaiserdados (Nota: omiti...
Eu tenho uma série temporal binária com 1 quando o carro não está em movimento e 0 quando o carro está em movimento. Quero fazer uma previsão para um horizonte de tempo de até 36 horas à frente e para cada hora. Minha primeira abordagem foi usar um Naive Bayes usando as seguintes entradas: t-24...
Estou tentando configurar um modelo de poisson inflado com zero no R e JAGS. Eu sou novo no JAGS e preciso de algumas orientações sobre como fazer isso. Eu tenho tentado com o seguinte onde y [i] é a variável observada model { for (i in 1:I) { y.null[i] <- 0 y.pois[i] ~ dpois(mu[i])...
Bloqueado . Esta pergunta e suas respostas estão bloqueadas porque a questão está fora de tópico, mas tem um significado histórico. No momento, não está aceitando novas respostas ou interações. Eu tenho um conjunto de dados ascii que consiste em três colunas, mas apenas...
Como aproximar a seguinte integral usando a simulação MC? ∫1−1∫1−1|x−y|dxdy∫−11∫−11|x−y|dxdy \int_{-1}^{1} \int_{-1}^{1} |x-y| \,\mathrm{d}x \,\mathrm{d}y Obrigado! Editar (em algum contexto): estou tentando aprender como usar a simulação para aproximar integrais e estou praticando alguma...
Estou tentando criar um SVM a partir de dados de treinamento em que um grupo é representado mais que o outro. No entanto, os grupos serão igualmente representados nos eventuais dados do teste. Portanto, eu gostaria de usar o class.weightsparâmetro da e1071interface do pacote R libsvmpara equilibrar...
Eu quero ajustar o modelo misto usando lme4, nlme, pacote de regressão baiano ou qualquer outro disponível. Modelo misto nas convenções de codificação Asreml-R Antes de entrar em detalhes, talvez desejemos ter detalhes sobre as convenções asreml-R, para aqueles que não estão familiarizados com...
Recentemente, fiz um curso introdutório sobre modelagem multinível. A maioria dos conjuntos de dados e exemplos que usamos foram das ciências sociais. Acabei de fazer um estágio de duas semanas em um departamento de bioestatística, onde eles querem que eu inicie um projeto referente à variação no...
Atualmente trabalhando no Octave, mas devido à documentação insuficiente, o progresso é muito lento. Qual idioma é fácil de aprender e usar, e bem documentado para resolver problemas de aprendizado de máquina? Estou procurando um protótipo em um pequeno conjunto de dados (milhares de exemplos),...
Alguém pode me dizer como fazer modelagem de recência, frequência e valor monetário (RFM) e modelagem de valor para cliente em R? Além disso, alguém pode me indicar alguma literatura sobre
A pergunta de ontem Determinar a precisão do modelo que estima a probabilidade de evento me deixou curioso sobre a pontuação de probabilidade. A pontuação Brier é uma medida de erro médio quadrático. O análogo significa a medida do desempenho do erro absoluto11N∑i = 1N( p r e di c t i o nEu- r e...
Desculpe se esta é uma pergunta nova; Estou tentando me ensinar estatística pela primeira vez. Acho que tenho o procedimento básico inativo, mas estou lutando para executá-lo com R. Então, eu estou tentando avaliar o significado dos coeficientes de regressão em uma regressão linear múltipla da...
Isso também foi solicitado na Ciência da Computação. Estou tentando calcular uma estimativa bayesiana de alguns coeficientes para uma regressão automática, com 11 amostras de dados: Yi=μ+α⋅Yi−1+ϵiYi=μ+α⋅Yi−1+ϵi Y_{i} = \mu + \alpha\cdot{}Y_{i-1} + \epsilon_{i} ondeϵiϵi\epsilon_{i} é normal de...
Eu sou muito novo em mínimos quadrados parciais (PLS) e tento entender a saída da função R plsr()no plspacote. Vamos simular dados e executar o PLS: library(pls) n <- 50 x1 <- rnorm(n); xx1 <- scale(x1) x2 <- rnorm(n); xx2 <- scale(x2) y <- x1 + x2 + rnorm(n,0,0.1); yy <-...
Estou executando o seguinte teste de raiz unitária (Dickey-Fuller) em uma série temporal usando a ur.df()função no urcapacote. O comando é: summary(ur.df(d.Aus, type = "drift", 6)) A saída é: ############################################### # Augmented Dickey-Fuller Test Unit Root Test #...
Em R, se eu definir.seed () e usar a função de amostra para randomizar uma lista, posso garantir que não gerarei a mesma permutação? ie ... set.seed(25) limit <- 3 myindex <- seq(0,limit) for (x in seq(1,factorial(limit))) { permutations <- sample(myindex) print(permutations) } Isso...