全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
5360 13
2011-10-21
如题
二维码

扫码加我 拉你入群

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

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

全部回复
2011-10-21 17:34:32
比如将以下过程重复1000次,计算z的均值
m=10;c=5;a=2;b=4
k=rgamma(1,a,b)
w=runif(m)
for(i in 1:m)
{
x[i]=(w[i]^(-1/k.star)-1)^(1/c)
}
x

S=sum(2*log(1+x^c))
y=10-2*k
y.star=1-(a+m)/(S+b)
y
y.star

z=(y-y.star)^2
z
大家帮忙看看,还是用FOR语句吗?
二维码

扫码加我 拉你入群

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

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

2011-10-21 17:48:28
看到这个貌似想起了当年学C语言的时候了!
二维码

扫码加我 拉你入群

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

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

2011-10-21 18:08:04
levyangel 发表于 2011-10-21 17:34
比如将以下过程重复1000次,计算z的均值
m=10;c=5;a=2;b=4
k=rgamma(1,a,b)
用while语句,设置一个记录循环次数的变量J,每次循环后J=J+1,while(J<=1000).
二维码

扫码加我 拉你入群

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

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

2011-10-21 22:08:31
levyangel 发表于 2011-10-21 17:34
比如将以下过程重复1000次,计算z的均值
m=10;c=5;a=2;b=4
k=rgamma(1,a,b)
可以外面再套一个for循环,从1到1000
二维码

扫码加我 拉你入群

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

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

2011-10-21 23:19:05
levyangel 发表于 2011-10-21 17:34
比如将以下过程重复1000次,计算z的均值
m=10;c=5;a=2;b=4
k=rgamma(1,a,b)
看了你的程序有几个小困惑:
1.k.star是什么?
2.你循环里面减了1又求了个^(1/c),后面又对x求了个^c还加了个1,这不抵消了吗?有什么意义吗?
你先把公式理清楚,然后再解释下大致目的吧,你这个程序也可以不用循环,而用向量形式来做的,循环的效率比较低,能向量化运算的话,要尽量避免循环。
二维码

扫码加我 拉你入群

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

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

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

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群