R コードを使用してシミュレーションを記述しています。ただし、ここでいくつかの問題に直面しています。シミュレーションを 10 回実行しています。関心のある変数は「I」で、最初は 10 です。
set.seed(12345)
Simulations<-10
Nsteps<-100
S<-50
I<-10
a<-2
singlesim <- function(alpha,beta,s,i){
a<- alpha*s*i
b<- beta*i
d<- runif(1)
ran <- max(a*d,beta*i*d)
z<-c(s,i)
j<- ifelse(c(ran>a&ran<=a+b, ran>a&ran<=a+b),c (s+1,i-1),z)
l<- ifelse(c(ran<=a,ran<=a),c(s-1,i+1),z)
x<- ifelse(c(ran<=a,ran<=a),l,j)
q<- ifelse(j==z&l==z,z,x)
g<- ifelse(c(ran==0,ran==0),z,q)
g}
onecompsim <- function(Nsteps,S,I,alpha,beta){
P<- array(0,c(Nsteps,2))
Z<- array(0,c(Nsteps,2))
Y<- array(0,c(Nsteps,2))
P[1,]<-c(S,I)
for(i in 2:Nsteps){
P[i,]<- singlesim(alpha,beta,P[i-1,1],P[i-1,2])}
P}
fisim<- function(Simulations,Nsteps,S,I,alpha,beta){
Q<-array(0,c(Nsteps,2,Simulations))
for(i in 1:Simulations){
Q[,,i] <- onecompsim(Nsteps,S,I,alpha,beta)}
Q}
res<-fisim(10,100,50,10,0.005,1)
print(res)
このプログラムの最後に、'S' と 'I' の両方の変数で I want の値を取得します。ただし、「I」を保存して「I ^ 2」を2乗し、「I ^ 2」の値を再度保存したいと思います。どうすれば書けるでしょうか？