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

В R, якщо пишу lm(a ~ b + c + b*c) це все-таки буде лінійною регресією? Як зробити інші види регресії в R? Буду вдячний за будь-яку рекомендацію щодо підручників чи навчальних посібників?

11 r  regression 

Заблокований . Це запитання та його відповіді заблоковано, оскільки це питання поза темою, але має історичне значення. Наразі не приймає нових відповідей чи взаємодій. Мені потрібно побудувати боксплот без будь-яких осей і додати його до поточного сюжету (крива ROC), але мені потрібно додати більше текстової інформації до boxplot: мітки для …

11 r  boxplot 

Я підходив до моделі з декількома незалежними змінними, одна з яких - відставання залежної змінної, використовуючи пакет dynlm. Якщо припустити, що я маю прогнози на 1 крок вперед для моїх незалежних змінних, як я можу отримати прогнози на 1 крок вперед для моїх залежних змінних? Ось приклад: library(dynlm) y<-arima.sim(model=list(ar=c(.9)),n=10) #Create …

11 r  dynamic-regression 

Ось, що я хочу зробити, але, здається, немає predictметоду для mlogit. Будь-які ідеї? library(mlogit) data("Fishing", package = "mlogit") Fish <- mlogit.data(Fishing, varying = c(2:9), shape = "wide", choice = "mode") Fish_fit<-Fish[-1,] Fish_test<-Fish[1,] m <- mlogit(mode ~price+ catch | income, data = Fish_fit) predict(m,newdata=Fish_test)

11 r  logistic  logit  multinomial 

Скажімо, у мене є дані: x1 <- rnorm(100,2,10) x2 <- rnorm(100,2,10) y <- x1+x2+x1*x2+rnorm(100,1,2) dat <- data.frame(y=y,x1=x1,x2=x2) res <- lm(y~x1*x2,data=dat) summary(res) Я хочу побудувати графік безперервної безперервної взаємодії таким чином, що x1 знаходиться на осі X, а x2 представлений 3 рядками, один з яких представляє x2 при Z-балі 0, один …

11 r  regression  ggplot2  interaction 

Я шукаю, як (візуально) пояснити просту лінійну кореляцію для студентів першого курсу. Класичним способом візуалізації було б дати графік розсіяння Y ~ X прямою регресійною лінією. Нещодавно мені прийшла ідея розширити цей тип графіки, додавши до сюжету ще 3 зображення, залишивши мене з: графік розкидання y ~ 1, потім y …

11 r  data-visualization  regression  correlation 

Отже, припускаючи, що є сенс перевірки припущення про нормальність для anova (див. 1 і 2 ) Як це можна перевірити на R? Я би сподівався зробити щось на кшталт: ## From Venables and Ripley (2002) p.165. utils::data(npk, package="MASS") npk.aovE <- aov(yield ~ N*P*K + Error(block), npk) residuals(npk.aovE) qqnorm(residuals(npk.aov)) Що не …

11 r  anova  normality-assumption  repeated-measures 

Кожен, хто слідкує за бейсболом, ймовірно, чув про нестандартне виконання MVP типу Жозе Баутіста в Торонто. За чотири роки тому він забивав приблизно 15 домашніх пробіжок за сезон. Минулого року він потрапив до 54, кількість перевершила лише 12 гравців в історії бейсболу. У 2010 році йому виплатили 2,4 мільйона, і …

11 r  regression  modeling 

У мене є такий кадр даних: case simulation temp plank oxygen 1 1 1 8 7 11 2 2 1 16 10 15 ... 17 17 2 26 12 17 18 18 2 15 8 12 19 19 2 28 11 21 20 20 2 24 6 14 Я хотів …

11 r 

Я хочу реалізувати (в R) наступну дуже просту динамічну лінійну модель, для якої у мене є 2 невідомі зміни параметрів часу (дисперсія помилки спостереження та дисперсія помилки стану ).ϵ1тϵт1\epsilon^1_tϵ2тϵт2\epsilon^2_t Yтθt + 1==θт+ ϵ1тθт+ ϵ2тYт=θт+ϵт1θт+1=θт+ϵт2 \begin{matrix} Y_t & = & \theta_t + \epsilon^1_t\\ \theta_{t+1} & = & \theta_{t}+\epsilon^2_t \end{matrix} Я хочу …

11 r  mcmc  dlm  particle-filter 

У мене є дані за два роки, які в основному виглядають так: Дата _ __ Насильство Y / N? _ Кількість пацієнтів 01.01.2008 _ ___ 0 __ _ __ _ ____ 11 1/1/2008 _ __ _ 0 _ __ _ __ _ __ 11 3/1/2008 _ ____ 1 __ _ …

11 r  mixed-model  autocorrelation  panel-data 

Приклад коду: (pc.cr <- princomp(USArrests)) summary(pc.cr) loadings(pc.cr) ## note that blank entries are small but not zero Я отримую різні результати від кожного, і я не впевнений, що розумію, в чому різниця. Ось вихід: > summary(pc.cr) Importance of components: Comp.1 Comp.2 Comp.3 Comp.4 Standard deviation 82.8908472 14.06956001 6.424204055 2.4578367034 Proportion …

11 r  pca 

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

11 r  correlation  categorical-data  data-visualization 

Зачинено. Це питання поза темою . Наразі відповіді не приймаються. Хочете вдосконалити це питання? Оновіть питання, щоб воно було тематичним для перехресної перевірки. Закрито 2 роки тому . Я хочу обчислити підсумок змінної у data.frame для кожної унікальної комбінації факторів у data.frame. Чи варто використовувати для цього plyr? Я все …

11 r  categorical-data  aggregation  plyr 

Ризикуючи поставити запитання специфічним для програмного забезпечення та з виправданням його всюдисущості та ідіосинкразії, я хочу запитати про функцію biplot()в R та, точніше, про обчислення та побудову графіку за замовчуванням, накладених червоними стрілками, відповідними до основних змінних. [Щоб мати сенс у деяких коментарях, спочатку розміщені сюжети мали надзвичайну проблему з …

11 r  pca  biplot