这个R输出的结果为什么为0？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

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

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

1. y1 <- rbinom(10,5,0.1)
   
2. y2 <- rbinom(10,5,0.4)
   
3. x <- c(y1,y2)
   
4. fir <- function(p,q,h,n){
   
5.   a1 <- 0
   
6.   a2 <- 0
   
7.   b1 <- 0
   
8.   b2 <- 0
   
9.   for (i in 0:n) {
   
10.     l <- x[i+1]
    
11.     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))
    
12.     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))
    
13.     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))
    
14.   }
    
15.   c(a1,a2,b1)
    
16. }
    
17. fir(0.1,0.4,0.5,10)

复制代码

R和一般程序语言不太一样，默认数组第一个为1位，没有0位

xwan288 发表于 2018-4-10 21:32
R和一般程序语言不太一样，默认数组第一个为1位，没有0位

谢谢！