Запитання з тегом «r»

Дорогі всі - я помітив щось дивне, чого я не можу пояснити, чи не так? Підсумовуючи: ручний підхід до обчислення довірчого інтервалу в моделі логістичної регресії та функції R confint()дають різні результати. Я пережив прикладну логістичну регресію Hosmer & Lemeshow (2-е видання). У 3-й главі є приклад обчислення коефіцієнта шансів …

34 r  regression  logistic  confidence-interval  profile-likelihood  correlation  mcmc  error  mixture  measurement  data-augmentation  r  logistic  goodness-of-fit  r  time-series  exponential  descriptive-statistics  average  expected-value  data-visualization  anova  teaching  hypothesis-testing  multivariate-analysis  r  r  mixed-model  clustering  categorical-data  unsupervised-learning  r  logistic  anova  binomial  estimation  variance  expected-value  r  r  anova  mixed-model  multiple-comparisons  repeated-measures  project-management  r  poisson-distribution  control-chart  project-management  regression  residuals  r  distributions  data-visualization  r  unbiased-estimator  kurtosis  expected-value  regression  spss  meta-analysis  r  censoring  regression  classification  data-mining  mixture 

Я працюю підручником зі статистикою під час навчання R і натрапив на камеру спотикання на наступному прикладі: Після перегляду ?quantileя спробував відтворити це в R за допомогою наступного: > nuclear <- c(7, 20, 16, 6, 58, 9, 20, 50, 23, 33, 8, 10, 15, 16, 104) > quantile(nuclear) 0% 25% …

33 r  quantiles 

Я новачок у моделюванні нейронних мереж, але мені вдалося створити нейронну мережу з усіма доступними точками даних, яка добре відповідає спостережуваним даним. Нейронна мережа була виконана в R за допомогою пакету nnet: require(nnet) ##33.8 is the highest value mynnet.fit <- nnet(DOC/33.80 ~ ., data = MyData, size = 6, decay …

33 r  neural-networks 

Я шукаю гарну книгу / підручник, щоб дізнатися про аналіз виживання. Мені також цікаві посилання на аналіз виживання в Р.

33 r  survival  references 

Я намагаюся зрозуміти стандартну помилку "кластеризації" і як виконати в R (це тривіально в Stata). У RI були невдалі, використовуючи plmабо написання власної функції. Я буду використовувати diamondsдані з ggplot2пакету. Я можу робити фіксовані ефекти з будь-якими фіктивними змінними > library(plyr) > library(ggplot2) > library(lmtest) > library(sandwich) > # with …

33 r  panel-data  standard-error  fixed-effects-model  clustered-standard-errors 

Я, як правило, використовую BIC, оскільки я розумію, що він оцінює парність сильніше, ніж AIC. Однак я вирішив зараз використовувати більш комплексний підхід і хотів би також використовувати AIC. Я знаю, що Raftery (1995) представив приємні вказівки щодо відмінностей BIC: 0-2 є слабким, 2-4 - позитивним доказом того, що одна …

32 r  model-selection  references  aic  bic 

З документації для anova(): Коли дається послідовність об'єктів, 'anova' перевіряє моделі один проти одного у визначеному порядку ... Що означає тестування моделей одна проти одної? І чому порядок має значення? Ось приклад з підручника GenABEL : > modelAdd = lm(qt~as.numeric(snp1)) > modelDom = lm(qt~I(as.numeric(snp1)>=2)) > modelRec = lm(qt~I(as.numeric(snp1)>=3)) anova(modelAdd, modelGen, …

32 r  anova 

Існує ряд надійних оцінювачів масштабу . Помітним прикладом є середнє абсолютне відхилення, яке відноситься до стандартного відхилення як . У байєсівській системі існує ряд способів чітко оцінити розташування приблизно нормального розподілу (скажімо, нормального зараження сторонніми людьми), наприклад, можна припустити, що дані поширюються як при розподілі, так і по розподілу Лапласа. …

32 r  bayesian  estimation  standard-deviation  robust 

Я вписую lm()модель до набору даних, що включає показники за фінансовий квартал (Q1, Q2, Q3, що робить Q4 за замовчуванням). Використовуючи lm(Y~., data = data) Я отримую NAяк коефіцієнт для Q3, і попередження про те, що одна змінна була виключена через особливості. Чи потрібно додати стовпчик Q4?

32 r  regression 

Ви можете мати дані у широкому або довгому форматі. Це досить важлива річ, оскільки корисні методи різні, залежно від формату. Я знаю, що вам потрібно працювати melt()і зcast() з пакетом переформатування, але, здається, деякі речі я не розумію. Чи може хтось дати мені короткий огляд, як ви це робите?

32 data-transformation  r 

Я хочу встановити багаторівневий GLMM з розподілом Пуассона (з надмірною дисперсією) за допомогою R. На даний момент я використовую lme4, але я помітив, що нещодавно quasipoissonсім’я була видалена. Я десь бачив, що можна моделювати адитивну наддисперсію для біноміальних розподілів, додаючи випадковий перехоплення з одним рівнем за спостереження. Чи стосується це …

32 r  mixed-model  poisson-distribution  lme4-nlme  overdispersion 

Я хотів би регресувати вектор B проти кожного зі стовпців у матриці А. Це тривіально, якщо відсутні дані, але якщо матриця A містить відсутні значення, то моя регресія проти A обмежена, щоб включати лише рядки, де всі значення присутні ( поведінка na.omit за замовчуванням ). Це дає неправильні результати для …

32 r  missing-data  linear-model 

Скажімо, у нас є наведена нижче статистика gender mean sd n f 1.666667 0.5773503 3 m 4.500000 0.5773503 4 Як ви проводите тест-випробування на двох зразках (щоб побачити, чи є значна різниця між засобами чоловіків і жінок у певній змінній), використовуючи подібну статистику, а не фактичні дані? Я не міг …

32 r  t-test 

Багато написано про вибір кольорових сліпих кольорів для карт, полігонів та затінених регіонів загалом (див., Наприклад, http://colorbrewer2.org ). Мені не вдалося знайти рекомендації щодо кольорів ліній та різної товщини ліній для графіків ліній. Цілі: легко розрізнити лінії, навіть коли вони переплітаються лінії легко відрізнити окремими особами з найпоширенішими формами кольорової …

31 r  data-visualization 

Які корисні вказівки щодо тестування параметрів (наприклад, глибина взаємодії, мінчік, частота вибірки тощо) з використанням ГБМ? Скажімо, у мене 70-100 функцій, кількість населення 200 000, і я маю намір перевірити взаємодію на глибині 3 і 4. Ясно, що мені потрібно зробити кілька тестувань, щоб побачити, яка комбінація параметрів найкраще вибірка. …

31 r  hypothesis-testing  cart  boosting