-
Notifications
You must be signed in to change notification settings - Fork 0
/
AucAf.R
58 lines (44 loc) · 1.34 KB
/
AucAf.R
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
rawdata1 <- read.csv(".csv",header=T,skip=40,nrows=2710)
rawdata2 <- read.csv(".csv",header=T,skip=40,nrows=2710)
rawdata3 <- read.csv(".csv",header=T,skip=40,nrows=2710)
name<-"Af"
library("signal")
fil_N <- 500
fs <- 10
fn <- fs/2
fc <- c(0.01,1)
fc_norm <- fc/fn
fir_filter <- fir1(fil_N, fc_norm, type="pass")
firProcessed1<-matrix(0,2710,22)
firProcessed2<-matrix(0,2710,22)
firProcessed3<-matrix(0,2710,22)
for(i in 1:22){
firProcessed1[,i] <- filtfilt(fir_filter,rawdata1[ ,i+1])
}
for(i in 1:22){
firProcessed2[,i] <- filtfilt(fir_filter,rawdata2[ ,i+1])
}
for(i in 1:22){
firProcessed3[,i] <- filtfilt(fir_filter,rawdata3[ ,i+1])
}
library("stats")
addAverage <- ( firProcessed1 + firProcessed2 + firProcessed3 )/3
MovAverage<-stats::filter(addAverage,rep(1,10),sides=1)/10
Average<-MovAverage[-1:-9,]
AUCfname<-paste("AUCf",name,sep="")
AUCsname<-paste("AUCs",name,sep="")
AUCf<-matrix(0,nrow=1, ncol=22)
AUCs<-matrix(0,nrow=1, ncol=22)
for(i in 1:22){
AUCf[1,i]<- sum(Average[601:1200,i])/10
}
for(i in 1:22){
AUCs[1,i]<- sum(Average[2101:2700,i])/10
}
assign(name,Average)
assign(AUCfname,AUCf)
assign(AUCsname,AUCs)
AUCfnameCsv<-paste("AucPA",name,".csv",sep="")
AUCsnameCsv<-paste("AucPO",name,".csv",sep="")
write.csv(AUCf,file=AUCfnameCsv)
write.csv(AUCs,file=AUCsnameCsv)