Вопросы по теме 'gam'
отдельные игровые участки в r
Я использовал обобщенную аддитивную модель для анализа некоторых данных и теперь хочу представить результат.
library(mgcv)
set.seed(2) ## simulate some data...
dat <- gamSim(1,n=400,dist="normal",scale=2)
b <-...
1288 просмотров
schedule
05.12.2022
Как вручную указать внешние узлы для более сглаживания в гамме (пакет mgcv)
Я подгоняю модели GAM к данным с помощью пакета mgcv в R. Некоторые из моих предикторов являются круговыми, поэтому я использую периодическое сглаживание. Я столкнулся с проблемой перекрестной проверки, когда мой набор данных удержания может...
2499 просмотров
schedule
02.01.2024
Проблема с функцией прогнозирования в R для одного биостека, а не для другого
Я занимаюсь моделированием распространения видов (моделирование экологических ниш), где я проецирую модель на текущие или будущие климатические растры (переменные Bioclim).
Когда я предсказываю текущие растры (object=bio_stack), все работает,...
267 просмотров
schedule
08.12.2022
ошибка на ggplot2 stat_smooth
Я получаю сообщение об ошибке, когда пытаюсь использовать ggplot2, но я не понимаю, почему. Я, наверное, где-то делаю небольшую ошибку...
> head(d)
TargetGroup2012 TargetGroup2000 bmi3 age3 PA_Score education3 asthma3 allasthma3...
3512 просмотров
schedule
11.11.2022
Как извлечь подогнанные значения GAM {mgcv} для каждой переменной в R?
Я ищу способ добавить предсказанные (реальные, не стандартизованные) значения каждой отдельной переменной в моей модели.
> model<-gam(LN_Brutto~s(agecont,by=Sex)+factor(Sex)+te(Month,Age)+s(Month,by=Sex),
data=bears)
Это краткое...
4532 просмотров
schedule
06.07.2022
Ошибка при использовании функции сглаживания micv в прогнозе .gam (mgcv)
При прогнозировании модели GAM с пакетом mgcv , построенным с использованием функции micv smooth (из пакета scam ) для одной из вариаций, я получаю следующую ошибку:
Error in predict.gam(mod_gam, newdata = dfGAMDev) :
number of items to...
132 просмотров
schedule
08.04.2022
Как изменить ось Y на графиках факторов MGCV GAM?
Я хотел бы настроить пределы y на графиках функций и факторов сглаживания GAM, чтобы было проще сравнивать эффекты двух разных моделей GAM для разных сезонов. При использовании параметра ylim в функции plot.gam я могу изменить ось Y только в...
1535 просмотров
schedule
03.06.2023
mgcv: predict.gam() дает разные результаты для type = terms и type = response
Я хочу оценить каждый компонент предиктора из модели GAM отдельно, используя опцию type="terms" . В качестве проверки работоспособности я сравнил результаты с оценкой общего прогноза с использованием опции type="response" .
Оказывается,...
1121 просмотров
schedule
16.01.2024
ошибки mgcv::gamm() и MuMIn::dredge()
Я пытался подогнать несколько GAM с помощью пакета mgcv внутри функции и грубо выбрать наиболее подходящую модель с помощью процедур выбора модели. Но моя функция запускает первую модель, а затем, похоже, снова не распознает входные данные dat ....
652 просмотров
schedule
08.10.2022
Пакет `gam`: обнаружен дополнительный сдвиг при наброске данных в `plot.gam`
Я пытаюсь установить GAM, используя пакет gam (я знаю, что mgcv более гибкий, но здесь мне нужно использовать gam ). У меня сейчас проблема в том, что модель выглядит хорошо, но по сравнению с исходными данными она как бы смещена по оси Y на...
324 просмотров
schedule
10.06.2022
Почему gam::step.gam возвращает NULL при прямом выборе?
У меня есть большая обобщенная аддитивная модель (GAM), сделанная из 10 000 наблюдений с ~ 100 переменными. Построение модели с прямым пошаговым выбором приводит к объекту класса "NULL". Почему это может быть и как решить?
library(gam)...
699 просмотров
schedule
12.05.2023
Порядок уровней факторов на термографике
Мне было интересно, можете ли вы изменить порядок уровней факторов на термплоте, поскольку в настоящее время они расположены в алфавитном порядке, а я требую, чтобы они были в хронологическом порядке.
Это был код, который я использовал для создания...
160 просмотров
schedule
27.02.2024
mgcv: как указать взаимодействие между Smooth и Factor?
В R я хотел бы подогнать игровую модель с категориальными переменными. Я думал, что смогу сделать это как с (cat — категориальная переменная).
lm(data = df, formula = y ~ x1*cat + x2 + x3);
Но я не могу делать такие вещи, как:
gam(data =...
6437 просмотров
schedule
18.04.2023
GAM в R: расчет вклада отдельных терминов
Скажем, у меня есть GAM в R, рассчитанный с помощью функции gam пакета mgcv . Я могу получить прогнозируемые значения модели для данного набора данных следующим образом:
# Load data
data("mtcars")
# Fit model
g <- gam(mpg ~ s(disp, k = 3)...
91 просмотров
schedule
02.11.2022
Удалите точки данных со слишком большим влиянием на игру.
Я поговорил со статистиком, чтобы спросить его, как я могу определить и удалить точки, которые имеют слишком большое влияние на гам (обобщенная аддитивная модель).
Он сказал мне, что я могу сделать это на основе матрицы влияние/проекция/шляпа. Я...
1930 просмотров
schedule
01.04.2023
Обобщенная аддитивная модель — Python
Я пытаюсь подобрать нелинейную модель, используя обобщенную аддитивную модель. Как определить количество используемых сплайнов. Есть ли особый способ выбрать количество сплайнов? Я использовал шлицевой фитинг 3-го порядка (кубический). Ниже...
4195 просмотров
schedule
06.05.2024
Извлечь объект модели GAM
Предположим, я создаю свою модель GAM следующим образом:
a <- runif(10)
b <- runif(10)
gm <- gam(a ~ ns(b, df=2))
plot(gm, all.terms=T, shade=T)
Тогда я получаю следующий график:
Мне нужен доступ к элементам gm , чтобы...
343 просмотров
schedule
17.07.2023
Могу ли я оценить изменяющийся во времени сезонный эффект в R с помощью GAMM?
Я хотел бы использовать обобщенную аддитивную модель для исследования данных временных рядов в R. Мои данные ежемесячные, и я хотел бы оценить сезонный эффект и эффект долгосрочного тренда. Я следил за некоторыми полезными сообщениями Гэвина Симпсона...
1095 просмотров
schedule
08.07.2022
Ошибка mgcv bam (): невозможно выделить вектор размером 99,6 ГБ
Я пытаюсь подогнать аддитивную смешанную модель с помощью bam (библиотека mgcv). В моем наборе данных содержится 10 ^ 6 наблюдений из лонгитюдного исследования роста 2,10 ^ 5 детей, находящихся в 300 медицинских центрах. Ищу уклон для каждого...
1313 просмотров
schedule
03.07.2022
Извлечение случайных эффектов из объекта gamlss
Я запустил бета-регрессию с моделью случайных перехватов, используя пакет gamlss в R. Вызов выглядит так:
uik.max.model <- gamlss(formula = max.opp.vote ~ n_commute + vote_commute + uik_dummy +
opp.gd + opp.m + pay +...
1353 просмотров
schedule
05.04.2023