Мне нужен пакет, который может дать мне уравнение для линейной модели 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 Тем...
Я знаю, что это, вероятно, основной вопрос ... Но я, кажется, не могу найти ответ. Я подгоняю GLM к семейству Пуассонов, а затем попытался взглянуть на прогнозы, однако смещение, похоже, принимается во внимание: model_glm=glm(cases~rhs(data$year,2003)+lhs(data$year,2003), offset=(log(population)),...
@whuber продемонстрировал, как моделировать многовариантные результаты ( , и ) для одной временной точки.у 2 у 3Y1y1y_1Y2y2y_2Y3y3y_3 Как известно, продольные данные часто встречаются в медицинских исследованиях. Мой вопрос заключается в том, как имитировать повторные измерения многовариантных...
Я использую Root Mean Squared Error(RMSE) для измерения точности значений, прогнозируемых с помощью модели. Я понимаю, что возвращаемое значение использует единицы измерения (а не процент). Тем не менее, я хотел бы процитировать мои значения в процентах. Подход, который я выбрал, заключается в...
Модель Брэдли – Терри – Люса (BTL) утверждает, что , где - вероятность того, что объект j будет оценен как «лучший», тяжелее, и т. д., чем объект i , и \ delta_i , и \ delta_j являются параметрами.pji=logit−1(δj−δi)pji=logit−1(δj−δi)p_{ji} = logit^{-1}(\delta_j -...
Мои данные - это временной ряд занятого населения, 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...
Я читаю A. Agresti (2007), Введение в категориальный анализ данных , 2-е. редакция, и я не уверен, правильно ли я понимаю этот параграф (с.106, 4.2.1) (хотя это должно быть легко): В Таблице 3.1, посвященной храпу и сердечным заболеваниям в предыдущей главе, 254 пациента сообщали о храпе каждую...
Я использую пакет randomForest в R (R версия 2.13.1, randomForest версия 4.6-2) для регрессии и заметил значительный сдвиг в моих результатах: ошибка прогнозирования зависит от значения переменной отклика. Высокие значения недооценены, а низкие значения переоценены. Сначала я подозревал, что это...
Следующие прививки взяты из этой статьи . Я новичок в начальной загрузке и пытаюсь реализовать параметрическую, полупараметрическую и непараметрическую загрузку начальной загрузки для линейной смешанной модели с R bootпакетом. Код R Вот мой Rкод: library(SASmixed) library(lme4) library(boot)...
Как я могу оценить 95% доверительные интервалы, используя профилирование для параметров, оцениваемых путем максимизации логарифмической функции правдоподобия с использованием optim в R? Я знаю, что могу асимптотически оценить ковариационную матрицу, инвертировав гессиан , но я обеспокоен тем, что...
Закрыто. Этот вопрос не по теме . В настоящее время он не принимает ответы. Хотите улучшить этот вопрос? Обновите вопрос, чтобы он соответствовал теме перекрестной проверки. Закрыто 2 года назад . Следующий код PredictNew <- predict (glm.fit, newdata = Predict, X1 =X1, Y1= Y1, type =...
У меня есть проблема, которая, я думаю, должна быть простой, но я не могу ее решить. Я смотрю на опыление семян, у меня есть растения (n = 36), которые цветут в кластерах, я выбираю 3 цветочных кластера от каждого растения и 6 семенных коробочек из каждого кластера (всего 18 семенных коробочек от...
Я пытаюсь написать сценарий R для имитации повторяющихся экспериментов интерпретации 95% доверительного интервала. Я обнаружил, что он переоценивает долю случаев, когда истинная популяционная ценность доли содержится в 95% ДИ выборки. Не большая разница - около 96% против 95%, но это, тем не менее,...
Я пытаюсь проанализировать влияние года на переменную logInd для конкретной группы лиц (у меня есть 3 группы). Самая простая модель: > fix1 = lm(logInd ~ 0 + Group + Year:Group, data = mydata) > summary(fix1) Call: lm(formula = logInd ~ 0 + Group + Year:Group, data = mydata) Residuals: Min 1Q...
Вот краткий вывод использованной мною модели Кокша (я использовал R, и вывод основан на наилучшей конечной модели, т.е. включены все значимые объясняющие переменные и их взаимодействия): coxph(formula = Y ~ LT + Food + Temp2 + LT:Food + LT:Temp2 + Food:Temp2 + LT:Food:Temp2) #...
Закрыто. Этот вопрос не по теме . В настоящее время он не принимает ответы. Хотите улучшить этот вопрос? Обновите вопрос, чтобы он соответствовал теме перекрестной проверки. Закрыто 2 года назад . R , кажется , чтобы иметь возможность вывод хороших краткие заговоры от bugsи jagsобъектов ,...
Я установил модель ARIMA (5,1,2), используя auto.arima()функцию в R, и, посмотрев порядок, мы можем сказать, что это не лучшая модель для прогнозирования. Если в рядах данных существуют выбросы, каков метод для подгонки модели к таким...
Я выставляю статистику по 5 IV (5 черт личности, экстраверсия, приятность, добросовестность, невротизм, открытость) против 3 DVs Отношение к PCT, Отношение к CBT, Отношение к PCT против CBT. Я также добавил возраст и пол, чтобы увидеть, какие еще есть эффекты. Я тестирую, чтобы понять, могут ли...
Предположим, у меня есть одна выборка частот из 4 возможных событий: Event1 - 5 E2 - 1 E3 - 0 E4 - 12 и у меня есть ожидаемые вероятности того, что мои события произойдут: p1 - 0.2 p2 - 0.1 p3 - 0.1 p4 - 0.6 С суммой наблюдаемых частот моих четырех событий (18) я могу рассчитать ожидаемые частоты...
Документация утверждает, что R gbm с distribution = "adaboost" может использоваться для задачи классификации 0-1. Рассмотрим следующий фрагмент кода: gbm_algorithm <- gbm(y ~ ., data = train_dataset, distribution = "adaboost", n.trees = 5000) gbm_predicted <- predict(gbm_algorithm,...