关于R - code，又晕又急！

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

set.seed(1)
p <- numeric(5200)
x <- (ifelse (rbinom (5200,1,0.5) == 1,1.005,0.995))  
for (i in 1:5200) {
    if (i==1) {
    p[i] <- 10*x[i]
}
    else {
    p[i] = x[i]*p[i-1]
}
}
plot(p, type="l")


s <- numeric()
n <- numeric()
share.v <- numeric()
j <- exp(0.07)/(260*20)-1
for (i in 1:40) {
   
    if (i==1) {
    s[i] <- 0
    n[i] <- 100
    share.v[i] <- p[i] * n[i-1]
}
   
    else {
         if (p[i] >= 10.1 && p[i] < 10.1 && share.v > 0) {
         s[i]=p[1]*n[i-1]*(1-0.001)
}
         else {
              if (p[i] <= 9.9 && p[i-1] < 9.9 && s[i-1] > 0） {
              n[i]=s[i-1]/p[i]
}
               else {
                    p[i] = p[i-1] * x[i]
                    s[i] = s[i-1] * (1+j)
}
}
}
}


到底哪错了呀？跪求正解！！！

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：code COD ODE seed see code

帮你顶顶，这块不熟悉