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

Используйте этот тег для любого * по теме * вопроса, который (a) включает `R` либо в качестве критической части вопроса, либо в ожидаемом ответе, & (b) не * просто * о том, как использовать` R`.

1
Логистическая регрессия для временных рядов
Я хотел бы использовать бинарную модель логистической регрессии в контексте потоковых данных (многомерных временных рядов), чтобы предсказать значение зависимой переменной данных (то есть строки), которые только что прибыли, учитывая прошлые наблюдения. Насколько я знаю, логистическая регрессия традиционно используется для посмертного анализа, где каждая зависимая переменная уже установлена ​​(либо путем проверки, …

1
Подгонка экспоненциальной модели к данным
Этот вопрос был перенесен из переполнения стека, потому что на него можно ответить по перекрестной проверке. Мигрировал 8 лет назад . У меня есть 2 переменные, оба из класса "числовой": > head(y) [1] 0.4651804 0.6185849 0.3766175 0.5489810 0.3695258 0.4002567 > head(x) [1] 59.32820 68.46436 80.76974 132.90824 216.75995 153.25551 Я построил …
21 r 

2
Как использовать веса в функции lm в R?
Locked . Этот вопрос и его ответы заблокированы, потому что вопрос не по теме, но имеет историческое значение. В настоящее время он не принимает новые ответы или взаимодействия. Может ли кто-нибудь предложить несколько советов о том, как использовать weightsаргумент в lmфункции R ? Скажем, например, вы пытались вписать модель в …
21 r  regression 

1
Генерация коррелированных биномиальных случайных величин
Мне было интересно, возможно ли генерировать коррелированные случайные биномиальные переменные, следуя подходу линейного преобразования? Ниже я попробовал что-то простое в R, и оно дает некоторую корреляцию. Но мне было интересно, есть ли принципиальный способ сделать это? X1 = rbinom(1e4, 6, .5) ; X2 = rbinom(1e4, 6, .5) ; X3 = …

2
Как работает метод обратного преобразования?
Как работает метод инверсии? Скажем , у меня случайную выборку X1,X2,...,XnX1,X2,...,XnX_1,X_2,...,X_n с плотностью f(x;θ)=1θx(1−θ)θf(x;θ)=1θx(1−θ)θf(x;\theta)={1\over \theta} x^{(1-\theta)\over \theta} более 0&lt;x&lt;10&lt;x&lt;10<x<1и, следовательно, cdfFX(x)=x1/θFX(x)=x1/θF_X(x)=x^{1/\theta}на(0,1)(0,1)(0,1). Тогда методом обращения я получаю распределениеXXXкакF−1X(u)=uθFX−1(u)=uθF_X^{-1}(u)=u^\theta. Так есть ли uθuθu^\theta распределение XXX ? Так работает метод инверсии? u&lt;-runif(n) x&lt;-u^(theta)

1
Почему квазипуассон в GLM не рассматривается как частный случай отрицательного бинома?
Я пытаюсь приспособить обобщенные линейные модели к некоторым наборам данных подсчета, которые могут быть или не быть перераспределены. Здесь применимы два канонических распределения: Пуассон и Отрицательный бином (Негбин) с EV и дисперсией.μμ\mu Вгп= μВaрпзнак равноμVar_P = \mu ВгNВ= μ + μ2θВaрNВзнак равноμ+μ2θVar_{NB} = \mu + \frac{\mu^2}{\theta} которые могут быть установлены …

4
Как проверить, является ли мой дистрибутив мультимодальным?
Когда я строю гистограмму моих данных, она имеет два пика: Означает ли это потенциальное мультимодальное распределение? Я запустил dip.testв R ( library(diptest)), и вывод: D = 0.0275, p-value = 0.7913 Я могу заключить, что мои данные имеют мультимодальное распределение? ДАННЫЕ 10346 13698 13894 19854 28066 26620 27066 16658 9221 13578 …

3
Как рассчитать качество посадки в glm (R)
У меня есть следующий результат от запуска функции GLM. Как я могу интерпретировать следующие значения: Нулевое отклонение Остаточное отклонение AIC Они имеют какое-то отношение к совершенству? Могу ли я рассчитать некоторую степень достоверности соответствия по этим результатам, таким как R-квадрат или любая другая мера? Call: glm(formula = tmpData$Y ~ tmpData$X1 …

1
Роль n.minobsinnode параметра GBM в R [закрыто]
Этот вопрос вряд ли поможет будущим посетителям; это относится только к небольшому географическому региону, конкретному моменту времени или необычайно узкой ситуации, которая обычно не применима к всемирной аудитории Интернета. Чтобы сделать этот вопрос более применимым, посетите справочный центр . Закрыто 7 лет назад . Я хотел знать, что означает параметр …
21 r  gbm 

3
Почему nls () выдаёт мне ошибку «матрица сингулярного градиента при начальных оценках параметров»?
У меня есть некоторые основные данные о сокращении выбросов и стоимости автомобиля: q24 &lt;- read.table(text = "reductions cost.per.car 50 45 55 55 60 62 65 70 70 80 75 90 80 100 85 200 90 375 95 600 ",header = TRUE, sep = "") Я знаю, что это экспоненциальная функция, …

4
Как спроецировать новый вектор на пространство PCA?
После выполнения анализа главных компонентов (PCA) я хочу спроецировать новый вектор на пространство PCA (т.е. найти его координаты в системе координат PCA). Я рассчитал PCA на языке R, используя prcomp. Теперь я должен быть в состоянии умножить свой вектор на матрицу вращения PCA. Должны ли главные компоненты в этой матрице …
21 r  pca  r  variance  heteroscedasticity  misspecification  distributions  time-series  data-visualization  modeling  histogram  kolmogorov-smirnov  negative-binomial  likelihood-ratio  econometrics  panel-data  categorical-data  scales  survey  distributions  pdf  histogram  correlation  algorithms  r  gpu  parallel-computing  approximation  mean  median  references  sample-size  normality-assumption  central-limit-theorem  rule-of-thumb  confidence-interval  estimation  mixed-model  psychometrics  random-effects-model  hypothesis-testing  sample-size  dataset  large-data  regression  standard-deviation  variance  approximation  hypothesis-testing  variance  central-limit-theorem  kernel-trick  kernel-smoothing  error  sampling  hypothesis-testing  normality-assumption  philosophical  confidence-interval  modeling  model-selection  experiment-design  hypothesis-testing  statistical-significance  power  asymptotics  information-retrieval  anova  multiple-comparisons  ancova  classification  clustering  factor-analysis  psychometrics  r  sampling  expectation-maximization  markov-process  r  data-visualization  correlation  regression  statistical-significance  degrees-of-freedom  experiment-design  r  regression  curve-fitting  change-point  loess  machine-learning  classification  self-study  monte-carlo  markov-process  references  mathematical-statistics  data-visualization  python  cart  boosting  regression  classification  robust  cart  survey  binomial  psychometrics  likert  psychology  asymptotics  multinomial 

2
Предупреждение «Модель не сходится» в lmer ()
В следующем наборе данных я хотел посмотреть, изменится ли ответ (эффект) в отношении сайтов, сезона, продолжительности и их взаимодействий. Некоторые онлайн-форумы по статистике предложили мне продолжить работу с линейными моделями со смешанными эффектами, но проблема в том, что, поскольку репликации рандомизированы на каждой станции, у меня мало шансов собрать образец …

4
Как создать произвольную ковариационную матрицу
Например, в R, MASS::mvrnorm()функция полезна для генерации данных, чтобы продемонстрировать различные вещи в статистике. Он принимает обязательный Sigmaаргумент, который представляет собой симметричную матрицу, определяющую ковариационную матрицу переменных. Как бы я создал симметричную матрицу с произвольными записями?n×nn×nn\times n

1
lme () и lmer () дают противоречивые результаты
Я работал с некоторыми данными, которые имеют некоторые проблемы с повторными измерениями. При этом я заметил очень различное поведение lme()и lmer()использование моих тестовых данных и хочу знать почему. Поддельный набор данных, который я создал, содержит измерения роста и веса для 10 предметов, взятых дважды каждый. Я настроил данные так, чтобы …

3
Найти способ симулировать случайные числа для этого распределения
Я пытаюсь написать программу на R, которая имитирует псевдослучайные числа из распределения с помощью кумулятивной функции распределения: F( х ) = 1 - опыт( - а х - бр + 1Икср + 1) ,х ≥ 0F(Икс)знак равно1-ехр⁡(-aИкс-бп+1Иксп+1),Икс≥0F(x)= 1-\exp \left(-ax-\frac{b}{p+1}x^{p+1}\right), \quad x \geq 0 гдеa , b &gt; 0 , p …

Используя наш сайт, вы подтверждаете, что прочитали и поняли нашу Политику в отношении файлов cookie и Политику конфиденциальности.
Licensed under cc by-sa 3.0 with attribution required.