y1 <- rbinom(10,5,0.1)
y2 <- rbinom(10,5,0.4)
x <- c(y1,y2)
fir <- function(p,q,h,n){
  a1 <- 0
  a2 <- 0
  b1 <- 0
  b2 <- 0
  for (i in 0:n) {
    l <- x[i]
    a1 <- a1 + (h*l*p^(l-1)*(1-p)^(n-l)-h*(n-l)*p^(l)*(1-p)^(n-l-1))/(h*((p^l))*p^(n-l)+(1-h)*(q^l)*(1-q)^(n-l))
    a2 <- a2 + ((1-h)*l*q^(l-1)*(1-q)^(n-l)-(1-h)*(n-l)*q^(l)*(1-q)^(n-l-1))/(h*((p^l))*p^(n-l)+(1-h)*(q^l)*(1-q)^(n-l))
    b1 <- b1 + ((p^l)*p^(n-l) - (q^l)*(1-q)^(n-l))/(h*(p^l)*p^(n-l)+(1-h)*(q^l)*(1-q)^(n-l))
  }
  c(a1,a2,b1)
}
fir(0.1,0.4,0.5,10)

这个结果为numeric0