res <- matrix( nrow=4, ncol=5 )

colnames(res) <- c( 'PBF', 'PsBF', 'logPBF', 'LSCV', 'LS')
rownames(res) <- c( 'm1','m2','m3','m4')


#-----------------------------------
x<-read.table('m1.txt',header=TRUE)

LPBF <- log(x[1,2])

PsBF <- 1/prod(x[2:26,2])

# -log(prod(x[2:26,2]))
LSCV <- -sum(log(x[2:26,2]))


LS <- log(prod(x[27:51,2]))

res[1,]<- c( exp(LPBF), PsBF, LPBF, LSCV, LS)
#-----------------------------------

#-----------------------------------
x<-read.table('m2.txt',header=TRUE)

LPBF <- log(x[1,2])

PsBF <- 1/prod(x[2:26,2])

# -log(prod(x[2:26,2]))
LSCV <- -sum(log(x[2:26,2]))


LS <- log(prod(x[27:51,2]))

res[2,]<- c( exp(LPBF), PsBF, LPBF, LSCV, LS)
#-----------------------------------
#-----------------------------------
x<-read.table('m3.txt',header=TRUE)

LPBF <- log(x[1,2])

PsBF <- 1/prod(x[2:26,2])

# -log(prod(x[2:26,2]))
LSCV <- -sum(log(x[2:26,2]))


LS <- log(prod(x[27:51,2]))

res[3,]<- c( exp(LPBF), PsBF, LPBF, LSCV, LS)
#-----------------------------------
#-----------------------------------
x<-read.table('m4.txt',header=TRUE)

LPBF <- log(x[1,2])

PsBF <- 1/prod(x[2:26,2])

# -log(prod(x[2:26,2]))
LSCV <- -sum(log(x[2:26,2]))


LS <- log(prod(x[27:51,2]))

res[4,]<- c( exp(LPBF), PsBF, LPBF, LSCV, LS)
#-----------------------------------


res1<- cbind( res[,3:5], log(res[,1]/res[4,1]), log(res[,2]/res[4,2]), res[,1]/sum(res[,1]), res[,2]/sum(res[,2]) )

colnames(res1) <- c('logPBF', 'LSCV', 'LS', 'logPBF', 'logPsBF', 'wPBF', 'wPsBF')

res2<- cbind( round( res1[,1:5],1) , round(res1[,6:7],4))
colnames(res2) <- c('logPBF', 'LSCV', 'LS', 'logPBF', 'logPsBF', 'wPBF', 'wPsBF')










# c2=1000
res <- matrix( nrow=4, ncol=5 )

colnames(res) <- c( 'PBF', 'PsBF', 'logPBF', 'LSCV', 'LS')
rownames(res) <- c( 'm1','m2','m3','m4')


#-----------------------------------
x<-read.table('c1000m1.txt',header=TRUE)

LPBF <- log(x[1,2])

PsBF <- 1/prod(x[2:26,2])

# -log(prod(x[2:26,2]))
LSCV <- -sum(log(x[2:26,2]))


LS <- log(prod(x[27:51,2]))

res[1,]<- c( exp(LPBF), PsBF, LPBF, LSCV, LS)
#-----------------------------------

#-----------------------------------
x<-read.table('c1000m2.txt',header=TRUE)

LPBF <- log(x[1,2])

PsBF <- 1/prod(x[2:26,2])

# -log(prod(x[2:26,2]))
LSCV <- -sum(log(x[2:26,2]))


LS <- log(prod(x[27:51,2]))

res[2,]<- c( exp(LPBF), PsBF, LPBF, LSCV, LS)
#-----------------------------------
#-----------------------------------
x<-read.table('c1000m3.txt',header=TRUE)

LPBF <- log(x[1,2])

PsBF <- 1/prod(x[2:26,2])

# -log(prod(x[2:26,2]))
LSCV <- -sum(log(x[2:26,2]))


LS <- log(prod(x[27:51,2]))

res[3,]<- c( exp(LPBF), PsBF, LPBF, LSCV, LS)
#-----------------------------------
#-----------------------------------
x<-read.table('c1000m4.txt',header=TRUE)

LPBF <- log(x[1,2])

PsBF <- 1/prod(x[2:26,2])

# -log(prod(x[2:26,2]))
LSCV <- -sum(log(x[2:26,2]))


LS <- log(prod(x[27:51,2]))

res[4,]<- c( exp(LPBF), PsBF, LPBF, LSCV, LS)
#-----------------------------------


res1<- cbind( res[,3:5], log(res[,1]/res[4,1]), log(res[,2]/res[4,2]), res[,1]/sum(res[,1]), res[,2]/sum(res[,2]) )

colnames(res1) <- c('logPBF', 'LSCV', 'LS', 'logPBF', 'logPsBF', 'wPBF', 'wPsBF')

res2<- cbind( round( res1[,1:5],1) , round(res1[,6:7],4))
colnames(res2) <- c('logPBF', 'LSCV', 'LS', 'logPBF', 'logPsBF', 'wPBF', 'wPsBF')