全部版块 我的主页
› 论坛 › 数据科学与人工智能 › 数据分析与数据科学 › R语言论坛
请问如何生成独立同分布的卡方分布或者指数分布数据
八神小弟
7811 9
2016-07-05
悬赏 10 个论坛币 已解决
比如有x1,x2,x3,...x10总共10个随机变量均服从同一卡方分布并且相互独立，并且对每个随机变量取样30，请问如何生成满足这些条件的数据？（即最后的数据是一个30*10的样本数据矩阵）

最佳答案

suzhzh 查看完整内容

If you use R, then it is quite easy. Like x
二维码

扫码加我 拉你入群

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

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

全部回复
suzhzh
2016-7-5 19:37:46
If you use R, then it is quite easy. Like
x<-matrix(nrow=30, ncol=10, byrow=FALSE)#define the structure to store the data
for (i in 1:10){x[,i]=rchisq(30, 1,1)}#x from Chisquare (1)
x#is what you wanted
二维码

扫码加我 拉你入群

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

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

suzhzh
2016-7-6 18:41:42
If you use sas, then PROC IML works as below:
PROC IML;
   call randseed(1234);
   x=j(30, 10);*30 row, 10 columns;
   call randgen (x, "CHISQ", 5 );**Chi-squre with freedom as 5, center parameter as 0;
   print x;
quit;

        
二维码

扫码加我 拉你入群

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

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

八神小弟
2016-7-6 19:44:40
suzhzh 发表于 2016-7-6 18:14
If you use R, then it is quite easy. Like
x
谢谢你的回复，看了之后我有如下问题：
（1）程序for (i in 1:30){x[,1]=rchisq(30, 1,1)}中for (i in 1:30)对于x好像没有起到循环作用
（2）这样生成的数据能满足两两独立吗？你的意思是对{x[,1]=rchisq(30, 1,1)}循环10次是吗？这样能得到这10次的数据是独立的吗？
二维码

扫码加我 拉你入群

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

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

suzhzh
2016-7-6 19:55:43
八神小弟 发表于 2016-7-6 19:44
谢谢你的回复，看了之后我有如下问题：
（1）程序for (i in 1:30){x[,1]=rchisq(30, 1,1)}中for (i in 1 ...
Typos, 1 should be i.

Yes, they are independent.
二维码

扫码加我 拉你入群

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

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

八神小弟
2016-7-6 20:37:50
suzhzh 发表于 2016-7-6 19:55
Typos, 1 should be i.

Yes, they are independent.
谢谢你的回复，问题如下
（1）for (i in 1:30){x[,1]=rchisq(30, 1,1)}是否应该改为for (i in 1:10){x[,i]=rchisq(30, 1,1)}，这样就会生成10列数据，每列数据容量为30
（2）为什么这样能保证两两独立？我不是很理解，能否给予详细的解释，谢谢。
二维码

扫码加我 拉你入群

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

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

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

微信QQ空间QQ微博
扫码加好友，拉您进群
各岗位、行业、专业交流群