mindreadtestshared.pl

%% Name, DOB, Date learned, psych appointment month=0 or 1, psych appointment day, thoughts count

sectest0 :-
meditators(A),
findall(_,(member(B,A),sectest(B)),_)
/*
sectest([first,last,dobd,dobm,doby,daylearned,monthlearned,yearlearned,0,1,16]),
sectest([first,last,dobd,dobm,doby,daylearned,monthlearned,yearlearned,0,1,16])
*/.

sectest1 :-
meditators2(A),
findall(_,(member(B,A),sectest(B)),_)
/*
sectest([first,last,dobd,dobm,doby,daylearned,monthlearned,yearlearned,1,0,16]),
sectest([first,last,dobd,dobm,doby,daylearned,monthlearned,yearlearned,1,0,16])
*/
.

find_time(H,M,S) :-
	trialy2_15("0",H11),
	trialy2_15("1",H12),
	trialy2_15("2",H13),
	H1L=[H11,H12,H13],
	sort(H1L,H1A),
	reverse(H1A,H1B),
	H1B=[[_,H1]|_Rest1],

	(H1="2"->(
		trialy2_30("0",H21),
		trialy2_30("1",H22),
		trialy2_30("2",H23),
		trialy2_30("3",H24),
		H2L=[H21,H22,H23,H24],
		sort(H2L,H2A),
		reverse(H2A,H2B),
		H2B=[[_,H2]|_Rest2]
	)
	;(	
		trialy2_30("0",H21),
		trialy2_30("1",H22),
		trialy2_30("2",H23),
		trialy2_30("3",H24),
		trialy2_30("4",H25),
		trialy2_30("5",H26),
		trialy2_30("6",H27),
		trialy2_30("7",H28),
		trialy2_30("8",H29),
		trialy2_30("9",H210),
	
		H2L=[H21,H22,H23,H24,H25,
		H26,H27,H28,H29,H210],
		sort(H2L,H2A),
		reverse(H2A,H2B),
		H2B=[[_,H2]|_Rest2]
	)),

	trialy2_15("0",M11),
	trialy2_15("1",M12),
	trialy2_15("2",M13),
	trialy2_15("3",M14),
	trialy2_15("4",M15),
	trialy2_15("5",M16),
	M1L=[M11,M12,M13,M14,M15,M16],
	sort(M1L,M1A),
	reverse(M1A,M1B),
	M1B=[[_,M1]|_Rest3],

	trialy2_30("0",M21),
	trialy2_30("1",M22),
	trialy2_30("2",M23),
	trialy2_30("3",M24),
	trialy2_30("4",M25),
	trialy2_30("5",M26),
	trialy2_30("6",M27),
	trialy2_30("7",M28),
	trialy2_30("8",M29),
	trialy2_30("9",M210),
	M2L=[M21,M22,M23,M24,M25,M26,M27,M28,M29,M210],
	sort(M2L,M2A),
	reverse(M2A,M2B),
	M2B=[[_,M2]|_Rest4],

	trialy2_15("0",S11),
	trialy2_15("1",S12),
	trialy2_15("2",S13),
	trialy2_15("3",S14),
	trialy2_15("4",S15),
	trialy2_15("5",S16),
	S1L=[S11,S12,S13,S14,S15,S16],
	sort(S1L,S1A),
	reverse(S1A,S1B),
	S1B=[[_,S1]|_Rest5],

	trialy2_30("0",S21),
	trialy2_30("1",S22),
	trialy2_30("2",S23),
	trialy2_30("3",S24),
	trialy2_30("4",S25),
	trialy2_30("5",S26),
	trialy2_30("6",S27),
	trialy2_30("7",S28),
	trialy2_30("8",S29),
	trialy2_30("9",S210),
	S2L=[S21,S22,S23,S24,S25,S26,S27,S28,S29,S210],
	sort(S2L,S2A),
	reverse(S2A,S2B),
	S2B=[[_,S2]|_Rest6],
	
	string_concat(H1,H2,H),
	string_concat(M1,M2,M),
	string_concat(S1,S2,S).

trialy2_6(Label,RA) :-
	%%writeln([testing,Label]),
	trialy1(R1),
	trialy1(R2),
	trialy1(R3),
	trialy1(R4),
	trialy1(R5),
	trialy1(R6), /**
	trialy1(R7),
	trialy1(R8),
	trialy1(R9),
	trialy1(R10),
	trialy1(R11),
	trialy1(R12),
	trialy1(R13),
	trialy1(R14),
	trialy1(R15),
	trialy1(R16),
	trialy1(R17),
	trialy1(R18),
	trialy1(R19),
	trialy1(R20),
	trialy1(R21),
	trialy1(R22),
	trialy1(R23),
	trialy1(R24),
	trialy1(R25),
	trialy1(R26),
	trialy1(R27),
	trialy1(R28),
	trialy1(R29),
	trialy1(R30), **/
	R=[R1,R2,R3,R4,R5,R6 /**,R7,R8,R9,R10,
	R11,R12,R13,R14,R15,R16,R17,R18,R19,R20,
	R21,R22,R23,R24,R25,R26,R27,R28,R29,R30 **/
	],
	%%(member(true,R)->(
	aggregate_all(count, member(true,R), Count),
	RA=[Count,Label].%%,writeln([Label,Count,"/10"]));true).

trialy2_15(Label,RA) :-
	%%writeln([testing,Label]),
	trialy1(R1),
	trialy1(R2),
	trialy1(R3),
	trialy1(R4),
	trialy1(R5),
	trialy1(R6),
	trialy1(R7),
	trialy1(R8),
	trialy1(R9),
	trialy1(R10),
	trialy1(R11),
	trialy1(R12),
	trialy1(R13),
	trialy1(R14),
	trialy1(R15),
	/**
	trialy1(R16),
	trialy1(R17),
	trialy1(R18),
	trialy1(R19),
	trialy1(R20),
	trialy1(R21),
	trialy1(R22),
	trialy1(R23),
	trialy1(R24),
	trialy1(R25),
	trialy1(R26),
	trialy1(R27),
	trialy1(R28),
	trialy1(R29),
	trialy1(R30),**/
	R=[R1,R2,R3,R4,R5,R6,R7,R8,R9,R10,
	R11,R12,R13,R14,R15 /**,R16,R17,R18,R19,R20,
	R21,R22,R23,R24,R25,R26,R27,R28,R29,R30
	**/],
	%%(member(true,R)->(
	aggregate_all(count, member(true,R), Count),
	RA=[Count,Label].%%,writeln([Label,Count,"/10"]));true).

trialy2_30(Label,RA) :-
	%%writeln([testing,Label]),
	trialy1(R1),
	trialy1(R2),
	trialy1(R3),
	trialy1(R4),
	trialy1(R5),
	trialy1(R6),
	trialy1(R7),
	trialy1(R8),
	trialy1(R9),
	trialy1(R10),
	trialy1(R11),
	trialy1(R12),
	trialy1(R13),
	trialy1(R14),
	trialy1(R15),
	trialy1(R16),
	trialy1(R17),
	trialy1(R18),
	trialy1(R19),
	trialy1(R20),
	trialy1(R21),
	trialy1(R22),
	trialy1(R23),
	trialy1(R24),
	trialy1(R25),
	trialy1(R26),
	trialy1(R27),
	trialy1(R28),
	trialy1(R29),
	trialy1(R30),
	R=[R1,R2,R3,R4,R5,R6,R7,R8,R9,R10,
	R11,R12,R13,R14,R15,R16,R17,R18,R19,R20,
	R21,R22,R23,R24,R25,R26,R27,R28,R29,R30],
	%%(member(true,R)->(
	aggregate_all(count, member(true,R), Count),
	RA=[Count,Label].%%,writeln([Label,Count,"/10"]));true).


trialy1(R1) :-
	%%control11(A1),
	trial0(A22), %% Control
	sum(A22,0,S22),
	mean(S22,A1),
	trial0(A21), %% Test 1
	sum(A21,0,S02),
	mean(S02,A2),
	(A1>A2->R1=true;R1=fail).

trial0(S3) :- N is 10, trial1(N,[],S),trial01(S,S3).
trial01(S1,S3) :-
	sort(S1,S),
	%%midpoint(S,MP),
	halves(S,H1,H2),
	midpoint(H1,Q1),
	midpoint(H2,Q3),
	IQR is Q3-Q1,
	sum(S,0,S02),
	mean(S02,Mean),
	furthestfrommean(S,Mean,V),
	D1 is 1.5*IQR,
	D2 is V-Mean,
	(D2>D1->(delete(S,V,S2),trial01(S2,S3));S=S3).
	
trial1(0,A,A) :- !.
trial1(N,A,B) :- mindreadtest(S), append(A,[S],A2),
	N1 is N-1,trial1(N1,A2,B).
	
midpoint(S,MP) :-
	length(S,L),
	A is mod(L,2),
	(A is 0->
		(M1 is L/2, M2 is M1+1,N1 is M1-1,N2 is M2-1,length(N11,N1),length(N21,N2),append(N11,[N12|_Rest1],S),append(N21,[N22|_Rest2],S),MP is (N12+N22)/2)
	;
		(L2 is L+1, M1 is L2/2, N1 is M1-1,length(N11,N1),append(N11,[MP|_Rest],S))).

halves(S,H1,H2) :-
	length(S,L),
	A is mod(L,2),
	(A is 0->
		(M1 is L/2,length(H1,M1),append(H1,H2,S))
	;
		(L2 is L-1,M1 is L2/2,length(H1,M1),append(H1,[_|H2],S))).

sum([],S,S):-!.
sum(S0,S1,S2) :-
	S0=[S3|S4],
	S5 is S1+S3,
	sum(S4,S5,S2).
	
mean(Sum,Mean) :-
	Mean is Sum/2.

furthestfrommean(S,Mean,V) :-
	absdiffmean(S,Mean,[],D),
	sort(D,D1),
	reverse(D1,[[_,V]|_Rest]).

absdiffmean([],_M,D,D) :- !.
absdiffmean(S,M,D1,D2) :-
	S=[S1|S2],
	S3 is abs(S1-M),
	append(D1,[[S3,S1]],D3),
	absdiffmean(S2,M,D3,D2).

mindreadtest(Sec) :-
	%% 250 br for characters to be br out with 10 br each from person to me - do when initial 250 br test done and doing 10 br test
	%%comment(fiftyastest),
	%%random(X),X1 is 10*X, X2 is floor(X1), (X2=<2 -> (
	%%texttobr,writeln(['true test']), %%); %% use breasonings breasoned out by computer for not by me, for job medicine for "me", at last time point
	%%true), %% leave last time point blank
	%%**texttobr2(640);true),%% make an A to detect reaction to gracious giving or blame of in following
	get_time(TimeStamp1),
	%%phrase_from_file(string(_String), 'file.txt'),
	texttobr2(2), %% 100 As for answer (must be br before this on same day)
	%% is gracious giving or blame
	get_time(TimeStamp2),
	%%comment(turnoffas),
   Sec is TimeStamp2 - TimeStamp1.

shell1(Command) :-
				(bash_command(Command,_)->
					true;
					(writeln(["Failed shell1 command: ",Command]),abort)
				),!.

bash_command(Command, Output) :-
        setup_call_cleanup(process_create(path(bash),
                ['-c', Command],
                [stdout(pipe(Out))]),
        read_string(Out, _, Output),
        close(Out)).