Вопросы по теме 'mgcv'
Как включить gam() в xyplot() из пакета Lattice?
Я пытаюсь включить обобщенные аддитивные модели gam() из пакета mgcv либо в функцию xyplot() , либо в функцию coplot() из пакета lattice в R.
Данные можно найти в http://statweb.stanford.edu/~tibs/ElemStatLearn/ , выбрав данные по озону....
433 просмотров
schedule
09.03.2023
mgcv доступ к функции Assessment.gam
Я использую библиотеку R mgcv
Функция gam возвращает объект со значением edf1
> g = gam(y~s(x, k = 10, bs = "cr"), data = data, scale = 1)
> g$edf1
(Intercept) s(x).1 s(x).2 s(x).3 s(x).4 s(x).5 s(x).6...
228 просмотров
schedule
17.07.2023
mgcv: Как сделать пошаговую регрессию с помощью модели ответа Твиди?
Есть ли у кого-нибудь идеи, как сделать пошаговую регрессию с Tweedie в R?
Я нашел пакет mgcv , который, по-видимому, рассматривает параметр мощности Tweedie как еще один параметр, который необходимо оценить. Кажется, это лучше, чем использование...
2496 просмотров
schedule
27.08.2022
Извлечение узлов, базиса, коэффициентов и прогнозов для P-сплайнов в адаптивном сглаживании.
Я использую пакет mgcv , чтобы подогнать полиномиальные сплайны к некоторым данным с помощью:
x.gam <- gam(cts ~ s(time, bs = "ad"), data = x.dd,
family = poisson(link = "log"))
Я пытаюсь извлечь...
7850 просмотров
schedule
06.04.2023
Пакет 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
Ошибка в нулевой завышенной модели (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
Экономия места при сохранении подогнанного объекта гам (mgcv :: gam and scam :: scam)
Я оцениваю очень простую модель на большом наборе данных. Формула выглядит как
scam::scam(formula = ratio ~ s(rate,bs="mpi"))
Эти модели затем используются для создания прогнозов для новых данных. Больше ничего о модели меня не волнует....
314 просмотров
schedule
15.09.2023
Как использовать сглаживающие сплайны в gam в пакете R mgcv
Вопрос в том, что это правильный способ указать узлы в сглаживающем сплайне в gam в mgcv?
Путаница заключается в том, что в виньетке говорится, что k - это размерность основы, используемой для представления гладкого члена.
(Ранее я думал, что в...
595 просмотров
schedule
14.07.2023
Есть ли способ получить значения predic.gam (, type = terms), которые НЕ центрированы?
Исходный вопрос:
Вызов predict.gam(..., type="terms") возвращает значения, центрированные по среднему значению. Есть ли способ получить необработанные значения прогнозируемых терминов (т. Е. Те, которые не были центрированы в среднем)?...
502 просмотров
schedule
05.12.2022
Как интерпретировать график случайных эффектов от mgcv
У меня есть несколько вопросов относительно использования случайного эффекта в GAM. Во-первых, как вы интерпретируете и передаете выходной график?
Я смоделировал огонь как случайный эффект в этом GAM, потому что это в значительной степени...
165 просмотров
schedule
11.07.2022
Функция mgcv::gam вызывает сообщение об ошибке Ошибка переформулирования (pav, env = p.env): неиспользуемый аргумент (env = p.env)
Я пытаюсь использовать пакет mgcv, но функция mgcv::gam возвращает только сообщение об ошибке
Tmean = c(1,3,10,5,7)
prec = c(300,500,700,600,500)
fit = mgcv::gam(Tmean ~ s(prec))
Сообщение об ошибке:
Ошибка в переформулировании (pav, env...
448 просмотров
schedule
17.08.2022
Циклический адаптивный сплайн в MGCV
Я хочу подогнать GAM к данным с циклическим предиктором и с наибольшей нестабильностью в определенной части цикла.
Есть ли способ сделать циклический сплайн (bs = 'cc' или 'cp') адаптивным? ... или, что то же самое, сделать адаптивный сплайн (bs =...
106 просмотров
schedule
20.12.2023