Вопросы по теме 'gam'
построение правильных значений для гамма-модели в r
У меня есть следующий пример (в формате, аналогичном моим данным):
set.seed(1)
RandData <- runif(8760*2)
Locations <- rep(c('UK','France'),each=8760)
Date = seq(from=as.POSIXct("1991-01-01 00:00"),
to=as.POSIXct("1991-12-31...
669 просмотров
schedule
21.11.2023
Как включить gam() в xyplot() из пакета Lattice?
Я пытаюсь включить обобщенные аддитивные модели gam() из пакета mgcv либо в функцию xyplot() , либо в функцию coplot() из пакета lattice в R.
Данные можно найти в http://statweb.stanford.edu/~tibs/ElemStatLearn/ , выбрав данные по озону....
433 просмотров
schedule
09.03.2023
Визуализация моделей GAM в R
Кто-нибудь знает, как очень хорошо визуализировать гладкую составляющую gam моделей в R? Мне бы очень хотелось визуализировать что-то вроде вывода функции visreg . Этот код ниже иллюстрирует мою проблему
library(gam)
f=function(v){exp(v)}...
818 просмотров
schedule
04.04.2022
mgcv: Как сделать пошаговую регрессию с помощью модели ответа Твиди?
Есть ли у кого-нибудь идеи, как сделать пошаговую регрессию с Tweedie в R?
Я нашел пакет mgcv , который, по-видимому, рассматривает параметр мощности Tweedie как еще один параметр, который необходимо оценить. Кажется, это лучше, чем использование...
2496 просмотров
schedule
27.08.2022
Использование CARET вместе с GAM (метод gamSpline) в регрессии Пуассона R
Я пытаюсь использовать пакет Caret для настройки параметра «df» игровой модели для моего когортного анализа.
Со следующими данными:
cohort = 1:60
age = 1:26
grid = data.frame(expand.grid(age = age, cohort = cohort))
size = data.frame(cohort =...
2204 просмотров
schedule
17.07.2022
Пакет Caret - GAM с перекрестной проверкой как с гладкими, так и с линейными предикторами
Я хотел бы перекрестно проверить модель GAM, используя каретку. Моя модель GAM имеет бинарную переменную результата, изотропное сглаживание пар координат широты и долготы, а затем линейные предикторы. Типичный синтаксис при использовании mgcv:...
6259 просмотров
schedule
27.03.2022
mgcv: коэффициент построения графика сглаживается
Я хочу построить сплайн-эффект параметра под названием «NO2» на массу тела при рождении, но мне нужны 4 графика для четырех квартилей. Мой текущий код дает только один график, не могли бы вы помочь мне разобраться в проблеме? Вы можете увидеть код в...
747 просмотров
schedule
15.07.2022
гамма-модели в R
Я пытаюсь адаптировать свою модель к gamm (которую я использовал в lmer ).
Моя бывшая формула была
b <- lmer(metric1~a+b+c+d+e+f+g+h+i+(1|X) + (1|Y), data = dataset)
Я пытаюсь адаптировать эту формулу к моделированию с помощью gamm ....
566 просмотров
schedule
11.07.2023
gam в mgcv R с большим количеством ковариат
Я хотел бы знать, есть ли другой способ написать функцию:
gam(VariableResponse ~ s(CovariateName1) + s(CovariateName2) + ... + s(CovariateName100),
family = gaussian(link = identity), data = MyData)
в пакете mgcv без ввода имени 100...
627 просмотров
schedule
24.06.2023
Матрица ковариации оцененных параметрических коэффициентов и оцененных параметров сглаживания в GAM (пакет: mgcv)?
Рассмотрим приведенный ниже код для соответствия обобщенной аддитивной модели, включающей два члена x0, который является линейным, и x1, который является нелинейным:
library(mgcv)
set.seed(2) ## simulate some data...
dat <-...
227 просмотров
schedule
12.06.2023
Как я могу получить данные за кривой ggplot2::geom_smooth?
взять этот пример:
ggplot(iris,aes(Sepal.Width,Sepal.Length)) + geom_smooth()
Как получить x и y , которые использовались для построения этой синей кривой?
Этот случай будет рассчитан loess автоматически, мой реальный случай будет...
139 просмотров
schedule
30.03.2024
Ошибка в нулевой завышенной модели (zip) с бамом
Я пытаюсь использовать bam для запуска следующей обобщенной аддитивной модели:
m <- bam(result ~ factor(city) + factor(year) + lnpopulation + s(lnincome), data=full_df, na.action=na.omit, family=ziP(theta = NULL, link = "identity",b=0))...
162 просмотров
schedule
04.06.2023
Как указать нелинейное взаимодействие двух факторных переменных в обобщенных аддитивных моделях [R]
У меня есть набор данных временного ряда, который содержит переменную результата, которая является непрерывной, и два факторных предиктора (один с 6 уровнями и один с 2 уровнями).
Я хотел бы смоделировать нелинейное взаимодействие двух факторных...
1301 просмотров
schedule
01.09.2022
Оптимизация степеней свободы в сплайн-регрессии
У меня есть два набора данных о времени экспрессии генов:
Во-первых, экспрессию генов измеряли в течение 14 моментов времени из 4 групп:
df1 <- structure(list(val = c(-0.1, -0.13, -0.4, -0.3, -0.3, -0.2, -0.24,...
342 просмотров
schedule
16.01.2023
Ошибка в функции gam в именах (x) ‹- значение: атрибут 'names' должен быть той же длины, что и вектор
Я использую пакет mgcv для моделирования концентрации загрязнения озоном в соответствии с некоторыми ковариатами окружающей среды. Модель принимает вид:
model1 <- gam(O3 ~ s(X, Y, bs = "tp", k = 10) + wd + s(date, bs = "cc", k = 100) +...
5115 просмотров
schedule
01.02.2023
Экономия места при сохранении подогнанного объекта гам (mgcv :: gam and scam :: scam)
Я оцениваю очень простую модель на большом наборе данных. Формула выглядит как
scam::scam(formula = ratio ~ s(rate,bs="mpi"))
Эти модели затем используются для создания прогнозов для новых данных. Больше ничего о модели меня не волнует....
314 просмотров
schedule
15.09.2023
Степени свободы сглаживающего сплайна соответствуют формуле GAM
Если сглаживающий сплайн — это просто натуральный кубический сплайн с узлами при каждом уникальном значении x_i . Тогда зачем gam::s() в R нужна степень свободы?
249 просмотров
schedule
15.12.2022
Несколько независимых переменных с использованием GAM в R?
Я новичок в R и пытаюсь использовать функцию mgcv::gam для моделирования своих данных, но у нее много пояснительных переменных. Я знаю, что с обычной моделью линейной регрессии вы можете просто использовать . для учета переменных, но работает ли...
547 просмотров
schedule
30.06.2023
Есть ли способ получить значения predic.gam (, type = terms), которые НЕ центрированы?
Исходный вопрос:
Вызов predict.gam(..., type="terms") возвращает значения, центрированные по среднему значению. Есть ли способ получить необработанные значения прогнозируемых терминов (т. Е. Те, которые не были центрированы в среднем)?...
502 просмотров
schedule
05.12.2022
Моделирует GAM в H2o от R
Я пытаюсь воспроизвести пример модели GAM в h2o documentation - GAM , однако я получаю следующую ошибку:
*Error: water.exceptions.H2OModelBuilderIllegalArgumentException: Illegal argument(s) for GAM model: GAM_model_R_1586448366888_1. Details:...
222 просмотров
schedule
16.04.2022