User Tools

Site Tools


lab:projects:19time_delay_embedding_dfa:r_code_to_process_results_from_mkfm6

R code related to Kalman filter methods

  • Process results

Click to display ⇲

Click to hide ⇱

## process results
allse<-NULL
allpar<-NULL
SS<-200
for (i in 1:500){
	filename<-paste('out', SS, i, '.txt', sep='')
	if (file.exists(filename)){
		temp<-readLines(filename)
		temp<-temp[substr(temp, 1, 4)==' nr ']
		temp<-substr(temp, 15, 71)
		if (length(temp)>0){
			temp<-read.table(textConnection(temp))
			allse<-rbind(allse, temp$V5)
			allpar<-rbind(allpar, temp$V1)
		}
	}
}

index<-(apply(allse, 1, min) > 0) & (apply(allse, 1, max) < 1)
allpar<-allpar[index, ]
allse<-allse[index, ]


cbind(apply(allpar, 2, mean), apply(allpar, 2, sd), apply(allse, 2, mean))

nrow(allpar)

Page Tools