figure 2.4について library(MASS) blue_ave<-c(1,0) red_ave<-c(0,1) blue<-mvrnorm(100,blue_ave,diag(2)) red<-mvrnorm(100,red_ave,diag(2)) cl<-factor(c(rep("blue",100),rep("red",100))) train<-rbind(blue,red) test<-rbind(mvrnorm(5000,blue_ave,…

検索してもあまりいいページ出てこなかったですが式は以下の通り library(deSolve)# 微分方程式の数値解法パッケージ #n:個体数 k:飽和数 s:食傷捕食数 r:再生産率 parameters <- c(k = 100, s = 1, r=0.1)# パラメータのセット initial <- c(n=5)# 初期条件…