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

Я просматриваю видео в бесплатном онлайн-курсе Эндрю Нг по машинному обучению в Стэнфорде. Он рассматривает Gradient Descent как алгоритм для решения функций линейной регрессии и записи в Octave для его выполнения. Предположительно я мог бы переписать эти функции в R, но мой вопрос в том, разве функция lm () уже …

14 r  regression  machine-learning  gradient-descent 

Может кто-нибудь сказать мне, как заставить R оценить точку разрыва в кусочно-линейной модели (как фиксированный или случайный параметр), когда мне также нужно оценить другие случайные эффекты? Ниже я привел пример с игрушкой, который соответствует регрессии хоккейной клюшки / сломанной палки со случайными отклонениями наклона и случайной дисперсией y-точки пересечения для …

14 r  mixed-model  lme4-nlme  change-point  piecewise-linear 

Закрыто. Этот вопрос не по теме . В настоящее время он не принимает ответы. Хотите улучшить этот вопрос? Обновите вопрос, чтобы он соответствовал теме перекрестной проверки. Закрыто 7 месяцев назад . Легко найти область вычисления пакета в ROC, но есть ли пакет, который вычисляет площадь под кривой точного возврата?

14 r  precision-recall 

Я собираюсь погрузиться в изучение R, и мой учебный проект повлечет за собой применение регрессии со смешанными или случайными эффектами к набору данных для разработки прогностического уравнения. Я разделяю озабоченность автора в этом посте Как выбрать библиотеку nlme или lme4 R для моделей со смешанными эффектами? задаваясь вопросом, является ли …

14 r  regression  random-effects-model 

У меня есть набор данных, который содержит, скажем, некоторые измерения для положения, скорости и ускорения. Все приходят от одного и того же «бега». Я мог бы построить линейную систему и подогнать полином для всех этих измерений. Но могу ли я сделать то же самое со сплайнами? Что такое способ «R»? …

14 r  model  fitting  splines  derivative 

Стандартной ошибкой является оценочное стандартное отклонение оценки для параметра . & thetas ; & thetasσ^( θ^)σ^(θ^)\hat \sigma(\hat\theta)θ^θ^\hat\thetaθθ\theta Почему расчетное стандартное отклонение от остатков называется «остаточной стандартной ошибкой» (например, при выводе функции R summary.lm), а не «остаточное стандартное отклонение»? Какую оценку параметров мы снабжаем здесь стандартной ошибкой? Рассматриваем ли мы каждый …

14 r  standard-error  residuals  terminology 

Я хотел бы провести прогноз на основе ARIMA-модели с несколькими временными рядами с несколькими экзогенными переменными. Поскольку я не настолько опытен в отношении статистики, которую ни RI хотят сохранить, это настолько просто, насколько это возможно (прогноз тренда на 3 месяца достаточно). У меня есть 1 зависимый временной ряд и 3-5 …

14 r  time-series  arima 

У меня есть некоторые данные, которые я использую loess. Я хотел бы найти точки перегиба сглаженной линии. Это возможно? Я уверен, что кто-то нашел причудливый метод, чтобы решить эту проблему ... Я имею в виду ... в конце концов, это R! Я в порядке с изменением функции сглаживания, которую я …

14 r  smoothing  loess 

Я пытаюсь понять, как вычислить оптимальную точку отсечения для кривой ROC (значение, при котором чувствительность и специфичность максимальны). Я использую набор данных aSAHиз пакета pROC. outcomeПеременная может быть объяснено двумя независимыми переменными: s100bи ndka. Используя синтаксис Epiпакета, я создал две модели: library(pROC) library(Epi) ROC(form=outcome~s100b, data=aSAH) ROC(form=outcome~ndka, data=aSAH) Вывод иллюстрируется на …

14 r  roc  sensitivity-analysis  sensitivity-specificity 

Я новичок в языке R. Я хотел бы знать, как имитировать модель множественной линейной регрессии, которая удовлетворяет всем четырем предположениям регрессии. хорошо спасибо. Допустим, я хочу смоделировать данные на основе этого набора данных: y<-c(18.73,14.52,17.43,14.54,13.44,24.39,13.34,22.71,12.68,19.32,30.16,27.09,25.40,26.05,33.49,35.62,26.07,36.78,34.95,43.67) x1<-c(610,950,720,840,980,530,680,540,890,730,670,770,880,1000,760,590,910,650,810,500) x2<-c(1,1,3,2,1,1,3,3,2,2,1,3,3,2,2,2,3,3,1,2) fit<-lm(y~x1+x2) summary(fit) тогда я получаю вывод: Call: lm(formula = y ~ x1 + x2) …

14 r  multiple-regression  simulation 

В настоящее время у меня проблема с пониманием синтаксиса R для подгонки GLM с использованием гамма-распределения. У меня есть набор данных, где каждая строка содержит 3 ко-вариации ( ), переменную ответа ( ) и параметр формы ( K ). Я хочу смоделировать масштаб гамма-распределения как линейную функцию от трех ковариат, …

14 r  generalized-linear-model  gamma-distribution  dglm 

Мне нужно объединить единицы в кластеров, чтобы минимизировать сумму квадратов внутри группы (WSS), но мне нужно убедиться, что каждый из кластеров содержит не менее единиц. Любая идея, если какая-либо из функций кластеризации R позволяет кластеризовать в кластеров с учетом ограничения минимального размера кластера? Кажется, kmeans () не предлагает опцию ограничения …

14 r  clustering 

Я пытался понять, что именно означает "coef" и "(exp) coef" coxph. Кажется, что "(exp) coef" - сравнение первой переменной в модели в соответствии с группой, назначенной в команде. Как функция coxph получает значения для "coef" и "(exp) coef"? Кроме того, как coxph определяет эти значения, когда проводится цензура?

14 r  survival  interpretation 

У меня есть вектор чисел, который я загрузил здесь (... / code / MyData.Rdata), используя dput. Я хотел бы получить bca ci, поэтому я написал этот код: my.mean <- function(dat, idx){ return (mean(dat[idx], na.rm = TRUE)) } boot.out<-boot(data=my.data, statistic = my.mean, R=1000) Но когда я запускаю следующее, я получаю это: …

14 r  bootstrap 

Основываясь на выборке из времен выживания, я хотел бы оценить вероятность времени выживания для некоторого конкретного , используя оценку Каплана-Мейера. Возможно ли сделать это в ? Пожалуйста, обратите внимание, что не обязательно время события.nnnttttttRttt

14 r  kaplan-meier