R 如何使用缺失值执行双向重复测量方差分析

对于我的数据集，我需要执行某种双因素重复测量ANOVA。我有一个被称为“治疗”的受试因素和一个被称为“频率”的受试因素，共有8个等级。我的问题是，我的大多数受试者对所有8个频率级别（缺失值）都没有反应，称为“阈值”。此外，我的两种治疗方法也不平衡（第一种治疗方法约23种，另一种治疗方法约21种）

你建议我试试什么r码？那代码会是什么样子？我一直在研究aov和Anova（汽车套件）功能。我还需要找出我的模型的公式。我在想

aov（阈值（治疗频率）+错误（受试者/（治疗频率））

但我一直收到类似“在aov（……）中，error（）模型是单数的”这样的错误消息

我这里的问题是，如果你只在错误术语中包含主题内因素，那么错误（主题/频率）或仅仅是错误（主题），或者我是否正确地包含了所有内容？此外，我是否应该为每只鸟的每个频率包含一行响应，即使我没有特定的值？我是否应该将NA放在缺少的值单元格中，或者删除该级别的整行数据

任何帮助都将不胜感激！我不熟悉更高级的统计和建模，所以请记住！如果我需要澄清或添加任何内容，请告诉我。谢谢