Фон
У меня есть данные полевого исследования, в котором есть четыре уровня обработки и шесть повторностей в каждом из двух блоков. (4x6x2 = 48 наблюдений)
Блоки находятся примерно в 1 миле друг от друга, а внутри блоков есть сетка из 42, 2 х 4 м участков и дорожка шириной 1 м; мое исследование использовало только 24 участка в каждом блоке.
Я хотел бы оценить оценку пространственной ковариации.
Вот пример анализа с использованием данных из одного блока без учета пространственной ковариации. В наборе данных plot
это идентификатор графика, x
это местоположение x и местоположение y
y каждого графика с графиком 1, центрированным на 0, 0. level
это уровень обработки и response
переменная отклика.
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)
Вопросов
- Как я могу вычислить ковариационную матрицу и включить ее в свою регрессию?
- Блоки сильно отличаются друг от друга, и есть сильные взаимодействия с блоком лечения. Уместно ли их анализировать отдельно?