Запитання з тегом «r»

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

20 r  regression  feature-selection  lasso  high-dimensional 

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

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

Нехай модель лінійної регресії, отримана функцією R, хотіла б знати, чи можливо її отримати за допомогою команди Середня квадратична помилка. У мене був наступний результат прикладу > lm <- lm(MuscleMAss~Age,data) > sm<-summary(lm) > sm Call: lm(formula = MuscleMAss ~ Age, data = data) Residuals: Min 1Q Median 3Q Max -16.1368 …

20 r  regression  error 

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

20 r  regression  t-test 

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

20 r  machine-learning  clustering  cart 

Привіт, я намагаюся знайти непараметричний еквівалент двосторонньої ANOVA (3x4 конструкції), яка здатна включати взаємодії. З мого читання в Зарі 1984 р. "Біостатистичного аналізу" це можливо за допомогою методу, викладеного в Шірера, Рея і Зайця (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 років, я хочу застосувати 2-річне вікно, тому що поведінка валют змінюється, і тому я хочу це відобразити. Однак у мене є основна проблема, тобто …

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 в Р. Я в основному використовую lmerкоманду, тому я буду ставити своє запитання через код, який використовує цей синтаксис. Я думаю, що може бути загальним простим питанням: чи добре порівнювати будь-які дві моделі, побудовані з lmerвикористанням коефіцієнтів ймовірності на основі однакових …

20 r  mixed-model  lme4-nlme  likelihood-ratio 

Ми іноді використовуємо графік спектральної щільності для аналізу періодичності у часових рядах. Зазвичай ми аналізуємо сюжет шляхом візуального огляду, а потім намагаємось зробити висновок про періодичність. Але чи розробив статистик будь-який тест, щоб перевірити, чи якісь шипи на ділянці статистично відрізняються від білого шуму? Чи розробили R-експерти який-небудь пакет для …

20 r  time-series  hypothesis-testing 

Я запустив багаточленну модель Logit у JMP і отримав результати, які включали AIC, а також значення c-квадрата p для кожної оцінки параметрів. Модель має один категоричний результат та 7 категоричних пояснювальних варіантів. Потім я підходив до того, що, як я думав, створив би ту саму модель в R, використовуючи multinomфункцію …

20 r  logistic  multinomial  logit  jmp