如何解决从 R 转换为 python 时列表索引出错
我正在尝试转换我的 R 代码
set.seed(2021)
startn <- 10
j <- 1
maxj <- 10^5
n <- startn
time <- 0
while (n[j] > 0 & j < maxj){
time[j+1] <- time[j] + rexp(1,rate=2*n[j]^2/5)
n[j+1] <- n[j] + sample(c(-1,1),1)
j <- j+1
}
plot(time,n,xlim=c(0,max(time)),ylim=c(0,max(n)))
进入蟒蛇
n_0 = 10
j = 0
maxj = 10^5
n = [n_0]
time = [0]
while n[j] > 0 & j < maxj:
time[j+1] = time[j] + np.random.exponential(2*n[j]**2/5)
n[j+1] = n[j] + np.random.choice([-1,1])
j = j+1
但由于错误列表分配索引超出范围而失败。有人能指出我的错误吗?
解决方法
我不确定这是否是您需要的,但这解决了索引问题。请尝试:
import numpy as np
n_0 = 10
j = 0
maxj = 10^5
n = [n_0]
time = [0]
while n[j] > 0 & j < maxj:
print(time,n)
time.append(np.random.exponential(2*n[j]**2/5))
n.append(np.random.choice([-1,1]))
j = j+1
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。