Предположим , у меня есть «экспертов», из которых я хотел бы, чтобы вызвать предварительное распределение по некоторой переменной X . Я хотел бы мотивировать их реальными деньгами . Идея состоит в том, чтобы выявить априоры, наблюдать n реализаций случайной величины X , а затем разделить некоторый...
Я выполнил трехсторонние повторные измерения ANOVA; какие специальные анализы действительны? Это полностью сбалансированный дизайн (2x2x2) с одним из факторов, имеющих повторное измерение внутри субъекта. Мне известны многовариантные подходы к повторным измерениям ANOVA в R, но мой первый инстинкт...
Существует много способов выполнения регуляризации - например, регуляризация на основе норм , и . Согласно Friedman Hastie & Tibsharani , лучший регуляризатор зависит от проблемы: а именно от природы истинной целевой функции, конкретной используемой основы, отношения сигнал / шум и размера...
Если я подгоняю свои данные к чему-то вроде lm(y~a*b), в синтаксисе R, где aэто двоичная переменная и bчисловая переменная, то a:bтермин взаимодействия - это разница между наклоном y~bat a= 0 и at a= 1. Теперь, скажем, отношения между yи bкриволинейные. Если я сейчас подхожу lm(y~a*poly(b,2)), то...
Коэффициент отношения правдоподобия (он же отклонение) и критерий несоответствия (или качества соответствия) довольно просто получить для модели логистической регрессии (подгонка с использованием функции) в R. Однако это может быть легко подсчитать количество клеток в конечном итоге достаточно...
Для линейной регрессии с несколькими группами (естественные группы определены априори) допустимо ли запускать две разные модели на одном наборе данных, чтобы ответить на следующие два вопроса? Есть ли в каждой группе ненулевой наклон и ненулевой перехват, и каковы параметры для каждой из них в...
Я пытаюсь определить, будут ли простые вероятности работать для моей проблемы или будет лучше использовать (и узнать о) более сложные методы, такие как логистическая регрессия. Переменная ответа в этой задаче является двоичным ответом (0, 1). У меня есть несколько переменных предикторов, которые...
Я пытаюсь повторить работу коллеги и перемещаю анализ из Stata в R. Модели, которые она использует, вызывают параметр «cluster» в функции nbreg для кластеризации стандартных ошибок. См. Http://repec.org/usug2007/crse.pdf для довольно полного описания того, что и почему этого параметра Мой вопрос:...
Я пытаюсь проанализировать отставание между временными рядами двух цен акций. В обычном анализе временных рядов мы можем выполнить Cross Correlaton, VECM (Причинность Грейнджера). Однако, как можно обращаться с тем же самым в нерегулярно расположенных временных рядах. Гипотеза состоит в том, что...
Я хочу смоделировать две разные переменные времени, некоторые из которых сильно коллинеарны в моих данных (возраст + группа = период). Делая это, я столкнулся с некоторыми проблемами lmerи взаимодействиями poly(), но это, вероятно, не ограничивалось lmer, я получил те же результаты с nlmeIIRC....
Изображение ниже из этой статьи в психологической науке . Коллега указал на две необычные вещи об этом: Согласно подписи, столбцы ошибок показывают «± 2,04 стандартных ошибок, 95% доверительный интервал». Я только видел, как ± 1,96 SE использовался для 95% -ного КИ, и я не могу найти ничего о том,...
Я делаю односторонний ANOVA (для каждого вида) с индивидуальными контрастами. [,1] [,2] [,3] [,4] 0.5 -1 0 0 0 5 1 -1 0 0 12.5 0 1 -1 0 25 0 0 1 -1 50 0 0 0 1 где я сравниваю интенсивность 0,5 против 5, 5 против 12,5 и так далее. Это данные, над которыми я работаю со следующими результатами...
Я заметил, что когда я пытаюсь найти критические значения для Манна-Уитни U, используя R, значения всегда 1 + критическое значение. Например, для критическое значение (двусторонний) равно 8, а для α = 0,05 , n = 12 , m = 8 критическое (двустороннее) значение равно 22 (проверьте таблицы ), но:α =...
В следующем примере у меня есть кадр данных, который состоит из временного ряда измерений температуры воды, зарегистрированных на 5 глубинах в океане, где каждое значение Tempсоответствует дате в DateTimeи глубине в Depth. set.seed(1) Temp <- rnorm(43800,sd=20) AirT <- rnorm(8760,sd=20) Depth...
Модели смещенных монет обычно имеют один параметр . Один из способов оценить θ по серии тиражей - использовать предварительное бета-тестирование и вычислить апостериорное распределение с биномиальной вероятностью.θ = P( Руководитель | θ )θ=P(Head|θ)\theta = P(\text{Head} | \theta)θθ\theta В моих...
Я снабдил логнормальную модель, используя R набором данных. Полученные параметры были: meanlog = 4.2991610 sdlog = 0.5511349 Я бы хотел перенести эту модель на Scipy, которой никогда раньше не пользовался. Используя Scipy, я смог получить форму и масштаб 1 и 3.1626716539637488e + 90 - очень разные...
Структура данных > str(data) 'data.frame': 6138 obs. of 10 variables: $ RT : int 484 391 422 516 563 531 406 500 516 578 ... $ ASCORE : num 5.1 4 3.8 2.6 2.7 6.5 4.9 2.9 2.6 7.2 ... $ HSCORE : num 6 2.1 7.9 1 6.9 8.9 8.2 3.6 1.7 8.6 ... $ MVMNT : Factor w/ 2 levels "_Withd","Appr": 2 2 1 1 2 1 2...
Я не получаю разницу между rfobject$importanceи importance(rfobject)в столбце MeanDecreaseAccuracy. Пример: > data("iris") > fit <- randomForest(Species~., data=iris, importance=TRUE) > fit$importance setosa versicolor virginica MeanDecreaseAccuracy MeanDecreaseGini Sepal.Length...
Мои вопросы касаются GAMs в пакете mgcv R. Из-за небольшого размера выборки я хочу определить ошибку прогнозирования, используя перекрестную проверку с пропуском. Это разумно? Есть ли пакет или код, как я могу это сделать? errorest()Функция в ipred пакете не работает. Простой тестовый набор данных:...
У меня есть эксперимент, который выполняется на сотнях компьютеров, распределенных по всему миру, который измеряет случаи определенных событий. Каждое событие зависит друг от друга, поэтому я могу расположить их в порядке возрастания, а затем рассчитать разницу во времени. События должны быть...