Stata The Stata listserver
[Date Prev][Date Next][Thread Prev][Thread Next][Date index][Thread index]

st: AR1


From   Dirk Nachbar <Dirk.Nachbar@dwp.gsi.gov.uk>
To   statalist <statalist@hsphsun2.harvard.edu>
Subject   st: AR1
Date   03 Oct 2005 09:43:00 +0100

Dear all

I am trying to simulate an AR(1) process but it proves to be hard in Stata. My idea was to declare a matrix and then do it in a loop, but it doesn't work. Maybe someone can help.

Dirk

------------------

clear

set matsize 800

set obs 700

gen time=_n

tsset time

gen y=invnorm(uniform())

mkmat y

local i = 1

while `i' < `obs' {

    matrix y[i+1,1]=1+0.5*y[i,1] +invnorm(uniform())

local i = `i' + 1

}

svmat y

arima y, arima(1,0,0)

---------------------------------------------------

Dirk Nachbar

Assistant Economist

Pensim2

Department for Work and Pensions

Level 4, The Adelphi

1-11 John Adam St

WC2N 6HT London

020 796 22705 ********************************************************************** This document is strictly confidential and is intended only for use by the addressee. If you are not the intended recipient, any disclosure, copying, distribution or other action taken in reliance of the information contained in this e-mail is strictly prohibited. Any views expressed by the sender of this message are not necessarily those of the Department for Work and Pensions. If you have received this transmission in error, please use the reply function to tell us and then permanently delete what you have received. Please note: Incoming and outgoing e-mail messages are routinely monitored for compliance with our policy on the use of electronic communications. **********************************************************************
The original of this email was scanned for viruses by the Government Secure Intranet (GSi) virus scanning service supplied exclusively by Energis in partnership with MessageLabs.

On leaving the GSi this email was certified virus-free




© Copyright 1996–2014 StataCorp LP   |   Terms of use   |   Privacy   |   Contact us   |   What's new   |   Site index