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

Використовуйте цей тег для будь-якого питання * на тему *, який (a) включає `R` як критичну частину запитання або очікувану відповідь, а (b) - не * просто * про те, як використовувати` R`.

2
Тест на коефіцієнт ймовірності в R
Припустимо, я збираюся зробити універсальну логістичну регресію на декількох незалежних змінних, як це: mod.a <- glm(x ~ a, data=z, family=binominal("logistic")) mod.b <- glm(x ~ b, data=z, family=binominal("logistic")) Я зробив порівняння моделі (тест коефіцієнта ймовірності), щоб побачити, чи модель краща за нульову модель за допомогою цієї команди 1-pchisq(mod.a$null.deviance-mod.a$deviance, mod.a$df.null-mod.a$df.residual) Потім я …
25 r  logistic  diagnostic 

2
Виявлення шаблонів обману на іспиті на багато запитань
ПИТАННЯ: У мене є двійкові дані щодо іспитових питань (правильні / неправильні). Деякі люди, можливо, мали попередній доступ до набору питань та їх правильних відповідей. Я не знаю, хто, скільки чи хто. Якби не було обману, припустимо, я б моделював вірогідність правильної відповіді для пунктуiii як, деявляє собою складність питання, …

2
Підмножина векторів часових рядів
У мене є часовий ряд, і я хочу його підмножити, зберігаючи його як часовий ряд, зберігаючи початок, кінець та частоту. Наприклад, скажімо, у мене є часовий ряд: > qs <- ts(101:110, start=c(2009, 2), frequency=4) > qs Qtr1 Qtr2 Qtr3 Qtr4 2009 101 102 103 2010 104 105 106 107 2011 …
25 r  time-series 

3
Щоденний аналіз часових рядів
Я намагаюся зробити аналіз часових рядів і я новачок у цій галузі. Я щодня перераховую подію 2006–2009 рр. І хочу приєднати до неї модель часових рядів. Ось прогрес, який я досяг: timeSeriesObj = ts(x,start=c(2006,1,1),frequency=365.25) plot.ts(timeSeriesObj) Отриманий сюжет я отримую: Щоб перевірити, чи є сезонність і тенденція в даних, чи ні, …

1
Вказання декількох (окремих) випадкових ефектів в lme [закрито]
Зачинено. Це питання поза темою . Наразі відповіді не приймаються. Хочете вдосконалити це питання? Оновіть питання, щоб воно було тематичним для Cross Valified. Закрито 6 місяців тому . Я працював у R-пакетах nlme та lme4 , намагаючись вказати моделі з декількома випадковими ефектами. Я виявив, що лише nlme дозволяє вказати …

3
Інтерпретація термінів взаємодії в логітній регресії з категоричними змінними
У мене є дані опитувального експерименту, в якому респонденти були випадковим чином віднесені до однієї з чотирьох груп: > summary(df$Group) Control Treatment1 Treatment2 Treatment3 59 63 62 66 Хоча три групи лікування дещо відрізняються залежно від застосованого подразника, головне розмежування, яке мене хвилює, - це контрольна та лікувальна групи. Тому …

3
Оцінка логістичної регресії та інтерпретації Хосмера-Лемешоу Goodness of Fit
Як ми всі знаємо, існує 2 методи оцінки логістичної регресійної моделі, і вони тестують дуже різні речі Прогнозова сила: Отримайте статистику, яка вимірює, наскільки добре ви можете передбачити залежну змінну на основі незалежних змінних. Добре відомі псевдо R ^ 2 - Макфадден (1974) і Кокс і Снелл (1989). Статистика придатності …

2
Байєсське ласо проти звичайного ласо
Для програми lasso доступні різні програми для впровадження . Я знаю, що багато обговорювалося про байєсівський підхід і частолістський підхід на різних форумах. Моє запитання дуже специфічне для лассо - Які відмінності чи переваги баснійського ласо проти звичайного ласо ? Ось два приклади реалізації в пакеті: # just example data …

4
Алгоритми виявлення аномалії часових рядів
Зараз я використовую AnomalyDetection Twitter в R: https://github.com/twitter/AnomalyDetection . Цей алгоритм забезпечує виявлення аномалії часових рядів для даних із сезонністю. Питання: чи існують подібні алгоритми, подібні до цього (контроль за сезонністю не має значення)? Я намагаюся оцінити якомога більше алгоритмів часових рядів на моїх даних, щоб я міг вибрати найкращий …

4
Це відповідний метод для перевірки сезонних наслідків даних про кількість самогубств?
У мене 17 років (з 1995 по 2011 рік) даних свідоцтва про смерть, пов’язаних із смертю від самогубства для штату в США. Існує багато міфологій про самогубства та місяці / пори року, багато чого суперечливе, а також про літературу, яку я ' Подивившись, я не розумію використаних методів або впевненості …

1
Як візуалізувати величезну таблицю з надзвичайними ситуаціями?
У мене є дві змінні: назва наркотиків (DN) та відповідні побічні явища (AE), які стоять у співвідношенні «багато до багатьох». Зафіксовано 33 556 найменувань препаратів та 9 556 побічних подій. Розмір вибірки становить близько 5,8 мільйона спостережень. Я хочу вивчити та зрозуміти зв'язок / зв’язок між DN та AE. Я …

2
Як включити термін взаємодії в GAM?
Наступний код оцінює схожість двох часових рядів: set.seed(10) RandData <- rnorm(8760*2) America <- rep(c('NewYork','Miami'),each=8760) Date = seq(from=as.POSIXct("1991-01-01 00:00"), to=as.POSIXct("1991-12-31 23:00"), length=8760) DatNew <- data.frame(Loc = America, Doy = as.numeric(format(Date,format = "%j")), Tod = as.numeric(format(Date,format = "%H")), Temp = RandData, DecTime = rep(seq(1, length(RandData)/2) / (length(RandData)/2), 2)) require(mgcv) mod1 <- gam(Temp …

3
Як зробити логістичну регресію в R, коли результат є дробовим (співвідношення двох рахунків)?
Я переглядаю статтю з наступним біологічним експериментом. Пристрій використовується для піддавання клітинам різної кількості напруги зсуву рідини. Оскільки сильні напруги зсуву прикладаються до клітин, більше їх починає відшаровуватися від субстрату. На кожному рівні напруги зсуву вони підраховують клітини, які залишаються приєднаними, і оскільки вони знають загальну кількість комірок, які були …

2
Чому lme та aov повертають різні результати для повторних заходів ANOVA в R?
Я намагаюся перейти від використання ezпакету до lmeповторних заходів ANOVA (оскільки, я сподіваюся, я зможу використовувати власні контрасти на с lme). Дотримуючись поради з цього допису в блозі, я зміг налаштувати ту саму модель, використовуючи як aov(так і ezза запитом) та lme. Однак, хоча у прикладі, наведеному в цій публікації, …

3
Яка підходяща модель для низькодисперсних даних підрахунку?
Я намагаюся моделювати дані підрахунку в R, які, мабуть, нерозподілені (параметр дисперсії ~ .40). Це, мабуть, причина glmз family = poissonабо негативною біноміальною ( glm.nb) моделлю не є істотною. Коли я дивлюсь на описи моїх даних, у мене немає типового перекосу даних підрахунку, і залишки в двох моїх експериментальних умовах …

Використовуючи наш веб-сайт, ви визнаєте, що прочитали та зрозуміли наші Політику щодо файлів cookie та Політику конфіденційності.
Licensed under cc by-sa 3.0 with attribution required.