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

st: toeplitz


From   David Airey <[email protected]>
To   [email protected]
Subject   st: toeplitz
Date   Tue, 20 May 2003 10:38:35 -0500

Can Stata estimate a model assuming the toeplitz variance-covariance structure (type=TOEP(1) in SAS Proc Mixed):

s^2 s1 s2 s3
s1 s^2 s1 s2
s2 s1 s^2 s1
s3 s2 s1 s^2

I have a data coming from a 13x13 matrix of "treatment conditions" where the row and column levels are the same, e.g.,

1 2 3 4 5
1
2
3
4
5

In SAS, we might use the following model and code:

u + Alpha_1(i) + Alpha_2(i) + M_{i}+ e_m(i) where
u = grand mean
Alpha = polygenic effects random (13 of them)
M = QTL effect fixed
e_m = random error

proc mixed data=a;
model trait = m1;
random alpha1-alpha13/type=TOEP(1);

What in Stata might come close to this model in terms of the command and covariance?

Thanks for any tips.

*
* For searches and help try:
* http://www.stata.com/support/faqs/res/findit.html
* http://www.stata.com/support/statalist/faq
* http://www.ats.ucla.edu/stat/stata/




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