R语言Circular包计算海龟孵化巢平均方向结果异常求助

我在处理濒危海龟孵化幼体的定向数据，需要为每个巢穴计算平均方向，但部分巢穴的计算结果存在异常（并非所有巢穴都有问题）。目前不确定问题出在以下哪个环节：

• Circular包的平均计算逻辑
• 按巢穴标签子集化数据的方式
• 绘图方向校正的zero=pi/2参数

操作流程 #

1. 从主数据表circdata中提取方位数据创建circular对象cdata（因Circular包不支持主表中的非数值型数据）
2. 绘制目标巢穴的数据点
3. 计算该巢穴的圆形平均方向并以箭头标注

代码片段 #

# 创建circular对象
cdata<- circular(circdata$Diff.360, type="directions", units="degrees", template="none",rotation = ("clock"),modulo = ("2pi") ,zero = pi/2)

# 绘制目标巢穴A401的数据点
plot(cdata [circdata$Nest.Tag == "A401"], cex=0.9, bin=180, stack=TRUE, sep=0.07, shrink=1.5)

# 计算并标注平均方向
df3mean<- mean.circular(cdata [circdata$Nest.Tag == "A401"], na.rm = TRUE) 
arrows.circular(df3mean,lwd = 3)

疑问点 #

通过主数据集circdata的巢穴标签对cdata进行子集化计算均值，这种方式是否会引发问题？

异常表现 #

图示中每个数据点代表一只幼体的定向，但Circular包计算出的圆形平均方向箭头与实际数据点位置明显不符。

内容的提问来源于stack exchange，提问作者a.wright