当前位置: 首页 > 知识库问答 >
问题:

如何计算循环数据的均方差

东门楚
2023-03-14

我遵循了此处提出的计算循环数据平均值的建议:

https://en.wikipedia.org/wiki/Mean_of_circular_quantities

但是我也想计算均方差。

#A vector of directional data (separated by 20 degrees each)
Dir2<-c(350,20,40)

#Degrees to Radians
D2R<-0.0174532925

#Radians to Degrees
Rad2<-Dir2 * D2R


Sin2<-sin(Rad2)
SinAvg<-mean(Sin2)

Cos2<-cos(Rad2)
CosAvg<-mean(Cos2)

RADAVG<-atan2(SinAvg, CosAvg)
DirAvg<-RADAVG * R2D

上面给出了平均值,但我不知道如何计算SD

我试着只取正弦和余弦的标准偏差的平均值,但我得到了不同的答案。

SinSD<-sd(Sin2)
CosSD<-sd(Cos2)
mean(CosSD, SinSD)

共有1个答案

亢建木
2023-03-14

您可以使用圆形软件包:

x <- circular(Rad2)
mean(x)
# Circular Data: 
# Type = angles 
# Units = radians 
# Template = none 
# Modulo = asis 
# Zero = 0 
# Rotation = counter 
# [1] 0.2928188 # The same as yours
sd(x)
# [1] 0.3615802

手动,

sqrt(-2 * log(sqrt(sum(Sin2)^2 + sum(Cos2)^2) / length(Rad2)))
# [1] 0.3615802

这可以从sd的源代码中看到。圆形

请参见此处和此处。

 类似资料:
  • 正如标题所示,我需要做的是使用for循环计算数组的平均值。我所说的平均值是指

  • 本文向大家介绍PHP计算加权平均数的方法,包括了PHP计算加权平均数的方法的使用技巧和注意事项,需要的朋友参考一下 本文实例讲述了PHP计算加权平均数的方法。分享给大家供大家参考。具体如下: 希望本文所述对大家的php程序设计有所帮助。

  • 问题内容: 任何人都知道如何计算这些列之一的平均值(在Linux上)? 例如:mean(第2栏) 问题答案: Awk: 读为: 对于每一行,将第2列添加到变量“总计”中。 在文件末尾,打印“总计”除以记录数。

  • 问题内容: 好的。所以我有这张桌子: |项目| 数量| 价格| |苹果| 2 | 2.00 | |橙色| 3 | 1.50 | |葡萄| 5 | 2.50 | 我想显示客户必须支付的 总计 。怎么做? 我真的不知道如何使用数组。谁能告诉我如何? 我的代码(某种) 使用此查询,价格显示在每一行中: 总计应显示在最后一行(总计= $ 21.00)。 谢谢你的帮助。 问题答案:

  • 我有一个简短的问题 我有一个数据帧,有许多测量列。我想计算具有相同(标题)名称的列的平均值。我使用下面的代码(在stackoverflow中找到)。。 如何计算数据帧中具有相同列名的列的平均值 作为示例数据... 结果是这样的。。。 这段代码告诉我具有相同(标题)名称的列的含义。 但是我也想要均方差。我试着用rowSds替换rowMeans,但是不起作用。 知道如何使用相同的代码来计算均方差吗??

  • 问题内容: 在下面的示例代码中,是否真的需要counter = 0,还是有更好,更多的Python方法来访问循环计数器?我看到了一些与循环计数器有关的PEP,但它们要么被延迟要么被拒绝(PEP 212 和PEP 281)。 这是我的问题的简化示例。在我的实际应用程序中,这是通过图形完成的,整个菜单必须每帧重新绘制一次。但这以易于复制的简单文本方式进行了演示。 也许我还应该补充一点,我正在使用Pyt