求外积的问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

想画一个gamma分布的极大似然值的3D图，

x <- rgamma(100,shape=3,scale=10)
sh<-sc <-seq(1,15)
like <- function(shx,scx) -sum(dgamma(x,shape=shx,scale=scx,log=T)) #外积函数，极大似然值
zz<-outer(sh,sc,like) #到这一步就不行了

persp(sh,sc,zz,theta=30,phi=15)


重新写了一段，可以运行
x <- rgamma(100,shape=3,scale=10)
sh<-sc <-seq(1,15,0.1)
zz <-matrix(rep(0,141*141),ncol=141)
#like <- function(shx,scx) sum(dgamma(x,shape=shx,scale=scx,log=T))
for (i in sh){
  for ( j in sc)
    zz[i,j]=sum(dgamma(x,shape=i,scale=j,log=T))

}

#zz<-outer(sh,sc)

persp(sh,sc,zz,theta=30,phi=30,expand=0.7,col="lightblue")

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：gamma 极大似然 3D图 GAM 似然值 问题

你自己已经搞定了，是吧？

没有，上面那段程序不行

在zz的上面加上
like <- Vectorize(like)

谢谢啦，我想再问一下，为什么下面这段程序能运行：

#做3D图

x <- y <- seq(-2*pi,2*pi,pi/15)
f <- function(x,y) sin(x)*sin(y)
z <- outer(x,y,f)
persp(x,y,z,theta=0,phi=30,expand=0.7,col="lightblue")

Vectorize到底起什么作用？

firelife 发表于 2012-4-15 08:33
谢谢啦，我想再问一下，为什么下面这段程序能运行：

#做3D图

因为sin已经是Vectorize过的函数了
比如sin(1:100)会得到100个结果

明白了，这个东西还是有点复杂的。谢谢啦