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

Я провожу небольшой эксперимент с регрессией LASSO в R, чтобы проверить, сможет ли она найти идеальную пару предикторов. Пара определяется следующим образом: f1 + f2 = исход Результатом здесь является предопределенный вектор, называемый «возраст». F1 и f2 создаются путем взятия половины вектора возраста и установки остальных значений в 0, например: …

20 r  regression  feature-selection  lasso  high-dimensional 

Меня интересует лучшее понимание дельта-метода для аппроксимации стандартных ошибок средних предельных эффектов регрессионной модели, включающей термин взаимодействия. Я посмотрел на связанные вопросы в рамках дельта-метода, но ни один из них не дал того, что я ищу. Рассмотрим следующие данные в качестве мотивирующего примера: set.seed(1) x1 <- rnorm(100) x2 <- rbinom(100,1,.5) …

20 r  regression  standard-error  effect-size  delta-method 

Пусть модель линейной регрессии, полученная функцией R lm, хотела бы знать, можно ли ее получить с помощью команды Mean Squared Error. У меня был СЛЕДУЮЩИЙ вывод примера > lm <- lm(MuscleMAss~Age,data) > sm<-summary(lm) > sm Call: lm(formula = MuscleMAss ~ Age, data = data) Residuals: Min 1Q Median 3Q Max …

20 r  regression  error 

У меня есть набор данных, состоящий из 5 функций: A, B, C, D, E. Все они являются числовыми значениями. Вместо кластеризации на основе плотности я хочу кластеризовать данные в виде дерева решений. Подход, который я имею в виду, выглядит примерно так: Алгоритм может делить данные на X исходных кластеров на …

20 r  machine-learning  clustering  cart 

В R, когда у меня есть (обобщенная) линейная модель ( lm, glm, gls, glmm, ...), как я могу проверить коэффициент (регрессионный наклон) в отношении любого другого значения , чем 0? В сводке модели результаты t-теста коэффициента автоматически сообщаются, но только для сравнения с 0. Я хочу сравнить его с другим …

20 r  regression  t-test 

Привет, я пытаюсь найти непараметрический эквивалент двухстороннего ANOVA (дизайн 3х4), который способен включать взаимодействия. Из моего прочтения в Zar 1984 г. «Биостатистический анализ» это возможно с использованием метода, предложенного в Scheirer, Ray и Hare (1976), однако, согласно другим публикациям в Интернете, было сделано заключение, что этот метод более не подходит …

20 r  anova  nonparametric  stata  interaction 

Я изучаю R и экспериментирую с анализом отклонений. Я бегал как kruskal.test(depVar ~ indepVar, data=df) и anova(lm(depVar ~ indepVar, data=dF)) Есть ли практическая разница между этими двумя тестами? Насколько я понимаю, они оба оценивают нулевую гипотезу о том, что популяции имеют одинаковое среднее значение.

20 r  anova  kruskal-wallis 

Я хочу реализовать алгоритм EM вручную , а затем сравнить его с результатами normalmixEMиз mixtoolsпакета. Конечно, я был бы счастлив, если бы они оба привели к одинаковым результатам. Основное упоминание - Джеффри МакЛахлан (2000), Модели конечных смесей . У меня плотность смеси двух гауссианов, в общем виде, логарифмическая вероятность определяется …

20 r  expectation-maximization  gaussian-mixture 

У меня есть 10-дневные данные ежедневных возвратов по 28 различным валютам. Я хотел бы извлечь первый основной компонент, но вместо того, чтобы использовать PCA в течение всех 10 лет, я хочу перенести двухлетнее окно, потому что поведение валют меняется, и поэтому я хочу это отразить. Однако у меня есть серьезная …

20 r  pca 

Я хотел бы понять, как генерировать интервалы прогнозирования для оценок логистической регрессии. Мне посоветовали следовать процедурам в Моделирующих двоичных данных Коллетта , 2-е издание, с.98-99. После реализации этой процедуры и сравнения ее с R predict.glm, я на самом деле думаю, что в этой книге показана процедура вычисления доверительных интервалов , …

20 r  regression  confidence-interval  logistic  prediction-interval 

Я немного новичок в изучении данных / машинного обучения / и т.д. и читали о нескольких способах объединения нескольких моделей и прогонов одной и той же модели для улучшения прогнозов. У меня сложилось впечатление, что после прочтения пары статей (которые часто интересны и хороши в теории и греческом языке, но …

20 r  random-forest  k-nearest-neighbour  ensemble 

Я использую библиотеку caretв R для тестирования различных процедур моделирования. trainControlОбъект позволяет указать метод повторной дискретизации. Эти методы описаны в документации разделе 2.3 , и включают в себя: boot, boot632, cv, LOOCV, LGOCV, repeatedcvи oob. Хотя некоторые из них легко вывести, не все эти методы четко определены. Какие процедуры соответствуют …

20 r  resampling  caret 

Я смотрел на моделирование смешанных эффектов с использованием пакета lme4 в R. Я в основном использую lmerкоманду, поэтому я задам свой вопрос через код, который использует этот синтаксис. Я предполагаю, что общий простой вопрос может состоять в том, можно ли сравнивать любые две модели, построенные с lmerиспользованием отношений правдоподобия на …

20 r  mixed-model  lme4-nlme  likelihood-ratio 

Иногда мы используем график спектральной плотности для анализа периодичности во временных рядах. Обычно мы анализируем сюжет путем визуального осмотра, а затем пытаемся сделать вывод о периодичности. Но разработали ли статистики какие-либо тесты для проверки того, являются ли какие-либо пики на графике статистически отличными от белого шума? Разработали ли R-специалисты какой-либо …

20 r  time-series  hypothesis-testing 

Я запустил полиномиальную модель логита в JMP и получил результаты, которые включали AIC, а также х-квадратные p-значения для каждой оценки параметра. Модель имеет один категорический результат и 7 категориальных объяснительных переменных. Затем я соответствую тому, что, как я думал, построю ту же модель в R, используя multinomфункцию из пакета nnet …

20 r  logistic  multinomial  logit  jmp