Hi,
I'm trying to carry out Cramer von Mises tests between pairs of vectors
belonging to a discrete distribution (concretely frequencies from 0 to 200).
However, the program crashes in the attempt. The problem seems to be that these
vectors only have positive integer numbers (+ zero). When I add a random very
small positive decimal to the non-decimal part everything works fine (files
prm1 & prpmr1). I attach two of these vectors so you can run the following
code. I've also thought to divide both vectors by a real constant such as pi.
Do you think these two approaches are acceptable?
setwd("")
require(CvM2SL2Test)
prm = scan('prm.txt')
prpmr = scan('prpmr.txt')
ct1 = cvmts.test(prm, prpmr) # here R crashes
ct1
cvmts.pval( ct1, length(prm), length(prpmr) )
Thank you for your help,
Santi199.09
193.59
199.99
173.89
179.99
200.89
198.89
200.09
186.59
171.79
118.79
44.19
155.79
200.49
201.29
199.99
201.09
201.19
200.39
200.59
201.19
200.79
201.09
201.29
200.79
200.99
201.19
200.99
201.49
201.39
200.99
199.99
201.29
201.39
201.19
164.89
200.79
200.29
201.49
200.99
198.99
198.29
200.09
201.09
194.29
189.49
170.29
173.99
23.99
200.39
200.49
200.79
151.19
201.09
131.79
25.39
0.29
0.69
1.39
170.49
200.59
201.89
200.89
201.19
201.29
200.99
197.89
200.89
185.09
166.69
172.59
185.99
201.59
201.59
184.49
200.99
201.99
200.19
200.69
201.19
201.89
197.29
201.29
200.39
201.69
200.39
140.99
200.69
200.49
201.69
197.39
198.79
201.09
200.39
201.09
200.79
201.39
200.89
201.29
201.39
201.49
201.29
150.39
178.29
199.29
201.49
200.69
201.29
200.69
182.29
200.99
201.19
200.99
201.59
136.69
200.59
200.89
201.89
200.39
201.49
184.29
185.49
187.09
124.29
187.19
64.19
65.19
201.39
201.09
201.39
201.89
175.39
200.19
200.79
200.99
200.19
201.19
201.19
201.19
201.59
200.69
200.79
201.29
201.69
193.49
197.99
118.19
115.89
95.29
92.89
201.29
200.99
201.09
201.59
201.49
193.99
179.59
198.39
139.09
36.09
0.99
1.49
141.39
174.89
169.19
128.59
3.59
18.19
0.89
0.69
0.49
1.39
59.89
1.69
135.09
186.99
185.49
190.19
19.79
5.79
1.19
1.19
1.49
102.09
200.79
201.39
201.09
193.19
201.49
180.59
0.69
189.09
98.79
123.69
119.09
148.19
179.69
185.29
199.09
195.59
197.99
200.09
172.19
198.59
199.99
164.69
201.39
200.59
198.39
175.69
200.59
192.09
143.49
142.39
196.79
191.69
200.99
200.79
200.89
193.39
193.39
201.59
201.39
200.39
200.89
200.79
195.39
200.49
200.99
200.69
201.49
201.49
200.59
201.29
200.59
200.79
200.79
201.69
200.99
201.39
201.09
201.39
200.49
201.49
201.09
201.29
200.39
197.39
201.59
196.69
201.89
200.79
21.89
10.49
86.09
12.49
5.09
25.19
76.19
200.59
196.19
74.69
77.79
8.29
177.19
183.89
173.79
150.69
41.69
186.69
200.59
200.89
197.09
5.09
0.89
3.09
5.69
1.69
1.59
40.89
201.59
200.79
201.39
201.09
201.39
201.39
200.29
140.59
72.79
201.49
201.39
200.39
201.39
192.59
200.19
200.79
201.09
200.39
201.49
162.29
200.29
201.79
186.49
200.89
200.69
191.69
192.49
183.49
112.09
8.69
1.09
200.39
171.89
21.19
201.09
201.29
201.09
200.59
195.69
187.69
191.59
181.79
94.59
0.69
201.89
195.99
200.69
201.49
200.99
193.79
191.69
200.69
200.99
199.89
200.79
200.49
200.59
200.99
192.09
166.59
186.49
178.29
201.29
200.79
192.59
4.59
0.89
1.69
200.89
168.99
200.89
200.69
196.29
200.69
200.79
201.29
201.19
201.09
200.79
201.09
200.89
201.49
201.79
201.09
198.69
201.59
200.39
201.09
200.79
200.69
201.09
201.19
165.69
200.59
162.89
145.49
124.29
168.89
201.29
201.49
201.49
200.69
201.89
94.39
2.09
95.99
0.49
0.69
1.09
171.19
18.59
0.89
1.39
1.39
189.09
200.39
27.49
201.59
200.89
67.49
184.99
160.89
53.39
1.09
1.19
0.89
0.39
1.59
201.09
200.49
166.99
119.59
95.99
64.49
0.29
1.19
1.29
2.39
106.89
42.39
201.19
200.99
200.69
200.09
201.09
198.19
201.29
200.89
200.09
200.29
200.99
93.69
177.09
196.39
129.79
143.49
19.89
64.39
200.89
201.09
153.49
-0.01
1.19
1.99
189.09
188.79
198.39
190.49
175.29
197.39
201.09
71.69
0.59
0.99
1.39
0.69
40.59
3.39
200.29
169.09
19.59
0.79
0.99
0.59
1.79
1.69
1.29
1.19
0.99
1.59
0.09
0.29
165.99
190.69
135.89
199.89
200.69
12.89
0.59
1.19
127.29
200.89
200.99
201.29
201.89
200.69
200.59
200.89
197.99
200.59
175.79
1.09
180.79
1.29
129.49
193.99
184.89
161.69
177.19
177.99
0.99
0.59
0.39
1.09
176.29
163.99
186.29
140.49
183.99
172.49
178.79
200.99
201.09
201.29
201.39
200.79
201.49
200.29
200.59
201.09
200.69
201.09
200.99
200.69
200.79
201.29
200.69
200.69
194.49
140.69
32.89
197.99
106.69
161.89
69.09
0.99
0.99
1.39
0.99
1.39
200.79
201.39
194.19
167.69
201.09
200.29
200.69
199.99
200.39
200.89
197.79
201.59
1.59
0.39
0.89
0.99
0.19
0.49
0.89
0.59
0.69
1.19
0.89
0.79
1.39
0.69
0.29
0.29
1.69
0.49
1.59
0.69
0.99
0.39
0.99
0.69
1.19
0.69
-0.01
0.09
0.69
0.49
1.59
1.09
0.69
1.49
1.59
1.29
0.99
0.39
0.59
1.89
1.09
0.19
0.49
0.49
1.19
0.49
0.59
0.99
1.29
1.09
1.49
0.89
0.99
1.69
0.69
1.59
1.29
0.69
0.09
0.99
0.39
0.29
0.39
1.79
0.59
0.19
1.49
1.29
1.29
0.99
1.29
1.09
1.19
1.49
1.39
1.29
1.39
1.39
1.59
1.89
0.39
1.09
1.69
1.39
0.79
0.99
0.39
0.99
0.79
1.89
1.09
1.09
0.69
1.09
0.19
1.19
0.69
0.79
1.29
1.29
0.79
1.19
1.29
1.29
0.89
0.99
0.99
0.59
1.39
0.09
1.39
0.29
0.49
1.09
1.09
1.19
1.09
-0.01
0.99
0.69
1.29
1.39
0.99
1.09
1.59
0.89
0.59
1.29
0.79
0.89
0.69
1.09
1.29
1.59
1.29
1.39
1.39
0.89
0.19
1.79
0.49
0.89
1.69
0.49
0.09
0.29
1.59
1.49
0.49
1.19
0.59
0.99
1.29
1.69
0.99
1.89
1.49
1.19
1.49
0.59
0.79
0.79
0.09
1.39
0.99
1.39
1.69
0.29
0.99
1.39
0.89
0.89
0.69
1.89
1.09
1.89
0.29
1.19
1.69
1.29
0.69
0.99
0.79
0.69
0.19
0.89
1.39
0.89
0.49
0.29
0.79
0.39
0.09
1.29
0.99
1.39
1.19
0.29
1.19
1.29
1.09
0.59
1.49
0.29
1.19
1.49
0.99
1.09
0.99
0.99
1.19
0.99
0.79
1.39
0.69
0.69
1.19
1.29
1.39
0.59
0.09
1.09
0.99
1.29
0.39
0.19
1.29
1.39
0.29
1.69
0.89
0.49
0.79
1.39
0.49
0.99
1.29
0.99
1.29
0.29
198.00045
193.00025
200.0006
173.00035
179.0004
200.00045
198.00055
199.0004
185.00035
171.00035
118.00015
44.00035
154.00025
200.0002
200.0008
199.00025
200.00035
200.00065
200.00055
200.00065
200.00065
200.0008
200.0007
200.0006
200.0003
200.0006
200.00015
200.00045
200.0005
200.00045
200.0004
200.00035
200.00055
200.0006
200.0009
163.001
200.00065
200.0007
200.00015
200.0003
198.0004
197.00065
200.00035
200.00045
193.0006
189.00045
170.0003
173.00085
23.0005
200.0007
200.0006
200.0007
150.0008
200.00085
131.00015
25.00065
0.00075
0.0008
0.00065
169.0005
200.00025
200.0004
200.00065
200.0004
200.0006
200.00085
197.0004
200.0005
184.00095
166.00025
172.00055
185.00075
200.0005
200.00055
184.0006
200.0001
200.00065
200.00085
200.00065
200.0007
200.0002
197.0005
200.00045
200.0005
200.0005
200.00025
140.0005
200.0004
200.00065
200.0006
197.00075
198.0007
200.0004
200.00065
200.0004
200.0005
200.0003
200.0001
200.00055
200.00035
200.0002
200.00045
150.0004
177.00065
198.0004
200.0006
200.0002
200.0006
200.0007
181.00055
200.0003
200.0003
200.0004
200.0002
136.00075
200.00085
200.00035
200.00065
200.0008
200.00055
183.00065
185.0007
186.0009
123.00035
186.00045
63.00035
64.00025
200.00045
200.0005
200.0006
200.0006
174.0007
200.0003
200.00045
200.0007
200.0007
200.0003
200.0006
200.00075
200.00015
200.00045
200.00065
200.0006
200.0004
193.00045
197.0007
117.0006
115.0004
94.00035
92.00045
200.0006
200.00095
200.0006
200.00065
200.00085
193.00075
179.0007
198.00025
138.00035
35.00065
0.0005
0.00065
141.00065
174.00055
168.0001
127.0007
3.00085
17.00045
0.0008
0.00065
0.0003
0.00055
59.00075
0.00035
134.00025
186.0004
184.00095
190.0005
18.0007
5.0005
0.00075
0.0003
0.00085
101.0006
200.0002
200.0006
200.0006
192.00065
200.00045
180.00055
0.0008
188.0007
98.00065
123.00055
118.0007
147.00065
179.00045
184.00075
198.00055
194.00035
197.00095
199.0007
171.00045
198.00095
199.0008
163.0004
200.00025
199.00045
198.00055
175.00055
200.00055
191.00025
142.0006
141.0007
196.00055
191.0003
200.00045
200.0002
200.0005
192.0004
193.00065
200.00065
200.0006
200.00025
200.00055
200.0005
194.00075
200.0006
200.00035
200.00055
200.00095
200.0006
200.00085
200.00055
200.0002
200.0005
200.0009
200.00045
200.0007
200.0005
200.00045
200.0006
200.0007
200.0006
200.00055
200.0006
200.0006
196.0005
200.00035
196.0006
200.00025
200.00035
21.00055
9.0009
85.0008
11.0006
4.0003
25.00055
76.0006
200.00045
195.0002
74.0008
77.00055
7.0008
176.00055
183.00085
173.0009
150.00035
41.00055
186.0003
200.00045
200.00065
196.0004
4.0007
0.00055
2.0004
4.00055
0.00075
0.0004
40.00075
200.0003
200.00065
200.00015
200.0004
200.00055
200.0007
200.0007
139.0006
72.00045
200.0005
200.0006
200.0006
200.0003
192.00075
200.0005
200.00055
200.0006
200.0007
200.0003
161.0003
200.0009
200.0008
186.00065
200.00035
200.0004
191.00055
191.0003
183.00075
111.00045
7.00065
0.00035
200.00025
171.0004
20.0007
200.00065
200.00075
200.00055
200.0009
194.0004
187.0003
191.0006
181.00045
94.00065
0.00025
200.0002
195.00075
199.00035
200.0006
200.0007
193.00045
191.00025
200.00085
200.00025
199.0005
200.0006
199.00045
200.0004
200.00055
192.00045
165.0007
186.00055
177.00085
200.0005
200.00065
192.00085
3.0006
0.0007
0.00055
200.00085
168.0005
200.001
200.0008
195.0006
199.00035
200.0007
200.0008
200.0006
200.0006
200.0003
200.0003
200.0009
200.0003
200.00045
200.0002
197.00075
200.00055
200.0008
200.0006
200.00055
200.00075
200.0005
200.00025
165.0007
199.00085
162.00075
145.00055
123.00055
168.00045
200.00065
200.0002
200.00025
200.0007
200.0003
93.00025
1.0005
95.0008
0.0004
0.00095
0.0003
171.00075
18.00035
0.0003
0.00065
0.0005
188.0006
200.00025
26.00085
200.00085
200.0008
67.0004
184.0005
160.00085
53.0005
0.00055
0.00045
0.0008
0.00015
0.00045
200.00055
200.0001
166.0006
119.0006
95.0006
63.0009
0.00055
0.0002
0.00035
1.00035
106.0007
41.0001
200.00065
200.00055
200.0006
200.00035
200.0007
198.00045
200.0008
200.00075
200.0006
200.00035
200.0006
93.0006
176.00045
195.0005
128.00065
142.00095
19.00065
63.0007
200.00045
200.0006
152.00045
0.0002
0.0005
1.0008
188.00055
188.00065
198.0005
189.0007
174.0008
196.0002
200.0003
70.0007
0.0002
0.0005
0.0002
0.00015
40.0004
2.0008
200.0007
168.00045
19.00025
0.001
0.00035
0.0009
0.0007
0.00065
0.00035
0.00065
0.0006
0.00025
0.00095
0.00065
165.0004
190.00035
135.00035
199.0003
200.0004
12.0001
0.0002
0.0004
126.00075
200.00055
200.0007
200.00055
200.0004
200.00035
200.0005
200.00095
197.00055
200.0006
175.0005
0.0006
179.00075
0.00045
129.0005
193.00065
184.0005
161.0008
176.0008
177.0005
0.00065
0.00035
0.00055
0.00025
176.00025
163.0005
186.0009
140.00025
183.00055
171.0006
178.0005
200.001
200.0004
200.0005
200.00025
200.0007
200.0003
200.00055
200.00025
200.0005
200.0005
200.0004
200.00045
200.00055
200.00055
200.0009
200.0002
200.00075
194.00075
139.00065
31.0005
197.0004
106.00055
161.00065
68.0006
0.0006
0.00035
0.0007
0.0005
0.00045
200.0007
200.0005
193.00025
167.0009
200.001
200.00065
200.00085
199.0001
200.00095
200.00065
196.00035
200.00055
______________________________________________
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.
