時系列規格化用,中央値時系列の作り方
apply(mat_food,"median",2) medi1_food<-NULL mean1_food<-apply(mat_food,2,"mean") for(i in 1:length(mat_food[,1])){ q<-(mat_food[i,]/mean1_food)[which(mean1_food>=100)] medi1_food[i]<-median(q) } #r0<-(mean_end+10^-3)/(mean_start+10^-3) mean_start_b<-NULL mean_end_b<-NULL