Perguntas com a marcação «r»

Estou trabalhando na estimativa da auto-regressão vetorial (VARs) e da função de resposta ao impulso (IRFs) com base nos dados do painel com 33 indivíduos em 77 trimestres. Como esse tipo de situação deve ser analisado? Quais algoritmos existem para esse fim? Eu preferiria realizar essas análises em R, portanto, …

9 r  econometrics  panel-data  var  impulse-response 

Estou interessado em encontrar um procedimento para simular dados que sejam consistentes com um modelo de mediação especificado. De acordo com a estrutura geral do modelo de equações estruturais lineares para testar modelos de mediação, primeiro descritos por Barron e Kenny (1986) e descritos em outros lugares como Judd, Yzerbyt …

9 r  regression  simulation  random-generation  mediation 

Em Dixon, Coles ( 1997 ), eles usaram a estimativa de máxima verossimilhança para os dois modelos independentes de Poisson modificados em (4.3) para modelar as pontuações no futebol. Estou tentando usar R para "reproduzir" os parâmetros alfa e beta, bem como os parâmetros de efeito doméstico (pág. 274, Tabela …

9 r  modeling  maximum-likelihood  games 

Estou analisando um conjunto de dados sobre comunidades entre marés. Os dados são cobertura percentual (de algas, cracas, mexilhões etc.) em quadrantes. Estou acostumado a pensar na análise de correspondência (CA) em termos de contagem de espécies e na análise de componentes principais (PCA) como algo mais útil para tendências …

9 r  pca  multivariate-analysis  correspondence-analysis 

Eu tenho os seguintes dados localizados aqui . Estou tentando calcular o intervalo de confiança de 95% sobre a pureza média quando a porcentagem de hidrocarbonetos é 1,0. Em R, insiro o seguinte. > predict(purity.lm, newdata=list(hydro=1.0), interval="confidence", level=.95) fit lwr upr 1 89.66431 87.51017 91.81845 No entanto, como posso obter …

9 r  regression  confidence-interval  prediction-interval 

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 …

9 r  hypothesis-testing  nonlinear-regression  goodness-of-fit  heteroscedasticity 

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 meus …

9 r  classification  unbalanced-classes  oversampling 

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 …

9 r  mixed-model  aic 

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 …

9 r  machine-learning  cross-validation  caret 

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 pontuações de três blocos …

9 r  meta-analysis  effect-size  meta-regression 

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 …

9 r  chi-squared  monte-carlo 

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 mensais por …

9 r  interpolation  fourier-transform 

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 …

9 r  random-forest  prediction  missing-data 

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 regressão (OLS) …

9 r  regression  standard-error  piecewise-linear 

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 …

9 r  multiple-regression  data-transformation  nonlinear-regression  gam