The RATS Software Forum

Code: Select all

*
* UhligAcceptLR(q,steps,constrained) returns a 1 (accept) or 0 (reject)
* for the weight vector <<q>> (applied to the impulse responses in the
* global RECT[SERIES] impulses). <<steps>> is the list of responses
* constrained. <<constrained>>  is a VECT[INTEGER] with the variable
* positions being constrained. To check for near-zeros for the first
* variable at steps 80, 100 and 120, use
*
* UhligAcceptLR(q,||80,100,120||,||1||)
*
function UhligAcceptLR q steps constrained
type integer   UhligAcceptLR
type vector    q
type vect[int] steps
type vect[int] constrained
*
local integer      i k
local vector       ik
local real         value
dec   rect[series] impulses
*
dofor k = steps
   compute ik=%xt(impulses,k)*q
   do i=1,%rows(constrained)
      compute value=ik(constrained(i))
      if value>.01.or.value<-.01 {
         compute UhligAcceptLR=0
         return
     }
   end do i
end do k
compute UhligAcceptLR=1
end