想画一个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")