Perguntas com a marcação «r»

Fiz uma pergunta sobre um modelo linear feito com R lm: "A regressão usou mínimos quadrados iterativos lineares ou não lineares?" Eu procurei um pouco e [acho que] entendi a diferença entre os dois, mas não consegui encontrar nenhuma evidência do uso de mínimos quadrados lineares por R lm(que é …

9 r  linear-model  lm 

Suponha que, em uma análise de regressão em R, eu tenha uma variável independente do tipo de fator com 3 níveis no meu conjunto de dados de trem. Mas no conjunto de dados de teste essa mesma variável de fator tem 5 níveis. Portanto, não posso prever os valores de …

9 r  regression  categorical-data  cross-validation  rare-events 

Meus dados brutos consistem em uma série temporal de 60 dias com uma tendência de queda. Os dados são semanais, portanto a frequência é definida como 7. Eu calculei a diferença dos dados que se parece com isso Quando executo gráficos de ACF e PACF sobre a diferença, pareço obter …

9 r  time-series  autocorrelation 

Ao selecionar um número apropriado de nós para um GAM, convém levar em consideração o número de dados e incrementos no eixo x. E se tivermos 100 incrementos no eixo x com 1000 pontos de dados em cada incremento. A informação aqui diz: Se eles não forem fornecidos, os nós …

9 r  gam  splines 

Eu tenho três variáveis ​​macroeconômicas (ICS - sentimento do consumidor, ER - taxa de emprego, DGO - pedido de bens duráveis) e executei testes de causalidade de Granger em R sobre eles. Eu realmente não sei como interpretar os resultados de um teste de Granger. Alguém poderia me ajudar com …

9 r  granger-causality 

Estou com problemas para obter resultados equivalentes entre um aovmodelo de medidas repetidas entre dentro e um lmermodelo misto. Meus dados e script têm a seguinte aparência data=read.csv("https://www.dropbox.com/s/zgle45tpyv5t781/fitness.csv?dl=1") data$id=factor(data$id) data id FITNESS TEST PULSE 1 1 pilates CYCLING 91 2 2 pilates CYCLING 82 3 3 pilates CYCLING 65 4 …

9 r  anova  mixed-model  repeated-measures  lme4-nlme 

Quero calcular as taxas de incidência para apresentar as taxas de risco, a fim de apresentar medidas de risco relativas e absolutas. Vi em outros estudos que essas taxas de incidência podem ser calculadas usando modelos de Poisson com tempo de acompanhamento no modelo como um deslocamento. Então, eu tentei …

9 r  survival  cox-model  poisson-regression  incidence-rate-ratio 

library(lme4) out <- glmer(cbind(incidence, size - incidence) ~ period + (1 | herd), data = cbpp, family = binomial, contrasts = list(period = "contr.sum")) summary(out) Fixed effects: Estimate Std. Error z value Pr(>|z|) (Intercept) -2.32337 0.22129 -10.499 < 2e-16 *** period1 0.92498 0.18330 5.046 4.51e-07 *** period2 -0.06698 0.22845 -0.293 …

9 r  binomial  logit  lme4-nlme 

Pelo que entendi, o argumento cp para a rpartfunção ajuda a pré-podar a árvore da mesma maneira que os argumentos minsplit ou minbucket. O que não entendo é como os valores de CP são calculados. Por exemplo df<-data.frame(x=c(1,2,3,3,3,4), y=as.factor(c(TRUE, TRUE, FALSE, TRUE, FALSE, FALSE)), method="class") mytree<-rpart(y ~ x, data = …

9 r  cart  rpart 

Estou lutando com pontos de projeção na análise discriminante linear (LDA). Muitos livros sobre métodos estatísticos multivariados ilustram a idéia do LDA com a figura abaixo. A descrição do problema é a seguinte. Primeiro, precisamos traçar o limite da decisão, adicionar uma linha perpendicular e, em seguida, plotar as projeções …

9 r  machine-learning  classification  discriminant-analysis 

Aqui estão os dados: > tires <- data.frame(Wear = c(17, 14, 12, 13, 14, 14, 12, 11, 13, 13, 10, 11, 13, 8, 9, 9), Brand = rep(LETTERS[1:4], 4), Car = as.character(as.roman(rep(1:4, each = 4)))) > tires Wear Brand Car 1 17 A I 2 14 B I 3 12 …

9 r  anova  dataset  degrees-of-freedom 

Eu faço uma regressão linear usando a função R lm: x = log(errors) plot(x,y) lm.result = lm(formula = y ~ x) abline(lm.result, col="blue") # showing the "fit" in blue mas não se encaixa bem. Infelizmente, não consigo entender o manual. Alguém pode me apontar na direção certa para ajustar isso …

9 r  regression 

É possível plotar um gráfico de dependência parcial para exibir a probabilidade da classe e estimar os efeitos de um preditor para um modelo GBM ? Algo semelhante ao partialPlotdo randomForestpacote. De acordo com este artigo , uma plotagem parcial é possível com o gbm. Agradeço antecipadamente por sua ajuda.

9 r  gbm 

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 …

9 r  pdf  kernel-smoothing  splines 

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 localização (mediana e …

9 r  nonparametric  wilcoxon-mann-whitney