求救，模拟纵向数据，各位大侠帮帮忙！

小弟想模拟纵向数据，有连续模拟和离散模拟两个（见图片中），有做过的吗，SAS或R均可，感激涕零……

不是很明白连续的要求，以下给出的是离散的
/*discrete*/
proc iml;
N=50;T=4;
bt1=0.5;bt0=0;
seed=j(T,N,0);
mu=j(T,N,0);
sig=i(T);/*不清楚方差矩阵的构造，这里以单位矩阵代替*/
Z=sig*rannor(seed)+mu;/*构造多元正态分布*/
X=ranuni(seed)*2-1;
Y=t(Z<=(bt0+bt1*X));
print Y;
quit;

各位大侠，大牛们帮帮忙，解决后论坛币会追加（大家都是抱着助人为乐的心，说这个可能不好，但只有这样才能表达感激之情）

ziyenano 发表于 2012-10-27 12:12
不是很明白连续的要求，以下给出的是离散的
/*discrete*/
proc iml;

你好，大侠，谢谢你抽时间回答我的问题，很感激，我改了一下那个连续的，麻烦你再看看……

proc iml;
N=50;T=10;
bt0=1;bt1=-1;
sig1=1;
sig2=2;
k=1:T;
seed=j(N,T,0);
X=sig1*rannor(seed)+j(N,1,1)*k;
mu=bt0+bt1*X;
Y=sig2*rannor(seed)+mu;
print Y;
quit;

ziyenano 发表于 2012-10-27 12:12
不是很明白连续的要求，以下给出的是离散的
/*discrete*/
proc iml;

你好，在哪里指定相关性，比如，这个是50个个体，每个个体有4个观测，想让这四个观测间的相关性符合如下的矩阵：
1.00  .     .    .

.70 1.00  .   .

.50  .60 1.00
.30  .40 0.50 1.00  

250603730 发表于 2012-10-27 15:33
你好，在哪里指定相关性，比如，这个是50个个体，每个个体有4个观测，想让这四个观测间的相关性符合如下的 ...

改变sig矩阵；
根据协方差阵
cov=
{1.00  .70    .50    .30,
  .70  1.00    .60   .40,
  .50    .60  1.00   .50,
  .30    .40    .50 1.00  };
计算sig矩阵，要求sig*sig'=cov，
进行cholesky分解，
sig=t(root(cov));
其余代码一样。

ziyenano 发表于 2012-10-27 14:35
proc iml;
N=50;T=10;
bt0=1;bt1=-1;

proc iml;
N=50;T=10;
bt0=1;bt1=-1;
sig1=1;
sig2=1;/*不好意思，上面写错了*/
k=1:T;
seed=j(N,T,0);
X=sig1*rannor(seed)+j(N,1,1)*k;
mu=bt0+bt1*X;
Y=sig2*rannor(seed)+mu;
print Y;
quit;

ziyenano 发表于 2012-10-28 16:38
proc iml;
N=50;T=10;
bt0=1;bt1=-1;

你好，再问两个问题：1.通过你的程序得到的y的相关性并不是指定的相关性，我在这儿看到了一段程序http://blog.sina.com.cn/s/blog_5d3b177c0100ffeb.html，麻烦你看一下，这个和你的有啥不同。2.我看的文献上说“模拟运行200次”，这个怎么实现，用循环语句吗？

250603730 发表于 2012-11-22 19:46
你好，再问两个问题：1.通过你的程序得到的y的相关性并不是指定的相关性，我在这儿看到了一段程序http:// ...

能不能把所有要求整合一下啊，时间长了，我也有点混了