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

Ортогональні многочлени в одновимірному наборі точок - це многочлени, які виробляють значення в цих точках таким чином, що його крапковий добуток і попарна кореляція дорівнюють нулю. R може виробляти ортогональні многочлени з функцією poly . Таку ж функцію має і варіант поліма, який виробляє ортогональні многочлени на множинні множинні точки. …

12 r  multiple-regression  polynomial  orthogonal 

Я хотів би отримати 95% довірчі інтервали для прогнозів нелінійної змішаної nlmeмоделі. Оскільки для цього не передбачено нічого стандартного nlme, мені було цікаво, чи правильно використовувати метод "інтервалів прогнозування населення", як це викладено в книзі книги Бен Болкер в контексті моделей, що відповідають максимальній вірогідності , заснованих на ідеї перекомпонування …

12 r  mixed-model  confidence-interval  lme4-nlme 

Я деякий час використовував лінійні моделі для виконання тестів на 2 вибірки, але зрозумів, що це може бути не зовсім коректно. Здається, що використання узагальненої лінійної моделі з двочленним сімейством + ідентифікаційним зв’язком дає саме непідкупчені 2-вибіркові пропорційні результати тесту. Однак використання лінійної моделі (або glm з гауссова сім'я) дає …

12 r  hypothesis-testing  generalized-linear-model  proportion 

Наприклад, у R, якщо ви викликаєте acf()функцію, вона за замовчуванням будує корелограму та малює довірчий інтервал 95%. Якщо ви зателефонували до коду plot(acf_object, ci.type="white"), ви побачите: qnorm((1 + ci)/2)/sqrt(x$n.used) як верхня межа для типу «білий шум». Чи може хтось пояснити теорію, що стоїть за цим методом? Чому ми отримуємо qnorm …

12 r  confidence-interval  autocorrelation 

Я намагаюся розібратися, як виявити кількість складів у корпусі аудіозаписів. Я думаю, що хороший проксі може бути піками у файлі хвиль. Ось що я спробував із файлом про мене, що розмовляв англійською мовою (мій фактичний випадок використання - у Kiswahili). Стенограма цього прикладу запису: "Це я намагаюся використовувати функцію таймера. …

12 r  signal-processing 

Я думав, що розумію це питання, але зараз я не такий впевнений і хотів би поговорити з іншими, перш ніж продовжувати. У мене є дві змінні, Xі Y. Yє співвідношенням, і воно не обмежене 0 і 1 і, як правило, нормально розподілене. Xє пропорцією, і вона обмежена 0 і 1 …

12 r  regression  data-transformation  nonlinear-regression  proportion 

Я зробив вимірювань двох змінних і . Вони обоє мають відомі невизначеності та пов'язані з ними. Я хочу знайти співвідношення між і . Як я можу це зробити?x y σ x σ y x yнnnхxxуyyσхσx\sigma_xσуσy\sigma_yхxxуyy EDIT : кожен має різні пов'язані з ним, і те ж саме, що і .σ …

12 r  regression  deming-regression 

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

12 r  count-data  poisson-regression  overdispersion  quasi-likelihood 

Мені хотілося б знати, як перетворити негативні значення Log(), оскільки у мене є гетерокедастичні дані. Я прочитав, що це працює з формулою, Log(x+1)але це не працює з моєю базою даних, і я продовжую отримувати NaNs в результаті. Наприклад, я отримую це Попереджувальне повідомлення (я не помістив повну базу даних, тому …

12 r  logarithm 

Я просто думав про акуратний (не обов'язково хороший) спосіб створення однієї розмірної оцінки щільності, і моє питання: Чи має цей метод оцінки щільності назву? Якщо ні, то це особливий випадок якогось іншого методу в літературі? Ось метод: У нас є вектор який, як ми вважаємо, виведений з якогось невідомого розподілу, …

12 r  estimation  nonparametric  pdf  kernel-smoothing 

TL; DR: lme4оптимізація виявляється лінійним по числу параметрів моделі за замовчуванням, і шлях повільніше , ніж аналогічна glmмодель з фіктивними змінними для груп. Чи можна щось зробити, щоб прискорити це? Я намагаюся вписати досить велику ієрархічну модель logit (~ 50k рядків, 100 стовпців, 50 груп). Встановлення нормальної моделі logit до …

12 r  mixed-model  optimization  lme4-nlme  numerics 

Як я читаю результати тесту Данна ? Зокрема, про що мені відповідають значення, наведені в таблиці нижче? У мене є непараметричні дані у 4 групах, і я спершу зробив тест Крускала-Уолліса, щоб підтвердити, що розподіли груп не відрізняються один від одного та сукупного набору даних. Потім я використав тест Данна, …

12 r  distributions  kruskal-wallis  dunn-test 

Скажімо, у мене є наступні чотири розмірні дані, де перші три можна вважати координатами, а останню можна вважати значеннями. c1, c2, c3, value 1, 2, 6, 0.456 34, 34, 12 0.27 12, 1, 66 0.95 Як краще візуалізувати вплив перших трьох координат на останнє значення? Я знаю три методи. Один …

12 r  data-visualization 

Я розробляю логістичну модель регресії на основі ретроспективних даних із національної бази травм травми голови у Великобританії. Основним результатом є смертність за 30 днів (позначена як міра "вижити"). Інші заходи з опублікованими доказами значного впливу на результати попередніх досліджень включають: Year - Year of procedure = 1994-2013 Age - Age …

12 r  logistic  generalized-linear-model  goodness-of-fit  regression-strategies 

R version 3.1.1 (2014-07-10) -- "Sock it to Me" > bl <- c(140, 138, 150, 148, 135) > fu <- c(138, 136, 148, 146, 133) > t.test(fu, bl, alternative = "two.sided", paired = TRUE) Error in t.test.default(fu, bl, alternative = "two.sided", paired = TRUE) : data are essentially constant Тоді …

12 r  t-test