Commit 8fec8d3f1c6c134359e9bdbda2a40cfdf466b124 - libpdl-stats-perl

Update upstream source from tag 'upstream/0.81' Update to upstream version '0.81' with Debian dir 4bad8075b948d1101ad96ac68cf431b00d78b04b Bas Couwenberg 2 years ago

10 changed file(s) with 834 addition(s) and 705 deletion(s). Raw diff Collapse all Expand all

-0

Changes less more

	0	0.81 2022-01-11
	1	- fix Kmeans NiceSlice for PDL 2.066
	2
0	3	0.80 2021-10-02
1	4	- add COV option to GLM::logistic - thanks David Myers for suggestion (in 2011)
2	5	- make all PLOT options default to off; only load PGPLOT when try to use

+214

-187

GENERATED/PDL/Stats/Basic.pm less more

2	2	#
3	3	package PDL::Stats::Basic;
4	4
5		our @EXPORT_OK = qw( binomial_test rtable which_id PDL::PP stdv PDL::PP stdv_unbiased PDL::PP var PDL::PP var_unbiased PDL::PP se PDL::PP ss PDL::PP skew PDL::PP skew_unbiased PDL::PP kurt PDL::PP kurt_unbiased PDL::PP cov PDL::PP cov_table PDL::PP corr PDL::PP corr_table PDL::PP t_corr PDL::PP n_pair PDL::PP corr_dev PDL::PP t_test PDL::PP t_test_nev PDL::PP t_test_paired );
	5	our @EXPORT_OK = qw(binomial_test rtable which_id stdv stdv_unbiased var var_unbiased se ss skew skew_unbiased kurt kurt_unbiased cov cov_table corr corr_table t_corr n_pair corr_dev t_test t_test_nev t_test_paired );
6	6	our %EXPORT_TAGS = (Func=>\@EXPORT_OK);
7	7
8	8	use PDL::Core;

19	19
20	20
21	21
	22
	23	#line 5 "stats_basic.pd"
	24
22	25	use PDL::LiteF;
23	26	use PDL::NiceSlice;
24	27	use Carp;

49	52	my $stdv = stdv( $data );
50	53
51	54	=cut
52
	55	#line 57 "Basic.pm"
53	56
54	57
55	58

63	66
64	67
65	68
	69	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
	70
66	71
67	72	=head2 stdv
68	73

86	91
87	92
88	93	=cut
89
90
91
92
93
94
	94	#line 96 "Basic.pm"
	95
	96
	97
	98	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
95	99	*stdv = \&PDL::stdv;
96
97
98
	100	#line 102 "Basic.pm"
	101
	102
	103
	104	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
99	105
100	106
101	107	=head2 stdv_unbiased

120	126
121	127
122	128	=cut
123
124
125
126
127
128
	129	#line 131 "Basic.pm"
	130
	131
	132
	133	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
129	134	*stdv_unbiased = \&PDL::stdv_unbiased;
130
131
132
	135	#line 137 "Basic.pm"
	136
	137
	138
	139	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
133	140
134	141
135	142	=head2 var

154	161
155	162
156	163	=cut
157
158
159
160
161
162
	164	#line 166 "Basic.pm"
	165
	166
	167
	168	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
163	169	*var = \&PDL::var;
164
165
166
	170	#line 172 "Basic.pm"
	171
	172
	173
	174	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
167	175
168	176
169	177	=head2 var_unbiased

188	196
189	197
190	198	=cut
191
192
193
194
195
196
	199	#line 201 "Basic.pm"
	200
	201
	202
	203	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
197	204	*var_unbiased = \&PDL::var_unbiased;
198
199
200
	205	#line 207 "Basic.pm"
	206
	207
	208
	209	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
201	210
202	211
203	212	=head2 se

229	238
230	239
231	240	=cut
232
233
234
235
236
237
	241	#line 243 "Basic.pm"
	242
	243
	244
	245	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
238	246	*se = \&PDL::se;
239
240
241
	247	#line 249 "Basic.pm"
	248
	249
	250
	251	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
242	252
243	253
244	254	=head2 ss

263	273
264	274
265	275	=cut
266
267
268
269
270
271
	276	#line 278 "Basic.pm"
	277
	278
	279
	280	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
272	281	*ss = \&PDL::ss;
273
274
275
	282	#line 284 "Basic.pm"
	283
	284
	285
	286	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
276	287
277	288
278	289	=head2 skew

297	308
298	309
299	310	=cut
300
301
302
303
304
305
	311	#line 313 "Basic.pm"
	312
	313
	314
	315	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
306	316	*skew = \&PDL::skew;
307
308
309
	317	#line 319 "Basic.pm"
	318
	319
	320
	321	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
310	322
311	323
312	324	=head2 skew_unbiased

331	343
332	344
333	345	=cut
334
335
336
337
338
339
	346	#line 348 "Basic.pm"
	347
	348
	349
	350	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
340	351	*skew_unbiased = \&PDL::skew_unbiased;
341
342
343
	352	#line 354 "Basic.pm"
	353
	354
	355
	356	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
344	357
345	358
346	359	=head2 kurt

365	378
366	379
367	380	=cut
368
369
370
371
372
373
	381	#line 383 "Basic.pm"
	382
	383
	384
	385	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
374	386	*kurt = \&PDL::kurt;
375
376
377
	387	#line 389 "Basic.pm"
	388
	389
	390
	391	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
378	392
379	393
380	394	=head2 kurt_unbiased

399	413
400	414
401	415	=cut
402
403
404
405
406
407
	416	#line 418 "Basic.pm"
	417
	418
	419
	420	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
408	421	*kurt_unbiased = \&PDL::kurt_unbiased;
409
410
411
	422	#line 424 "Basic.pm"
	423
	424
	425
	426	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
412	427
413	428
414	429	=head2 cov

433	448
434	449
435	450	=cut
436
437
438
439
440
441
	451	#line 453 "Basic.pm"
	452
	453
	454
	455	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
442	456	*cov = \&PDL::cov;
443
444
445
	457	#line 459 "Basic.pm"
	458
	459
	460
	461	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
446	462
447	463
448	464	=head2 cov_table

499	515
500	516
501	517	=cut
502
503
504
505
506
507
	518	#line 520 "Basic.pm"
	519
	520
	521
	522	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
508	523	*cov_table = \&PDL::cov_table;
509
510
511
	524	#line 526 "Basic.pm"
	525
	526
	527
	528	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
512	529
513	530
514	531	=head2 corr

555	572
556	573
557	574	=cut
558
559
560
561
562
563
	575	#line 577 "Basic.pm"
	576
	577
	578
	579	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
564	580	*corr = \&PDL::corr;
565
566
567
	581	#line 583 "Basic.pm"
	582
	583
	584
	585	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
568	586
569	587
570	588	=head2 corr_table

617	635
618	636
619	637	=cut
620
621
622
623
624
625
	638	#line 640 "Basic.pm"
	639
	640
	641
	642	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
626	643	*corr_table = \&PDL::corr_table;
627
628
629
	644	#line 646 "Basic.pm"
	645
	646
	647
	648	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
630	649
631	650
632	651	=head2 t_corr

661	680
662	681
663	682	=cut
664
665
666
667
668
669
	683	#line 685 "Basic.pm"
	684
	685
	686
	687	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
670	688	*t_corr = \&PDL::t_corr;
671
672
673
	689	#line 691 "Basic.pm"
	690
	691
	692
	693	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
674	694
675	695
676	696	=head2 n_pair

695	715
696	716
697	717	=cut
698
699
700
701
702
703
	718	#line 720 "Basic.pm"
	719
	720
	721
	722	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
704	723	*n_pair = \&PDL::n_pair;
705
706
707
	724	#line 726 "Basic.pm"
	725
	726
	727
	728	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
708	729
709	730
710	731	=head2 corr_dev

733	754
734	755
735	756	=cut
736
737
738
739
740
741
	757	#line 759 "Basic.pm"
	758
	759
	760
	761	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
742	762	*corr_dev = \&PDL::corr_dev;
743
744
745
	763	#line 765 "Basic.pm"
	764
	765
	766
	767	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
746	768
747	769
748	770	=head2 t_test

775	797
776	798
777	799	=cut
778
779
780
781
782
783
	800	#line 802 "Basic.pm"
	801
	802
	803
	804	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
784	805	*t_test = \&PDL::t_test;
785
786
787
	806	#line 808 "Basic.pm"
	807
	808
	809
	810	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
788	811
789	812
790	813	=head2 t_test_nev

813	836
814	837
815	838	=cut
816
817
818
819
820
821
	839	#line 841 "Basic.pm"
	840
	841
	842
	843	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
822	844	*t_test_nev = \&PDL::t_test_nev;
823
824
825
	845	#line 847 "Basic.pm"
	846
	847
	848
	849	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
826	850
827	851
828	852	=head2 t_test_paired

847	871
848	872
849	873	=cut
850
851
852
853
854
855
	874	#line 876 "Basic.pm"
	875
	876
	877
	878	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
856	879	*t_test_paired = \&PDL::t_test_paired;
857
858
859
860
861		#line 1252 "Basic/stats_basic.pd"
	880	#line 882 "Basic.pm"
	881
	882
	883
	884	#line 1252 "stats_basic.pd"
	885
	886	#line 1252 "stats_basic.pd"
862	887
863	888	=head2 binomial_test
864	889

1292	1317	All rights reserved. There is no warranty. You are allowed to redistribute this software / documentation as described in the file COPYING in the PDL distribution.
1293	1318
1294	1319	=cut
1295
1296
1297		;
	1320	#line 1687 "stats_basic.pd"
	1321	#line 1323 "Basic.pm"
	1322
	1323
	1324
1298	1325
1299	1326
1300	1327

+266

-231

GENERATED/PDL/Stats/Distr.pm less more

2	2	#
3	3	package PDL::Stats::Distr;
4	4
5		our @EXPORT_OK = qw( PDL::PP mme_beta PDL::PP pdf_beta PDL::PP mme_binomial PDL::PP pmf_binomial PDL::PP mle_exp PDL::PP pdf_exp PDL::PP mme_gamma PDL::PP pdf_gamma PDL::PP mle_gaussian PDL::PP pdf_gaussian PDL::PP mle_geo PDL::PP pmf_geo PDL::PP mle_geosh PDL::PP pmf_geosh PDL::PP mle_lognormal PDL::PP mme_lognormal PDL::PP pdf_lognormal PDL::PP mme_nbd PDL::PP pmf_nbd PDL::PP mme_pareto PDL::PP pdf_pareto PDL::PP mle_poisson PDL::PP pmf_poisson PDL::PP pmf_poisson_stirling PDL::PP _pmf_poisson_factorial );
	5	our @EXPORT_OK = qw(mme_beta pdf_beta mme_binomial pmf_binomial mle_exp pdf_exp mme_gamma pdf_gamma mle_gaussian pdf_gaussian mle_geo pmf_geo mle_geosh pmf_geosh mle_lognormal mme_lognormal pdf_lognormal mme_nbd pmf_nbd mme_pareto pdf_pareto mle_poisson pmf_poisson pmf_poisson_stirling _pmf_poisson_factorial );
6	6	our %EXPORT_TAGS = (Func=>\@EXPORT_OK);
7	7
8	8	use PDL::Core;

19	19
20	20
21	21
	22
	23	#line 4 "distr.pd"
	24
22	25	use strict;
23	26	use warnings;
24	27

67	70	$data->plot_distr( 'gaussian', 'lognormal' );
68	71
69	72	=cut
70
	73	#line 75 "Distr.pm"
71	74
72	75
73	76

81	84
82	85
83	86
	87	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
	88
84	89
85	90	=head2 mme_beta
86	91

108	113
109	114
110	115	=cut
111
112
113
114
115
116
	116	#line 118 "Distr.pm"
	117
	118
	119
	120	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
117	121	*mme_beta = \&PDL::mme_beta;
118
119
120
	122	#line 124 "Distr.pm"
	123
	124
	125
	126	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
121	127
122	128
123	129	=head2 pdf_beta

142	148
143	149
144	150	=cut
145
146
147
148
149
150
	151	#line 153 "Distr.pm"
	152
	153
	154
	155	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
151	156	*pdf_beta = \&PDL::pdf_beta;
152
153
154
	157	#line 159 "Distr.pm"
	158
	159
	160
	161	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
155	162
156	163
157	164	=head2 mme_binomial

180	187
181	188
182	189	=cut
183
184
185
186
187
188
	190	#line 192 "Distr.pm"
	191
	192
	193
	194	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
189	195	*mme_binomial = \&PDL::mme_binomial;
190
191
192
	196	#line 198 "Distr.pm"
	197
	198
	199
	200	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
193	201
194	202
195	203	=head2 pmf_binomial

214	222
215	223
216	224	=cut
217
218
219
220
221
222
	225	#line 227 "Distr.pm"
	226
	227
	228
	229	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
223	230	*pmf_binomial = \&PDL::pmf_binomial;
224
225
226
	231	#line 233 "Distr.pm"
	232
	233
	234
	235	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
227	236
228	237
229	238	=head2 mle_exp

252	261
253	262
254	263	=cut
255
256
257
258
259
260
	264	#line 266 "Distr.pm"
	265
	266
	267
	268	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
261	269	*mle_exp = \&PDL::mle_exp;
262
263
264
	270	#line 272 "Distr.pm"
	271
	272
	273
	274	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
265	275
266	276
267	277	=head2 pdf_exp

286	296
287	297
288	298	=cut
289
290
291
292
293
294
	299	#line 301 "Distr.pm"
	300
	301
	302
	303	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
295	304	*pdf_exp = \&PDL::pdf_exp;
296
297
298
	305	#line 307 "Distr.pm"
	306
	307
	308
	309	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
299	310
300	311
301	312	=head2 mme_gamma

324	335
325	336
326	337	=cut
327
328
329
330
331
332
	338	#line 340 "Distr.pm"
	339
	340
	341
	342	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
333	343	*mme_gamma = \&PDL::mme_gamma;
334
335
336
	344	#line 346 "Distr.pm"
	345
	346
	347
	348	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
337	349
338	350
339	351	=head2 pdf_gamma

358	370
359	371
360	372	=cut
361
362
363
364
365
366
	373	#line 375 "Distr.pm"
	374
	375
	376
	377	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
367	378	*pdf_gamma = \&PDL::pdf_gamma;
368
369
370
	379	#line 381 "Distr.pm"
	380
	381
	382
	383	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
371	384
372	385
373	386	=head2 mle_gaussian

396	409
397	410
398	411	=cut
399
400
401
402
403
404
	412	#line 414 "Distr.pm"
	413
	414
	415
	416	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
405	417	*mle_gaussian = \&PDL::mle_gaussian;
406
407
408
	418	#line 420 "Distr.pm"
	419
	420
	421
	422	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
409	423
410	424
411	425	=head2 pdf_gaussian

430	444
431	445
432	446	=cut
433
434
435
436
437
438
	447	#line 449 "Distr.pm"
	448
	449
	450
	451	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
439	452	*pdf_gaussian = \&PDL::pdf_gaussian;
440
441
442
	453	#line 455 "Distr.pm"
	454
	455
	456
	457	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
443	458
444	459
445	460	=head2 mle_geo

464	479
465	480
466	481	=cut
467
468
469
470
471
472
	482	#line 484 "Distr.pm"
	483
	484
	485
	486	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
473	487	*mle_geo = \&PDL::mle_geo;
474
475
476
	488	#line 490 "Distr.pm"
	489
	490
	491
	492	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
477	493
478	494
479	495	=head2 pmf_geo

498	514
499	515
500	516	=cut
501
502
503
504
505
506
	517	#line 519 "Distr.pm"
	518
	519
	520
	521	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
507	522	*pmf_geo = \&PDL::pmf_geo;
508
509
510
	523	#line 525 "Distr.pm"
	524
	525
	526
	527	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
511	528
512	529
513	530	=head2 mle_geosh

532	549
533	550
534	551	=cut
535
536
537
538
539
540
	552	#line 554 "Distr.pm"
	553
	554
	555
	556	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
541	557	*mle_geosh = \&PDL::mle_geosh;
542
543
544
	558	#line 560 "Distr.pm"
	559
	560
	561
	562	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
545	563
546	564
547	565	=head2 pmf_geosh

566	584
567	585
568	586	=cut
569
570
571
572
573
574
	587	#line 589 "Distr.pm"
	588
	589
	590
	591	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
575	592	*pmf_geosh = \&PDL::pmf_geosh;
576
577
578
	593	#line 595 "Distr.pm"
	594
	595
	596
	597	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
579	598
580	599
581	600	=head2 mle_lognormal

604	623
605	624
606	625	=cut
607
608
609
610
611
612
	626	#line 628 "Distr.pm"
	627
	628
	629
	630	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
613	631	*mle_lognormal = \&PDL::mle_lognormal;
614
615
616
	632	#line 634 "Distr.pm"
	633
	634
	635
	636	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
617	637
618	638
619	639	=head2 mme_lognormal

642	662
643	663
644	664	=cut
645
646
647
648
649
650
	665	#line 667 "Distr.pm"
	666
	667
	668
	669	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
651	670	*mme_lognormal = \&PDL::mme_lognormal;
652
653
654
	671	#line 673 "Distr.pm"
	672
	673
	674
	675	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
655	676
656	677
657	678	=head2 pdf_lognormal

676	697
677	698
678	699	=cut
679
680
681
682
683
684
	700	#line 702 "Distr.pm"
	701
	702
	703
	704	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
685	705	*pdf_lognormal = \&PDL::pdf_lognormal;
686
687
688
	706	#line 708 "Distr.pm"
	707
	708
	709
	710	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
689	711
690	712
691	713	=head2 mme_nbd

714	736
715	737
716	738	=cut
717
718
719
720
721
722
	739	#line 741 "Distr.pm"
	740
	741
	742
	743	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
723	744	*mme_nbd = \&PDL::mme_nbd;
724
725
726
	745	#line 747 "Distr.pm"
	746
	747
	748
	749	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
727	750
728	751
729	752	=head2 pmf_nbd

748	771
749	772
750	773	=cut
751
752
753
754
755
756
	774	#line 776 "Distr.pm"
	775
	776
	777
	778	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
757	779	*pmf_nbd = \&PDL::pmf_nbd;
758
759
760
	780	#line 782 "Distr.pm"
	781
	782
	783
	784	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
761	785
762	786
763	787	=head2 mme_pareto

786	810
787	811
788	812	=cut
789
790
791
792
793
794
	813	#line 815 "Distr.pm"
	814
	815
	816
	817	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
795	818	*mme_pareto = \&PDL::mme_pareto;
796
797
798
	819	#line 821 "Distr.pm"
	820
	821
	822
	823	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
799	824
800	825
801	826	=head2 pdf_pareto

820	845
821	846
822	847	=cut
823
824
825
826
827
828
	848	#line 850 "Distr.pm"
	849
	850
	851
	852	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
829	853	*pdf_pareto = \&PDL::pdf_pareto;
830
831
832
	854	#line 856 "Distr.pm"
	855
	856
	857
	858	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
833	859
834	860
835	861	=head2 mle_poisson

858	884
859	885
860	886	=cut
861
862
863
864
865
866
	887	#line 889 "Distr.pm"
	888
	889
	890
	891	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
867	892	*mle_poisson = \&PDL::mle_poisson;
868
869
870
	893	#line 895 "Distr.pm"
	894
	895
	896
	897	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
871	898
872	899
873	900	=head2 pmf_poisson

892	919
893	920
894	921	=cut
895
896
897
898
899
900
	922	#line 924 "Distr.pm"
	923
	924
	925
	926	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
901	927	*pmf_poisson = \&PDL::pmf_poisson;
902
903
904
	928	#line 930 "Distr.pm"
	929
	930
	931
	932	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
905	933
906	934
907	935	=head2 pmf_poisson_stirling

926	954
927	955
928	956	=cut
929
930
931
932
933
934
	957	#line 959 "Distr.pm"
	958
	959
	960
	961	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
935	962	*pmf_poisson_stirling = \&PDL::pmf_poisson_stirling;
936
937
938
939
940		#line 1139 "Distr/distr.pd"
	963	#line 965 "Distr.pm"
	964
	965
	966
	967	#line 1139 "distr.pd"
	968
	969	#line 1139 "distr.pd"
941	970
942	971	=head2 pmf_poisson_factorial
943	972

962	991	return _pmf_poisson_factorial(@_);
963	992	}
964	993	}
965
966
967
968
	994	#line 1165 "distr.pd"
	995	#line 997 "Distr.pm"
	996
	997
	998
	999	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
969	1000	*_pmf_poisson_factorial = \&PDL::_pmf_poisson_factorial;
970
971
972
973		;
974
975
976		#line 1201 "Distr/distr.pd"
	1001	#line 1003 "Distr.pm"
	1002
	1003
	1004
	1005
	1006
	1007	#line 1201 "distr.pd"
	1008
	1009	#line 1201 "distr.pd"
977	1010
978	1011	=head2 plot_distr
979	1012

1117	1150	All rights reserved. There is no warranty. You are allowed to redistribute this software / documentation as described in the file COPYING in the PDL distribution.
1118	1151
1119	1152	=cut
	1153	#line 1346 "distr.pd"
	1154	#line 1156 "Distr.pm"
1120	1155
1121	1156
1122	1157

+124

-106

GENERATED/PDL/Stats/GLM.pm less more

2	2	#
3	3	package PDL::Stats::GLM;
4	4
5		our @EXPORT_OK = qw( ols_t anova anova_rptd dummy_code effect_code effect_code_w interaction_code ols ols_rptd r2_change logistic pca pca_sorti plot_means plot_residuals plot_screes PDL::PP fill_m PDL::PP fill_rand PDL::PP dev_m PDL::PP stddz PDL::PP sse PDL::PP mse PDL::PP rmse PDL::PP pred_logistic PDL::PP d0 PDL::PP dm PDL::PP dvrs );
	5	our @EXPORT_OK = qw(ols_t anova anova_rptd dummy_code effect_code effect_code_w interaction_code ols ols_rptd r2_change logistic pca pca_sorti plot_means plot_residuals plot_screes fill_m fill_rand dev_m stddz sse mse rmse pred_logistic d0 dm dvrs );
6	6	our %EXPORT_TAGS = (Func=>\@EXPORT_OK);
7	7
8	8	use PDL::Core;

18	18
19	19
20	20
	21
	22
	23	#line 4 "glm.pd"
21	24
22	25	use strict;
23	26	use warnings;

66	69	print "$_\t$m{$_}\n" for (sort keys %m);
67	70
68	71	=cut
69
	72	#line 74 "GLM.pm"
70	73
71	74
72	75

79	82
80	83
81	84
	85
	86	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
82	87
83	88
84	89	=head2 fill_m

116	121
117	122
118	123	=cut
119
120
121
122
123
124
	124	#line 126 "GLM.pm"
	125
	126
	127
	128	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
125	129	*fill_m = \&PDL::fill_m;
126
127
128
	130	#line 132 "GLM.pm"
	131
	132
	133
	134	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
129	135
130	136
131	137	=head2 fill_rand

164	170
165	171
166	172	=cut
167
168
169
170
171
172
	173	#line 175 "GLM.pm"
	174
	175
	176
	177	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
173	178	*fill_rand = \&PDL::fill_rand;
174
175
176
	179	#line 181 "GLM.pm"
	180
	181
	182
	183	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
177	184
178	185
179	186	=head2 dev_m

198	205
199	206
200	207	=cut
201
202
203
204
205
206
	208	#line 210 "GLM.pm"
	209
	210
	211
	212	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
207	213	*dev_m = \&PDL::dev_m;
208
209
210
	214	#line 216 "GLM.pm"
	215
	216
	217
	218	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
211	219
212	220
213	221	=head2 stddz

231	239
232	240
233	241	=cut
234
235
236
237
238
239
	242	#line 244 "GLM.pm"
	243
	244
	245
	246	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
240	247	*stddz = \&PDL::stddz;
241
242
243
	248	#line 250 "GLM.pm"
	249
	250
	251
	252	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
244	253
245	254
246	255	=head2 sse

265	274
266	275
267	276	=cut
268
269
270
271
272
273
	277	#line 279 "GLM.pm"
	278
	279
	280
	281	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
274	282	*sse = \&PDL::sse;
275
276
277
	283	#line 285 "GLM.pm"
	284
	285
	286
	287	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
278	288
279	289
280	290	=head2 mse

299	309
300	310
301	311	=cut
302
303
304
305
306
307
	312	#line 314 "GLM.pm"
	313
	314
	315
	316	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
308	317	*mse = \&PDL::mse;
309
310
311
	318	#line 320 "GLM.pm"
	319
	320
	321
	322	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
312	323
313	324
314	325	=head2 rmse

333	344
334	345
335	346	=cut
336
337
338
339
340
341
	347	#line 349 "GLM.pm"
	348
	349
	350
	351	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
342	352	*rmse = \&PDL::rmse;
343
344
345
	353	#line 355 "GLM.pm"
	354
	355
	356
	357	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
346	358
347	359
348	360	=head2 pred_logistic

373	385
374	386
375	387	=cut
376
377
378
379
380
381
	388	#line 390 "GLM.pm"
	389
	390
	391
	392	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
382	393	*pred_logistic = \&PDL::pred_logistic;
383
384
385
	394	#line 396 "GLM.pm"
	395
	396
	397
	398	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
386	399
387	400
388	401	=head2 d0

410	423
411	424
412	425	=cut
413
414
415
416
417
418
	426	#line 428 "GLM.pm"
	427
	428
	429
	430	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
419	431	*d0 = \&PDL::d0;
420
421
422
	432	#line 434 "GLM.pm"
	433
	434
	435
	436	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
423	437
424	438
425	439	=head2 dm

450	464
451	465
452	466	=cut
453
454
455
456
457
458
	467	#line 469 "GLM.pm"
	468
	469
	470
	471	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
459	472	*dm = \&PDL::dm;
460
461
462
	473	#line 475 "GLM.pm"
	474
	475
	476
	477	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
463	478
464	479
465	480	=head2 dvrs

484	499
485	500
486	501	=cut
487
488
489
490
491
492
	502	#line 504 "GLM.pm"
	503
	504
	505
	506	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
493	507	*dvrs = \&PDL::dvrs;
494
495
496
497
498		#line 594 "GLM/glm.pd"
	508	#line 510 "GLM.pm"
	509
	510
	511
	512	#line 594 "glm.pd"
	513
	514	#line 594 "glm.pd"
499	515
500	516	# my tmp var for PDL 2.007 slice upate
501	517	my $_tmp;

2557	2573	All rights reserved. There is no warranty. You are allowed to redistribute this software / documentation as described in the file COPYING in the PDL distribution.
2558	2574
2559	2575	=cut
2560
2561
2562		;
	2576	#line 2657 "glm.pd"
	2577	#line 2579 "GLM.pm"
	2578
	2579
	2580
2563	2581
2564	2582
2565	2583

+73

-55

GENERATED/PDL/Stats/Kmeans.pm less more

2	2	#
3	3	package PDL::Stats::Kmeans;
4	4
5		our @EXPORT_OK = qw( random_cluster iv_cluster PDL::PP _random_cluster PDL::PP which_cluster PDL::PP assign PDL::PP centroid PDL::PP _d_p2l );
	5	our @EXPORT_OK = qw(random_cluster iv_cluster _random_cluster which_cluster assign centroid _d_p2l );
6	6	our %EXPORT_TAGS = (Func=>\@EXPORT_OK);
7	7
8	8	use PDL::Core;

18	18
19	19
20	20
	21
	22
	23	#line 4 "kmeans.pd"
21	24
22	25	use Carp;
23	26	use PDL::LiteF;

76	79	for (0 .. $k{cluster}->dim(1)-1);
77	80
78	81	=cut
79
	82	#line 84 "Kmeans.pm"
80	83
81	84
82	85

90	93
91	94
92	95
93		#line 75 "Kmeans/kmeans.pd"
	96	#line 75 "kmeans.pd"
	97
	98	#line 75 "kmeans.pd"
94	99
95	100	# my tmp var for PDL 2.007 slice upate
96	101	my $_tmp;

121	126	} while (PDL::any $cluster->sumover == 0 );
122	127	return $cluster;
123	128	}
124
125
126
127
	129	#line 107 "kmeans.pd"
	130	#line 132 "Kmeans.pm"
	131
	132
	133
	134	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
128	135	*_random_cluster = \&PDL::_random_cluster;
129
130
131
	136	#line 138 "Kmeans.pm"
	137
	138
	139
	140	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
132	141
133	142
134	143	=head2 which_cluster

163	172
164	173
165	174	=cut
166
167
168
169
170
171
	175	#line 177 "Kmeans.pm"
	176
	177
	178
	179	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
172	180	*which_cluster = \&PDL::which_cluster;
173
174
175
	181	#line 183 "Kmeans.pm"
	182
	183
	184
	185	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
176	186
177	187
178	188	=head2 assign

223	233
224	234
225	235	=cut
226
227
228
229
230
231
	236	#line 238 "Kmeans.pm"
	237
	238
	239
	240	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
232	241	*assign = \&PDL::assign;
233
234
235
	242	#line 244 "Kmeans.pm"
	243
	244
	245
	246	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
236	247
237	248
238	249	=head2 centroid

292	303
293	304
294	305	=cut
295
296
297
298
299
300
	306	#line 308 "Kmeans.pm"
	307
	308
	309
	310	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
301	311	*centroid = \&PDL::centroid;
302
303
304
305
306		#line 432 "Kmeans/kmeans.pd"
	312	#line 314 "Kmeans.pm"
	313
	314
	315
	316	#line 432 "kmeans.pd"
	317
	318	#line 432 "kmeans.pd"
307	319
308	320	sub _scree_ind {
309	321	# use as scree cutoff the point with max distance to the line formed

318	330	croak "1D pdl only please";
319	331
320	332	my $a = zeroes 2, $self->nelem;
321		($_tmp = $a((0), )) .= sequence $self->nelem;
322		($_tmp = $a((1), )) .= $self;
323
324		my $d = _d_point2line( $a, $a( ,(0)), $a( ,(-1)) );
	333	($_tmp = $a->slice('(0)')) .= sequence $self->nelem;
	334	($_tmp = $a->slice('(1)')) .= $self;
	335
	336	my $d = _d_point2line( $a, $a->slice(':,(0)'), $a->slice(':,(-1)') );
325	337
326	338	return $d->maximum_ind;
327	339	}

336	348
337	349	return _d_p2l( $self->mv(0,-1)->dog, $p1->mv(0,-1)->dog, $p2->mv(0,-1)->dog );
338	350	}
339
340
341
342
	351	#line 466 "kmeans.pd"
	352	#line 354 "Kmeans.pm"
	353
	354
	355
	356	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
343	357	*_d_p2l = \&PDL::_d_p2l;
344
345
346
347
348		#line 495 "Kmeans/kmeans.pd"
	358	#line 360 "Kmeans.pm"
	359
	360
	361
	362	#line 495 "kmeans.pd"
	363
	364	#line 495 "kmeans.pd"
349	365
350	366	=head2 kmeans
351	367

709	725	All rights reserved. There is no warranty. You are allowed to redistribute this software / documentation as described in the file COPYING in the PDL distribution.
710	726
711	727	=cut
712
713
714		;
	728	#line 860 "kmeans.pd"
	729	#line 731 "Kmeans.pm"
	730
	731
	732
715	733
716	734
717	735

+147

-119

GENERATED/PDL/Stats/TS.pm less more

2	2	#
3	3	package PDL::Stats::TS;
4	4
5		our @EXPORT_OK = qw(PDL::PP _acf PDL::PP _acvf PDL::PP diff PDL::PP inte PDL::PP dseason PDL::PP _fill_ma PDL::PP filter_exp PDL::PP filter_ma PDL::PP mae PDL::PP mape PDL::PP wmape PDL::PP portmanteau PDL::PP _pred_ar );
	5	our @EXPORT_OK = qw(_acf _acvf diff inte dseason _fill_ma filter_exp filter_ma mae mape wmape portmanteau _pred_ar );
6	6	our %EXPORT_TAGS = (Func=>\@EXPORT_OK);
7	7
8	8	use PDL::Core;

18	18
19	19
20	20
	21
	22
	23	#line 2 "ts.pd"
21	24
22	25	=encoding utf8
23	26

48	51	use PDL::Stats::Kmeans;
49	52
50	53	my $DEV = ($^O =~ /win/i)? '/png' : '/xs';
51
	54	#line 56 "TS.pm"
52	55
53	56
54	57

62	65
63	66
64	67
65
	68	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
66	69	*_acf = \&PDL::_acf;
67
68
69
70
71
	70	#line 72 "TS.pm"
	71
	72
	73
	74	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
72	75	*_acvf = \&PDL::_acvf;
73
74
75
76
77		#line 113 "TS/ts.pd"
	76	#line 78 "TS.pm"
	77
	78
	79
	80	#line 113 "ts.pd"
	81
	82	#line 113 "ts.pd"
78	83
79	84	=head2 acf
80	85

144	149	$h \|\|= $self->dim(0) - 1;
145	150	return $self->_acvf($h+1);
146	151	}
147
148
	152	#line 184 "ts.pd"
	153	#line 155 "TS.pm"
	154
	155
	156
	157	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
149	158
150	159
151	160	=head2 diff

168	177
169	178
170	179	=cut
171
172
173
174
175
176
	180	#line 182 "TS.pm"
	181
	182
	183
	184	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
177	185	*diff = \&PDL::diff;
178
179
180
	186	#line 188 "TS.pm"
	187
	188
	189
	190	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
181	191
182	192
183	193	=head2 inte

200	210
201	211
202	212	=cut
203
204
205
206
207
208
	213	#line 215 "TS.pm"
	214
	215
	216
	217	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
209	218	*inte = \&PDL::inte;
210
211
212
	219	#line 221 "TS.pm"
	220
	221
	222
	223	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
213	224
214	225
215	226	=head2 dseason

233	244
234	245
235	246	=cut
236
237
238
239
240
241
	247	#line 249 "TS.pm"
	248
	249
	250
	251	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
242	252	*dseason = \&PDL::dseason;
243
244
245
246
247		#line 363 "TS/ts.pd"
	253	#line 255 "TS.pm"
	254
	255
	256
	257	#line 363 "ts.pd"
	258
	259	#line 363 "ts.pd"
248	260
249	261	=head2 fill_ma
250	262

276	288
277	289	return $x_filled;
278	290	}
279
280
281
282
	291	#line 396 "ts.pd"
	292	#line 294 "TS.pm"
	293
	294
	295
	296	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
283	297	*_fill_ma = \&PDL::_fill_ma;
284
285
286
	298	#line 300 "TS.pm"
	299
	300
	301
	302	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
287	303
288	304
289	305	=head2 filter_exp

309	325
310	326
311	327	=cut
312
313
314
315
316
317
	328	#line 330 "TS.pm"
	329
	330
	331
	332	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
318	333	*filter_exp = \&PDL::filter_exp;
319
320
321
	334	#line 336 "TS.pm"
	335
	336
	337
	338	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
322	339
323	340
324	341	=head2 filter_ma

342	359
343	360
344	361	=cut
345
346
347
348
349
350
	362	#line 364 "TS.pm"
	363
	364
	365
	366	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
351	367	*filter_ma = \&PDL::filter_ma;
352
353
354
	368	#line 370 "TS.pm"
	369
	370
	371
	372	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
355	373
356	374
357	375	=head2 mae

381	399
382	400
383	401	=cut
384
385
386
387
388
389
	402	#line 404 "TS.pm"
	403
	404
	405
	406	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
390	407	*mae = \&PDL::mae;
391
392
393
	408	#line 410 "TS.pm"
	409
	410
	411
	412	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
394	413
395	414
396	415	=head2 mape

420	439
421	440
422	441	=cut
423
424
425
426
427
428
	442	#line 444 "TS.pm"
	443
	444
	445
	446	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
429	447	*mape = \&PDL::mape;
430
431
432
	448	#line 450 "TS.pm"
	449
	450
	451
	452	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
433	453
434	454
435	455	=head2 wmape

459	479
460	480
461	481	=cut
462
463
464
465
466
467
	482	#line 484 "TS.pm"
	483
	484
	485
	486	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
468	487	*wmape = \&PDL::wmape;
469
470
471
	488	#line 490 "TS.pm"
	489
	490
	491
	492	#line 1059 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
472	493
473	494
474	495	=head2 portmanteau

510	531
511	532
512	533	=cut
513
514
515
516
517
518
	534	#line 536 "TS.pm"
	535
	536
	537
	538	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
519	539	*portmanteau = \&PDL::portmanteau;
520
521
522
523
524		#line 701 "TS/ts.pd"
	540	#line 542 "TS.pm"
	541
	542
	543
	544	#line 701 "ts.pd"
	545
	546	#line 701 "ts.pd"
525	547
526	548	=head2 pred_ar
527	549

578	600	return $x->append($ext($b->dim(0) : -1));
579	601	}
580	602	}
581
582
583
584
	603	#line 759 "ts.pd"
	604	#line 606 "TS.pm"
	605
	606
	607
	608	#line 1061 "/home/osboxes/.perlbrew/libs/perl-5.32.0@normal/lib/perl5/x86_64-linux/PDL/PP.pm"
585	609	*_pred_ar = \&PDL::_pred_ar;
586
587
588
589
590		#line 791 "TS/ts.pd"
	610	#line 612 "TS.pm"
	611
	612
	613
	614	#line 791 "ts.pd"
	615
	616	#line 791 "ts.pd"
591	617
592	618	=head2 season_m
593	619

813	839	All rights reserved. There is no warranty. You are allowed to redistribute this software / documentation as described in the file COPYING in the PDL distribution.
814	840
815	841	=cut
816
817
818		;
	842	#line 1018 "ts.pd"
	843	#line 845 "TS.pm"
	844
	845
	846
819	847
820	848
821	849

-4

Kmeans/kmeans.pd less more

442	442	croak "1D pdl only please";
443	443
444	444	my $a = zeroes 2, $self->nelem;
445		($_tmp = $a((0), )) .= sequence $self->nelem;
446		($_tmp = $a((1), )) .= $self;
447
448		my $d = _d_point2line( $a, $a( ,(0)), $a( ,(-1)) );
	445	($_tmp = $a->slice('(0)')) .= sequence $self->nelem;
	446	($_tmp = $a->slice('(1)')) .= $self;
	447
	448	my $d = _d_point2line( $a, $a->slice(':,(0)'), $a->slice(':,(-1)') );
449	449
450	450	return $d->maximum_ind;
451	451	}

-1

META.json less more

58	58	"web" : "https://github.com/PDLPorters/PDL-Stats"
59	59	}
60	60	},
61		"version" : "0.80",
	61	"version" : "0.81",
62	62	"x_serialization_backend" : "JSON::PP version 4.04"
63	63	}

-1

META.yml less more

24	24	PDL: '2.057'
25	25	resources:
26	26	repository: git://github.com/PDLPorters/PDL-Stats
27		version: '0.80'
	27	version: '0.81'
28	28	x_serialization_backend: 'CPAN::Meta::YAML version 0.018'

-1

Stats.pm less more

2	2	use strict;
3	3	use warnings;
4	4
5		our $VERSION = '0.80';
	5	our $VERSION = '0.81';
6	6
7	7	sub import {
8	8	my $pkg = (caller())[0];