R语言程序出现错误问题？

> ralcond<-function(d, b, t0, a1, b1)
+ {
+  grid<-seq(-10,10,.001)
+   
+   p<-numeric(length=length(grid))
+   for(i in 1:length(p))
+ {
+   p[i]<- sum( (grid[i]*(n-t0)*log(exp (b))+(grid[i]*(n-t0)-1) *log(exp(d)))-log(exp(gamma(grid[i]*(n-t0)))))+(a1-1)*log(exp(al))-al*b1
+  }
+   
+  draw<-sample(p, size = 1, prob = exp(1-p/max(p)))
+  return(draw)
+ }

```r
Error in sample.int(length(x), size, replace, prob) : 机率矢量中有NA
此外: There were 50 or more warnings (use warnings() to see the first 50)

警告信息：
1: In gamma(grid[i] * (n - t0)) : 产生了NaNs
2: In p[i] <- sum((grid[i] * (n - t0) * log(exp(b)) + (grid[i] *  ... :
  被替换的项目不是替换值长度的倍数
3: In p[i] <- sum((grid[i] * (n - t0) * log(exp(b)) + (grid[i] *  ... :
  被替换的项目不是替换值长度的倍数
4: In p[i] <- sum((grid[i] * (n - t0) * log(exp(b)) + (grid[i] *  ... :
  被替换的项目不是替换值长度的倍数
5: In p[i] <- sum((grid[i] * (n - t0) * log(exp(b)) + (grid[i] *  ... :
  被替换的项目不是替换值长度的倍数
6: In p[i] <- sum((grid[i] * (n - t0) * log(exp(b)) + (grid[i] *  ... :
  被替换的项目不是替换值长度的倍数

程序是一篇文章的框架，我把自己的内容放进去了。上面的程序只是一部分，运行的时候出现这两个错误，一直不知道怎么解决，尤其是第二个机率矢量中有NA，在百度里查不到。。希望有专家能帮我分析一下问题的原因

你好，我是有问必答小助手，非常抱歉，本次您提出的有问必答问题，技术专家团超时未为您做出解答

本次提问扣除的有问必答次数，将会以问答VIP体验卡（1次有问必答机会、商城购买实体图书享受95折优惠）的形式为您补发到账户。

​​​​因为有问必答VIP体验卡有效期仅有1天，您在需要使用的时候【私信】联系我，我会为您补发。