title 'conjugacy example';
data conjugacy;
input task	lnuy	diffyes;
cards;
1	0.5	0.05
1	0.55	0.2
1	0.5	0
1	0.5	0
1	0.3	0.2
1	0.5	0.2
1	0.2	0
1	0.15	0.1
1	0.35	0.1
1	0.5	0.1
1	0.5	0
1	0.1	0.1
1	0.45	0.1
1	0.45	0.1
1	0.4	0.2
1	0.15	0.1
1	0.5	0.05
1	0.2	0
1	0.48	0.05
1	0.55	0.1
1	0.5	0.4
1	0.5	0.4
1	0.65	0.2
1	0.5	0.8
1	0.25	0.3
1	0.4	0.3
1	0.3	0.3
1	0.5	0
1	0.45	0.2
1	0.4	0.4
1	0.5	0.1
1	0.4	0.2
1	0.5	0.6
1	0.35	0.1
1	0.3	0.1
1	0.4	0.15
1	0.5	0.1
1	0.14	0
1	0.5	0.1
1	0.5	0.2
1	0.3	0.2
1	0.08	0.05
1	0.48	0.05
1	0.5	0.14
1	0.4	0.2
1	0.13	0.05
1	0.5	0.5
1	0.38	0.04
1	0.5	0.5
1	0.4	0.2
1	0.5	0
1	0.5	0.4
1	0.4	0.2
1	0.45	0
1	0.5	0.1
1	0.35	0.2
1	0.52	0.05
1	0.6	0.6
1	0.5	0.4
1	0.3	0.1
1	0.5	0.2
1	0.5	0.2
1	0.65	0.2
1	0.15	0.15
1	0.4	0.2
1	0.55	0.6
1	0.35	0.4
1	0.5	0.5
1	0.5	0.15
1	0.5	0.15
1	0.5	0.2
1	0.5	0.1
1	0.02	0.01
1	0.38	0.05
1	0.5	0.2
1	0.5	0.98
1	0.5	0.1
1	0.3	0.1
1	0.5	0.8
1	0.55	0.2
1	0.3	0.2
1	0.5	0.1
1	0.35	0.3
1	0.5	0.2
1	0.55	0.1
1	0.25	0.1
1	0.5	0.4
1	0.45	0.4
1	0.55	0.2
1	0.75	0.5
1	0.75	0.8
1	0.5	0.1
1	0.5	0.1
1	0.35	0.2
1	0.5	0.2
1	0.5	0.2
1	0.3	0.3
1	0.45	0.2
1	0.5	0.15
1	0.25	0.1
1	0.4	0.5
1	0.5	0.08
1	0.23	0.35
1	0.18	0.05
1	0.57	0.18
1	0.5	0.8
1	0.5	0.6
1	0.35	0.25
1	0.48	0.05
1	0.7	0.6
1	0.5	0.2
1	0.2	0.09
1	0.4	0.2
1	0.5	0.15
1	0.53	0.21
1	0.38	0.25
1	0.5	0.05
1	0.5	0.01
1	0.25	0.2
1	0.45	0.1
1	0.45	0.1
1	0.35	0.7
1	0.2	0.1
1	0.5	0.72
1	0.3	0.2
1	0.3	0.1
1	0.4	0.1
1	0.45	0.4
1	0.5	0.3
1	0.5	0.2
1	0.5	1
1	0.5	0.1
1	0.5	0.2
1	0.45	0.3
1	0.35	0.2
1	0.5	0.08
1	0.3	0.2
1	0.4	0.2
1	0.5	0.2
1	0.5	0.1
1	0.3	0.15
1	0.5	0
1	0.15	0.1
1	0.4	0.2
1	0.5	0.1
1	0.5	0.17
1	0.5	0
1	0.3	0.15
1	0.5	0.08
1	0.45	0.1
1	0.5	0.1
1	0.35	0.2
1	0.3	0.1
1	0.5	0.01
1	0.33	0.15
1	0.5	0.1
1	0.07	0.13
1	0.4	0.1
1	0.5	0.1
1	0.45	0.1
1	0.93	0.86
1	0.45	0.1
1	0.3	0.4
1	0.45	0.2
1	0.5	0.2
1	0.5	0.05
1	0.14	0
1	0.5	0.1
1	0.5	0.1
1	0.5	0.01
1	0.25	0.1
-1	0.33	0.4
-1	0.5	0.4
-1	0.25	0.45
-1	0.28	0.4
-1	0.45	0.7
-1	0.3	0.3
-1	0.28	0.35
-1	0.5	0.35
-1	0.45	0.6
-1	0.43	0.55
-1	0.33	0.3
-1	0.25	0.3
-1	0.5	0.9
-1	0.23	0.1
-1	0.55	0.4
-1	0.5	0.8
-1	0.35	0.4
-1	0.45	0.7
-1	0.43	0.75
-1	0.45	0.7
-1	0.45	0.45
-1	0.48	0.6
-1	0.5	0.5
-1	0.48	0.19
-1	0.48	0.55
-1	0.45	0.7
-1	0.53	0.55
-1	0.3	0.4
-1	0.5	0.55
-1	0.4	0.5
-1	0.5	0.6
-1	0.58	0.65
-1	0.45	0.55
-1	0.53	0.55
-1	0.48	0.29
-1	0.5	0.5
-1	0.53	0.65
-1	0.5	0.1
-1	0.5	0.3
-1	0.5	0.7
-1	0.53	0.35
-1	0.35	0.5
-1	0.35	0.55
-1	0.48	0.6
-1	0.48	0.55
-1	0.55	0.45
-1	0.45	0.2
-1	0.4	0.55
-1	0.48	0.25
-1	0.33	0.55
-1	0.35	0.6
-1	0.5	0.05
-1	0.53	0.6
-1	0.43	0.15
-1	0.5	0.4
-1	0.5	0.05
-1	0.6	0.4
-1	0.53	0.35
-1	0.45	0.55
-1	0.58	0.55
-1	0.6	0.5
-1	0.48	0.2
-1	0.48	0.2
-1	0.58	0.4
-1	0.63	0.6
-1	0.53	0.15
-1	0.23	0.24
-1	0.5	0.04
-1	0.5	0.14
-1	0.1	0.15
-1	0.63	0.4
;
run;
;
proc nlmixed data = conjugacy tech = trureg hess cov itdetails;
*Starting values;
parms B0 = -0.9,  B1 = 1.6, B2 = 0.1, B3 = 0.1, G0 = -2.5, Q = -0.8, Q1 = 0.5;
title 'lnuy on task and diffyes (example 2 code), mmixture model';
*task is whether the task was Boeing or Sunday;
*diffyes is the difference between upper and lower P(Yes);
*Composition submodel;
C = 1/(1+EXP(Q + Q1*task));
*Location submodel;
B = B0 + B1*diffyes + B2*task + B3*diffyes*task;
M1 = EXP(B)/(1+EXP(B));
M2 = 1/2;
*Dispersion submodel;
PHI1 = EXP(-G0);
PHI2 = EXP(8);
ll = log(C*PDF('BETA',lnuy, M1*PHI1, PHI1 - M1*PHI1,0,1)+(1-C)* PDF('BETA',lnuy, M2*PHI2, PHI2 - M2*PHI2,0,1));
model lnuy ~ general(ll);
;
run;