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

12
Какой язык программирования вы рекомендуете для создания прототипа проблемы машинного обучения?

В настоящее время работает в Octave, но из-за плохой документации прогресс очень медленный. Какой язык прост в изучении и использовании и хорошо документирован для решения проблем машинного обучения? Я ищу прототип на небольшом наборе данных (тысячи примеров), поэтому скорость не важна....

12
RFM и моделирование ценности жизни клиента в R

Кто-нибудь может сказать мне, как сделать моделирование недавних событий, частоты и денежной стоимости (RFM) и моделирование потребительской ценности в R? Кроме того, кто-нибудь может направить мне немного литературы об...

12
Наименование средней абсолютной ошибки, аналогичной шкале Бриера?

Вчерашний вопрос « Определить точность модели, которая оценивает вероятность события» , заинтересовал меня оценкой вероятности. Оценка Бриера - это мера среднего квадрата ошибки. Показывает ли аналогичная средняя абсолютная погрешность показатели эффективности есть имя тоже?11NΣя = 1N( Р г е дя с т...

12
Регрессия частичных наименьших квадратов в R: почему PLS на стандартизированных данных не эквивалентно максимизации корреляции?

Я очень новичок в частичных наименьших квадратах (PLS) и пытаюсь понять вывод функции R plsr()в plsпакете. Давайте смоделируем данные и запустим PLS: library(pls) n <- 50 x1 <- rnorm(n); xx1 <- scale(x1) x2 <- rnorm(n); xx2 <- scale(x2) y <- x1 + x2 + rnorm(n,0,0.1); yy <-...

12
Интерпретация результатов Rs ur.df (модульный тест Дикки-Фуллера)

Я выполняю следующий модульный корневой тест (Dickey-Fuller) для временного ряда, используя ur.df()функцию в urcaпакете. Команда: summary(ur.df(d.Aus, type = "drift", 6)) Выход: ############################################### # Augmented Dickey-Fuller Test Unit Root Test #...

12
Как выполнить повторную выборку в R, не повторяя перестановок?

Если в R установить set.seed (), а затем использовать функцию примера для рандомизации списка, могу ли я гарантировать, что не сгенерирую такую ​​же перестановку? то есть ... set.seed(25) limit <- 3 myindex <- seq(0,limit) for (x in seq(1,factorial(limit))) { permutations <-...

12
Как получить область эллипса из двумерных нормальных распределенных данных?

У меня есть данные, которые выглядят так: Я попытался применить нормальное распределение (оценка плотности ядра работает лучше, но мне не нужна такая большая точность), и это работает довольно хорошо. Плотность графика составляет эллипс. Мне нужно получить эту функцию эллипса, чтобы решить,...

12
Оценка распределения по данным

У меня есть образец данных , полученных в Rпути rnorm(50,0,1), поэтому данные , очевидно , берет на себя нормальное распределении. Однако Rне «знает» эту информацию о распределении данных. Существует ли метод R, позволяющий оценить, из какого распределения исходит моя выборка? Если нет, я...

12
Бутстрап, Монте-Карло

Мне задали следующий вопрос в рамках домашней работы: Разработайте и внедрите имитационное исследование, чтобы изучить производительность начальной загрузки для получения 95% доверительных интервалов по среднему значению однофакторной выборки данных. Ваша реализация может быть в R или SAS....

12
На что это указывает, когда корреляция Спирмена на определенную величину меньше Пирсона?

У меня есть куча связанных наборов данных. Корреляции Пирсона между их парами обычно определенно больше, чем корреляции Спирмена. Это говорит о том, что любая корреляция является линейной, но можно ожидать, что даже если бы Пирсон и Спирмен были одинаковыми. Что это означает, когда существует...

12
Как мне подогнать нелинейную модель смешанных эффектов для данных повторных измерений, используя nlmer ()?

Я пытаюсь проанализировать данные повторных измерений и пытаюсь заставить их работать R. Мои данные по существу следующие, у меня есть две группы лечения. Каждый субъект в каждой группе тестируется каждый день и получает оценку (процент правильных на тесте). Данные в длинном формате: Time Percent...

12
Какой смысл в неинформативных априорах?

Почему даже есть неинформативные приоры? Они не предоставляют информацию о . Так зачем их использовать? Почему бы не использовать только информационные приоры? Например, предположим, что θ ∈ [ 0 , 1 ] . Тогда является ли θ ∼ U ( 0 , 1 ) неинформативным априорным для θ ?θθ\thetaθ∈[0,1]θ∈[0,1] \theta...

12
Логистическая квантильная регрессия - как лучше всего передать результаты

В предыдущем посте я задавался вопросом, как справиться с оценками EQ-5D . Недавно я наткнулся на логистическую квантильную регрессию, предложенную Bottai и McKeown, которая представляет элегантный способ справиться с ограниченными результатами. Формула проста: л о гя т ( у) = Л о г( у- ум я нYм а...

12
Различия между PROC Mixed и lme / lmer в R - степени свободы

Примечание: этот вопрос является репостом, так как мой предыдущий вопрос пришлось удалить по юридическим причинам. Сравнивая PROC MIXED из SAS с функцией lmeиз nlmeпакета в R, я наткнулся на некоторые довольно запутанные различия. Более конкретно, степени свободы в разных тестах различаются между...

12
Как генерировать прогнозы с помощью rjags?

Я использовал rjags для запуска MCMC на модели, указанной на языке JAGS. Есть ли хороший способ извлечь эту модель и выполнить с ней предсказания (используя апостериорные распределения моих параметров)? Я могу переопределить модель в R и подключить режимы своих параметров для постеров; Мне просто...

12
Как симулировать функциональные данные?

Я пытаюсь проверить различные подходы анализа функциональных данных. В идеале я хотел бы протестировать панель подходов, которые у меня есть, на смоделированных функциональных данных. Я попытался сгенерировать смоделированный FD, используя подход, основанный на суммировании гауссовских шумов (код...

12
Как я могу оптимизировать вычислительную эффективность при многократной подгонке сложной модели к большому набору данных?

У меня проблемы с производительностью при использовании MCMCglmmпакета в R для запуска модели смешанных эффектов. Код выглядит так: MC1<-MCMCglmm(bull~1,random=~school,data=dt,family="categorical" , prior=list(R=list(V=1,fix=1), G=list(G1=list(V=1, nu=0))) , slice=T, nitt=iter, ,burnin=burn,...

12
Построение кривой вероятности для логит-модели с несколькими предикторами

У меня есть следующая функция вероятности: Prob=11+e−zProb=11+e−z\text{Prob} = \frac{1}{1 + e^{-z}} где z=B0+B1X1+⋯+BnXn.z=B0+B1X1+⋯+BnXn.z = B_0 + B_1X_1 + \dots + B_nX_n. Моя модель выглядит Pr(Y=1)=11+exp(−[−3.92+0.014×(bid)])Pr(Y=1)=11+exp⁡(−[−3.92+0.014×(bid)])\Pr(Y=1) = \frac{1}{1 +...

12
Разница между поисковым и подтверждающим факторным анализом при определении независимости конструкции

Исследователи часто используют две меры, которые имеют очень похожие предметы, и утверждают, что они измеряют разные вещи (например, «я всегда волнуюсь, когда я рядом с машинами»; «я боюсь машин»). Назовем гипотетические меры «Мера страха перед автомобилем» и «Беспокойство от автомобильной шкалы»....