Запитання з тегом «r»

Я знаю, що це, мабуть, основне питання ... Але я, здається, не знаходжу відповіді. Я підходив до ГМ з родиною Пуассона, а потім намагався ознайомитись з прогнозами, однак зміщення, здається, враховується: 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") Я отримую випадки не ставки ... Я також спробував model_glm=glm(cases~rhs(data$year,2003)+lhs(data$year,2003)+ offset(log(population)), …

9 r  generalized-linear-model  prediction  offset 

Мої дані - це часовий ряд зайнятого населення, 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 

@whuber продемонстрував, як імітувати багатоваріантні результати ( , та ) за один момент часу.у1y1y_1у2y2y_2у3y3y_3 Як ми знаємо, поздовжні дані часто зустрічаються в медичних дослідженнях. Моє запитання - як імітувати багаторазові результати в R? Наприклад, ми неодноразово вимірюємо , і у 5 різних часових точках для двох різних груп лікування.у1y1y_1у2y2y_2у3y3y_3

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 

Я шукаю простий приклад коду, як запустити фільтр частинок в Р. Здається, що пакет пам’яті підтримує біт математики простору стану, але приклади трохи складні, щоб програматично наслідувати простого розробника OO, такого як я, зокрема як завантажити спостережувані дані в помпезний об’єкт. Приклади тут: http://cran.r-project.org/web/packages/pomp/vignettes/intro_to_pomp.pdf Скажімо, у мене є вхідний файл …

9 r 

Зачинено. Це питання поза темою . Наразі відповіді не приймаються. Хочете вдосконалити це питання? Оновіть питання, щоб воно було тематичним для перехресної перевірки. Закрито 9 місяців тому . Проблема, згадана в цьому питанні, виправлена ​​у версії 1.7.3 glmnet пакета R. У мене виникають проблеми з керуванням glmnet з сім'єю = …

9 r  multinomial  glmnet 

Я створив власну трохи вдосконалену версію термплота, яку я використовую в цьому прикладі, ви можете знайти її тут . Я раніше розміщував повідомлення про SO, але чим більше я думаю про це, я вважаю, що це, ймовірно, більше пов'язане з інтерпретацією моделі небезпеки Кокса, ніж з фактичним кодуванням. Проблема Коли …

9 r  survival  cox-model 

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

9 r  regression  python 

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

9 r  regression  categorical-data  categorical-encoding 

Я намагаюся візуалізувати деякі дані споживачів, у яких є 4 категорії. Користувачі можуть перемикатися між різними категоріями. Я хотів би візуалізувати останні три-чотири вимикачі для кожної людини. Отже, ми б почали з сюжету зі стовпцем з 4 складеними пропорціями. Після цього у нас було б 16, оскільки кожна категорія розбивається …

9 r  data-visualization  sankey-diagram