source('../../Data/get_visual_fusion_data.R')
DF = get_visual_fusion_data()

nv=DF[DF$Treatment=='NV', 'Time']
vv=DF[DF$Treatment=='VV', 'Time']

## The raw data:
##
par(mfrow=c(1, 2))
qqnorm(nv, main='NV')
qqline(nv)
grid()
qqnorm(vv, main='VV')
qqline(vv)
grid()
par(mfrow=c(1, 1))

## Taking logarithms of the data:
##
log_nv = log(nv)
log_vv = log(vv)
par(mfrow=c(1, 2))
qqnorm(log_nv, main='log(NV)')
qqline(log_nv)
grid()
qqnorm(log_vv, main='log(VV)')
qqline(log_vv)
grid()
par(mfrow=c(1, 1))

res = t.test(log_nv, log_vv)
print(res)