Hi,
I have the binned data (observed and generated from model) that I would like
to
test using the anderson-darling goodness of fit test. But I'm not sure which
package in R to use.
I tried ad.test(...) but it does not recognise the test by Vito Ricci in
FITTING
DISTRIBUTIONS WITH R
> ad.test(hist_hume_beec[,1],hist_hume_beec[,2])
Error: could not find function "ad.test"
I also tried stat.ad from surv2.ks {surv2sample} but not really sure how to use
it.
Here are the data:
rel_obs rel_gen
[1,] 0.000000000 0.0000000000
[2,] 0.052576236 0.0282122554
[3,] 0.085173502 0.0752419321
[4,] 0.080967403 0.0986985317
[5,] 0.091482650 0.1060457970
[6,] 0.092534175 0.1035785945
[7,] 0.073606730 0.0956528954
[8,] 0.077812829 0.0851286202
[9,] 0.078864353 0.0738023237
[10,] 0.060988433 0.0627435684
[11,] 0.057833859 0.0525393532
[12,] 0.037854890 0.0434658038
[13,] 0.044164038 0.0356059397
[14,] 0.031545741 0.0289287619
[15,] 0.030494217 0.0233411924
[16,] 0.026288118 0.0187212879
[17,] 0.027339642 0.0149387493
[18,] 0.014721346 0.0118669624
[19,] 0.008412198 0.0093895074
[20,] 0.010515247 0.0074031461
[21,] 0.001051525 0.0058186459
[22,] 0.005257624 0.0045603385
[23,] 0.003154574 0.0035649996
[24,] 0.003154574 0.0027804220
[25,] 0.003154574 0.0021639094
[26,] 0.001051525 0.0016808219
[27,] 0.000000000 0.0013032428
[28,] 0.000000000 0.0010087536
[29,] 0.000000000 0.0007790357
[30,] 0.000000000 0.0005975463
[31,] 0.000000000 0.0004469915
Thank you for any help given.
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
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.
