Tôi đã đo hai thông số (Carbon hữu cơ hòa tan DOC = y và xả = x). Khi hai biến này được vẽ với nhau, chúng ta sẽ có một vòng lặp trễ (xem ví dụ mã và hình ảnh).
Bây giờ, để phân tích sâu hơn, tôi muốn xác định khu vực của vòng lặp cuồng loạn này. Tôi đã tìm ra rằng điều này có thể được thực hiện bằng phương pháp phóng phi tiêu Monte Carlo. Phương pháp này nói rằng diện tích của một khu vực không xác định tỷ lệ với diện tích của một hình chữ nhật đã biết nhân với số lần truy cập trong trường bên trong (vòng lặp).
Vấn đề của tôi bây giờ là, làm thế nào để giải quyết vấn đề bên trong / bên ngoài bằng R. Làm cách nào tôi có thể vẽ một hình chữ nhật với một khu vực đã biết và làm thế nào tôi có thể vượt trội các lần truy cập ngẫu nhiên bên trong và bên ngoài vòng lặp cuồng loạn?
Xin lưu ý rằng tôi mở cho bất kỳ phương pháp nào khác ...
Tôi googled, và tìm kiếm các trang web thống kê khác nhau nhưng không thể tìm thấy câu trả lời. Bất kỳ trợ giúp trực tiếp hoặc liên kết đến các trang web / bài viết khác được đánh giá rất cao.
Data <- read.table("http://dl.dropbox.com/u/2108381/DOC_Q_hystersis.txt", sep = ";",
header = T)
head(Data)
plot(Data$Q, Data$DOC, type = "o", xlab = "Discharge (m3 s-1)", ylab = "DOC (mg C l-1)",
main = "Hystersis loop of the C/Q relationship")