時系列規格化用,中央値時系列の作り方

時系列規格化用,中央値時系列の作り方

     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