Описание тега vgam
1
ответ
Tobit: Ошибка в lm.fit(X.vlm, y = z.vlm, ...): NA/NaN/Inf в "x" (смоделированные данные)
Я провожу исследование Монте-Карло. У меня есть линейная модель с гетероскедастичностью и левой цензурой зависимой переменной на 0. Среднее значение цензуры составляет 25,9. Я получаю ошибку Error in lm.fit(X.vlm, y = z.vlm, ...) : NA/NaN/Inf in 'x'…
02 окт '15 в 15:56
1
ответ
Полиномиальный / условный логит с сгруппированными данными в R
Я ищу способ подгонки условного логита к сгруппированным данным в R. У меня есть полиномиальный выбор подсчетов с переменными отклика и индивидуально. К сожалению, расширенный набор данных будет слишком большим, чтобы с ним работать. Таким образом, …
04 июн '16 в 10:59
1
ответ
Связь функции "идентичности" с гамма-распределением в пакете VGAM R
Я использую VGAM пакет, чтобы сделать регрессию между двумя выборками, которые распределены по гамме. Можно ли выбрать функцию идентификации для ссылки? Если да, то как? Я использую эти строки кода, и это работает: fit <- vglm(Y ~ X , family = ga…
08 янв '15 в 22:31
1
ответ
Переполнение функции логита VGAM
VGAM версия 0.93 > logit(1000, inverse=T) [1] 0 # it should be 1 Проблема здесь: exp(1000 - log1p(exp(1000))) Вот log1p(exp(1000)) становится Inf Таким образом, численный метод, который он использует, не обрабатывает большие числа, по сравнению с…
25 апр '14 в 13:51
1
ответ
Как я могу получить logLik или AIC VGLM подходят
Я не мог получить их обычными методами. > x <- rnorm(1000, 10, 1) > y <-rgamma(1000, shape=1, rate=x) > myfit <- vglm(y~x, family=gammaff(link="reciprocal")) > summary(myfit) Call: vglm(formula = y ~ x, family = gammaff(link = "…
25 июл '14 в 13:33
1
ответ
График процентили VGAM в R
Я запускаю следующий код из файлов справки VGAM: library(VGAM) fit4 <- vgam(BMI ~ s(age, df = c(4, 2)), lms.bcn(zero = 1), data = bmi.nz, trace = TRUE) qtplot(fit4, percentiles = c(5,50,90,99), main = "Quantiles", las = 1, xlim = c(15, 90), ylab …
18 дек '14 в 17:59
0
ответов
Таблица вывода итоговых данных латекса или HTML для объектов регрессии vglm (VGAM)
Я пытаюсь получить латексный или html-вывод результатов регрессии модели VGAM (в приведенном ниже примере это обобщенный порядковый логит). Но пакеты, которые я знаю для этой цели, не работают с объектом vglm. Здесь вы можете увидеть небольшой приме…
04 май '15 в 18:25
1
ответ
Пакет V VGAM: после добавления второй пояснительной
Почему подгонка ухудшается после добавления второй объясняющей переменной? require("VGAM") df = data.frame(x = c(1,2,3,4,5,6,7,8,9,10), y = c(1,4,8,15,25,36,48,65,80,105), z = c(0,0,0,1,100,400,900,1600,1800,200) ) vgt1 = vgam(y~s(x, df=2), data=df,…
05 фев '14 в 01:03
1
ответ
Сравните модель подходит для двух GAM
У меня есть матрица Expr со строками, представляющими переменные и образцы столбцов. У меня есть категорический вектор под названием groups (содержит "A","B" или "C") Я хочу проверить, какая из переменных "Expr" может быть объяснена тем фактом, что …
21 май '14 в 17:06
1
ответ
Нужна помощь, я застрял при установке (VGAM) и как использовать mlogit. Я не понимаю из этого примера PDF
Мне нужна помощь для моих проблем. У меня 2 проблемы. Я не могу установить библиотеку (VGAM) на RStudio. Есть идеи для другого регрессионного логистического порядкового пакета или есть решение для моей проблемы? Застрял на первом шаге, когда использ…
11 апр '16 в 19:53
1
ответ
s() в mgcv не работает при загрузке пакета VGAM
Это x <- rnorm(100) y <- rnorm(100) gam(y ~ s(x)) ## Family: gaussian ## Link function: identity ## Formula: ## y ~ s(x) ## Estimated degrees of freedom: ## 1 total = 2 ## GCV score: 0.8116283 ломается, когда VGAM пакет загружен: library(VGAM)…
19 дек '13 в 23:45
1
ответ
vglm в пакете VGAM для порядковых переменных: что обозначают буквы L и Q?
Я подгоняю модель кумулятивного логита с использованием порядкового выражения, коэффициенты, которые я получаю, следующие: coefficients(cum_mod_podds_ord,matrix=TRUE) logit(P[Y>=2]) logit(P[Y>=3]) (Intercept) -0.4106068 -1.2504301 depr_ord.L 1…
04 июн '14 в 14:44
1
ответ
Оценка максимального правдоподобия бета-нормы в R
Я хочу оценить параметры бета-нормального распределения. Я использовал пакет maxLik library(VGAM) library(maxLik) alfa=2;beta=3;mu=0;sigma=1 n=100 x=rbetanorm(n,alfa,beta,mu,sigma) logLikFun=function(w){ alfa=w[1] beta=w[2] mu=w[3] sigma=w[4] ll={-n…
09 окт '17 в 22:48
1
ответ
R: Несовместимые ошибки измерения vglm функция в VGAM
TL;DR Я запускаю регрессии Tobit с пакетом VGAM в R - Вот игрушечный набор данных, который постоянно дает мне ошибку, которую я не смог диагностировать: library(data.table) library(VGAM) > sessionInfo()$otherPkgs $VGAM Package: VGAM Version: 0.9-…
19 дек '14 в 16:44
1
ответ
Процентили от VGAM
Я использую следующий пример со страниц справки пакета VGAM library(VGAM) fit4 <- vgam(BMI ~ s(age, df = c(4, 2)), lms.bcn(zero = 1), data = bmi.nz, trace = TRUE) qtplot(fit4, percentiles = c(5,50,90,99), main = "Quantiles", las = 1, xlim = c(15,…
18 дек '14 в 12:55
2
ответа
Ошибка добавления взаимодействий в пользовательскую модель VGAM::vglm в карете
Я построил пользовательскую модель с caret с помощью vglm() от VGAM, Он отлично работает с простыми эффектами, но когда я пытаюсь добавить взаимодействия, он не работает с object 'x1:x2' not found сообщение об ошибке, где x1 а также x2 переменные пр…
14 мар '16 в 17:53
0
ответов
gamma2 ссылка на функцию выбора VGAM пакета R
Я использую пакет VGAM в R с gamma2 семейная функция (два параметра: му и форма). Я колеблюсь в выборе функции ссылки. На самом деле, у меня есть два следующих подбора: fit1 <- vglm(y ~ x, family = gamma2(lmu ="identitylink", lshape="identitylink…
22 сен '15 в 18:30
1
ответ
Как использовать надежный se и кластерный se с моделью vglm tobit?
Я пытаюсь перенести модель тобита из Stata в R. Команды Stata для робастных было бы просто добавить ,vce(robust) к модели. И для кластеризации это было бы ,vce(cluster idvar), Воспроизводимый Stata пример: use http://www.ats.ucla.edu/stat/stata/dae/…
06 фев '14 в 16:30
1
ответ
Ошибка с VGAM? vglm family=posnegbinomial => "Ошибка в if (take.half.step) {: пропущенное значение там, где требуется TRUE/FALSE"
У меня есть некоторые реальные данные, которые я боюсь, что-то противное. По сути, это положительное отрицательное биномиальное распределение (без нуля). Тем не менее, есть некоторые выбросы, которые, по-видимому, приводят к ошибочным вычислениям (м…
14 фев '16 в 22:07
0
ответов
Построение полиномиальной логистической регрессии VGLM с 95% ДИ
Я хотел бы построить прогнозируемые значения полиномиальной логистической регрессии, полученные из vglm() функция в пакете VGAM. Важно, чтобы я использовал VGAM, потому что я пытаюсь повторить анализ коллеги, проведенный в Stata, который я достиг с …
28 июн '16 в 05:17