No Quick-R de Robert Kabacoff, tenho # Bootstrap 95% CI for regression coefficients library(boot) # function to obtain regression weights bs <- function(formula, data, indices) { d <- data[indices,] # allows boot to select sample fit <- lm(formula, data=d) return(coef(fit)) } #...
Estou fazendo análise de séries temporais usando R. Tenho que decompor meus dados em tendência, componente sazonal e aleatório. Eu tenho dados semanais por 3 anos. Eu encontrei duas funções em R - stl()e decompose(). Eu li que isso stl()não é bom para decomposição multiplicativa. Alguém pode me...
Sou relativamente novo em R e estou tentando ajustar um modelo aos dados que consistem em uma coluna categórica e uma coluna numérica (número inteiro). A variável dependente é um número contínuo. Os dados têm o seguinte formato: predCateg, predIntNum, ResponseVar Os dados são mais ou menos...
Ouvi falar um pouco sobre o uso de redes neurais para prever séries temporais. Como posso comparar, qual método para prever minhas séries temporais (dados diários de varejo) é melhor: auto.arima (x), ets (x) ou nnetar (x). Eu posso comparar auto.arima com ets por AIC ou BIC. Mas como posso...
A fórmula que você precisa especificar para treinar um modelo multinível (usando lmerda lme4 Rbiblioteca) sempre me entende. Eu li inúmeros livros e tutoriais, mas nunca entendi direito. Então, aqui está um exemplo deste tutorial que eu gostaria de ver formulado em uma equação. Estamos tentando...
Eu tenho oito variáveis independentes e uma dependente. Eu corri uma matriz de correlação, e 5 deles têm uma baixa correlação com o DV. Em seguida, executei uma regressão múltipla passo a passo para ver se algum / todos os IVs podem prever o DV. A regressão mostrou que apenas dois IVs podem...
A figura abaixo (Figura 1 da p. 646 deste artigo ) compara os valores observados com os valores esperados na distribuição de Poisson. Em seguida, ele executa um teste qui-quadrado para verificar se os valores observados diferem dos valores esperados na distribuição de Poisson. Usando R, como é...
Eu tenho algumas experiências com modelagem de séries temporais, na forma de modelos ARIMA simples e assim por diante. Agora eu tenho alguns dados que exibem clustering de volatilidade e gostaria de tentar começar ajustando um modelo GARCH (1,1) nos dados. Eu tenho uma série de dados e acho que...
Eu estou analisando alguns dados comportamentais usando lme4em R, principalmente após excelentes tutoriais do Bodo de Inverno , mas eu não entendo se eu estou lidando com interações corretamente. Pior, ninguém mais envolvido nesta pesquisa usa modelos mistos, então estou um pouco à deriva quando se...
Estou trabalhando em um algoritmo que precisa calcular o tamanho de um conjunto gerado pelas interseções de pelo menos 2 conjuntos. Mais especificamente: z= | UMA0 0∩ … ∩ An|z=|A0∩…∩An| z = \left |A_0 \cap \ldots \cap A_n \right | Os conjuntos que são cruzados são gerados por consultas SQL e, em...
Foi-me perguntado por um revisor se as correlações de Pearson (valores-r) apresentadas em uma tabela podem ser comparadas entre si, de modo que se possa afirmar que uma é "mais forte" que a outra (além de apenas observar os valores-reais) . Como você faria isso? Eu encontrei esse...
Minha situação: tamanho pequeno da amostra: 116 variável de resultado binário lista longa de variáveis explicativas: 44 variáveis explicativas não vieram do alto da minha cabeça; sua escolha foi baseada na literatura. a maioria dos casos na amostra e a maioria das variáveis têm valores...
Sou relativamente novo em estatística e R. Gostaria de conhecer o processo para determinar os parâmetros ARIMA para meu conjunto de dados. Você pode me ajudar a descobrir o mesmo usando R e teoricamente (se possível)? Os dados variam de 12 de janeiro a 14 de março e retratam as vendas mensais....
Estou trabalhando com alguns dados do mundo real e os modelos de regressão estão produzindo alguns resultados contra-intuitivos. Normalmente confio nas estatísticas, mas, na realidade, algumas dessas coisas não podem ser verdadeiras. O principal problema que estou vendo é que um aumento em uma...
Estou usando ccfpara encontrar uma correlação entre duas séries temporais. Estou recebendo uma trama que se parece com isso: Observe que estou interessado principalmente na correlação para o atraso = 0. Questões: Interprete corretamente que existe uma correlação cruzada para o atraso = 0, pois...
Estou tendo problemas para entender como a varImpfunção funciona para um modelo randomForest com o caretpacote. No exemplo abaixo, o recurso var3 recebe zero importância usando a varImpfunção de cursor , mas o modelo final randomForest subjacente tem importância diferente de zero para o recurso...
Toda a literatura de Modelagem de Distribuição de Espécies sugere que, ao prever a presença / ausência de uma espécie usando um modelo que gera probabilidades (por exemplo, RandomForests), é importante escolher a probabilidade de limiar pela qual classificar uma espécie como presença ou ausência é...
Seja uma amostra aleatória de uma distribuição para . Ou seja,X1,...,XnX1,...,Xn X_1,
Como posso obter efeitos aleatórios agrupados para lmer após imputação múltipla? Estou usando mouses para atribuir múltiplos quadros de dados. E lme4 para um modelo misto com interceptação aleatória e inclinação aleatória. O pool de pool funciona bem, exceto que ele não combina os efeitos...
Espero sinceramente que tenha formulado esta pergunta de forma que ela possa ser respondida definitivamente - caso contrário, informe-me e tentarei novamente! Eu também acho que devo usar R para essas análises. Eu tenho várias medidas plant performance (Ys)que eu suspeito que foram influenciadas...