$SIZES LIM6=500 LTH=-200
$PROBLEM 
$INPUT ID STUDYID TIME DV TEST ITEM REP MMAX OCC MDV 
; Explanation of the columns in the dataset
; ID = ID number
; STUDYID = study number
; TIME = time in days since start of study
; DV = patients response to ADAS-cog assessment item
;		binary items:  0 = succeeded task, 1 = failed task
;		word recall/recognition: X = number of words not recalled/recognized
;		rater assessed items: 0 - 5 = category (no impairment - severely impaired)
; TEST = ADAS-cog subcomponent (see specification document for coding)
; ITEM = subcomponent item (see specification document for coding)
; REP = repetition of the word recall test
; MMAX = maximal DV value for current item 
; OCC = occassion
; MDV = missing variable flag

$DATA Simulated_adascog_irt.csv IGNORE=@
$ABBR DECLARE INTEGER IOCC
$ABBR DECLARE ADASARR(8)
$ABBR DECLARE TESTSUM(8)
$PRED 

;-------------------------------------------------------------------------------------- 
;---------------change of cognitive disability over time-------------------------------
	BASELINE=THETA(1)+ETA(1)
	SLP=THETA(2)+ETA(2)
	;ALZ = cognitive disability
	ALZ=BASELINE+SLP*TIME/365.242

;---------------constants to select model type----------------------------------------
	MODEL=0
	BINOMIAL=1
	OC=2
	BINOMCENSOR=3

	;set number of repetitions for word based items
	N=MMAX
;--------------assignment of item parameters------------------------------------------
	
	;commands subcomponent
	IF(TEST.EQ.1.AND.ITEM.EQ.1) THEN
		MODEL=BINOMIAL
		UPP=1
		DIS=THETA(3)
		DIF=THETA(4)
		GUE=THETA(5)
	ENDIF
	 
	IF(TEST.EQ.1.AND.ITEM.EQ.2) THEN
		MODEL=BINOMIAL
		UPP=1
		DIS=THETA(6)
		DIF=THETA(7)
		GUE=THETA(8)
	ENDIF
	 
	IF(TEST.EQ.1.AND.ITEM.EQ.3) THEN
		MODEL=BINOMIAL
		UPP=1
		DIS=THETA(9)
		DIF=THETA(10)
		GUE=THETA(11)
	ENDIF
	 
	IF(TEST.EQ.1.AND.ITEM.EQ.4) THEN
		MODEL=BINOMIAL
		UPP=1
		DIS=THETA(12)
		DIF=THETA(13)
		GUE=THETA(14)
	ENDIF
	 
	IF(TEST.EQ.1.AND.ITEM.EQ.5) THEN
		MODEL=BINOMIAL
		UPP=1
		DIS=THETA(15)
		DIF=THETA(16)
		GUE=THETA(17)
	ENDIF
	
	;comprehension subcomponent
	IF(TEST.EQ.2) THEN
		MODEL=OC
		DIS=THETA(18)
		DIF1=THETA(19)
		DIF2=THETA(20)
		DIF3=THETA(21)
		DIF4=THETA(22)
		DIF5=THETA(23) 
	ENDIF
	
	;construction subcomponent
	IF(TEST.EQ.3.AND.ITEM.EQ.1) THEN
		MODEL=BINOMIAL
		UPP=1
		DIS=THETA(24)
		DIF=THETA(25)
		GUE=THETA(26)
	ENDIF
	 
	IF(TEST.EQ.3.AND.ITEM.EQ.2) THEN
		MODEL=BINOMIAL
		UPP=1
		DIS=THETA(27)
		DIF=THETA(28)
		GUE=THETA(29)
	ENDIF
	 
	IF(TEST.EQ.3.AND.ITEM.EQ.3) THEN
		MODEL=BINOMIAL
		UPP=1
		DIS=THETA(30)
		DIF=THETA(31)
		GUE=THETA(32)
	ENDIF
	 
	IF(TEST.EQ.3.AND.ITEM.EQ.4) THEN
		MODEL=BINOMIAL
		UPP=1
		DIS=THETA(33)
		DIF=THETA(34)
		GUE=THETA(35)
	ENDIF
	
	;ideational praxis subcomponent
	IF(TEST.EQ.6.AND.ITEM.EQ.1) THEN
		MODEL=BINOMIAL
		UPP=1
		DIS=THETA(36)
		DIF=THETA(37)
		GUE=THETA(38)
	ENDIF
	 
	IF(TEST.EQ.6.AND.ITEM.EQ.2) THEN
		MODEL=BINOMIAL
		UPP=1
		DIS=THETA(39)
		DIF=THETA(40)
		GUE=THETA(41)
	ENDIF
	 
	IF(TEST.EQ.6.AND.ITEM.EQ.3) THEN
		MODEL=BINOMIAL
		UPP=1
		DIS=THETA(42)
		DIF=THETA(43)
		GUE=THETA(44)
	ENDIF
	 
	IF(TEST.EQ.6.AND.ITEM.EQ.4) THEN
		MODEL=BINOMIAL
		UPP=1
		DIS=THETA(45)
		DIF=THETA(46)
		GUE=THETA(47)
	ENDIF
	 
	IF(TEST.EQ.6.AND.ITEM.EQ.5) THEN
		MODEL=BINOMIAL
		UPP=1
		DIS=THETA(48)
		DIF=THETA(49)
		GUE=THETA(50)
	ENDIF
	
	;naming objects & fingers subcomponent
	IF(TEST.EQ.7.AND.ITEM.EQ.1) THEN
		MODEL=BINOMIAL
		UPP=1
		DIS=THETA(51)
		DIF=THETA(52)
		GUE=THETA(53)
	ENDIF
	 
	IF(TEST.EQ.7.AND.ITEM.EQ.2) THEN
		MODEL=BINOMIAL
		UPP=1
		DIS=THETA(54)
		DIF=THETA(55)
		GUE=THETA(56)
	ENDIF
	 
	IF(TEST.EQ.7.AND.ITEM.EQ.3) THEN
		MODEL=BINOMIAL
		UPP=1
		DIS=THETA(57)
		DIF=THETA(58)
		GUE=THETA(59)
	ENDIF
	 
	IF(TEST.EQ.7.AND.ITEM.EQ.4) THEN
		MODEL=BINOMIAL
		UPP=1
		DIS=THETA(60)
		DIF=THETA(61)
		GUE=THETA(62)
	ENDIF
	 
	IF(TEST.EQ.7.AND.ITEM.EQ.5) THEN
		MODEL=BINOMIAL
		UPP=1
		DIS=THETA(63)
		DIF=THETA(64)
		GUE=THETA(65)
	ENDIF
	 
	IF(TEST.EQ.7.AND.ITEM.EQ.6) THEN
		MODEL=BINOMIAL
		UPP=1
		DIS=THETA(66)
		DIF=THETA(67)
		GUE=THETA(68)
	ENDIF
	 
	IF(TEST.EQ.7.AND.ITEM.EQ.7) THEN
		MODEL=BINOMIAL
		UPP=1
		DIS=THETA(69)
		DIF=THETA(70)
		GUE=THETA(71)
	ENDIF
	 
	IF(TEST.EQ.7.AND.ITEM.EQ.8) THEN
		MODEL=BINOMIAL
		UPP=1
		DIS=THETA(72)
		DIF=THETA(73)
		GUE=THETA(74)
	ENDIF
	 
	IF(TEST.EQ.7.AND.ITEM.EQ.9) THEN
		MODEL=BINOMIAL
		UPP=1
		DIS=THETA(75)
		DIF=THETA(76)
		GUE=THETA(77)
	ENDIF
	 
	IF(TEST.EQ.7.AND.ITEM.EQ.10) THEN
		MODEL=BINOMIAL
		UPP=1
		DIS=THETA(78)
		DIF=THETA(79)
		GUE=THETA(80)
	ENDIF
	 
	IF(TEST.EQ.7.AND.ITEM.EQ.11) THEN
		MODEL=BINOMIAL
		UPP=1
		DIS=THETA(81)
		DIF=THETA(82)
		GUE=THETA(83)
	ENDIF
	 
	IF(TEST.EQ.7.AND.ITEM.EQ.12) THEN
		MODEL=BINOMIAL
		UPP=1
		DIS=THETA(84)
		DIF=THETA(85)
		GUE=THETA(86)
	ENDIF
	 
	IF(TEST.EQ.7.AND.ITEM.EQ.13) THEN
		MODEL=BINOMIAL
		UPP=1
		DIS=THETA(87)
		DIF=THETA(88)
		GUE=THETA(89)
	ENDIF
	 
	IF(TEST.EQ.7.AND.ITEM.EQ.14) THEN
		MODEL=BINOMIAL
		UPP=1
		DIS=THETA(90)
		DIF=THETA(91)
		GUE=THETA(92)
	ENDIF
	 
	IF(TEST.EQ.7.AND.ITEM.EQ.15) THEN
		MODEL=BINOMIAL
		UPP=1
		DIS=THETA(93)
		DIF=THETA(94)
		GUE=THETA(95)
	ENDIF
	 
	IF(TEST.EQ.7.AND.ITEM.EQ.16) THEN
		MODEL=BINOMIAL
		UPP=1
		DIS=THETA(96)
		DIF=THETA(97)
		GUE=THETA(98)
	ENDIF
	 
	IF(TEST.EQ.7.AND.ITEM.EQ.17) THEN
		MODEL=BINOMIAL
		UPP=1
		DIS=THETA(99)
		DIF=THETA(100)
		GUE=THETA(101)
	ENDIF
	
	;orientation subcomponent
	IF(TEST.EQ.9.AND.ITEM.EQ.1) THEN
		MODEL=BINOMIAL
		UPP=1
		DIS=THETA(102)
		DIF=THETA(103)
		GUE=THETA(104)
	ENDIF
	 
	IF(TEST.EQ.9.AND.ITEM.EQ.2) THEN
		MODEL=BINOMIAL
		UPP=1
		DIS=THETA(105)
		DIF=THETA(106)
		GUE=THETA(107)
	ENDIF
	 
	IF(TEST.EQ.9.AND.ITEM.EQ.3) THEN
		MODEL=BINOMIAL
		UPP=1
		DIS=THETA(108)
		DIF=THETA(109)
		GUE=THETA(110)
	ENDIF
	 
	IF(TEST.EQ.9.AND.ITEM.EQ.4) THEN
		MODEL=BINOMIAL
		UPP=1
		DIS=THETA(111)
		DIF=THETA(112)
		GUE=THETA(113)
	ENDIF
	 
	IF(TEST.EQ.9.AND.ITEM.EQ.5) THEN
		MODEL=BINOMIAL
		UPP=1
		DIS=THETA(114)
		DIF=THETA(115)
		GUE=THETA(116)
	ENDIF
	 
	IF(TEST.EQ.9.AND.ITEM.EQ.6) THEN
		MODEL=BINOMIAL
		UPP=1
		DIS=THETA(117)
		DIF=THETA(118)
		GUE=THETA(119)
	ENDIF
	 
	IF(TEST.EQ.9.AND.ITEM.EQ.7) THEN
		MODEL=BINOMIAL
		UPP=1
		DIS=THETA(120)
		DIF=THETA(121)
		GUE=THETA(122)
	ENDIF
	 
	IF(TEST.EQ.9.AND.ITEM.EQ.8) THEN
		MODEL=BINOMIAL
		UPP=1
		DIS=THETA(123)
		DIF=THETA(124)
		GUE=THETA(125)
	ENDIF
	
	;spoken language subcomponent
	IF(TEST.EQ.10) THEN
		MODEL=OC
		DIS=THETA(126)
		DIF1=THETA(127)
		DIF2=THETA(128)
		DIF3=THETA(129)
		DIF4=THETA(130)
		DIF5=THETA(131) 
	ENDIF
	
	;remembering subcomponent
	IF(TEST.EQ.11) THEN
		MODEL=OC
		DIS=THETA(132)
		DIF1=THETA(133)
		DIF2=THETA(134)
		DIF3=THETA(135)
		DIF4=THETA(136)
		DIF5=THETA(137) 
	ENDIF
	
	;word finding subcomponent
	IF(TEST.EQ.12) THEN
		MODEL=OC
		DIS=THETA(138)
		DIF1=THETA(139)
		DIF2=THETA(140)
		DIF3=THETA(141)
		DIF4=THETA(142)
		DIF5=THETA(143) 
	ENDIF
	
	
	;word recall subcomponent
	IF(TEST.EQ.13.AND.REP.EQ.1) THEN
		MODEL=BINOMIAL
		UPP=1
		DIS=THETA(144)
		DIF=THETA(145)
		GUE=THETA(146)
	ENDIF
	

	IF(TEST.EQ.13.AND.REP.EQ.2) THEN
		MODEL=BINOMIAL
		UPP=1
		DIS=THETA(147)
		DIF=THETA(148)
		GUE=THETA(149)
	ENDIF
	

	IF(TEST.EQ.13.AND.REP.EQ.3) THEN
		MODEL=BINOMIAL
		UPP=1
		DIS=THETA(150)
		DIF=THETA(151)
		GUE=THETA(152)
	ENDIF
			
	;word recognition subcomponent
	RDV=DV
	IF(TEST.EQ.14) THEN
		UPP=THETA(153)
		DIS=THETA(154)
		DIF=THETA(155)
		GUE=THETA(156)
		MODEL=BINOMIAL
		N=24
	ENDIF

	
;---------------------------binary and binomial data model----------------------------- 
	IF(MODEL.EQ.BINOMIAL.OR.MODEL.EQ.BINOMCENSOR) P=GUE+(UPP-GUE)*EXP(DIS*(ALZ-DIF))/(1+EXP(DIS*(ALZ-DIF)))
	IF(P.LT.1E-16) P=1E-16
	IF(P.GT.(1-1E-16))P=1-1E-16
	IF(MODEL.EQ.BINOMIAL.AND.RDV.EQ.0) Y=-2*N*LOG(1-P)
	IF(MODEL.EQ.BINOMIAL.AND.RDV.EQ.N) Y=-2*N*LOG(P)
	IF(MODEL.EQ.BINOMIAL.AND.RDV.GT.0.AND.RDV.LT.N) THEN
"LNFAC=LGAMMA(N+1)
"LKFAC=LGAMMA(RDV+1)
"LNKFAC=LGAMMA(N-RDV+1)
	ENDIF
	IF(MODEL.EQ.BINOMIAL.AND.RDV.GT.0.AND.RDV.LT.N) Y=-2*(RDV*LOG(P)+(N-RDV)*LOG(1-P)+LNFAC-LKFAC-LNKFAC)
	
;---------------------------ordered categorical data model----------------------------- 	
	IF(MODEL.EQ.OC) THEN
		DIF1=DIF1 
		DIF2=DIF1+DIF2
		DIF3=DIF2+DIF3
		DIF4=DIF3+DIF4
		DIF5=DIF4+DIF5
		
		PGE1=EXP(DIS*(ALZ-DIF1))/(1+EXP(DIS*(ALZ-DIF1)))
		PGE2=EXP(DIS*(ALZ-DIF2))/(1+EXP(DIS*(ALZ-DIF2)))
		PGE3=EXP(DIS*(ALZ-DIF3))/(1+EXP(DIS*(ALZ-DIF3)))
		PGE4=EXP(DIS*(ALZ-DIF4))/(1+EXP(DIS*(ALZ-DIF4)))
		PGE5=EXP(DIS*(ALZ-DIF5))/(1+EXP(DIS*(ALZ-DIF5)))
  
		P0=1-PGE1
		P1=PGE1-PGE2
		P2=PGE2-PGE3
		P3=PGE3-PGE4
		P4=PGE4-PGE5
		P5=PGE5 
	ENDIF

	IF(MODEL.EQ.OC.AND.DV.EQ.0) P=P0
	IF(MODEL.EQ.OC.AND.DV.EQ.1) P=P1
	IF(MODEL.EQ.OC.AND.DV.EQ.2) P=P2
	IF(MODEL.EQ.OC.AND.DV.EQ.3) P=P3 
	IF(MODEL.EQ.OC.AND.DV.EQ.4) P=P4
	IF(MODEL.EQ.OC.AND.DV.EQ.5) P=P5
	
	IF(P.LT.1E-16) P=1E-16	
	IF(P.GT.(1-1E-16)) P=1-1E-16
	IF(MODEL.EQ.OC) Y=-2*LOG(P)
;---------------------------simulation of item responses--------------------------- 	
	IF(ICALL.EQ.4) THEN
		ISIM=IREP
		SDV=0
		IF(MODEL.EQ.OC) THEN
			CALL RANDOM (2,R)
			SDV=0
			IF(R.LT.PGE1) SDV=1
			IF(R.LT.PGE2) SDV=2
			IF(R.LT.PGE3) SDV=3
			IF(R.LT.PGE4) SDV=4
			IF(R.LT.PGE5) SDV=5
		ENDIF
		IF(MODEL.EQ.BINOMIAL) THEN
			I=0
			SDV=0
			DO WHILE (I.LT.N)
				CALL RANDOM (2,R)
			IF(P.GT.R) SDV=SDV+1
			I=I+1
			END DO
		ENDIF
		DV=SDV
	ENDIF
;---------------------------calculation of ADAS-cog score---------------------------
	IF(NEWIND.NE.2) THEN
		ADASARR(1)=0
		ADASARR(2)=0
		ADASARR(3)=0
		ADASARR(4)=0
		ADASARR(5)=0
		ADASARR(6)=0
		ADASARR(7)=0
		ADASARR(8)=0

		TESTSUM(1)=0
		TESTSUM(2)=0
		TESTSUM(3)=0
		TESTSUM(4)=0
		TESTSUM(5)=0
		TESTSUM(6)=0
		TESTSUM(7)=0
		TESTSUM(8)=0
	ENDIF
	IOCC=0
	"IOCC=INT(OCC)
	IF(TEST.NE.0.AND.TEST.NE.7.AND.TEST.NE.13.AND.TEST.NE.14) THEN
		ADASARR(IOCC)=ADASARR(IOCC) + DV
	ENDIF

	IF(ITEM.EQ.1.AND.REP.EQ.1) TESTSUM(IOCC) = 0

	TESTSUM(IOCC) = TESTSUM(IOCC) + DV
	 
	IF(TEST.EQ.7.AND.ITEM.EQ.17) THEN
		IF(TESTSUM(IOCC).GT.2) ADASARR(IOCC)=ADASARR(IOCC) + 1
		IF(TESTSUM(IOCC).GT.5) ADASARR(IOCC)=ADASARR(IOCC) + 1
		IF(TESTSUM(IOCC).GT.8) ADASARR(IOCC)=ADASARR(IOCC) + 1
		IF(TESTSUM(IOCC).GT.11) ADASARR(IOCC)=ADASARR(IOCC) + 1
		IF(TESTSUM(IOCC).GT.14) ADASARR(IOCC)=ADASARR(IOCC) + 1
	ENDIF

	IF(TEST.EQ.13.AND.REP.EQ.3) THEN
		ADASARR(IOCC)=ADASARR(IOCC) + TESTSUM(IOCC)/3
	ENDIF

	IF(TEST.EQ.14)THEN
		IF(DV.GT.12) THEN
			ADASARR(IOCC) = ADASARR(IOCC) + 12
		ELSE
			ADASARR(IOCC) = ADASARR(IOCC) + DV
		ENDIF
	ENDIF
	ADASCOG=ADASARR(IOCC)
;---------------------------estimation task -------------------------------------------
$EST METHOD=1 LAPLACE -2LL PRINT=1 MSFO=msf MCETA=4 MAXEVAL=0
$TABLE ID STUDYID TIME DV TEST ITEM REP MMAX OCC MDV ADASCOG
NOAPPEND ONEHEADER NOPRINT FILE=all_data.csv FORMAT=,1PE11.4
;---------------------------random effects for slope and baseline---------------------- 
$OMEGA BLOCK(2) 4.57E-01	;baseline variance
1.37E-01 1.27E-01			;slope variance
;---------------------------fixed effects for slope and baseline-----------------------
$THETA (-INF,7.67E-01,INF)		;baseline
$THETA (-INF,3.27E-01,INF)		;slope
;---------------------------item parameters -------------------------------------------
$THETA 1.34063 FIX;T1I1DIS
$THETA 4.68236 FIX;T1I1DIF
$THETA 0.00350057 FIX;T1I1GUE
$THETA 2.24119 FIX;T1I2DIS
$THETA 3.05678 FIX;T1I2DIF
$THETA 0.00498419 FIX;T1I2GUE
$THETA 1.9185 FIX;T1I3DIS
$THETA 2.21697 FIX;T1I3DIF
$THETA 0.00849066 FIX;T1I3GUE
$THETA 2.06947 FIX;T1I4DIS
$THETA 1.96115 FIX;T1I4DIF
$THETA 0.006864 FIX;T1I4GUE
$THETA 1.55198 FIX;T1I5DIS
$THETA 1.59181 FIX;T1I5DIF
$THETA 0.0813223 FIX;T1I5GUE
$THETA 1.53738 FIX;T2I1DIS
$THETA 1.16952 FIX;T2I1DIFOC
$THETA 1.03948 FIX;T2I2DIFOC
$THETA 0.799086 FIX;T2I3DIFOC
$THETA 1.57975 FIX;T2I4DIFOC
$THETA 1.37699 FIX;T2I5DIFOC
$THETA 1.20071 FIX;T3I1DIS
$THETA 4.39975 FIX;T3I1DIF
$THETA 0.00120455 FIX;T3I1GUE
$THETA 1.6495 FIX;T3I2DIS
$THETA 2.22316 FIX;T3I2DIF
$THETA 0.00149816 FIX;T3I2GUE
$THETA 1.09901 FIX;T3I3DIS
$THETA 2.34689 FIX;T3I3DIF
$THETA 2.07857e-05 FIX;T3I3GUE
$THETA 1.40713 FIX;T3I4DIS
$THETA 0.168809 FIX;T3I4DIF
$THETA 0.212824 FIX;T3I4GUE
$THETA 2.38741 FIX;T6I1DIS
$THETA 2.64408 FIX;T6I1DIF
$THETA 0.00211518 FIX;T6I1GUE
$THETA 2.43999 FIX;T6I2DIS
$THETA 2.65469 FIX;T6I2DIF
$THETA 0.00058608 FIX;T6I2GUE
$THETA 1.69972 FIX;T6I3DIS
$THETA 2.32036 FIX;T6I3DIF
$THETA 0.0181654 FIX;T6I3GUE
$THETA 2.44295 FIX;T6I4DIS
$THETA 1.34097 FIX;T6I4DIF
$THETA 0.0119545 FIX;T6I4GUE
$THETA 1.46274 FIX;T6I5DIS
$THETA 2.18147 FIX;T6I5DIF
$THETA 0.0057591 FIX;T6I5GUE
$THETA 1.81336 FIX;T7I1DIS
$THETA 2.56472 FIX;T7I1DIF
$THETA 2.14515e-10 FIX;T7I1GUE
$THETA 1.5894 FIX;T7I2DIS
$THETA 1.26058 FIX;T7I2DIF
$THETA 0.0484782 FIX;T7I2GUE
$THETA 1.04474 FIX;T7I3DIS
$THETA 0.232459 FIX;T7I3DIF
$THETA 0.0258602 FIX;T7I3GUE
$THETA 1.05361 FIX;T7I4DIS
$THETA 1.57735 FIX;T7I4DIF
$THETA 2.8266e-07 FIX;T7I4GUE
$THETA 1.19563 FIX;T7I5DIS
$THETA 1.52464 FIX;T7I5DIF
$THETA 0.116568 FIX;T7I5GUE
$THETA 1.94142 FIX;T7I6DIS
$THETA 3.49457 FIX;T7I6DIF
$THETA 0.00237255 FIX;T7I6GUE
$THETA 2.3115 FIX;T7I7DIS
$THETA 2.93832 FIX;T7I7DIF
$THETA 0.00370298 FIX;T7I7GUE
$THETA 3.35186 FIX;T7I8DIS
$THETA 3.1271 FIX;T7I8DIF
$THETA 0.0730606 FIX;T7I8GUE
$THETA 1.89587 FIX;T7I9DIS
$THETA 2.36795 FIX;T7I9DIF
$THETA 2.05765e-09 FIX;T7I9GUE
$THETA 1.6237 FIX;T7I10DIS
$THETA 1.9243 FIX;T7I10DIF
$THETA 0.0158487 FIX;T7I10GUE
$THETA 1.42877 FIX;T7I11DIS
$THETA 2.82543 FIX;T7I11DIF
$THETA 0.00510635 FIX;T7I11GUE
$THETA 2.30901 FIX;T7I12DIS
$THETA 2.10143 FIX;T7I12DIF
$THETA 0.000364909 FIX;T7I12GUE
$THETA 1.98626 FIX;T7I13DIS
$THETA 2.9834 FIX;T7I13DIF
$THETA 0.00195687 FIX;T7I13GUE
$THETA 2.21178 FIX;T7I14DIS
$THETA 1.57519 FIX;T7I14DIF
$THETA 0.0044263 FIX;T7I14GUE
$THETA 2.25038 FIX;T7I15DIS
$THETA 1.25083 FIX;T7I15DIF
$THETA 2.94378e-09 FIX;T7I15GUE
$THETA 0.855837 FIX;T7I16DIS
$THETA 2.06231 FIX;T7I16DIF
$THETA 0.0127098 FIX;T7I16GUE
$THETA 0.479465 FIX;T7I17DIS
$THETA 2.05567 FIX;T7I17DIF
$THETA 9.60659e-07 FIX;T7I17GUE
$THETA 1.54909 FIX;T9I1DIS
$THETA 3.88607 FIX;T9I1DIF
$THETA 9.23772e-09 FIX;T9I1GUE
$THETA 2.12315 FIX;T9I2DIS
$THETA 0.987855 FIX;T9I2DIF
$THETA 1.67205e-09 FIX;T9I2GUE
$THETA 1.84366 FIX;T9I3DIS
$THETA 0.901343 FIX;T9I3DIF
$THETA 0.00709225 FIX;T9I3GUE
$THETA 1.61835 FIX;T9I4DIS
$THETA 1.31643 FIX;T9I4DIF
$THETA 1.143e-07 FIX;T9I4GUE
$THETA 2.56301 FIX;T9I5DIS
$THETA 0.612652 FIX;T9I5DIF
$THETA 0.00950114 FIX;T9I5GUE
$THETA 2.30238 FIX;T9I6DIS
$THETA 0.969427 FIX;T9I6DIF
$THETA 0.0138761 FIX;T9I6GUE
$THETA 2.55149 FIX;T9I7DIS
$THETA 1.07487 FIX;T9I7DIF
$THETA 1.20167e-07 FIX;T9I7GUE
$THETA 1.18255 FIX;T9I8DIS
$THETA 1.44985 FIX;T9I8DIF
$THETA 0.00767535 FIX;T9I8GUE
$THETA 1.50851 FIX;T10I1DIS
$THETA 1.03234 FIX;T10I1DIFOC
$THETA 0.696581 FIX;T10I2DIFOC
$THETA 0.598863 FIX;T10I3DIFOC
$THETA 0.508038 FIX;T10I4DIFOC
$THETA 0.440452 FIX;T10I5DIFOC
$THETA 1.5971 FIX;T11I1DIS
$THETA 1.57006 FIX;T11I1DIFOC
$THETA 0.896699 FIX;T11I2DIFOC
$THETA 0.852933 FIX;T11I3DIFOC
$THETA 1.00695 FIX;T11I4DIFOC
$THETA 1.4643 FIX;T11I5DIFOC
$THETA 1.25435 FIX;T12I1DIS
$THETA 0.539365 FIX;T12I1DIFOC
$THETA 1.40176 FIX;T12I2DIFOC
$THETA 1.05259 FIX;T12I3DIFOC
$THETA 1.05351 FIX;T12I4DIFOC
$THETA 179.341 FIX;T12I5DIFOC
$THETA 0.733913 FIX;T13I1DIS
$THETA -0.629249 FIX;T13I1DIF
$THETA 0.127288 FIX;T13I1GUE
$THETA 0.792436 FIX;T13I1R2DIS
$THETA 0.0256737 FIX;T13I1R2DIF
$THETA 5.75612e-09 FIX;T13I1R2GUE
$THETA 0.852149 FIX;T13I1R3DIS
$THETA 0.395876 FIX;T13I1R3DIF
$THETA 1.3128e-09 FIX;T13I1R3GUE
$THETA 0.433285 FIX;T14I1UPP
$THETA 1.07606 FIX;T14I1DIS
$THETA -0.0626467 FIX;T14I1DIF
$THETA 1.18921e-07 FIX;T14I1GUE