sspace solution
!loglik = -1000000
table results
results (1,1) = "loglik"
results (1,2)= "result"
results (1,3) = "AIC"
results (1,4) = "Schwarz"
results (1,5) = "Hannan-Quinn"
!row = 1
for !i = 1 to 1000 step 1
		smpl @all
		series start = (rnd)
		
		for !j = 1 to 9 step 1
			!startvar = start (!j)
			param c(!j) !startvar
		next

smpl 06/1/1996 12/01/2013
sspace ss1
ss1.append @ename e1
ss1.append @ename e2
ss1.append @ename e3
ss1.append @ename e4
ss1.append @ename e5
ss1.append @evar var(e1) = c(1)^2
ss1.append @evar var(e2) = c(2)^2
ss1.append @evar var(e3) = c(3)^2
ss1.append @evar var(e4) = c(4)^2
ss1.append @evar var(e5) = c(5)^2
ss1.append @signal logy=cyc + pot + e1
ss1.append @state cyc = c(11)*cyc(-1) + e3
ss1.append @state pot = pot (-1) + c(7) + e4
ss1.append @signal enf = sv1 + c(8)*sv2 + c(10)*cyc + e2
ss1.append @state sv1= sv1(-1) + e5
ss1.append @state sv2 = sv1(-1)
ss1.ml(m=1000, c=1e-6)
freeze(out1)ss1.output
%st1 = @left(out1(6,1),11)

if %st1="Convergence" then
	!row = !row+1
	results(!row,1) = ss1.@logl
	results(!row,2) = out1(6,1)
	results(!row,3) = ss1.@aic
	results(!row,4) = ss1.@sc
	results(!row,5) = ss1.@hq

	for !k=1 to 9 step 1
			results(!row, 9+!k)= start(!k)
	next
	
	for !k=1 to 9 step 1
	results(!row,5+9+!k)=C(!k)
	next
	
	if ss1.@logl>!loglik then
		!loglik = ss1.@logl
		delete solution
		copy ss1 solution	
		endif
endif

delete out1
delete ss1

next
stop