-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathactib modified latent change model_ibs ipqn_bayes_free b.inp
95 lines (72 loc) · 2.98 KB
/
actib modified latent change model_ibs ipqn_bayes_free b.inp
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
Data:
File is "ACTIB long data std.dat" ;
Variable:
Names are
tcbt ibs0 ibs1 ibs2 ibs3 ipqn0 ipqn1 ipqn2 ipqn3 ipqn0s ipqn1s ipqn2s
ipqn3s ibs0s ibs1s ibs2s ibs3s;
USEVARIABLES ARE ibs0s ibs1s ibs2s ibs3s ipqn0s ipqn1s ipqn2s ipqn3s tcbt ;
Missing are all (99) ;
ANALYSIS:
estimator = bayes;
FBITERATIONS = 40000;
CHAINS = 3;
POINT = mean;
MODEL:
fm0 by ipqn0s ipqn1s ipqn2s ipqn3s (estm) ;
fm1 by ipqn1s ipqn2s ipqn3s (estm) ;
fm2 by ipqn2s ipqn3s (estm) ;
fm3 by ipqn3s@1;
ipqn0s ipqn1s ipqn2s ipqn3s (varipq);
fy0 by ibs0s ibs1s ibs2s ibs3s (esty) ;
fy1 by ibs1s ibs2s ibs3s (esty) ;
fy2 by ibs2s ibs3s (esty) ;
fy3 by ibs3s@1;
ibs0s ibs1s ibs2s ibs3s (varibs);
fm1 on tcbt*0 (a11) ;
fm2 on tcbt*0 (a12) ;
fm3 on tcbt@0 ;
fy1 on tcbt*0 (c11) ;
fy2 on tcbt*0 (c12) ;
fy3 on tcbt@0 ;
fy1 on fm1*0 (b1) ;
fy2 on fm2*0 (b2) ;
fy3 on fm3*0 (b3) ;
fm0 with fy0*0 ;
fm1 with fy1@0 ;
fm2 with fy2@0 ;
fm3 with fy3@0 ;
fy1 with fy2@0 ;
fy2 with fy3@0 ;
fy1 with fy3@0 ;
fm1 with fm2@0 ;
fm2 with fm3@0 ;
fm1 with fm3@0 ;
ipqn0s with ibs0s (cor) ;
ibs1s with ipqn1s (cor) ;
ibs2s with ipqn2s (cor) ;
ibs3s with ipqn3s (cor) ;
MODEL INDIRECT: ibs1s ind tcbt ;
ibs2s ind tcbt ;
MODEL CONSTRAINT:
NEW(TP1, INT1, TP21, TP22, INT2, TP31, TP32, INT3);
!FIRST TIME POINT
TP1 = a11*b1*1;
INT1 = TP1;
!SECOND TIME POINT
TP21 = a11*b1*esty;
TP22 = a12*b2*1;
INT2 = TP21 + TP22;
!THIRD TIME POINT
TP31 = a11*b1*esty;
TP32 = a12*b2*esty;
INT3 = TP31 + TP32;
MODEL PRIORS:
a11 ~ N (-0.516, 0.0576);
a12 ~ N (-0.497, 0.065536);
b1 ~ N (0.969, 0.1764);
b2 ~ N (1.150, 34.035556);
b3 ~ N (0.732, 99.820081);
c11 ~ N (-0.097, 0.085849);
c12 ~ N (0.052, 1.149184);
OUTPUT: CINTERVAL tech4 tech8 ;
PLOT: TYPE = PLOT2;