Вопросы с тегом «r»

Например, в R, если вы вызываете acf()функцию, она строит коррелограмму по умолчанию и рисует 95% доверительный интервал. Глядя на код, если вы звоните plot(acf_object, ci.type="white"), вы видите: qnorm((1 + ci)/2)/sqrt(x$n.used) в качестве верхнего предела для типа белого шума. Кто-нибудь может объяснить теорию, стоящую за этим методом? Почему мы получаем qnorm …

12 r  confidence-interval  autocorrelation 

Я пытаюсь выяснить, как определить количество слогов в корпусе аудиозаписей. Я думаю, что хороший прокси может быть пиками в волновом файле. Вот то, что я попробовал с файлом моего разговора на английском языке (мой фактический случай использования на кисуахили). Стенограмма записи этого примера: «Я пытаюсь использовать функцию таймера. Я смотрю …

12 r  signal-processing 

Я сделал измерений двух переменных и . Они оба имеют известные неопределенности и связанные с ними. Я хочу найти связь между и . Как мне это сделать?x y σ x σ y x ynnnxxxyyyσxσx\sigma_xσyσy\sigma_yxxxyyy РЕДАКТИРОВАТЬ : каждый имеет различные связанные с ним, и то же самое с .σ x , …

12 r  regression  deming-regression 

Я думал, что понял эту проблему, но теперь я не так уверен, и я хотел бы проверить с другими, прежде чем продолжить. У меня есть две переменные, Xи Y. Yявляется отношением, и оно не ограничено 0 и 1 и обычно нормально распределено. Xявляется пропорцией, и он ограничен 0 и 1 …

12 r  regression  data-transformation  nonlinear-regression  proportion 

При моделировании данных подсчета претензий в страховой среде я начал с Пуассона, но затем заметил чрезмерную дисперсию. Квази-Пуассон лучше моделировал большее отношение средней дисперсии, чем основной Пуассон, но я заметил, что коэффициенты были идентичны как в модели Пуассона, так и в модели Квази-Пуассона. Если это не ошибка, почему это происходит? …

12 r  count-data  poisson-regression  overdispersion  quasi-likelihood 

Я хотел бы знать, как преобразовать отрицательные значения в Log(), так как у меня есть гетероскедастические данные. Я прочитал, что это работает с формулой, Log(x+1)но это не работает с моей базой данных, и я продолжаю получать NaN в результате. Например, я получаю это предупреждение (я не поместил свою полную базу …

12 r  logarithm 

Я просто подумал о аккуратном (не обязательно хорошем) способе создания одномерных оценок плотности, и мой вопрос: У этого метода оценки плотности есть имя? Если нет, то является ли это частным случаем какого-либо другого метода в литературе? Вот метод: Мы имеем вектор который мы предполагаем, взят из некоторого неизвестного распределения, которое …

12 r  estimation  nonparametric  pdf  kernel-smoothing 

TL; DR: lme4оптимизация кажется линейной по количеству параметров модели по умолчанию и намного медленнее, чем эквивалентная glmмодель с фиктивными переменными для групп. Что я могу сделать, чтобы ускорить это? Я пытаюсь соответствовать довольно большой иерархической модели логита (~ 50 тыс. Строк, 100 столбцов, 50 групп). Подгонка нормальной модели логита к …

12 r  mixed-model  optimization  lme4-nlme  numerics 

Как мне прочитать результаты теста Данна ? В частности, что говорят значения в таблице ниже? У меня есть непараметрические данные в 4 группах, и я сначала выполнил тест Крускала-Уоллиса, чтобы подтвердить, что распределения групп отличались друг от друга и совокупного набора данных. Затем я использовал тест Данна, чтобы увидеть, какие …

12 r  distributions  kruskal-wallis  dunn-test 

Скажем, у меня есть следующие четырехмерные данные, где первые три можно рассматривать как координаты, а последние можно рассматривать как значения. c1, c2, c3, value 1, 2, 6, 0.456 34, 34, 12 0.27 12, 1, 66 0.95 Как лучше визуализировать влияние первых трех координат на последнее значение? Я знаю о трех …

12 r  data-visualization 

Я разрабатывал модель логистической регрессии на основе ретроспективных данных из национальной базы данных о травмах головы в Великобритании. Ключевым результатом является 30-дневная смертность (обозначается как «выживаемая» мера). Другие меры с опубликованным доказательством существенного влияния на результат в предыдущих исследованиях включают в себя: Year - Year of procedure = 1994-2013 Age …

12 r  logistic  generalized-linear-model  goodness-of-fit  regression-strategies 

R version 3.1.1 (2014-07-10) -- "Sock it to Me" > bl <- c(140, 138, 150, 148, 135) > fu <- c(138, 136, 148, 146, 133) > t.test(fu, bl, alternative = "two.sided", paired = TRUE) Error in t.test.default(fu, bl, alternative = "two.sided", paired = TRUE) : data are essentially constant Затем …

12 r  t-test 

Я - пользователь, более знакомый с R, и пытался оценить случайные уклоны (коэффициенты отбора) примерно для 35 особей в течение 5 лет для четырех переменных среды обитания. Переменная ответа - является ли место «использованным» (1) или «доступным» (0) местом обитания («использование» ниже). Я использую Windows 64-битный компьютер. В версии R …

12 r  binomial  sas  random-effects-model  lme4-nlme 

Резюме: существует ли статистическая теория, поддерживающая использование распределения (со степенями свободы, основанными на остаточном отклонении) для тестов коэффициентов логистической регрессии, а не стандартного нормального распределения?Ttt Некоторое время назад я обнаружил, что при подборе модели логистической регрессии в SAS PROC GLIMMIX при настройках по умолчанию коэффициенты логистической регрессии тестируются с использованием …

12 r  logistic  mathematical-statistics  sas  degrees-of-freedom 

Я пытаюсь разобраться с пакетом посредничества в R, используя виньетку для пакета. Я изо всех сил пытаюсь понять вывод mediate()функции. require("mediation") require("sandwich") data("framing") med.fit <- lm(emo ~ treat + age + educ + gender + income, data = framing) out.fit <- glm(cong_mesg ~ emo + treat + age + educ …

12 r  mediation