Câu hỏi được gắn thẻ «r»

Tôi đã được hỏi một câu hỏi liên quan đến một mô hình tuyến tính được thực hiện với R's lm: "Hồi quy có sử dụng bình phương tối thiểu lặp tuyến tính hay không tuyến tính không?" Tôi đã tìm kiếm một chút và [nghĩ rằng tôi] hiểu sự …

9 r  linear-model  lm 

Giả sử trong phân tích hồi quy trong R, tôi có một biến độc lập loại yếu tố với 3 cấp độ trong tập dữ liệu đào tạo của tôi. Nhưng trong tập dữ liệu thử nghiệm, biến nhân tố tương tự có 5 cấp độ. Do đó, tôi không …

9 r  regression  categorical-data  cross-validation  rare-events 

Dữ liệu thô của tôi bao gồm chuỗi thời gian 60 ngày với xu hướng giảm. Dữ liệu là hàng tuần nên tần số được đặt thành 7. Tôi đã tính toán sự khác biệt của dữ liệu trông như thế này Khi tôi chạy các âm mưu ACF và …

9 r  time-series  autocorrelation 

Khi chọn số lượng nút thắt thích hợp cho GAM, người ta có thể muốn tính đến số lượng dữ liệu và gia số trên trục x. Điều gì xảy ra nếu chúng ta có 100 gia số trên trục x với 1000 điểm dữ liệu ở mỗi mức tăng. …

9 r  gam  splines 

Tôi có ba biến kinh tế vĩ mô (ICS - tâm lý người tiêu dùng, ER - tỷ lệ việc làm, DGO - đơn hàng lâu bền) và đã thực hiện các thử nghiệm quan hệ nhân quả Granger trong R về chúng. Tôi thực sự không biết làm thế …

9 r  granger-causality 

Tôi gặp một số khó khăn khi có được kết quả tương đương giữa aovmô hình đo lặp lại giữa bên trong và lmermô hình hỗn hợp. Dữ liệu và tập lệnh của tôi trông như sau data=read.csv("https://www.dropbox.com/s/zgle45tpyv5t781/fitness.csv?dl=1") data$id=factor(data$id) data id FITNESS TEST PULSE 1 1 pilates CYCLING 91 2 …

9 r  anova  mixed-model  repeated-measures  lme4-nlme 

Tôi muốn tính toán tỷ lệ mới mắc để trình bày theo tỷ lệ rủi ro để trình bày các biện pháp rủi ro tương đối và tuyệt đối. Tôi đã thấy trong các nghiên cứu khác rằng tỷ lệ mới mắc như vậy có thể được tính bằng cách …

9 r  survival  cox-model  poisson-regression  incidence-rate-ratio 

library(lme4) out <- glmer(cbind(incidence, size - incidence) ~ period + (1 | herd), data = cbpp, family = binomial, contrasts = list(period = "contr.sum")) summary(out) Fixed effects: Estimate Std. Error z value Pr(>|z|) (Intercept) -2.32337 0.22129 -10.499 < 2e-16 *** period1 0.92498 0.18330 5.046 4.51e-07 *** period2 -0.06698 0.22845 -0.293 …

9 r  binomial  logit  lme4-nlme 

Theo những gì tôi hiểu, đối số cp cho rparthàm giúp cắt tỉa trước cây theo cách tương tự như đối số minsplit hoặc minbucket. Điều tôi không hiểu là cách tính giá trị CP. Ví dụ df<-data.frame(x=c(1,2,3,3,3,4), y=as.factor(c(TRUE, TRUE, FALSE, TRUE, FALSE, FALSE)), method="class") mytree<-rpart(y ~ x, data = …

9 r  cart  rpart 

Tôi đang vật lộn với các điểm chiếu trong phân tích phân biệt tuyến tính (LDA). Nhiều cuốn sách về các phương pháp thống kê đa biến minh họa ý tưởng của LDA với hình dưới đây. Mô tả vấn đề như sau. Trước tiên, chúng ta cần vẽ ranh …

9 r  machine-learning  classification  discriminant-analysis 

Đây là dữ liệu: > tires <- data.frame(Wear = c(17, 14, 12, 13, 14, 14, 12, 11, 13, 13, 10, 11, 13, 8, 9, 9), Brand = rep(LETTERS[1:4], 4), Car = as.character(as.roman(rep(1:4, each = 4)))) > tires Wear Brand Car 1 17 A I 2 14 B I 3 12 …

9 r  anova  dataset  degrees-of-freedom 

Tôi thực hiện hồi quy tuyến tính bằng hàm R lm: x = log(errors) plot(x,y) lm.result = lm(formula = y ~ x) abline(lm.result, col="blue") # showing the "fit" in blue nhưng nó không phù hợp Thật không may, tôi không thể hiểu ý nghĩa của hướng dẫn. Ai đó có …

9 r  regression 

Có thể vẽ sơ đồ phụ thuộc một phần để hiển thị xác suất lớp và ước tính tác động của bộ dự đoán cho mô hình GBM không? Một cái gì đó tương tự partialPlottừ randomForestgói. Theo bài viết này , một phần cốt truyện có thể thực hiện …

9 r  gbm 

Tôi có dữ liệu cho các vụ tai nạn xe máy theo giờ trong ngày. Như bạn mong đợi, chúng ở mức cao vào giữa ngày và cao điểm vào giờ cao điểm. Geom_d mật độ mặc định của ggplot2 làm mịn nó ra Một tập hợp con của dữ …

9 r  pdf  kernel-smoothing  splines 

Nói rằng chúng tôi có các dữ liệu sau: set.seed(123) data <- data.frame(x = c(rnorm(50, 1, 1), rnorm(50, 5, 2)), y = c(rep('A', 50), rep('B', 50))) Mà mang lại boxplot sau ( boxplot(data$x ~ data$y)): Bây giờ hãy nói rằng tôi muốn kiểm tra xem hai mẫu có cùng tham …

9 r  nonparametric  wilcoxon-mann-whitney