rgb-noise: hvdvm initial analysis

gistfile1.r

hvm.all <- hvdvm$new(has.RGGB.pattern = TRUE)
hvm.all$digest(photo.conds.file = 'ISO100/crops/full-set.csv',
               crop.files.path  = 'ISO100/crops',
               min.raw          = 4,
               max.raw          = c(16379, 15774, 15774 ,16379))

# Fit a quadratic model
hvm.all$fit.model(model.name   = 'weighted',
                  model.family = 'lmrob',
                  degree       = 2,
                  weights      = 1/mean^2)

# Plot the data close to the axes origin
hvm.all$plot(model.name = 'weighted',
             with= ~ channel != 'Green Avg',
             xlim=c(0,50), ylim=c(0, 30))

# Plot the whole data and fitted model
hvm.all$plot(model.name = 'weighted',
             with = ~ channel != 'Green Avg')