我有一个气候模型运行数据集。它们目前存储在如下列表中：
$ ensemble :List of 25
..$ run_name :List of 2
.. ..$ variable: num [1:72, 1:36, 1:12, 1:40] 255 256 256 257 257 ...
variable是特定模型输出，例如'表面温度'，尺寸为[lat, long, month, year]（不要问我为什么输出不仅仅是按月...）
这不一定是存储这些数据的最佳方式，我想知道是否有一种R-ish方法可以使操作更容易。特别是，我想看看每个变量的年平均值，包括所有在内部和整体的运行（即每个整体/变量一个图，25行），以及时间序列上每个集合的统计数据（即移动PDF） ），以后可能会更复杂。
理想情况下，我想避免for循环，而是使用*apply函数。我一直在尝试使用这种结构，但仍然需要在lapply()调用中组合两个函数，这不起作用。