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

Một trong những vấn đề tôi luôn gặp phải với các mô hình hỗn hợp là tìm ra trực quan hóa dữ liệu - loại có thể kết thúc trên giấy hoặc áp phích - một khi có kết quả. Ngay bây giờ, tôi đang làm việc trên một mô …

15 r  data-visualization  mixed-model 

Tôi đang tham gia một khóa học sau đại học về Thống kê ứng dụng sử dụng sách giáo khoa sau (để cho bạn cảm nhận về mức độ của tài liệu được bảo hiểm): Khái niệm và phương pháp thống kê , của GK Bhattacharyya và RA Johnson. Giáo …

15 r  sas  java 

Tôi đang thực hiện hồi quy tỷ lệ rủi ro Cox trong R bằng cách sử dụng coxph, bao gồm nhiều biến. Phần dư Martingale trông rất tuyệt, và phần dư Schoenfeld là tuyệt vời cho tất cả các biến. Có ba biến số mà phần dư Schoenfeld không bằng …

15 r  model  cox-model  hazard  schoenfeld-residuals 

Giả sử chúng ta có X 2 ~ unif ( n , 0 , 1 ) ,X1∼unif(n,0,1),X1∼unif(n,0,1),X_1 \sim \textrm{unif}(n,0,1), X2∼unif(n,0,1),X2∼unif(n,0,1),X_2 \sim \textrm{unif}(n,0,1), trong đó là mẫu ngẫu nhiên đồng nhất có kích thước n vàunif(n,0,1)unif(n,0,1)\textrm{unif}(n,0,1) Y=X1,Y=X1,Y=X_1, Z=0.4X1+1−0.4−−−−−−√X2.Z=0.4X1+1−0.4X2.Z = 0.4 X_1 + \sqrt{1 - 0.4}X_2. Khi đó tương quan giữa và …

15 r  correlation  random-generation  uniform 

Tôi đang sử dụng rlm trong gói R MASS để hồi quy mô hình tuyến tính đa biến. Nó hoạt động tốt đối với một số mẫu nhưng tôi nhận được hệ số gần đúng cho một mô hình cụ thể: Call: rlm(formula = Y ~ X1 + X2 + …

15 r  multiple-regression  robust 

Thông qua trang web này, gần đây tôi đã phát hiện ra Sơ đồ Sankey, một cách tuyệt vời để hình dung những gì đang xảy ra trong biểu đồ dòng truyền thống . Dưới đây là một ví dụ điển hình về Sơ đồ Sankey của George M. Whitesides …

15 r  data-visualization  sankey-diagram 

Tôi có một câu hỏi liên quan đến việc sử dụng biến nhóm trong mô hình phi tuyến tính. Vì hàm nls () không cho phép các biến nhân tố, tôi đã phải vật lộn để tìm hiểu xem liệu người ta có thể kiểm tra ảnh hưởng của một …

15 r  mixed-model  nls 

Đã khóa . Câu hỏi này và câu trả lời của nó bị khóa vì câu hỏi không có chủ đề nhưng có ý nghĩa lịch sử. Nó hiện không chấp nhận câu trả lời hoặc tương tác mới. Tôi đang gặp vấn đề trong khi thực hiện một số …

15 r 

Làm thế nào để tính toán thống kê R bình phương ( ) trong R cho và / hoặc đầu ra hàm? Ví dụ cho dữ liệu này:r2r2r^2loesspredict cars.lo <- loess(dist ~ speed, cars) cars.lp <- predict(cars.lo, data.frame(speed = seq(5, 30, 1)), se = TRUE) cars.lpcó hai mảng fitcho mô …

15 r  r-squared  loess 

CARET sẽ tự động sử dụng lưới điều chỉnh được chỉ định trước để xây dựng các mô hình khác nhau trước khi chọn một mô hình cuối cùng, sau đó đào tạo mô hình cuối cùng trên dữ liệu đào tạo đầy đủ. Tôi có thể cung cấp lưới …

15 r  caret 

Tôi đang phân tích dữ liệu từ một thí nghiệm giai thừa không cân bằng cả với SASvà R. Cả hai SASvà Rcung cấp tổng bình phương loại I tương tự nhưng tổng bình phương loại III của chúng khác nhau. Dưới đây là SASvà Rmã và đầu ra. DATA …

15 r  anova  sas  sums-of-squares 

Lâu lắm rồi tôi mới nhìn phân vùng cây. Lần trước tôi đã làm điều này, tôi thích tiệc tùng trong R (được tạo bởi Hothorn). Ý tưởng về suy luận có điều kiện thông qua lấy mẫu có ý nghĩa với tôi. Nhưng rpart cũng đã kháng cáo. Trong …

15 r  cart  rpart  partitioning 

Đã có một số câu hỏi hay và bộ câu trả lời về sách giới thiệu hoặc phương pháp tiếp cận để học R, ví dụ ở đây và đây . Nhưng tôi có một vấn đề hơi khác - cách tốt nhất để chạy một phiên dài một giờ …

15 r  teaching 

Đã khóa . Câu hỏi này và câu trả lời của nó bị khóa vì câu hỏi không có chủ đề nhưng có ý nghĩa lịch sử. Nó hiện không chấp nhận câu trả lời hoặc tương tác mới. Tôi có biểu đồ như thế này: Mã R để tạo …

15 r  data-visualization  ggplot2 

Giả sử tôi có một thí nghiệm nhân đôi 2x2 đơn giản mà tôi muốn thực hiện ANOVA. Như thế này, ví dụ: d <- data.frame(a=factor(sample(c('a1','a2'), 100, rep=T)), b=factor(sample(c('b1','b2'), 100, rep=T))); d$y <- as.numeric(d$a)*rnorm(100, mean=.75, sd=1) + as.numeric(d$b)*rnorm(100, mean=1.2, sd=1) + as.numeric(d$a)*as.numeric(d$b)*rnorm(100, mean=.5, sd=1) + rnorm(100); Trong trường hợp không …

15 r  anova  contrasts