library(MASS)

ntw.lrn <- read.table("../lrn/ntw_lrn.dat",header=T,sep=",")
ntw.val <- read.table("../val/ntw_val.dat",header=T,sep=",")
ntw.tst <- read.table("../tst/ntw_tst.dat",header=T,sep=",")

ntw.names <- names(ntw.lrn)

y.lrn   <- ntw.lrn$connection.type
f23.lrn <- ntw.lrn$count
f29.lrn <- ntw.lrn$same.srv.rate
f.lrn   <- cbind(f23.lrn,f29.lrn)

y.val   <- ntw.val$connection.type
f23.val <- ntw.val$count
f29.val <- ntw.val$same.srv.rate
f.val   <- cbind(f23.val,f29.val)

y.tst   <- ntw.tst$connection.type
f23.tst <- ntw.tst$count
f29.tst <- ntw.tst$same.srv.rate
f.tst   <- cbind(f23.tst,f29.tst)

rm(ntw.lrn,ntw.val,ntw.tst)

f23.lrn.neptune <- f23.lrn[y.lrn=="neptune"]
f23.lrn.normal <- f23.lrn[y.lrn=="normal"]
f23.lrn.smurf <- f23.lrn[y.lrn=="smurf"]

f29.lrn.neptune <- f29.lrn[y.lrn=="neptune"]
f29.lrn.normal <- f29.lrn[y.lrn=="normal"]
f29.lrn.smurf <- f29.lrn[y.lrn=="smurf"]

f.lrn.neptune <- f.lrn[y.lrn=="neptune",]
f.lrn.normal <- f.lrn[y.lrn=="normal",]
f.lrn.smurf <- f.lrn[y.lrn=="smurf",]

m.lrn.neptune <- cbind(mean(f23.lrn.neptune),mean(f29.lrn.neptune))
m.lrn.normal <- cbind(mean(f23.lrn.normal),mean(f29.lrn.normal))
m.lrn.smurf <- cbind(mean(f23.lrn.smurf),mean(f29.lrn.smurf))

s.lrn.neptune <- cbind(sd(f23.lrn.neptune),sd(f29.lrn.neptune))
s.lrn.normal <- cbind(sd(f23.lrn.normal),sd(f29.lrn.normal))
s.lrn.smurf <- cbind(sd(f23.lrn.smurf),sd(f29.lrn.smurf))

c.lrn.neptune <- cor(f23.lrn.neptune,f29.lrn.neptune)
c.lrn.normal <- cor(f23.lrn.normal,f29.lrn.normal)
c.lrn.smurf <- cor(f23.lrn.smurf,f29.lrn.smurf)

v.lrn.neptune <- var(f.lrn.neptune)
v.lrn.normal <- var(f.lrn.normal)
v.lrn.smurf <- var(f.lrn.smurf)

v.lrn <- (v.lrn.neptune + v.lrn.normal + v.lrn.smurf)/3
s.lrn <- c(sqrt(v.lrn[1,1]),sqrt(v.lrn[2,2]))
c.lrn <- v.lrn[1,2]/(s.lrn[1]*s.lrn[2])

print(" "); print(s.lrn) ; print(c.lrn)
print(" "); print(m.lrn.neptune) ; print(s.lrn.neptune); print(c.lrn.neptune)
print(" "); print(m.lrn.normal) ; print(s.lrn.normal); print(c.lrn.normal)
print(" "); print(m.lrn.smurf) ; print(s.lrn.smurf); print(c.lrn.smurf)