У меня есть следующие простые векторы X и Y: > X [1] 1.000 0.063 0.031 0.012 0.005 0.000 > Y [1] 1.000 1.000 1.000 0.961 0.884 0.000 > > plot(X,Y) Я хочу сделать регрессию с использованием журнала X. Чтобы избежать получения журнала (0), я стараюсь положить +1 или +0,1 или +0,00001 или...
Функция R cv.glm (library: boot) вычисляет предполагаемую K-кратную ошибку прогнозирования перекрестной проверки для обобщенных линейных моделей и возвращает дельту. Имеет ли смысл использовать эту функцию для регрессии лассо (library: glmnet) и, если да, то как ее можно выполнить? Библиотека...
Я хочу сравнить 2 типовых средства для 1-минутного возврата. Я предполагаю, что они распределены по Лапласу (уже проверены), и я разделил результаты на 2 группы. Как я могу проверить, значительно ли они отличаются? Я думаю, что не могу относиться к ним как к нормальному распределению, потому что,...
В lmerфункции lme4in in Rесть вызов для построения модельной матрицы случайных эффектов , как описано здесь , стр. 7 - 9.ZZZ Вычисление влечет за собой произведения ХатриРао и / или Кронекера двух матриц, и . ZZZJiJiJ_iXiXiX_i Матрица представляет собой глоток: «Матрица индикаторов группирующих...
Согласно документации removeSparseTermsфункции из tmпакета, это то, что влечет за собой редкость: A term-document matrix where those terms from x are removed which have at least a sparse percentage of empty (i.e., terms occurring 0 times in a document) elements. I.e., the resulting matrix contains...
Каковы преимущества выражения модели ARMA как модели пространства состояний и прогнозирования с использованием фильтра Калмана? Эта методология, например, используется в реализации SARIMAX для python-statsmodels:...
Существует ли пакет R (или даже базовая функция R), который реализует метод Фишера или Стоуффера для объединения p-значений? Кодирование этого должно быть почти тривиальным, но я бы предпочел использовать (и процитировать) пакет. Пример кода в этом вопросе: метод Фишера для комбинирования...
Я хочу оценить среднее значение функции f, т. где и - независимые случайные величины. У меня есть образцы f, но не iid: есть образцы для и для каждого есть выборки из :X Y Y 1 , Y 2 , … Y n Y i n i X X i , 1 , X i , 2 , … , X i , n...
Если таковые имеются, между подгонкой линии к нескольким отдельным «экспериментам», затем усреднением подгонки или усреднением данных из отдельных экспериментов, затем подгонкой усредненных данных. Позвольте мне уточнить: Я выполняю компьютерное моделирование, которое генерирует кривую, показанную...
Задачи классификации с нелинейными границами не могут быть решены простым персептроном . Следующий код R предназначен для иллюстративных целей и основан на этом примере в Python): nonlin <- function(x, deriv = F) { if (deriv) x*(1-x) else 1/(1+exp(-x)) } X <- matrix(c(-3,1, -2,1, -1,1, 0,1,...
Как я могу получить p,d and qзначения в ARIMA(p,d,q)модели по оценкам auto.arima(mytimeseries)? arima_model <- auto.arima (mytimeseries, ic = 'bic') Если мы посмотрим на вывод arima_model $ арма мы получили, [1] 1 0 0 0 1 2 0 Какое значение имеют цифры в приведенной выше...
контекст Я пытаюсь понять, как R coxph () принимает и обрабатывает повторяющиеся записи для субъектов (или пациента / клиента, если вы предпочитаете). Некоторые называют это длинным форматом, другие называют это «повторными измерениями». Посмотрите, например, набор данных, который включает столбец...
Доказать / Не подтвердить E[1A|Ft]=0 or 1 a.s. ⇒E[1A|Fs]=E[1A|Ft] a.s.E[1A|Ft]=0 or 1 a.s. ⇒E[1A|Fs]=E[1A|Ft] a.s.E[1_A | \mathscr{F_t}] = 0 \ \text{or} \ 1 \ \text{a.s.} \ \Rightarrow E[1_A | \mathscr{F_{s}}] = E[1_A | \mathscr{F_t}] \ \text{a.s.} Учитывая , отфильтрованный вероятностное...
Я занимаюсь исследованиями в области образовательных игр, и некоторые из моих текущих проектов включают использование данных из BoardGameGeek (BGG) и VideoGameGeek (VGG) для изучения взаимосвязей между элементами дизайна игр (т. Е. «Набор во Второй мировой войне», «включает в себя бросание...
Загрузите пакет, необходимый. library(ggplot2) library(MASS) Генерация 10000 номеров, приспособленных к гамма-распределению. x <- round(rgamma(100000,shape = 2,rate = 0.2),1) x <- x[which(x>0)] Нарисуйте функцию плотности вероятности, предположим, что мы не знаем, к какому распределению x...
При выполнении байесовской линейной регрессии необходимо назначить априор для наклона и точки пересечения . Поскольку является параметром местоположения, имеет смысл назначить унифицированный априор; однако, мне кажется, что сродни параметру масштаба, и кажется неестественным назначать униформу до...
Как бы вы назвали этот тип сюжета, и возможно ли его создать в R? РЕДАКТИРОВАТЬ: большое спасибо всем - очень полезно. Лучший титул на данный момент: квантовые сюжеты для скрипки!...
Ситуация: две птицы (самец и самка) защищают свои яйца в гнезде от злоумышленника. Каждая птица может использовать для защиты либо атаку, либо угрозу, а также присутствовать или отсутствовать. Из данных возникает закономерность, что поведение может быть взаимодополняющим - атаки мужчин, в то время...
Являются ли «изученные остатки» и «стандартизированные остатки» одинаковыми в регрессионных моделях? Я построил модель линейной регрессии в R и хотел построить график Studentized Остатки v / s, подобранные значения, но не нашел автоматизированный способ сделать это в R. Предположим, у меня есть...
Я попытался реализовать численную оценку дивергенции Кульбака-Лейблера для двух выборок. Для отладки реализации возьмем образцы из двух нормальных распределений N(0,1)N(0,1)\mathcal N (0,1) и N(1,2)N(1,2)\mathcal N (1,2) . Для простой оценки я сгенерировал две гистограммы и попытался численно...