enquanto passo para a previsão com os modelos ARIMA, estou tentando entender como posso melhorar uma previsão com base no ajuste do ARIMA à sazonalidade e à deriva. Meus dados são as seguintes séries temporais (mais de 3 anos, com clara tendência ascendente e sazonalidade visível, o que parece...
Portanto, estou trabalhando com modelos de regressão logística em R. Embora ainda seja novo em estatística, sinto que já compreendi um pouco os modelos de regressão, mas ainda há algo que me incomoda: Observando a imagem vinculada, você vê o resumo R impresso para um modelo de exemplo que eu...
Eu gostaria de gerar dados aleatórios a partir de uma distribuição normal restrita usando R. Por exemplo, eu posso querer simular uma variável de uma distribuição normal com mean=3, sd= 2e quaisquer valores maiores que 5 são reamostrados da mesma distribuição normal. Assim, para a função geral,...
Em R, se eu chamar a lm()função da seguinte maneira: lm.1 = lm(response ~ var1 + var2 + var1 * var2) summary(lm.1) Isso me dá um modelo linear da variável resposta com var1, var2e a interação entre eles. No entanto, como exatamente interpretamos numericamente o termo de interação? A...
Modelos penalizados podem ser usados para estimar modelos em que o número de parâmetros é igual ou até maior que o tamanho da amostra. Essa situação pode surgir em modelos log-lineares de grandes tabelas esparsas de dados categóricos ou de contagem. Nessas configurações, geralmente também é...
Estou procurando a equação correta para calcular a covariância de amostra imparcial ponderada. As fontes da Internet são bastante raras nesse tema e todas usam equações diferentes. A equação mais provável que encontrei é esta: qj k= ∑Ni = 1WEu( ∑Ni = 1WEu)2- ∑Ni = 1W2Eu∑Ni = 1WEu( xeu j- x¯j) (...
EDIT: Desde que fiz este post, segui com um post adicional aqui . Resumo do texto abaixo: Estou trabalhando em um modelo e tentei regressão linear, transformações de Box Cox e GAM, mas não progredi muito Usando R, Eu estou trabalhando atualmente em um modelo para prever o sucesso de jogadores...
A situação Eu tenho um conjunto de dados com um dependente uma variável independente . Quero ajustar uma regressão linear por partes contínua com pontos de interrupção conhecidos / fixos que ocorrem em . Os breakpoins são conhecidos sem incerteza, então não quero calculá-los. Então eu ajuste uma...
Eu tenho um bom randomForestmodelo de classificação que gostaria de usar em um aplicativo que prevê a classe de um novo caso. O novo caso tem inevitavelmente valores ausentes. Prever não funcionará como tal para NAs. Como devo fazer isso então? data(iris) # create first the new case with missing...
Gostaria de entender o uso da simulação de Monte Carlo na chisq.test()função em R. Eu tenho uma variável qualitativa que tem 128 níveis / classes. O tamanho da minha amostra é 26 (não pude provar mais "indivíduos"). Então, obviamente, terei alguns níveis com 0 "indivíduos". Mas o fato é que eu...
fundo Em um artigo de Epstein (1991): Ao obter valores climatológicos diários a partir de médias mensais , são dadas a formulação e um algoritmo para calcular a interpolação de Fourier para valores periódicos e com espaçamento uniforme. No artigo, o objetivo é obter valores diários de médias...
Eu não sou especialista, então me perdoe se parte da terminologia for um pouco desajeitada. É um prazer fornecer mais informações quando necessário. Eu tenho dois vetores de 50 valores numéricos emparelhados em R. Eu quero executar um teste de randomização ou permutação bicaudal para determinar se...
Estou conduzindo uma meta-análise dos tamanhos de efeito d em R usando o pacote metafor. d representa diferenças nos escores de memória entre pacientes e saudáveis. No entanto, alguns estudos relatam apenas subescores da medida de interesse d (por exemplo, várias pontuações de memória diferentes ou...
Talvez essa seja uma pergunta tola, mas, ao gerar um modelo com sinal de intercalação e usar algo como LOOCVou (ainda mais direto ao ponto) LGOCV, qual é o benefício de dividir dados em conjuntos de trens e testes, se é basicamente isso que a etapa de validação cruzada faz assim mesmo? Eu li...
Eu tenho modelos como este: require(nlme) set.seed(123) n <- 100 k <- 5 cat <- as.factor(rep(1:k, n)) cat_i <- 1:k # intercept per kategorie x <- rep(1:n, each = k) sigma <- 0.2 alpha <- 0.001 y <- cat_i[cat] + alpha * x + rnorm(n*k, 0, sigma) plot(x, y) m1 <- lm(y ~...
Estou pensando em um problema que é prever o log (gasto) de um cliente usando regressão linear. Estou pensando em quais recursos usar como entrada e me perguntando se seria bom usar o percentil de uma variável como entradas. Por exemplo, eu poderia usar a receita das empresas como insumo. O que...
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...
Estou tentando usar o SMOTE para corrigir o desequilíbrio no meu problema de classificação de várias classes. Embora o SMOTE funcione perfeitamente no conjunto de dados da íris conforme o documento de ajuda do SMOTE, ele não funciona em um conjunto de dados semelhante. Aqui está a aparência dos...
Eu tenho dados para falhas de veículos a motor por hora do dia. Como seria de esperar, eles são altos no meio do dia e atingem o pico na hora do rush. O padrão geom_density do ggplot2 suaviza muito bem Um subconjunto dos dados, para falhas relacionadas ao consumo de bebida, é alto no final do dia...
Eu tenho um exemplo de conjunto de dados da seguinte maneira: Volume <- seq(1,20,0.1) var1 <- 100 x2 <- 1000000 x3 <- 30 x4 = sqrt(x2/pi) H = x3 - Volume r = (x4*H)/(H + Volume) Power = (var1*x2)/(100*(pi*Volume/3)*(x4*x4 + x4*r + r*r)) Power <- jitter(Power, factor = 1, amount...