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

Мои данные - это временной ряд занятого населения, L, и временной интервал, год. n.auto=auto.arima(log(L),xreg=year) summary(n.auto) Series: log(L) ARIMA(2,0,2) with non-zero mean Coefficients: ar1 ar2 ma1 ma2 intercept year 1.9122 -0.9567 -0.3082 0.0254 -3.5904 0.0074 s.e. NaN NaN NaN NaN 1.6058 0.0008 sigma^2 estimated as 1.503e-06: log likelihood=107.55 AIC=-201.1 AICc=-192.49 BIC=-193.79 …

9 r  regression  arima 

Я знаю, что это, вероятно, основной вопрос ... Но я, кажется, не могу найти ответ. Я подгоняю GLM к семейству Пуассонов, а затем попытался взглянуть на прогнозы, однако смещение, похоже, принимается во внимание: model_glm=glm(cases~rhs(data$year,2003)+lhs(data$year,2003), offset=(log(population)), data=data, subset=28:36, family=poisson()) predict (model_glm, type="response") Я получаю случаи, а не ставки ... Я пробовал …

9 r  generalized-linear-model  prediction  offset 

@whuber продемонстрировал, как моделировать многовариантные результаты ( , и ) для одной временной точки.у 2 у 3Y1y1y_1Y2y2y_2Y3y3y_3 Как известно, продольные данные часто встречаются в медицинских исследованиях. Мой вопрос заключается в том, как имитировать повторные измерения многовариантных результатов в R? Например, мы неоднократно измеряем , и в 5 различных временных точках …

9 r  repeated-measures  simulation 

Мне нужен пакет, который может дать мне уравнение для линейной модели SVM. В настоящее время я использую E1071 так: library(e1071) m = svm(data, labels, type='C', kernel='linear', cost=cost, probability=FALSE, scale=scale) w = t(m$coefs) %*% data[m$index,] #Weight vector b = -model$rho #Offset Тем не менее, я не уверен, как e1071::svm()выбирать положительные и …

9 r  svm  e1071 

У меня есть проблемы в использовании cor()и cor.test()функции. У меня просто есть две матрицы (только числовые значения и одинаковое количество строк и столбцов), и я хочу получить номер корреляции и соответствующее значение p. Когда я использую, cor(matrix1, matrix2)я получаю коэффициенты корреляции для всех ячеек. Я просто хочу одно число как …

9 r  correlation 

Я использую пакет kernlab в R для создания SVM для классификации некоторых данных. SVM работает хорошо, поскольку он обеспечивает «предсказания» с приличной точностью, однако мой список входных переменных больше, чем мне бы хотелось, и я не уверен относительно относительной важности различных переменных. Я хотел бы реализовать генетический алгоритм, чтобы выбрать …

9 r  machine-learning  svm  genetic-algorithms 

Скажем, у вас есть данные о выживании, как это: 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, status) ~ group, data = obs, rho = 0) …

9 r  survival  wilcoxon-signed-rank 

В качестве дополнительного хобби я изучал прогнозирование временных рядов (в частности, с использованием R). По моим данным, у меня есть количество посещений в день, за каждый день, уходящий почти на 4 года. В этих данных есть несколько четких закономерностей: Понедельник-пятница имеет много посещений (самый высокий в понедельник / вторник), но …

9 r  time-series  seasonality  forecasting 

Закрыто. Этот вопрос не по теме . В настоящее время он не принимает ответы. Хотите улучшить этот вопрос? Обновите вопрос, чтобы он соответствовал теме перекрестной проверки. Закрыто 3 года назад . Мне очень интересно узнать, как создавать повторяющиеся отчеты из моего кода R и визуализации ggplot2. Я понимаю, что LaTeX …

9 r 

Я ищу простой пример кода для того, как запустить Particle Filter в R. Пакет pomp, кажется, поддерживает математический бит пространства состояний, но примеры являются немного сложными, чтобы следовать программно для простого разработчика OO, такого как я, особенно как загрузить наблюдаемые данные в объект помпы. Примеры здесь: http://cran.r-project.org/web/packages/pomp/vignettes/intro_to_pomp.pdf Допустим, у меня …

9 r 

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

9 r  survival  cox-model 

Закрыто. Этот вопрос не по теме . В настоящее время он не принимает ответы. Хотите улучшить этот вопрос? Обновите вопрос, чтобы он соответствовал теме перекрестной проверки. Закрыто 9 месяцев назад . Проблема, упомянутая в этом вопросе, исправлена ​​в версии 1.7.3 пакета R glmnet. У меня возникли некоторые проблемы при запуске …

9 r  multinomial  glmnet 

У меня есть данные временного ряда, где измеряемая переменная представляет собой дискретные положительные целые числа (числа). Я хочу проверить, есть ли тенденция со временем (или нет). Независимая переменная (x) находится в диапазоне 0-500, а зависимая переменная (y) находится в диапазоне 0-8. Я думал, что я отвечу на это, подгоняя регрессию …

9 r  regression  python 

Если у меня есть переменная с 4 уровнями, теоретически мне нужно использовать 3 фиктивные переменные. На практике, как это на самом деле осуществляется? Я использую 0-3, я использую 1-3 и оставляю 4 пустыми? Какие-либо предложения? ПРИМЕЧАНИЕ: я собираюсь работать в R. ОБНОВЛЕНИЕ: Что случилось бы, если бы я только использовал …

9 r  regression  categorical-data  categorical-encoding 

Я пытаюсь визуализировать некоторые потребительские данные, которые имеют 4 категории. Пользователи могут свободно переключаться между различными категориями. Я хотел бы визуализировать последние три или четыре переключателя для каждого человека. Таким образом, мы начнем с графика с колонкой с 4 сложенными пропорциями. После этого у нас будет 16, поскольку каждая категория …

9 r  data-visualization  sankey-diagram