数据拟合检验属于包装柯西或冯米塞斯

Question

我的工作中有错误数据，我需要测试它属于冯·米塞斯或包装的柯西发行版。我用 watson.test 来表示 von Mises 和 mwraped cauchy test 来表示 wraped cauchy distribution。 但是在我的工作中，这个按摩师出现在R中的问题找不到功能"mwrapped"，"watson"。 例如，我们测试ε

epsilon<-c( 2.49723058,-0.3950603,-0.5350603 ,-1.899642  , 1.18661
, -1.773387,-2.93964, 2.16494 , -0.8042236  , -0.2996419,
0.98035807 , -0.9119328  ,2.238904   ,-1.333387 ,-0.13193,
0.443486, -0.5242236  ,2.368067  ,0.400358  ,-1.748805,
-1.437969,0.1497405  ,2.12661  , 1.248067 ,0.055776,
0.0126489 , 1.48807 ,  0.69432  ,1.37745 ,-0.78881,
0.11578 ,1.97723 , 1.76494  ,  2.109521  ,0.1849397,
-1.39796863 ,-1.754841  ,1.525159  ,-1.936514  ,-1.579642,
-0.276514 , 2.16494,  1.100358 ,-1.18255  ,0.772868,
1.460577 ,-0.6994228 ,-0.125678,  -1.127351  ,-0.9550603,
-1.127351,  0.3220314    , 1.08515885  ,  0.74057718 ,-0.05797,
-1.139642,-0.231933,-1.567351 ,   0.49723058 ,   2.42974053,
-0.08713198 , -1.151714 ,-0.9917137   ,-1.268805 ,-0.25506 ,
-0.25506,  1.49723058  ,-2.557969 ,-2.407132   , -2.43171366,
-0.367132 ,-2.041096  ,-2.23506  ,1.60494  ,  1.60494,
2.746613,-0.207132 ,-1.14255  ,1.39578 ,-0.7950603,
0.2818122, -0.8904786  ,  0.49723058 ,0.7126489 ,-1.1691,
-1.29339 ,-2.14277,-0.5458969  ,-0.3950603 , 1.60494,
1.60494,2.820358, 0.5434855,1.30036 ,-1.299642,
-0.1150603)
watson.test(epsilon, alpha=0, dist= "vonmises")
mwrappedg("cauchy",epsilon,starts=c(1,1),K=10,method="BFGS")

Answer 1

像watson.test和mwrappedg这样的函数不适用于 base R，这就是它显示这些错误的原因。对于watson.test，您必须使用circularR 包，对于mwrappedg，您必须使用 R 包Wrapped。首先，您必须安装它们，然后您可以使用

install.packages(c("circular","Wrapped")
library(circular)
watson.test(x, alpha=0.05, dist="vonmises")
library(Wrapped)
mwrappedg("cauchy", epsilon, starts=c(1,1), K=10, method="BFGS")