Nested / SplitModel - Rep LặpMeasure / MixedModel ANOVA: mức độ lồng và kịch bản trong R

Tập dữ liệu của tôi có các biến sau:

• Điều trị (4 loại - cố định)
• Vị trí (8 vị trí cố định)
• Vị trí trong Vị trí (3 vị trí trên mỗi vị trí - cố định)
• Các mẫu được lấy ở mỗi vị trí (3 mẫu cho mỗi vị trí ngẫu nhiên)
  • Thời gian (hai lần lấy mẫu - cố định)
  • Tỷ lệ khoáng hóa (kết quả phân tích mẫu đã lấy)

Hai địa điểm được sử dụng để kiểm tra mỗi lần điều trị (nghĩa là 4 địa điểm điều trị, 2 địa điểm mỗi lần điều trị, tổng cộng 8 địa điểm).

Tôi muốn thực hiện một biện pháp lặp lại phân tách (/ lồng nhau?) (/ Mô hình hỗn hợp?) ANOVA trong Rviệc sử dụng các biến trên.

Q.1. Liệu âm thanh này có phù hợp?

Mục tiêu của tôi là để xem liệu có ảnh hưởng của 1) vị trí, 2) điều trị, 3) thời gian và 4) tương tác của tất cả (ví dụ: pos * Treat *, pos * time, Treat * time, pos * Treat * time) trên tỷ lệ khoáng hóa.

Q 2. Vị trí được lồng trong điều trị? Là mẫu lồng nhau ở vị trí?

Câu 3. Các yếu tố giữa và trong là gì?

Câu 4. Chủ đề / cốt truyện là gì? - Đó có phải là vị trí hoặc vị trí hoặc mẫu hoặc tỷ lệ?

Câu hỏi 5. Làm thế nào tôi có thể đặt thời gian như các biện pháp lặp lại trong công thức R của mình?

Câu 6. Bạn sẽ sử dụng aov, lme hoặc ezANOVA?

Câu 7. Làm cách nào để mã hóa các biến độc lập riêng biệt và các tương tác của chúng thành một công thức R thích hợp?

Tôi thực sự đã cố gắng để tìm ra điều này trong nhiều ngày và dường như tôi không thể tìm thấy một câu trả lời có ý nghĩa ...

Vấn đề rắc rối! Là vị trí cố định hay ngẫu nhiên? Là vị trí cố định hay ngẫu nhiên? Tôi cho rằng mẫu đó là ngẫu nhiên.

• Vì điều trị được chỉ định cho vị trí, vị trí là đơn vị lấy mẫu. Về cơ bản, việc so sánh giữa các phương pháp điều trị được thực hiện ở cấp độ đó. .$n=8$
• Đơn vị đo lường là quan sát bạn thực hiện trên "mẫu" của mình tại một thời điểm nhất định.
• Vị trí không được lồng trong điều trị. Việc điều trị được áp dụng cho vị trí.
• Vị trí được lồng bên trong vị trí.
• Mẫu được lồng vào vị trí bên trong.
• Thời gian được lồng bên trong Mẫu.
• Thời gian được vượt qua với điều trị.

Bạn có 3 cấp độ lồng nhau (thời gian trong mẫu, mẫu trong vị trí, vị trí trong vị trí).

Nếu vị trí, vị trí và mẫu là ngẫu nhiên, tôi nghĩ công thức R sẽ như thế này:

 Y ~ Treatment * Time +(1|location|position|sample)

Bạn có 1 hàng trong khung dữ liệu của mình cho mỗi lần quan sát mẫu mỗi lần - với mã hóa phù hợp cho tất cả các đặc điểm thiết kế của bạn.

Nó sẽ hoạt động để kết hợp các biện pháp lặp đi lặp lại thành một số điểm như trung bình hoặc sự khác biệt của chúng? Điều đó có thể làm cho mô hình dễ giải thích hơn.