你不需要使用一个主题之间的因素
ranova
。你的设计有两个试因素,这是很好。
我在代码中看到两个小问题。首先,你用“条件”试设计的名字,其中一个变量的名字在设计。可能不是一个好主意。下面,我改变了设计的名称
withinDesign
。其次,您需要定义
条件
和
一边
作为分类变量。我添加了两行,和有一个清洁运行代码的一些随机的测试数据。见下文。
%测试数据
MoS_AP_side1_CWS_MSZP =兰德(20日1);
MoS_AP_side2_CWS_MSZP =兰德(20日1);
MoS_AP_side1_FWS_MSZP =兰德(20日1);
MoS_AP_side2_FWS_MSZP =兰德(20日1);
MoS_AP_side1_CWS2_MSZP =兰德(20日1);
MoS_AP_side2_CWS2_MSZP =兰德(20日1);
t_APMoS_MSZP =表(MoS_AP_side1_CWS_MSZP MoS_AP_side2_CWS_MSZP MoS_AP_side1_FWS_MSZP,…
MoS_AP_side2_FWS_MSZP、MoS_AP_side1_CWS2_MSZP MoS_AP_side2_CWS2_MSZP…
,“VariableNames”,{“CWSside1”,“CWSside2”,“FWSside1”,“FWSside2”,“CWS2side1”,“CWS2side2”});
withinDesign =表([1 1 2 2 3 3],[1 2 1 2 1 2]”,“VariableNames”,{“条件”“一边”});
withinDesign。条件= categorical(withinDesign.Condition);
withinDesign。一边= categorical(withinDesign.Side);
rm_APMoS_MSZP = fitrm (t_APMoS_MSZP,“CWSside1-CWS2side2 ~ 1”,“WithinDesign”,withinDesign);
[b1_APMoS, A_APMoS C_APMoS D_APMoS] = ranova (rm_APMoS_MSZP,“WithinModel”,“条件*边”);
