Я анализирую результаты сельскохозяйственного эксперимента следующим образом:
Сюжеты располагаются в рандомизированном полном блочном дизайне с 3 повторениями. В течение десятилетий ими управляли либо обычным, либо органическим образом (M) , и культуры (C) чередуются между кукурузой или томатом, которые в течение каждого года содержат 4 комбинации управления урожаем * , Я взял три образца почвы с каждого участка (сс) . Я делал это в течение двух лет (Y) . Места отбора проб на участке не были одинаковыми за два года.
Блокировка устраняет сильную изменчивость текстуры почвы.
Подводя итог:
- 3 блока (B)
- 2 системы управления (M)
- 2 урожая (C)
- 3 подвыборки (сс)
Я делаю миллион различных анализов на этих почвах, и меня интересует сравнение изменчивость из-за управления с изменчивостью из-за других факторов (включая блок) или их взаимодействия с управлением для различных испытаний почвы.
В некоторых местах я запутался
Должен ли он быть сплит-сюжетом во времени? Я измеряю одни и те же графики два года подряд. Но это в разных местах, и они довольно большие.
Каковы соответствующие условия ошибки? Смотрите ниже, что я имею сейчас. Я думаю, что Блок дает изменчивость из-за блока, Year*Crop*Management*Block дает изменчивость, связанную с конкретным c графиком? Должен ли я включить какие-либо другие взаимодействия с блоком? Они все немного интересны, но они очень быстро меняют мои степени свободы.
Я получаю довольно разные результаты, если использую Type3 или reml в качестве метода оценки. Я использую тип 3, потому что в моей текущей модели параметры ковариации обычно равны 0 или меньше, а также потому, что мне интересно видеть Type3. Это хорошее оправдание?
Пожалуйста, смотрите мой текущий код SAS ниже
Proc mixed data=soil covtest method=type3;
class Block Year Mgt Crop;
model Nitrogen=Year|Mgt|Crop/outp=pred ddfm=satterth;
random Block Year*Mgt*Crop*Block;
lsmeans Year Mgt Crop
Year*Mgt Year*Crop Mgt*Crop Year*Mgt*Crop
/ slice= (Year Mgt Crop) diff adjust=tukey cl;
ods exclude diffs;
ods output diffs=dif lsmeans=lsm;
title 'FDA';
Run;