Às 18:04 de 08/01/2023, Nick Wray escreveu:
Hello I'm using the Rao Spacing test from the circular package on a series
of (circular) data vectors but I can't find a way of simply recording
whether the null is rejected or not for each vector in the series - the
test gives me an individual result but unlike with other tests in R (say
MannKendall or pettitt etc) I can't find a way of delving into the guts of
the output to extract either1)both the test statistic and critical value
or 2)the status of ejection or not of the null eg as T/F
I''ve tried various things including str(rao) and rao[[1]]
as.character(rao) etc
rao<-rao.spacing.test(pi_data,alpha=0.05) ## my input
Rao's Spacing Test of Uniformity
Test Statistic = 169.6952
Level 0.05 critical value = 155.19
Reject null hypothesis of uniformity-- can't break this down
Des anyone know how I could do this?
Thanks, Nick Wray
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
Hello,
This test returns a non-standard object, meaning, not an object of class
"htest" and its print method creates and invisibly returns an object
with the value you want.
library(circular)
#>
#> Attaching package: 'circular'
#> The following objects are masked from 'package:stats':
#>
#> sd, var
# data
x <- circular(runif(200, 0, 2*pi))
# run the test
test_result <- rao.spacing.test(x, alpha=0.1)
# save the print method result
print_result <- print(test_result)
#>
#> Rao's Spacing Test of Uniformity
#>
#> Test Statistic = 135.2864
#> Level 0.1 critical value = 140.06
#> Do not reject null hypothesis of uniformity
print_result$accepted
#> [1] TRUE
If you run
getAnywhere("print.rao.spacing.test")
you will see that the printing is made by a series of cat instructions
so to get the TRUE/FALSE above you will need to print the result like above.
______________________________________________
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
