【R言語】統計解析フリーソフトＲ第５章【GNU R】

【R言語】統計解析フリーソフトＲ第５章【GNU R】at MATH

【R言語】統計解析フリーソフトＲ第５章【GNU R】 - 暇つぶし2ch481:下に貼ってみましたこれで分かりますでしょうか t<-c(1:100) pre<- c(1:100) k1=0.3 k2=0.1 graphplot<-function(x){ for(i in 1:x){ pre[i]<- k1/(k1-k2)*(exp(1)^(-k2*t[i])-exp(1)^(-k1*t[i]))+ ceiling(min(1,max(0,t[i]-7)))*k1/(k1-k2)*(exp(1)^(-k2*(t[i]-7))-exp(1)^(-k1*(t[i]-7)))+ ceiling(min(1,max(0,t[i]-21)))*k1/(k1-k2)*(exp(1)^(-k2*(t[i]-21))-exp(1)^(-k1*(t[i]-21)))+ ceiling(min(1,max(0,t[i]-35)))*k1/(k1-k2)*(exp(1)^(-k2*(t[i]-35))-exp(1)^(-k1*(t[i]-35))) par(new=T) plot(t[i],pre[i],xlim=c(0,100),ylim=c(-5,5)) } } graphplot(100) 続きます

次ページ

続きを表示

1を表示