title 'segale example';
data segale;
input partit	elicit	pyes;
cards;
-1	-1	0.949
-1	-1	0.5
-1	-1	0.151
-1	-1	0.35
-1	-1	0.5
-1	-1	0.151
-1	-1	0.141
-1	-1	0.151
-1	-1	0.144
-1	-1	0.091
-1	-1	0.081
-1	-1	0.301
-1	-1	0.301
-1	-1	0.144
-1	-1	0.151
-1	-1	0.141
-1	-1	0.251
-1	-1	0.45
-1	-1	0.101
-1	-1	0.251
-1	-1	0.161
-1	-1	0.151
-1	-1	0.141
-1	-1	0.051
-1	-1	0.144
-1	-1	0.151
-1	-1	0.141
-1	-1	0.101
-1	-1	0.144
-1	-1	0.201
-1	-1	0.051
-1	-1	0.5
-1	-1	0.101
-1	-1	0.35
-1	-1	0.4
-1	-1	0.201
-1	-1	0.101
-1	-1	0.151
-1	-1	0.144
-1	-1	0.201
-1	-1	0.5
-1	-1	0.201
-1	-1	0.301
-1	-1	0.141
-1	-1	0.151
-1	-1	0.141
-1	-1	0.251
-1	-1	0.5
-1	-1	0.101
-1	-1	0.151
-1	-1	0.38
-1	-1	0.301
-1	-1	0.121
-1	-1	0.101
-1	-1	0.144
-1	-1	0.101
-1	-1	0.121
-1	-1	0.5
-1	-1	0.301
-1	-1	0.101
-1	-1	0.151
-1	-1	0.101
-1	-1	0.5
-1	-1	0.131
-1	-1	0.141
-1	-1	0.151
-1	-1	0.251
-1	-1	0.141
-1	-1	0.5
-1	-1	0.5
-1	-1	0.141
-1	-1	0.141
-1	-1	0.699
-1	-1	0.5
-1	-1	0.021
-1	-1	0.151
-1	-1	0.101
-1	-1	0.5
-1	-1	0.201
-1	-1	0.5
-1	-1	0.301
-1	-1	0.5
-1	-1	0.141
-1	-1	0.5
-1	-1	0.5
-1	-1	0.141
-1	-1	0.121
-1	-1	0.699
-1	-1	0.5
-1	-1	0.051
-1	1	0.176
-1	1	0.201
-1	1	0.5
-1	1	0.5
-1	1	0.201
-1	1	0.5
-1	1	0.201
-1	1	0.151
-1	1	0.151
-1	1	0.151
-1	1	0.5
-1	1	0.101
-1	1	0.151
-1	1	0.151
-1	1	0.201
-1	1	0.151
-1	1	0.076
-1	1	0.201
-1	1	0.176
-1	1	0.151
-1	1	0.5
-1	1	0.301
-1	1	0.5
-1	1	0.5
-1	1	0.251
-1	1	0.251
-1	1	0.251
-1	1	0.5
-1	1	0.4
-1	1	0.4
-1	1	0.151
-1	1	0.5
-1	1	0.5
-1	1	0.151
-1	1	0.5
-1	1	0.226
-1	1	0.151
-1	1	0.141
-1	1	0.151
-1	1	0.201
-1	1	0.201
-1	1	0.076
-1	1	0.475
-1	1	0.071
-1	1	0.5
-1	1	0.126
-1	1	0.251
-1	1	0.141
-1	1	0.251
-1	1	0.201
-1	1	0.5
-1	1	0.5
-1	1	0.5
-1	1	0.5
-1	1	0.5
-1	1	0.4
-1	1	0.126
-1	1	0.4
-1	1	0.5
-1	1	0.151
-1	1	0.5
-1	1	0.5
-1	1	0.5
-1	1	0.076
-1	1	0.201
-1	1	0.6
-1	1	0.5
-1	1	0.5
-1	1	0.126
-1	1	0.126
-1	1	0.301
-1	1	0.151
-1	1	0.016
-1	1	0.126
-1	1	0.201
-1	1	0.5
-1	1	0.151
-1	1	0.051
-1	1	0.5
-1	1	0.5
-1	1	0.301
-1	1	0.201
-1	1	0.35
-1	1	0.5
-1	1	0.849
-1	1	0.151
-1	1	0.5
-1	1	0.301
-1	1	0.5
1	-1	0.181
1	-1	0.131
1	-1	0.151
1	-1	0.141
1	-1	0.201
1	-1	0.144
1	-1	0.4
1	-1	0.281
1	-1	0.051
1	-1	0.144
1	-1	0.5
1	-1	0.144
1	-1	0.021
1	-1	0.141
1	-1	0.144
1	-1	0.146
1	-1	0.141
1	-1	0.101
1	-1	0.141
1	-1	0.011
1	-1	0.131
1	-1	0.144
1	-1	0.141
1	-1	0.141
1	-1	0.151
1	-1	0.144
1	-1	0.144
1	-1	0.151
1	-1	0.151
1	-1	0.151
1	-1	0.141
1	-1	0.151
1	-1	0.301
1	-1	0.201
1	-1	0.161
1	-1	0.301
1	-1	0.101
1	-1	0.151
1	-1	0.4
1	-1	0.301
1	-1	0.101
1	-1	0.101
1	-1	0.5
1	-1	0.151
1	-1	0.101
1	-1	0.151
1	-1	0.141
1	-1	0.051
1	-1	0.201
1	-1	0.031
1	-1	0.101
1	-1	0.151
1	-1	0.051
1	-1	0.141
1	-1	0.201
1	-1	0.151
1	-1	0.141
1	-1	0.051
1	-1	0.031
1	-1	0.141
1	-1	0.201
1	-1	0.141
1	-1	0.141
1	-1	0.699
1	-1	0.141
1	-1	0.141
1	-1	0.144
1	-1	0.151
1	-1	0.301
1	-1	0.141
1	-1	0.201
1	-1	0.151
1	-1	0.051
1	-1	0.101
1	-1	0.141
1	-1	0.101
1	-1	0.141
1	-1	0.201
1	-1	0.6
1	-1	0.151
1	-1	0.101
1	-1	0.151
1	-1	0.011
1	-1	0.141
1	1	0.5
1	1	0.5
1	1	0.151
1	1	0.151
1	1	0.5
1	1	0.201
1	1	0.4
1	1	0.251
1	1	0.301
1	1	0.176
1	1	0.251
1	1	0.35
1	1	0.161
1	1	0.226
1	1	0.176
1	1	0.211
1	1	0.5
1	1	0.5
1	1	0.276
1	1	0.126
1	1	0.4
1	1	0.101
1	1	0.096
1	1	0.301
1	1	0.176
1	1	0.176
1	1	0.326
1	1	0.126
1	1	0.146
1	1	0.201
1	1	0.051
1	1	0.151
1	1	0.35
1	1	0.151
1	1	0.5
1	1	0.301
1	1	0.251
1	1	0.151
1	1	0.5
1	1	0.251
1	1	0.201
1	1	0.5
1	1	0.151
1	1	0.201
1	1	0.251
1	1	0.201
1	1	0.061
1	1	0.4
1	1	0.5
1	1	0.5
1	1	0.251
1	1	0.326
1	1	0.141
1	1	0.151
1	1	0.201
1	1	0.151
1	1	0.216
1	1	0.141
1	1	0.276
1	1	0.161
1	1	0.151
1	1	0.151
1	1	0.4
1	1	0.151
1	1	0.146
1	1	0.276
1	1	0.151
1	1	0.066
1	1	0.151
1	1	0.151
1	1	0.151
1	1	0.57
1	1	0.251
1	1	0.301
1	1	0.5
1	1	0.101
1	1	0.126
1	1	0.141
1	1	0.151
1	1	0.151
1	1	0.146
1	1	0.151
;
run;
;
proc nlmixed data = segale tech = trureg hess cov itdetails;
*Starting values;
parms B01 = -1.7, G01 = -3.4, Q = -1.2, Q1 = -0.55, Q2 = 0.47;
title 'pyes on partit and elicit (example 1 code), mmixture model';
*partit is the partition prime condition;
*elicit is whether the judged probabilities were imprecise or precise;
*Composition submodel;
C = 1/(1+EXP(Q + Q1*partit + Q2*elicit));
*Location submodel;
M1 = EXP(B01)/(1+EXP(B01));
M2 = 1/2;
*Dispersion submodel;
PHI1 = EXP(-G01);
PHI2 = EXP(8);
ll = log(C*PDF('BETA',pyes, M1*PHI1, PHI1 - M1*PHI1,0,1)+(1-C)* PDF('BETA',pyes, M2*PHI2, PHI2 - M2*PHI2,0,1));
model pyes ~ general(ll);
;
run;