Вопросы с тегом «r»

Я рассчитал автокорреляцию по данным временного ряда о характере движения рыбы по ее положениям: X ( x.ts) и Y ( y.ts). Используя R, я запустил следующие функции и создал следующие графики: acf(x.ts,100) acf(y.ts,100) Мой вопрос, как мне интерпретировать эти сюжеты? Какая информация необходима, чтобы сообщить о каком-либо шаблоне? Я занимался …

13 r  time-series  autocorrelation 

Я хотел бы использовать вменение для замены отсутствующих значений в моем наборе данных при определенных ограничениях. Например, я бы хотел, чтобы вмененная переменная x1была больше или равна сумме двух других моих переменных, скажем, x2и x3. Я также хочу x3быть вмененным либо 0или, >= 14и я хочу x2быть вмененным либо 0или …

13 r  spss  missing-data  multiple-imputation 

Из моих результатов видно, что GLM Gamma отвечает большинству допущений, но стоит ли это значительного улучшения по сравнению с лог-преобразованным LM? Большая часть литературы, которую я нашел, посвящена пуассоновским или биномиальным GLM. Я нашел статью ОЦЕНКА ОБОБЩЕНИЙ ОБОБЩЕННОЙ ЛИНЕЙНОЙ МОДЕЛИ С ИСПОЛЬЗОВАНИЕМ Рандомизации очень полезной, но в ней отсутствуют реальные …

13 r  generalized-linear-model  model-selection  gamma-distribution 

Я хотел бы сравнить модели, выбранные с ребристой, лассо и эластичной сеткой. На рисунке ниже показаны коэффициенты пути, используя все 3 метода: гребень (рис. A, альфа = 0), лассо (рис. B; альфа = 1) и эластичная сетка (рис. C; альфа = 0,5). Оптимальное решение зависит от выбранного значения лямбда, которое …

13 r  multiple-regression  modeling  regularization 

Я хочу вычислить предсказанную функцию выживания для модели пропорциональных рисков Кокса с ненадежными терминами [используя пакет выживания]. Похоже, что когда в модели присутствуют слабые члены, предсказанная функция выживания не может быть вычислена. ## Example require(survival) data(rats) ## Create fake weight set.seed(90989) rats$weight<-runif(nrow(rats),0.2,0.9) ## Cox model with gamma frailty on litter …

13 r  survival  cox-model  frailty 

Я ищу какой-нибудь открытый исходный код или существующую библиотеку, которую я могу использовать. Насколько я знаю, пакет glmnet не очень легко расширяется, чтобы охватить неотрицательный случай. Я могу ошибаться, Любой, у кого есть идеи, высоко ценится. Под неотрицательным я подразумеваю, что все коэффициенты должны быть положительными (> 0).

13 r  lasso 

Лучше ли различать ряды (если это необходимо) перед использованием Arima ИЛИ лучше использовать параметр d в Arima? Я был удивлен тем, насколько разные подходящие значения зависят от того, какой маршрут выбран с той же моделью и данными. Или я что-то делаю неправильно? install.packages("forecast") library(forecast) wineindT<-window(wineind, start=c(1987,1), end=c(1994,8)) wineindT_diff <-diff(wineindT) #coefficients …

13 r  time-series  arima 

У меня есть набор данных с двумя перекрывающимися классами, семь точек в каждом классе, точки находятся в двухмерном пространстве. В R, и я бегу svmиз e1071пакета, чтобы построить разделяющую гиперплоскость для этих классов. Я использую следующую команду: svm(x, y, scale = FALSE, type = 'C-classification', kernel = 'linear', cost = …

13 r  machine-learning  svm  libsvm  e1071 

Я хочу приспособить смешанную модель, используя lme4, nlme, пакет байсовой регрессии или любой доступный. Смешанная модель в соглашениях о кодировании Asreml-R прежде чем углубляться в детали, мы можем захотеть узнать подробности о соглашениях asreml-R для тех, кто не знаком с кодами ASREML. y = Xτ + Zu + e ........................(1) …

13 r 

Закрыто. Этот вопрос не по теме . В настоящее время он не принимает ответы. Хотите улучшить этот вопрос? Обновите вопрос, чтобы он соответствовал теме перекрестной проверки. Закрыто 4 года назад . Есть ли пакет R, который я могу использовать, чтобы выяснить, существуют ли отношения между переменными? Обычно, когда я ищу …

13 r  data-visualization  correlation  eda 

Я хочу применить PCA к набору данных, который состоит из переменных смешанного типа (непрерывных и двоичных). Чтобы проиллюстрировать процедуру, я вставил минимальный воспроизводимый пример в R ниже. # Generate synthetic dataset set.seed(12345) n <- 100 x1 <- rnorm(n) x2 <- runif(n, -2, 2) x3 <- x1 + x2 + rnorm(n) …

13 r  pca 

Я пытаюсь выяснить, как R вычисляет автокорреляцию lag-k (очевидно, это та же формула, что используется Minitab и SAS), так что я могу сравнить ее с использованием функции CORREL в Excel, примененной к серии, и ее версии с k-lagged. R и Excel (используя CORREL) дают немного разные значения автокорреляции. Мне также …

13 r  sas  autocorrelation  excel 

Я работаю над множественной логистической регрессией в R, используя glm. Переменные предиктора являются непрерывными и категориальными. Выдержка из резюме модели показывает следующее: Coefficients: Estimate Std. Error z value Pr(>|z|) (Intercept) 2.451e+00 2.439e+00 1.005 0.3150 Age 5.747e-02 3.466e-02 1.658 0.0973 . BMI -7.750e-02 7.090e-02 -1.093 0.2743 ... --- Signif. codes: 0 …

13 r  logistic  interpretation  p-value 

У меня проблемы с выяснением, что plot.stlименно означают диапазоны . Я нашел пост Гэвина по этому вопросу и прочитал документацию, я понимаю, что они говорят об относительной величине разложенных компонентов, но все же я не совсем уверен, как они работают. Например: данные: крошечный столбик, без шкалы сезонный: полный столбик, с …

13 r  time-series 

Как определить, значимы ли корреляции при разных лагах, полученные из взаимной корреляции (функция ccf) двух временных рядов?

13 r  statistical-significance  cross-correlation