Я хочу использовать библиотечную функцию PoweR pvalueM C p-значение для проверки соответствия (GOF), но я не понимаю параметр null. law.index . Может быть, кто-то должен был использовать что-то подобное? Сталкивались ли вы с чем-нибудь подобным? Я прочитал статью Пьера Лафе де Мишо «PoweR: воспроизводимый инструмент исследования, облегчающий исследования по моделированию мощности Монте-Карло для тестов на пригодность для испытаний в R». Я знаю, что значение параметра null.law.index получается функцией getindex . Функция getindex сообщает:
Index Law Nbparams Default1 Default2 Default3 Default4
1 1 Laplace(mu,b) 2 0.000000 1 NA NA
2 2 Normal(mu,sigma) 2 0.000000 1 NA NA
3 3 Cauchy(mu,sigma) 2 0.000000 1 NA NA
4 4 Logistic(mu,sigma) 2 0.000000 1 NA NA
5 5 Gamma(shape,rate) 2 2.000000 1 NA NA
6 6 Beta(a,b) 2 1.000000 1 NA NA
7 7 Uniform(a,b) 2 0.000000 1 NA NA
8 8 Student-t(df) 1 1.000000 NA NA NA
9 9 Chi-squared(df) 1 1.000000 NA NA NA
10 10 Lognormal(logmean,logsd) 2 0.000000 1 NA NA
11 11 Weibull(shape,scale) 2 1.000000 1 NA NA
12 12 ShiftedExp(l,rate) 2 0.000000 1 NA NA
13 13 U^{j+1} 1 1.000000 NA NA NA
14 14 AveUnif(k,a,b) 3 2.000000 0 1.0 NA
15 15 UUnif(j) 1 1.000000 NA NA NA
16 16 VUnif(j) 1 1.000000 NA NA NA
17 17 JSU(mu,sigma,nu,tau) 4 0.000000 1 0.0 5e-01
18 18 Tukey(l) 1 1.000000 NA NA NA
19 19 LoConN(p,m) 2 0.200000 3 NA NA
20 20 JSB(g,d) 2 0.000000 1 NA NA
21 21 SkewN(xi,omega,alpha) 3 0.000000 1 0.0 NA
22 22 ScConN(p,d) 2 0.200000 2 NA NA
23 23 GP(mu,sigma,xi) 3 0.000000 1 0.0 NA
24 24 GED(mu,sigma,p) 3 0.000000 1 1.0 NA
25 25 Stable(alpha,beta,c,mu) 4 1.000000 0 1.0 0e+00
26 26 Gumbel(mu,sigma) 2 1.000000 1 NA NA
27 27 Frechet(mu,sigma,alpha) 3 0.000000 1 1.0 NA
28 28 GEV(mu,sigma,xi) 3 0.000000 1 0.0 NA
29 29 GArcSine(alpha) 1 0.500000 NA NA NA
30 30 FoldN(mu,sigma) 2 0.000000 1 NA NA
31 31 MixN(p,m,d) 3 0.500000 0 1.0 NA
32 32 TruncN(a,b) 2 0.000000 1 NA NA
33 33 Nout(a) 1 1.000000 NA NA NA
34 34 GEP(t1,t2,t3,crit) 4 0.500000 0 1.0 1e-06
35 35 Exponential(lambda) 1 1.000000 NA NA NA
36 36 ALaplace(mu,b,k) 3 0.000000 1 2.0 NA
37 37 NIG(alpha,beta,delta,mu) 4 1.000000 0 1.0 0e+00
Index Stat Alter Nbparams
1 1 K-S 3 0
2 2 AD^* 3 0
3 3 Z_C 3 0
4 4 Z_A 3 0
5 5 P_S 3 0
6 6 K^2 3 0
7 7 JB 3 0
8 8 DH 3 0
9 9 RJB 3 0
10 10 T_{Lmom} 3 0
11 11 T_{Lmom}^{(1)} 3 0
12 12 T_{Lmom}^{(2)} 3 0
13 13 T_{Lmom}^{(3)} 3 0
14 14 BM_{3-4} 3 0
15 15 BM_{3-6} 3 0
16 16 T_{MC-LR} 3 0
17 17 T_w 0,1,2 0
18 18 T_{MC-LR}-T_w 3 0
19 19 T_{S,5} 3 0
20 20 T_{K,5} 3 0
21 21 W 4 0
22 22 W' 4 0
23 23 tilde{W} 4 0
24 24 D 0,1,2 0
25 25 r 4 0
26 26 CS 3 0
27 27 Q 0,1,2 0
28 28 Q-Q* 0,1,2 0
29 29 BCMR 3 0
30 30 beta_3^2 3 0
31 31 T^*(alpha) 1 1
32 32 I_n 3 0
33 33 R_{sJ} 3 0
34 34 Q* 0,1,2 0
35 35 R_n 3 0
36 36 X_{APD} 3 0
37 37 Z_{EPD} 0,1,2 0
Цель моего исследования - проверить правильность гипотез соответствия с помощью различных тестов. Я попытался сделать это, используя:
statistikos <- statcompute(stat.index =3, data = X, level = 0.05, alter =3)
stat <- statistikos$statistic
pval <-statistikos$pvalue
Но с тестами 3-5, 18-20, 23-32, 34, 35 pval = NA. Затем я попытался использовать:
pval<-pvalueMC(X, stat.index=3, null.law.index =2, M=1000, alter =3)
Мой главный вопрос: каково значение параметра добротности параметра null.law.index 3-5, 18-20, 23-32, 34, 35 ( Нумерация библиотеки PoweR) тесты?