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

Чи є спосіб у R (вбудована функція) обчислити матрицю переходу ланцюга Маркова з набору спостережень? Наприклад, взявши набір даних, як описано нижче, та обчислити матрицю переходу першого порядку? dat<-data.frame(replicate(20,sample(c("A", "B", "C","D"), size = 100, replace=TRUE)))

29 r  markov-process 

Я досить новачок у статистиці, і мені потрібна ваша допомога. У мене невеликий зразок: H4U 0.269 0.357 0.2 0.221 0.275 0.277 0.253 0.127 0.246 Я провів тест Шапіро-Вілка за допомогою R: shapiro.test(precisionH4U$H4U) і я отримав такий результат: W = 0.9502, p-value = 0.6921 Тепер, якщо я припускаю, що рівень значущості …

29 r  distributions  interpretation  goodness-of-fit  normality-assumption 

Я використовую функцію auto.arima () у пакеті прогнозів, щоб підходити до моделей ARMAX з різними коваріатами. Однак мені часто доводиться вибирати велику кількість змінних і, як правило, закінчують остаточну модель, яка працює з їх підмножиною. Мені не подобаються спеціальні методи для вибору змінних, тому що я людина і підданий упередженості, …

29 r  time-series  lasso  regularization  elastic-net 

Я завжди використовую lm()в R для виконання лінійної регресії на . Ця функція повертає коефіцієнт такий, щоуyyхxxββ\betaу= βх .y=βx.y = \beta x. Сьогодні я дізнався про загальні найменші квадрати і цю princomp()функцію (аналіз основних компонентів, PCA) можна використовувати для її виконання. Це має бути добре для мене (точніше). Я зробив …

29 r  pca  least-squares  deming-regression  total-least-squares 

Ми можемо використовувати lm()для прогнозування значення, але нам все ж потрібне рівняння формули результату в деяких випадках. Наприклад, додайте рівняння до графіків.

29 r  regression  lm 

Зачинено. Це питання поза темою . Наразі відповіді не приймаються. Хочете вдосконалити це питання? Оновіть питання, щоб воно було тематичним для перехресної перевірки. Закрито 2 роки тому . Я використовую caret, щоб запустити перехрещений випадковий ліс над набором даних. Змінна Y - фактор. У моєму наборі даних немає NaN, Inf …

29 r  random-forest  caret  regression  prediction  fitting  social-science  poisson-distribution  distributions  characteristic-function  bayesian  prior  regression  normal-distribution  interaction  nonparametric  skewness  svm  standard-deviation  standard-error  regression-coefficients  igraph  natural-language  word2vec  word-embeddings  regression  machine-learning  sampling  r  regression  machine-learning  random-forest  ensemble  sampling  unbiased-estimator  proof  estimators  mse  probability  conditional-probability  bayes  anova  missing-data  neural-networks  recommender-system  r  confidence-interval  sample  multiple-imputation  r  time-series  forecasting  mase 

Я використовую пакет robustbase, щоб виконати оцінку glm . Однак коли я це роблю, я отримую таку помилку: Error in solve.default(crossprod(X, DiagB * X)/nobs, EEq) : system is computationally singular: reciprocal condition number = 1.66807e-16 Що це означає / вказує? І як я можу це налагодити? PS. Якщо вам потрібно …

29 r  generalized-linear-model  robust 

Я використовую моделі R (3.1.1) та ARIMA для прогнозування. Я хотів би знати, яким повинен бути параметр "частота", який призначається у ts()функції , якщо я використовую дані часових рядів, які є: розділяється на хвилини і поширюється на 180 днів (1440 хвилин / день) розділена на секунди і поширюється на 180 …

28 r  time-series  arima  multiple-seasonalities  mape 

Я використовую glmnet для розрахунку оцінок регресії хребта. Я отримав певні результати, які змусили мене підозріти в тому, що glmnet справді робить те, що, на мою думку, це робить. Щоб перевірити це, я написав простий R-скрипт, де я порівнюю результат регресії хребта, виконаного рішенням, і той, який є в glmnet, …

28 r  ridge-regression  glmnet 

Я хотів би отримати p-значення та розмір ефекту незалежної категоріальної змінної (з декількома рівнями) - тобто "загальний", а не для кожного рівня окремо, як це нормальний вихід lme4у R. Це просто так про що люди повідомляють під час роботи програми ANOVA. Як я можу це отримати?

28 r  hypothesis-testing  anova  mixed-model  lme4-nlme 

[Початкова назва "Вимірювання подібності для ієрархічних дерев кластеризації" згодом @ttnphns змінила, щоб краще відобразити тему] Я виконую ряд ієрархічних кластерних аналізів у кадрі даних пацієнтів (наприклад, подібний до http://www.biomedcentral.com/1471-2105/5/126/figure/F1?highres=y ) Я експериментую з різними мірами відстані , різною вагою параметрів та різними ієрархічними методами , щоб зрозуміти їх вплив на …

28 r  clustering  distance-functions  similarities  dendrogram 

Заблокований . Це питання та його відповіді заблоковано, оскільки це питання поза темою, але має історичне значення. Наразі не приймає нових відповідей чи взаємодій. Я розумію, що це питання досить широке, але мені цікаво, що має бути вирішальним моментом у вирішенні питання про створення (чи ні) нового пакету для Р. …

28 r  software 

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

28 r  mixed-model  interpretation  panel-data 

Як приклад, розглянемо ChickWeightнабір даних у Р. Дисперсія очевидно зростає з часом, тому якщо я використовую просту лінійну регресію, наприклад: m <- lm(weight ~ Time*Diet, data=ChickWeight) Мої запитання: Які аспекти моделі будуть сумнівними? Чи обмежуються проблеми екстраполяцією поза Timeдіапазоном? Наскільки толерантною є лінійна регресія до порушення цього припущення (тобто, якою …

28 r  regression  heteroscedasticity  assumptions 

Я використовую пакет "boot" для обчислення приблизного двостороннього завантаженого p-значення, але результат занадто далекий від p-значення використання t.test. Я не можу зрозуміти, що я зробив не так у своєму R-коді. Може хтось, будь ласка, підкаже мені про це time = c(14,18,11,13,18,17,21,9,16,17,14,15, 12,12,14,13,6,18,14,16,10,7,15,10) group=c(rep(1:2, each=12)) sleep = data.frame(time, group) require(boot) diff …

28 r  hypothesis-testing  p-value  bootstrap  permutation-test