```
library(tidyverse)
library(circular)
#>
#> Attaching package: 'circular'
#> The following objects are masked from 'package:stats':
#>
#> sd, var
x_prex_wtt <- data.frame(c(57,88,90,100,25,90,46,77,38,35,29,23,27,13,18,55,33,55 ))
y_prex_wtt <- data.frame(c(2,24,210,180,145,1,88,30,170,177,16,44,315,300,320,189,15,177))
xCircular_wtt <- circular(x_prex_wtt, units = "degrees", template = "geographics")
yCircular_wtt <- circular(y_prex_wtt, units = "degrees", template = "geographics")
watson.two.test(xCircular_wtt,yCircular_wtt)
#>
#> Watson's Two-Sample Test of Homogeneity
#>
#> Test Statistic: 0.3839
#> 0.001 < P-value < 0.01
#>
```

Instead of getting a range for where the p-value is, I need a specific p-value.

Any suggestions?

Thanks,

Bro