Мне нужно сделать вывод о положительном параметре . Чтобы подчеркнуть положительность, я репараметризовал . Используя процедуру MLE, я вычислил точечную оценку и нашел для . Свойство инвариантности MLE напрямую дает мне точечную оценку для p , но я не уверен, как вычислить se для p . Заранее...
Я запустил ANOVA, обнаружив, например, взаимодействие между полом и классом, по которому я хочу узнать, в каких классах мальчики и девочки отличаются, но во многих случаях я нахожу (скорректированные) значения p 0 и 1. Как / почему это возможно? Не кажется правильным ... as.factor(gender) 1 16 16.2...
Locked . Этот вопрос и его ответы заблокированы, потому что вопрос не по теме, но имеет историческое значение. В настоящее время он не принимает новые ответы или взаимодействия. Функции, с которыми я знаком, включают масштабирование из базы R, масштабирование из ARM. Возможно, лучшим способом было...
Насколько я понимаю, мы не должны позволять тому же набору данных, который мы анализируем, управлять / определять, как выглядят предыдущие распределения в байесовском анализе. В частности, неуместно определять предыдущие распределения для байесовского анализа, основанные на сводной статистике из...
Проверяя, были ли две переменные коррелированными, я заметил, что применение корреляции Пирсона позволило получить числа, равные 0,1, что указывает на отсутствие корреляции. Что я могу сделать, чтобы усилить это утверждение? Набор данных (подмножество из-за ограничений публикации), на который я...
Я изучаю книгу «Вводные временные ряды с R» Каупертвейта и Меткалфа. На странице 36 написано, что строки: . Я прочитал здесь R форум, что строки в . - 1 / n ± 2 / n--√-1/N±2/N-1/n \pm 2/\sqrt{n}± 1,96 / н--√±1,96/N\pm 1.96/\sqrt{n} Я запустил следующий код: b = c(3,1,4,1) acf(b) и я вижу, что...
Википедия предполагает, что одним из способов оценки надежности между оценками является использование модели случайных эффектов для вычисления внутриклассовой корреляции . Пример внутриклассовой корреляции говорит о взгляде на...
У меня есть вектор x <- c(1,2,3,4,5,5,5,6,6,6,6, 7,7,7,7,7,7,7,7,7,7,7,7,7,7,7,7,7,7,7, 7,7,7,7,7,7,7,7,8,8,8,8,9,9,9,10) (мой фактический вектор имеет длину> 10000), и я хотел бы найти интервалы, в которых лежит 90% плотности. Является ли quantile(x, probs=c(0.05,0.95), type=5)наиболее...
Я пытаюсь визуализировать некоторые потребительские данные, которые имеют 4 категории. Пользователи могут свободно переключаться между различными категориями. Я хотел бы визуализировать последние три или четыре переключателя для каждого человека. Таким образом, мы начнем с графика с колонкой с 4...
Если у меня есть переменная с 4 уровнями, теоретически мне нужно использовать 3 фиктивные переменные. На практике, как это на самом деле осуществляется? Я использую 0-3, я использую 1-3 и оставляю 4 пустыми? Какие-либо предложения? ПРИМЕЧАНИЕ: я собираюсь работать в R. ОБНОВЛЕНИЕ: Что случилось бы,...
У меня есть данные временного ряда, где измеряемая переменная представляет собой дискретные положительные целые числа (числа). Я хочу проверить, есть ли тенденция со временем (или нет). Независимая переменная (x) находится в диапазоне 0-500, а зависимая переменная (y) находится в диапазоне 0-8. Я...
Закрыто. Этот вопрос не по теме . В настоящее время он не принимает ответы. Хотите улучшить этот вопрос? Обновите вопрос, чтобы он соответствовал теме перекрестной проверки. Закрыто 9 месяцев назад . Проблема, упомянутая в этом вопросе, исправлена в версии 1.7.3 пакета R glmnet. У меня возникли...
Я создал свою немного улучшенную версию termplot, которую я использую в этом примере, вы можете найти ее здесь . Ранее я писал о SO, но чем больше я об этом думаю, тем больше думаю, что это скорее связано с интерпретацией модели пропорциональных рисков Кокса, чем с фактическим кодированием....
Я ищу простой пример кода для того, как запустить Particle Filter в R. Пакет pomp, кажется, поддерживает математический бит пространства состояний, но примеры являются немного сложными, чтобы следовать программно для простого разработчика OO, такого как я, особенно как загрузить наблюдаемые данные...
Существует ли общий принцип о том, следует ли вычислять корреляцию Пирсона для двух случайных величин X и Y перед выполнением их лог-преобразования или после? Есть ли процедура для проверки, которая более подходит? Они дают одинаковые, но разные значения, поскольку логарифмическое преобразование...
Закрыто. Этот вопрос не по теме . В настоящее время он не принимает ответы. Хотите улучшить этот вопрос? Обновите вопрос, чтобы он соответствовал теме перекрестной проверки. Закрыто 3 года назад . Мне очень интересно узнать, как создавать повторяющиеся отчеты из моего кода R и визуализации...
В качестве дополнительного хобби я изучал прогнозирование временных рядов (в частности, с использованием R). По моим данным, у меня есть количество посещений в день, за каждый день, уходящий почти на 4 года. В этих данных есть несколько четких закономерностей: Понедельник-пятница имеет много...
Скажем, у вас есть данные о выживании, как это: obs <- data.frame( time = c(floor(runif(100) * 30), floor((runif(100)^2) * 30)), status = c(rbinom(100, 1, 0.2), rbinom(100, 1, 0.7)), group = gl(2,100) ) Чтобы выполнить стандартный тест ранга журнала, можно использовать survdiff(Surv(time,...
Я использую пакет kernlab в R для создания SVM для классификации некоторых данных. SVM работает хорошо, поскольку он обеспечивает «предсказания» с приличной точностью, однако мой список входных переменных больше, чем мне бы хотелось, и я не уверен относительно относительной важности различных...
У меня есть проблемы в использовании cor()и cor.test()функции. У меня просто есть две матрицы (только числовые значения и одинаковое количество строк и столбцов), и я хочу получить номер корреляции и соответствующее значение p. Когда я использую, cor(matrix1, matrix2)я получаю коэффициенты...