move examples upwards.

packages/CLPBN/examples/School/README

@@ -0,0 +1,53 @@
+
+This is a version of the school database, based on the PRM School example.
+
+There are four main files:
+
+school_128.yap: a school with 128 professors, 256 courses and 4096 students.
+school_64.yap: medium size school
+school_32.yap: small school
+
+schema.yap: the schema
+tables: CPTs
+
+
+
+
+
+professor_ability(p0,X).
+
+professor_popularity(p0,X).
+
+professor_ability(p0,X), professor_popularity(p0,h).
+
+professor_ability(p0,h), professor_popularity(p0,X).
+
+registration_grade(r0,X).
+
+registration_grade(r0,X), registration_course(r0,C), course_difficulty(C,h).
+
+registration_grade(r0,X), registration_course(r0,C), course_difficulty(C,h), registration_student(r0,S), student_intelligence(S,h).
+
+registration_grade(r0,X), registration_course(r0,C), course_difficulty(C,l), registration_student(r0,S), student_intelligence(S,h).
+
+registration_satisfaction(r0,X).
+
+registration_satisfaction(r0,X), registration_student(r0,S), student_intelligence(S,h).
+
+registration_satisfaction(r0,X), registration_grade(r0,a).
+
+registration_satisfaction(r0,X), registration_grade(r0,d).
+
+registration_satisfaction(r0,h), registration_grade(r0,X).
+
+course_rating(c0,X).
+
+course_rating(c0,h), course_difficulty(c0,X).
+
+course_difficulty(c0,X).
+
+student_ranking(s0,X).
+
+student_ranking(s0,X), student_intelligence(s0,h).
+
+

packages/CLPBN/examples/School/evidence_128.yap

@@ -0,0 +1,17 @@
+
+:- [school_128].
+
+professor_popularity(p0,h) :- {}.
+professor_popularity(p3,h) :- {}.
+professor_popularity(p5,l) :- {}.
+professor_popularity(p45,h) :- {}.
+professor_popularity(p15,m) :- {}.
+
+course_rating(c0, h) :- {}.
+course_rating(c1, m) :- {}.
+course_rating(c2, l) :- {}.
+course_rating(c3, h) :- {}.
+course_rating(c4, m) :- {}.
+course_rating(c5, l) :- {}.
+course_rating(c62, m) :- {}.
+

packages/CLPBN/examples/School/parlearn.yap

@@ -0,0 +1,44 @@
+
+:- [pos:sample32].
+
+:- ['~/Yap/work/CLPBN/clpbn/examples/School/school_32'].
+
+% These libraries provide same functionality.
+:- [library('clpbn/learning/mle')].
+%:- [library('clpbn/learning/bnt_parms')].
+
+:- [library(matrix)].
+
+main :-
+	findall(X,goal(X),L),
+	learn_parameters(L,CPTs),
+	write_cpts(CPTs).
+
+goal(professor_ability(P,V)) :-
+	pos:professor_ability(P,V),
+	p(pa, M), random < M.
+goal(professor_popularity(P,V)) :-
+	pos:professor_popularity(P,V),
+	p(pp, M), random < M.
+goal(registration_grade(P,V)) :-
+	pos:registration_grade(P,V),
+	p(rg, M), random < M.
+goal(student_intelligence(P,V)) :-
+	pos:student_intelligence(P,V),
+	p(si, M), random < M.
+goal(course_difficulty(P,V)) :-
+	pos:course_difficulty(P,V),
+	p(cd, M), random < M.
+goal(registration_satisfaction(P,V)) :-
+	pos:registration_satisfaction(P,V),
+	p(rs, M), random < M.
+
+% sampling parameter
+p(_, 1.0).
+
+write_cpts([]).
+write_cpts([CPT|CPTs]) :-
+	matrix_to_list(CPT,L),
+	format('CPT=~w~n',[L]),
+	write_cpts(CPTs).
+	

packages/CLPBN/examples/School/schema.yap

@@ -0,0 +1,71 @@
+
+/* base file for school database. Supposed to be called from school_*.yap */
+
+professor_key(Key) :-
+	professor(Key).
+
+professor_ability(Key,Abi) :-
+	abi_table(Key, AbiDist),
+	{ Abi = ability(Key) with p([h,m,l], AbiDist)   }.
+
+professor_popularity(Key, Pop) :-
+	professor_ability(Key, Abi),
+	pop_table(Key,PopTable),
+	{ Pop = popularity(Key) with
+		p([h,m,l], PopTable,[Abi]) }.
+
+registration_key(Key) :-
+	registration(Key, _, _).
+
+registration_course(Key, CKey) :-
+	registration(Key, CKey, _).
+
+registration_student(Key, SKey) :-
+	registration(Key, _, SKey).
+
+registration_grade(Key, Grade) :-
+	registration(Key, CKey, SKey),
+	course_difficulty(CKey, Dif),
+	student_intelligence(SKey, Int),
+	grade_table(Int, Dif, Table),
+	{ Grade = grade(Key) with Table }.
+
+% registration_satisfaction(r0, h) :- {}.
+registration_satisfaction(Key, Sat) :-
+	registration_course(Key, CKey),
+	course_professor(CKey, PKey),
+	professor_ability(PKey, Abi),
+	registration_grade(Key, Grade),
+	satisfaction_table(Abi, Grade, Table),
+	{ Sat = satisfaction(Key) with Table }.
+
+course_key(Key) :-
+	course(Key,_).
+	
+course_professor(Key, PKey) :-
+	course(Key, PKey).
+	
+course_rating(CKey, Rat) :-
+	setof(Sat, RKey^(registration_course(RKey,CKey), registration_satisfaction(RKey,Sat)), Sats),
+	{ Rat =  rating(CKey) with avg([h,m,l],Sats) }.
+
+course_difficulty(Key, Dif) :-
+	dif_table(Key, Dist),
+	{ Dif = difficulty(Key) with p([h,m,l], Dist) }.
+
+student_key(Key) :-
+	student(Key).
+
+student_intelligence(Key, Int) :-
+	int_table(Key, IDist, Domain),
+	{ Int = intelligence(Key) with p(Domain, IDist) }.
+
+student_ranking(Key, Rank) :-
+	setof(Grade, CKey^(registration_student(CKey,Key),
+			 registration_grade(CKey, Grade)), Grades),
+	{ Rank = ranking(Key) with avg([a,b,c,d],Grades) }.
+
+:- ensure_loaded(tables).
+
+
+

packages/CLPBN/examples/School/tables.yap

@@ -0,0 +1,45 @@
+
+int_table(_, [0.5,
+	      0.4,
+	      0.1],[h, m, l]).
+
+grade_table(I, D, 
+	   /* h h   h m  h l  m h  m m  m l  l h  l m l l */
+	    p([a,b,c,d],
+		[ 0.2,  0.7, 0.85, 0.1, 0.2, 0.5, 0.01, 0.05,0.1 ,
+	          0.6, 0.25, 0.12, 0.3, 0.6,0.35,0.04, 0.15, 0.4 ,
+                  0.15,0.04, 0.02, 0.4,0.15,0.12, 0.5, 0.6, 0.4,
+	          0.05,0.01, 0.01, 0.2,0.05,0.03, 0.45, 0.2, 0.1 ], [I,D])).
+    
+	
+/*
+  A: professor's ability;
+  B: student's grade (for course registration).
+*/
+satisfaction_table(A, G,
+       /* h a  h b  h c  h d  m a  m b  m c  m d  l a  l b  l c  l d */
+	p([h,m,l],
+/*h*/     [0.98, 0.9,0.8 , 0.6, 0.9, 0.4, 0.2, 0.01, 0.5, 0.2, 0.01, 0.01,
+/*m*/	   0.01, 0.09,0.15, 0.3,0.05, 0.4, 0.3, 0.04,0.35, 0.3, 0.09, 0.01 ,
+/*l*/	   0.01, 0.01,0.05, 0.1,0.05, 0.2, 0.5, 0.95,0.15, 0.5, 0.9, 0.98], [A,G])).
+
+
+% The idea is quite simple:
+%  hs = h -> r = ( 0.9, 0.1, 0)
+%  hs = m -> r = ( 0.2, 0.6, 0.2)
+%  hs = l -> r = ( 0, 0.1, 0.9)
+%
+% add all and divide on the number of elements on the table!
+%
+rating_prob_table([0.9,0.05,0.01,
+		   0.09,0.9,0.09,
+		   0.01,0.05,0.9]).
+
+abi_table( _, [0.50, 0.40,  0.10]).
+
+
+pop_table(_, [0.9, 0.2, 0.01,
+              0.09, 0.6, 0.09,
+	      0.01, 0.2, 0.9]).
+
+dif_table( _, [0.25, 0.50, 0.25]).