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

13
Почему R's lm () возвращает оценки коэффициентов, отличные от моего учебника?

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

13
Интеграция эмпирического CDF

У меня есть эмпирическое распределение . Я рассчитываю это следующим образомG(x)G(x)G(x) x <- seq(0, 1000, 0.1) g <- ecdf(var1) G <- g(x) Я обозначаю , т. - это pdf, а - это cdf.h Gh(x)=dG/dxh(x)=dG/dxh(x) = dG/dxhhhGGG Теперь я хочу решить уравнение для верхнего предела интегрирования...

13
Как проверить, следует ли распределение степенному закону?

У меня есть данные о том, сколько пользователей публикуют сколько вопросов. Например, [UserCount, QuestionCount] [2, 100] [9, 10] [3, 80] ... ... Это означает, что 2 пользователя разместили по 100 вопросов, 9 пользователей - по 10 вопросов и т. Д. Итак, как я могу определить, UserCount,...

13
С помощью пакета каретки можно ли получить матрицы путаницы для конкретных пороговых значений?

Я получил модель логистической регрессии (через train) для бинарного ответа, и я получил логистическую матрицу спутанности через confusionMatrixв caret. Это дает мне путаницу в логистической модели, хотя я не уверен, какой порог используется для ее получения. Как получить матрицу путаницы для...

13
Расчет

Я читал о расчете значений в смешанных моделях и после прочтения FAQ по R-sig, других постов на этом форуме (я бы связал несколько, но мне не хватает репутации) и нескольких других ссылок, которые я понимаю, используя Значения в контексте смешанных моделей сложны.R 2р2R2R^2р2R2R^2 Однако недавно я...

13
Пример противоречивой оценки максимального правдоподобия

Я читаю комментарий к статье, и автор заявляет, что иногда, даже если оценки (найденные по ОД или максимальному квазиликтуализму) могут быть непоследовательными, сила отношения правдоподобия или теста отношения квази-правдоподобия все еще может сходиться к 1, поскольку число наблюдаемых данных...

13
Неверный вывод, когда наблюдения не являются независимыми

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

13
R: проверить нормальность остатков линейной модели - какие остатки использовать

Я хотел бы сделать W-тест Шапиро Уилка и тест Колмогорова-Смирнова на невязках линейной модели, чтобы проверить на нормальность. Мне было просто интересно, какие остатки следует использовать для этого - необработанные остатки, остатки Пирсона, студентизированные остатки или стандартизированные...

13
Извлечение откосов для случаев из модели смешанных эффектов (lme4)

Я хотел бы извлечь наклоны для каждого человека в модели смешанного эффекта, как обрисовано в общих чертах в следующем параграфе Модели смешанных эффектов использовались для характеристики отдельных путей изменения показателей когнитивной сводки, включая термины для возраста, пола и года обучения в...

13
Тест на коинтеграцию между двумя временными рядами с использованием двухэтапного метода Энгла – Грейнджера

Я пытаюсь проверить коинтеграцию между двумя временными рядами. Обе серии имеют еженедельные данные, охватывающие ~ 3 года. Я пытаюсь сделать двухэтапный метод Энгла-Грейнджера. Мой порядок действий следующий. Протестируйте каждый временной ряд для корневого модуля с помощью Augmented...

13
Аддитивная ошибка или мультипликативная ошибка?

Я относительно новичок в статистике и был бы признателен за помощь в понимании этого вопроса. В моей области есть широко используемая модель вида: пT= Pо( VT)αпTзнак равнопо(ВT)αP_t = P_o(V_t)^\alpha Когда люди подгоняют модель к данным, они обычно линеаризуют ее и соответствуют следующим журнал(...

13
Процедура и методы анализа временных рядов с использованием R

Я работаю над небольшим проектом, в котором мы пытаемся прогнозировать цены на товары (нефть, алюминий, олово и т. Д.) На ближайшие 6 месяцев. У меня есть 12 таких переменных для прогнозирования, и у меня есть данные за апрель 2008 года - май 2013 года. Как я должен идти о предсказании? Я сделал...

13
Смещенная оценка для регрессии, достигающая лучших результатов, чем объективная оценка в модели Error In Variables

Я работаю над некоторыми синтетическими данными для модели Error In Variable для некоторых исследований. В настоящее время у меня есть одна независимая переменная, и я предполагаю, что знаю дисперсию для истинного значения зависимой переменной. Таким образом, с помощью этой информации я могу...

13
Значение оси Y на графике частичной зависимости Random Forest

Я использую RandomForestпакет R и не понимаю, как интерпретировать значения оси Y на графиках их частичной зависимости. Справочные документы утверждают, что график представляет собой «графическое изображение предельного влияния переменной на вероятность класса». Тем не менее, я все еще не понимаю,...

13
Какова интерпретация ковариации коэффициентов регрессии?

Функция lm в R может выводить оценочную ковариацию коэффициентов регрессии. Что эта информация дает нам? Можем ли мы теперь лучше интерпретировать модель или диагностировать проблемы, которые могут присутствовать в...

13
Интерпретация случайного отклонения эффекта в блеске

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

13
Настройка гиперпараметра в регрессии Гаусса

Я пытаюсь настроить гиперпараметры алгоритма гауссовой регрессии, который я реализовал. Я просто хочу максимизировать предельное правдоподобие, определяемое формулой где K - ковариационная матрица с элементы K_ {ij} = k (x_i, x_j) = b ^ {- 1} \ exp (- \ frac {1} {2} (x_i-x_j) ^ TM (x_i-x_j)) + a ^...

13
Формула для 95% доверительного интервала для

Я гуглил и искал по stats.stackexchange, но не могу найти формулу для расчета 95% доверительного интервала для значения для линейной регрессии. Кто-нибудь может это предоставить?р2R2R^2 Еще лучше, скажем, я выполнил линейную регрессию ниже в R. Как бы я вычислил 95% доверительный интервал для...

13
Какова цель нормализации строк

Я понимаю причину нормализации столбцов, поскольку она приводит к одинаковому взвешиванию объектов, даже если они не измеряются в одном и том же масштабе - однако часто в литературе ближайшего соседа столбцы и строки нормализуются. Что такое нормализация строк для / почему нормализация строк? В...

13
Объясните, как `eigen` помогает инвертировать матрицу

Мой вопрос относится к вычислительной технике, используемой в geoR:::.negloglik.GRFили geoR:::solve.geoR. В линейной смешанной модели: Y=Xβ+Zb+eY=Xβ+Zb+e Y=X\beta+Zb+e где ββ\beta и bbb - фиксированные и случайные эффекты соответственно. Также Σ=cov(Y)Σ=cov(Y)\Sigma=\text{cov}(Y) При оценке...