Francisco Coelho / sasp

Download as
Browse Code »

Commit 5bcfee4d5d465618db31bea4d2e408aa27e94e89

Authored by Francisco Coelho
1 parent e3447353
Exists in master

Added README

Showing 4 changed files with 202 additions and 191 deletions   Show diff stats
README.md
  Diff comments   View file @5bcfee4
... ... @@ -0,0 +1,9 @@
  1 +# SASP - Stochastic ASP
  2 +
  3 +> This repository contains code and data used in the paper _An Algebraic Approach to Stochastic ASP_
  4 +
  5 +
  6 +There are two [julia](https://julialang.org/) files:
  7 +
  8 +- `sampling.jl` generates (and saves) a sample following the process described in the _Model+Noise_ item in subsection 4.1 of the paper.
  9 +- `sample_analysis.jl` analyses, summarizes and outputs `LaTeX` rows of the samples above.
0 10 \ No newline at end of file
... ...
data/model_sample-[0.0;0.3;0.2].csv
  Diff comments   View file @5bcfee4
1 1 event
2 2 A
3 3 A
4   -A
5   -A
6   -A
7   -A
8   -A
  4 +ac
9 5 ab
10 6 A
11   -A
12   -A
13   -A
14   -A
15 7 ac
16   -A
17   -A
  8 +ac
18 9 ac
19 10 A
20 11 A
21 12 A
22   -ab
23 13 A
24 14 A
25 15 A
26 16 A
27 17 A
28 18 A
29   -ac
30 19 A
31 20 A
32 21 A
33   -ac
34 22 A
35 23 A
36   -ac
37   -ac
38 24 A
39 25 ac
40 26 A
41 27 A
42   -ac
43 28 A
44 29 A
45 30 A
  31 +ab
  32 +A
46 33 A
47 34 A
48 35 A
... ... @@ -50,7 +37,6 @@ A
50 37 A
51 38 A
52 39 ac
53   -ab
54 40 A
55 41 A
56 42 A
... ... @@ -59,8 +45,13 @@ A
59 45 ac
60 46 A
61 47 A
  48 +ab
  49 +ab
62 50 A
63 51 ac
  52 +ac
  53 +A
  54 +ab
64 55 A
65 56 A
66 57 A
... ... @@ -68,45 +59,50 @@ ac
68 59 A
69 60 A
70 61 A
71   -ac
72 62 A
73 63 A
74 64 A
75 65 A
76 66 ac
  67 +ac
77 68 A
78 69 A
79 70 A
80   -ab
81   -A
82 71 A
  72 +ac
83 73 A
84   -ab
85 74 A
86   -ab
  75 +ac
87 76 A
  77 +ac
88 78 A
89 79 A
90 80 ab
91 81 A
92 82 A
  83 +ac
93 84 A
94 85 ac
95 86 ac
  87 +ac
96 88 A
97   -ab
98 89 A
99 90 A
100   -ac
101 91 A
102 92 A
  93 +ac
103 94 A
104 95 A
  96 +ab
  97 +ac
  98 +ac
  99 +ac
105 100 ac
106 101 A
107 102 A
108 103 A
109 104 A
  105 +ac
110 106 A
111 107 A
112 108 ac
... ... @@ -115,14 +111,10 @@ ac
115 111 A
116 112 A
117 113 ac
118   -ac
119 114 A
120   -ac
121   -ac
122 115 A
123 116 A
124 117 A
125   -ac
126 118 A
127 119 A
128 120 A
... ... @@ -130,18 +122,22 @@ ac
130 122 A
131 123 ac
132 124 A
  125 +A
  126 +A
  127 +A
  128 +A
  129 +A
133 130 ab
134   -ac
135 131 A
136 132 ac
  133 +ac
  134 +A
137 135 A
138 136 A
139   -ab
140 137 A
141 138 ac
142 139 A
143 140 ac
144   -ab
145 141 ac
146 142 A
147 143 A
... ... @@ -149,18 +145,20 @@ ac
149 145 A
150 146 A
151 147 A
  148 +ac
152 149 A
153 150 A
154 151 A
155 152 ac
  153 +ac
156 154 A
157 155 A
158 156 A
159 157 A
160 158 A
161   -ac
162 159 A
163 160 A
  161 +ab
164 162 ac
165 163 A
166 164 A
... ... @@ -174,22 +172,18 @@ A
174 172 A
175 173 A
176 174 A
177   -A
178   -A
179 175 ac
180 176 A
181 177 A
182   -ac
183   -A
184 178 A
  179 +ac
185 180 A
186 181 ac
  182 +A
187 183 ac
188 184 ac
189 185 A
190 186 A
191   -ac
192   -A
193 187 A
194 188 A
195 189 A
... ... @@ -197,21 +191,16 @@ A
197 191 A
198 192 A
199 193 A
200   -ab
201 194 A
  195 +ac
202 196 A
203 197 A
204 198 A
205   -ac
206   -ac
207 199 A
208 200 ab
209 201 A
210   -A
211 202 ac
212   -ab
213   -A
214   -A
  203 +ac
215 204 A
216 205 A
217 206 ac
... ... @@ -220,15 +209,11 @@ ac
220 209 A
221 210 ac
222 211 A
223   -A
224 212 ac
225 213 A
226 214 A
227   -ab
228   -ac
229 215 A
230   -ac
231   -ac
  216 +A
232 217 A
233 218 A
234 219 ac
... ... @@ -238,26 +223,27 @@ ac
238 223 A
239 224 A
240 225 A
241   -ac
242 226 A
  227 +ac
243 228 A
244 229 A
245 230 A
  231 +ac
246 232 A
247 233 A
248 234 A
249 235 A
250 236 A
251 237 A
  238 +ab
252 239 A
253 240 A
254 241 A
255 242 A
256   -ab
  243 +ac
257 244 ac
258 245 A
259 246 A
260   -ac
261 247 A
262 248 A
263 249 A
... ... @@ -265,41 +251,53 @@ A
265 251 A
266 252 ab
267 253 ac
268   -A
269 254 ac
270 255 A
271 256 ac
272   -A
273 257 ac
274 258 A
275 259 A
276   -ab
277 260 A
278 261 A
  262 +A
  263 +A
  264 +A
  265 +A
  266 +ac
279 267 ac
280 268 ac
281 269 A
282 270 A
  271 +A
283 272 ac
284 273 A
285 274 A
286 275 A
287 276 A
288 277 ac
  278 +ac
  279 +A
  280 +ab
  281 +ab
  282 +A
289 283 A
290 284 A
  285 +A
  286 +ab
  287 +ab
  288 +ac
291 289 ab
292 290 ac
293 291 A
294 292 A
295   -ac
296 293 A
297 294 A
298 295 A
299 296 ab
  297 +ab
  298 +A
300 299 A
301 300 ac
302   -ab
303 301 A
304 302 A
305 303 ac
... ... @@ -307,12 +305,19 @@ A
307 305 A
308 306 A
309 307 A
  308 +ab
  309 +A
  310 +ac
  311 +A
310 312 A
311 313 A
312 314 ac
  315 +A
313 316 ac
314 317 A
  318 +ab
315 319 A
  320 +ac
316 321 A
317 322 ac
318 323 A
... ... @@ -325,30 +330,33 @@ A
325 330 A
326 331 A
327 332 A
  333 +ab
  334 +ac
328 335 A
  336 +ac
329 337 A
330 338 A
331 339 A
332 340 A
333 341 A
334 342 A
335   -ac
336 343 A
337   -ac
338 344 A
339 345 ac
  346 +A
340 347 ac
341 348 ac
342 349 A
343 350 A
  351 +A
  352 +ac
  353 +A
344 354 ac
345 355 A
346 356 A
347 357 A
348 358 A
349 359 ac
350   -ac
351   -A
352 360 A
353 361 A
354 362 A
... ... @@ -359,125 +367,110 @@ A
359 367 ab
360 368 A
361 369 A
362   -ab
363 370 ac
364   -ab
  371 +A
  372 +A
365 373 ac
  374 +ab
366 375 ac
367 376 A
368 377 A
369   -ab
370 378 ac
371 379 A
372 380 ac
373 381 A
  382 +ac
374 383 A
375 384 A
376 385 A
377 386 ac
378   -ac
379   -A
380 387 A
381 388 ac
382 389 A
383 390 A
384 391 A
385 392 A
386   -ab
387   -A
388   -A
389   -ac
390 393 ac
391 394 A
392 395 A
393 396 A
394   -A
395   -A
396 397 ac
397 398 A
398 399 A
399 400 A
400 401 A
401 402 A
402   -A
403   -A
404 403 ac
405 404 A
406 405 A
407 406 ac
408   -ac
409   -A
410 407 A
411 408 A
412 409 A
413   -ac
414   -ac
415 410 A
416   -ac
417 411 A
418 412 ac
419 413 A
420 414 A
421 415 A
422   -A
423 416 ac
424   -A
425   -A
426   -A
427   -A
428 417 ac
  418 +ab
429 419 A
430 420 A
  421 +ac
  422 +ac
431 423 A
432 424 A
433 425 A
434 426 A
  427 +ac
  428 +ab
435 429 A
436 430 A
437 431 A
438   -ac
439 432 A
440 433 A
441 434 A
442 435 A
443 436 A
  437 +ab
444 438 A
445 439 A
446 440 ac
447 441 ac
  442 +ac
  443 +A
448 444 A
449 445 ac
450 446 A
451 447 A
452 448 A
453 449 A
454   -ac
455 450 A
456 451 ac
457 452 A
458 453 ac
  454 +ab
459 455 A
460 456 A
461 457 A
462 458 A
463   -ab
464 459 A
465   -ac
466 460 A
467 461 A
468   -ab
469 462 A
470   -ac
471 463 A
472 464 ac
473 465 ac
474   -ac
475 466 A
476 467 A
477 468 A
478 469 A
479 470 A
480 471 A
  472 +ab
  473 +ab
481 474 A
482 475 A
483 476 A
... ... @@ -487,9 +480,6 @@ A
487 480 A
488 481 A
489 482 A
490   -ac
491   -ab
492   -ab
493 483 A
494 484 A
495 485 A
... ... @@ -500,31 +490,37 @@ ac
500 490 A
501 491 A
502 492 A
  493 +ac
503 494 A
504 495 A
505 496 A
  497 +ac
  498 +ac
506 499 A
  500 +ab
507 501 A
508 502 A
509 503 A
510   -A
511   -ab
512   -A
513 504 ac
514 505 A
  506 +ac
515 507 A
516 508 A
517 509 A
518 510 A
519 511 A
  512 +ac
520 513 A
  514 +ab
521 515 A
522 516 ac
523 517 A
524 518 A
525 519 A
526 520 A
  521 +ac
527 522 A
  523 +ac
528 524 A
529 525 A
530 526 ac
... ... @@ -532,15 +528,16 @@ A
532 528 A
533 529 A
534 530 A
  531 +ab
535 532 ac
536 533 A
537 534 ac
538   -ac
539   -A
540 535 A
541 536 A
542 537 A
  538 +ac
543 539 A
  540 +ac
544 541 A
545 542 A
546 543 A
... ... @@ -548,21 +545,16 @@ A
548 545 ab
549 546 A
550 547 A
551   -ab
552 548 A
553   -ac
554 549 A
  550 +ac
  551 +ac
555 552 A
556 553 A
557   -ab
558 554 A
559   -ac
560   -ab
561 555 A
562 556 A
563   -ac
564 557 A
565   -ab
566 558 A
567 559 ac
568 560 A
... ... @@ -570,11 +562,11 @@ A
570 562 A
571 563 A
572 564 A
573   -ac
574 565 A
575   -ac
576 566 ab
577 567 ac
  568 +ac
  569 +A
578 570 A
579 571 A
580 572 A
... ... @@ -584,74 +576,68 @@ A
584 576 ac
585 577 A
586 578 ac
587   -ab
588 579 A
589   -ab
590   -ac
591 580 A
592 581 ac
593 582 ac
594 583 ac
595 584 A
596   -A
597   -ab
598   -A
599   -ab
600   -A
601   -A
602 585 ac
603 586 A
604 587 ac
605 588 A
  589 +A
606 590 ac
607 591 A
608 592 A
  593 +ac
609 594 A
610 595 A
  596 +ac
611 597 A
612 598 A
613 599 A
614 600 ac
615 601 A
  602 +ac
616 603 A
617 604 ac
618 605 ac
  606 +A
619 607 ab
620 608 A
621 609 A
622 610 A
623 611 A
624 612 ac
625   -ac
626 613 A
627 614 A
628 615 A
629 616 A
630   -ab
631   -A
632 617 A
  618 +ac
633 619 A
634 620 A
635 621 A
636 622 A
637 623 A
638   -ac
639   -ac
640 624 A
  625 +ab
641 626 A
642 627 A
643 628 A
  629 +ac
644 630 A
645 631 A
646 632 A
  633 +ac
647 634 A
648   -ab
649 635 A
650 636 ac
651 637 A
652 638 A
653 639 A
654   -ac
  640 +A
655 641 A
656 642 ac
657 643 ac
... ... @@ -661,9 +647,6 @@ A
661 647 A
662 648 A
663 649 A
664   -ab
665   -ac
666   -ac
667 650 A
668 651 A
669 652 A
... ... @@ -671,13 +654,19 @@ A
671 654 A
672 655 A
673 656 A
674   -ab
  657 +ac
  658 +A
  659 +A
  660 +A
  661 +A
675 662 A
676 663 A
677 664 ac
678   -ac
679 665 A
680 666 A
  667 +A
  668 +ab
  669 +ac
681 670 ac
682 671 A
683 672 A
... ... @@ -689,14 +678,15 @@ ac
689 678 A
690 679 A
691 680 ac
  681 +ac
692 682 A
693 683 A
  684 +ac
694 685 A
695   -ab
696 686 A
697   -ab
698 687 A
699 688 A
  689 +ac
700 690 A
701 691 A
702 692 A
... ... @@ -705,6 +695,10 @@ A
705 695 ac
706 696 A
707 697 A
  698 +ac
  699 +ac
  700 +A
  701 +A
708 702 A
709 703 A
710 704 A
... ... @@ -712,27 +706,21 @@ A
712 706 A
713 707 A
714 708 A
715   -ac
716 709 A
717   -ac
718 710 A
719 711 A
720   -ac
721   -ac
722 712 ac
723 713 A
724 714 A
725 715 A
726   -ac
727 716 A
728 717 ac
729 718 A
730 719 A
731 720 A
732 721 A
733   -ab
734   -ac
735   -ac
  722 +A
  723 +A
736 724 ac
737 725 A
738 726 A
... ... @@ -740,19 +728,24 @@ A
740 728 A
741 729 A
742 730 A
  731 +ac
  732 +A
  733 +A
  734 +ab
  735 +A
743 736 A
744 737 A
745 738 ac
746 739 ac
747 740 A
748   -A
749 741 ac
750 742 A
751 743 A
752 744 A
753   -ab
  745 +A
754 746 A
755 747 ac
  748 +ac
756 749 A
757 750 ac
758 751 A
... ... @@ -763,6 +756,8 @@ A
763 756 A
764 757 A
765 758 ac
  759 +ac
  760 +A
766 761 A
767 762 A
768 763 A
... ... @@ -770,34 +765,37 @@ A
770 765 ab
771 766 ab
772 767 A
773   -ac
774 768 A
775 769 A
776 770 A
777 771 A
778 772 A
  773 +ab
779 774 A
780 775 A
  776 +ac
  777 +ac
781 778 A
782 779 A
783   -ac
784 780 A
785 781 A
786 782 ac
787 783 A
788 784 A
789 785 A
790   -ac
791   -A
  786 +ab
792 787 A
793 788 ac
794 789 ac
795 790 A
796   -ac
797 791 A
798 792 A
  793 +ab
  794 +ac
  795 +ac
799 796 A
800 797 A
  798 +ac
801 799 A
802 800 A
803 801 A
... ... @@ -808,28 +806,26 @@ A
808 806 A
809 807 A
810 808 A
  809 +ac
  810 +ab
811 811 A
812 812 A
813 813 A
  814 +ac
814 815 A
815 816 A
816   -ac
817 817 A
818 818 A
819 819 A
820 820 ac
821   -ac
  821 +ab
  822 +A
822 823 A
823 824 A
824   -ac
825 825 A
826   -ab
827 826 A
828   -ac
829 827 A
830 828 A
831   -ac
832   -ac
833 829 A
834 830 A
835 831 ab
... ... @@ -837,65 +833,57 @@ A
837 833 A
838 834 A
839 835 A
840   -ac
841 836 A
842   -ac
843   -ac
844 837 A
845 838 A
846   -ab
847 839 A
848 840 A
849 841 ac
850 842 A
851 843 A
852 844 A
  845 +A
  846 +A
853 847 ac
854 848 A
855 849 A
856 850 A
857 851 A
858 852 A
  853 +ac
859 854 A
860 855 A
861 856 A
862 857 A
863 858 ac
  859 +ab
864 860 A
865 861 A
  862 +ac
866 863 A
867 864 A
868   -ac
869 865 A
870 866 A
871 867 ab
872   -ab
873   -ac
874   -A
875 868 A
876 869 ac
877 870 A
878 871 A
879 872 ac
880   -ac
881 873 A
882   -ac
883 874 A
884 875 A
885   -ac
886 876 ab
887 877 A
  878 +ac
  879 +A
888 880 A
889 881 A
890   -ac
891   -ac
892   -ac
893 882 A
  883 +ab
894 884 A
895 885 A
896 886 ac
897   -ac
898   -ac
899 887 A
900 888 A
901 889 A
... ... @@ -908,62 +896,68 @@ A
908 896 A
909 897 A
910 898 A
911   -ac
  899 +ab
912 900 A
  901 +ab
913 902 ac
914 903 A
  904 +ab
915 905 A
916 906 A
917 907 ac
918 908 A
919   -A
  909 +ac
920 910 A
921 911 ac
  912 +ac
922 913 A
923 914 A
924 915 A
925   -ac
  916 +ab
926 917 A
927 918 A
928 919 A
929 920 ac
930   -ac
931   -ac
932   -ab
933 921 A
934   -ac
  922 +A
  923 +A
935 924 A
936 925 A
937 926 A
938 927 A
939 928 A
940 929 ac
941   -ab
942 930 A
943 931 A
  932 +ac
944 933 A
  934 +ac
945 935 A
946 936 ac
947 937 A
948 938 A
949 939 ac
950 940 A
951   -ab
952 941 A
953 942 A
  943 +A
  944 +ab
  945 +A
954 946 ac
955 947 A
956 948 A
  949 +ab
  950 +A
  951 +A
957 952 A
958 953 A
959 954 ac
960 955 A
961 956 A
962   -ac
963 957 A
964 958 A
965 959 ac
966   -ac
  960 +A
967 961 ac
968 962 A
969 963 A
... ... @@ -971,11 +965,14 @@ A
971 965 A
972 966 A
973 967 A
  968 +ab
  969 +A
974 970 A
975 971 A
976 972 ac
977 973 ac
978 974 ac
  975 +A
979 976 ac
980 977 A
981 978 A
... ... @@ -983,19 +980,22 @@ A
983 980 A
984 981 A
985 982 ac
  983 +ac
986 984 A
987 985 A
  986 +ac
  987 +A
  988 +ac
988 989 A
989 990 A
990 991 A
991   -ac
992 992 A
993   -ac
994 993 A
995   -ac
996 994 A
997 995 A
998 996 A
999 997 A
  998 +ac
  999 +A
1000 1000 A
1001 1001 A
... ...
sample_analysis.jl
  Diff comments   View file @5bcfee4
... ... @@ -63,11 +63,13 @@ core_of = Dict{String,String}(
63 63 "Δ" => "A,ab,ac",
64 64 )
65 65  
66   -α = 0.1
  66 +
  67 +α = 0.0
67 68 β = 0.3
68   -γ = 0.8
  69 +γ = 0.2
69 70  
70   -data = CSV.read("model_sample-[$α;$β;$γ].csv", DataFrame)
  71 +
  72 +data = CSV.read("data/model_sample-[$α;$β;$γ].csv", DataFrame)
71 73 transform!(data, :event => ByRow(e -> get(core_of, e, "z")) => :core)
72 74 groups = groupby(data, :core)
73 75 totals = combine(groups, :core => length => :count)
... ...
sampling.jl
  Diff comments   View file @5bcfee4
... ... @@ -74,9 +74,9 @@ consistent = filter(is_consistent, omega)
74 74 consistent = join.(consistent)
75 75  
76 76  
77   -α = 0.1
  77 +α = 0.0
78 78 β = 0.3
79   -γ = 0.5
  79 +γ = 0.2
80 80  
81 81 events = model_sample(1000;
82 82 α = α,
... ... @@ -85,7 +85,7 @@ events = model_sample(1000;
85 85  
86 86 println( events )
87 87  
88   -open("model_sample-[$α;$β;$γ].csv", "w") do io
  88 +open("data/model_sample-[$α;$β;$γ].csv", "w") do io
89 89 writedlm(io, [ "event" ])
90 90 writedlm(io, events)
91 91 end
92 92 \ No newline at end of file
... ...