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

250603730

3894

收藏 2012-10-26

悬赏 200 个论坛币已解决

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

连续模拟.jpg

原图尺寸 65.64 KB

连续模拟要求

离散模拟.jpg

原图尺寸 65.14 KB

离散模拟的要求

连续.jpg

原图尺寸 57.11 KB

最佳答案

ziyenano 查看完整内容

不是很明白连续的要求，以下给出的是离散的 /*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

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

全部回复

ziyenano

2012-10-26 16:04:07

不是很明白连续的要求，以下给出的是离散的
/*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;

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

250603730

2012-10-26 16:07:30

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

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

250603730

2012-10-27 14:10:44

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

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

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

ziyenano

2012-10-27 14:35:55

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;

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

250603730

2012-10-27 15:33:46

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

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

点击查看更多内容…

ziyenano

2012-10-27 21:34:47

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-28 16:38:03

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;

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

250603730

2012-11-22 19:46:31

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次”，这个怎么实现，用循环语句吗？

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

ziyenano

2012-11-23 18:14:05

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

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

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

250603730

2012-11-24 14:44:54

ziyenano 发表于 2012-11-23 18:14
能不能把所有要求整合一下啊，时间长了，我也有点混了

把论坛币给你吧，谢谢你以前的解答。这个挺复杂的，我也没太看明白，这样在帖子上也不一定说得明白。谢谢……

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

jingju11

2012-11-26 10:06:35

你最后问题是如何解决的呢？
比如E[Yij] = MUij =B*Xij,这和Yij～N(BXij,1)等价吗？也就是说每个Y具备一个可能不同的分布？若是如此，如何控制其中的相关系数呢？你的具体问题应该远比你推荐的帖子复杂。
谢谢，京剧

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

250603730

2012-11-26 11:07:45

jingju11 发表于 2012-11-26 10:06
你最后问题是如何解决的呢？
比如E[Yij] = MUij =B*Xij,这和Yij～N(BXij,1)等价吗？也就是说每个Y具备一个 ...

你好，这个问题没解决呐。我是看的一篇大论文上用模拟的纵向数据来比较两种模型。那上面有具体的要求，我这个帖子中就是截的图。大侠，你有时间的话，我把文献传给你看看怎样？数据模拟方面你做的多不多，你是什么专业的？

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

栏目导航

最佳答案

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

分享

扫码加好友，拉您进群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群