我收集了有关船舶海鸟干扰的数据。我在船上带着测距双筒望远镜和角度板。对于我调查的每只鸟,我都有一个相对于船舶航向的起始距离和方位。我还有鸟做出反应(或在某些情况下没有反应)的距离和方位。
我想制作一个两面板图,其中一个显示起始距离和轴承位置,另一个显示终止距离和轴承。理想情况下,第二个图将进行颜色编码(或 pch 编码)以显示不同的反应类型。
我的数据是这种格式
date_id dist bear act
550 40711_027 200 30 f
551 40711_028 500 45 n
552 40711_028 450 60 n
553 40711_028 400 75 n
554 40711_028 371 80 f
555 40711_029 200 5 f
556 40711_030 200 10 d
557 40711_031 400 30 n
558 40711_031 350 30 d
Run Code Online (Sandbox Code Playgroud)
这是您可以使用的格式的数据
id <- c(1,2,2,2,2,3,4,5,5)
dist <- c(200,500,450,400,371,200,200,400,350)
bear <- c(30,45,60,75,80,5,10,30,30)
act <- c("f","n","n","n","f","f","d","n","d")
dat <- data.frame(id, dist, bear, act)
Run Code Online (Sandbox Code Playgroud)
正如您所看到的,有些 id 是重复的,有些则只有一行。我想在一个地块上绘制第一个距离和熊,在另一个地块上绘制最后一个距离和熊(每个 ID)。对于只有一次观察的鸟类来说,这些可能是相同的。最好根据“act”列对第二个图中的点进行颜色编码。此外,方位没有左侧或右侧的指定,所以我可以接受所有点都位于中线的一侧或另一侧,但如果您知道将它们随机放置在中心线的左侧或右侧会很酷。理想情况下,绘图看起来像这样。

更新:遵循@jbaums 的建议,使用他在此处 …
我在R中找到了plotrix包但是还不知道如何在R中做这个简单的圆.基本上,我怎样才能做半径为1和0的极坐标:360度角度,产生一个圆?
例
$$ r\COS \左(\压裂{2\PI} {3} \左(\压裂{3\THETA} {2\PI} - \左\ lfloor \压裂{3\THETA} {2\PI}\right\rfloor\right) - \frac {\ pi} {3}\right)= 1 $$
也许相关