全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
3884 7
2014-06-24
悬赏 30 个论坛币 已解决
问题如上,
现有一个i*j矩阵h,以h的每一个元素为积分上限,求积分,并且结果也存在一个i*j的矩阵中
假设i为100,j为50
for (i in 1:100){
  for (j in 1:50){

Ef <- integrate(dnorm,-Inf, h[i,j], mean=0, sd=2)

j <- j+1}
i <- i+1
}
问题,是不是Ef变成一个向量了?如何让其结果为 i*j 的矩阵呢?问题很简单,求大神指导菜鸟

最佳答案

二维码

扫码加我 拉你入群

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

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

全部回复
2014-6-24 15:36:12
Ef<-matrix(0,nrow=500,ncol=50)
for (i in 1:100){
  for (j in 1:50){
    Ef[i,j]<-integrate(dnorm,-Inf, h[i,j])
}
}
二维码

扫码加我 拉你入群

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

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

2014-6-25 18:15:11
integrate(dnorm,-Inf, h[i,j], mean=0, sd=2)表达式好像不对
二维码

扫码加我 拉你入群

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

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

2014-6-25 18:15:45
可能你要自己写被积函数表达式
二维码

扫码加我 拉你入群

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

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

2014-6-26 11:26:07
lww1993 发表于 2014-6-25 18:15
integrate(dnorm,-Inf, h, mean=0, sd=2)表达式好像不对
这个积分表达式还是对的,是对正太分布密度函数求积分,这个正太分布函数均值为0,方差是2,之所以这么定义,是后面还有计算方差会一改变。这个积分可以算出来结果,就是单个结果而已。明白了对Ef矩阵的定义要在循环前,我定义在循环内了!晚上回家尝试运行下!马上给分、现在在车上,今天一天在外面!非常感谢!
二维码

扫码加我 拉你入群

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

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

2014-6-27 11:38:31
运行了一下,打印出来的ef矩阵是0矩阵啊,数值都是0,为什么计算数值代不进去呢?
二维码

扫码加我 拉你入群

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

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

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

说点什么

分享

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