Lý lịch
Tôi có dữ liệu từ một nghiên cứu thực địa, trong đó có bốn cấp độ điều trị và sáu lần lặp lại trong mỗi hai khối. (4x6x2 = 48 quan sát)
Các khối cách nhau khoảng 1 dặm, và trong các khối, có một lưới gồm các ô 42, 2m x 4m và một lối đi rộng 1m; nghiên cứu của tôi chỉ sử dụng 24 lô trong mỗi khối.
Tôi muốn đánh giá hiệp phương sai không gian.
Dưới đây là một phân tích ví dụ sử dụng dữ liệu từ một khối duy nhất, mà không tính đến hiệp phương sai không gian. Trong tập dữ liệu, plot
là id của biểu đồ, x
là vị trí x và vị trí y
y của mỗi ô với ô 1 tập trung vào 0, 0. level
là mức xử lý và response
là biến trả lời.
layout <- structure(list(plot = c(1L, 3L, 5L, 7L, 8L, 11L, 12L, 15L, 16L,
17L, 18L, 22L, 23L, 26L, 28L, 30L, 31L, 32L, 35L, 36L, 37L, 39L,
40L, 42L), level = c(0L, 10L, 1L, 4L, 10L, 0L, 4L, 10L, 0L, 4L,
0L, 1L, 0L, 10L, 1L, 10L, 4L, 4L, 1L, 1L, 1L, 0L, 10L, 4L), response = c(5.93,
5.16, 5.42, 5.11, 5.46, 5.44, 5.78, 5.44, 5.15, 5.16, 5.17, 5.82,
5.75, 4.48, 5.25, 5.49, 4.74, 4.09, 5.93, 5.91, 5.15, 4.5, 4.82,
5.84), x = c(0, 0, 0, 3, 3, 3, 3, 6, 6, 6, 6, 9, 9, 12, 12, 12,
15, 15, 15, 15, 18, 18, 18, 18), y = c(0, 10, 20, 0, 5, 20, 25,
10, 15, 20, 25, 15, 20, 0, 15, 25, 0, 5, 20, 25, 0, 10, 20,
25)), .Names = c("plot", "level", "response", "x", "y"), row.names = c(NA,
-24L), class = "data.frame")
model <- lm(response ~ level, data = layout)
summary(model)
Câu hỏi
- Làm thế nào tôi có thể tính toán một ma trận hiệp phương sai và đưa nó vào hồi quy của mình?
- Các khối rất khác nhau, và có sự tương tác mạnh mẽ * tương tác khối. Có thích hợp để phân tích chúng một cách riêng biệt?