#!/usr/bin/env perl

use strict;
use warnings;

use Data::Table;
use Data::Table::Excel qw(tables2xlsx);
use Statistics::Lite qw(mean stddev);

my $dt = Data::Table::fromCSV("data.csv");
my $melt = $dt->melt(['Month', 'Zone', 'Replicate']);

my $cast_mean = $melt->cast(
    ['Month', 'Zone'],
    'variable',
    Data::Table::STRING,
    'value',
    \&mean
); 

my $cast_stddev = $melt->cast(
    ['Month', 'Zone'],
    'variable',
    Data::Table::STRING,
    'value',
    \&stddev
); 

tables2xlsx("descriptive_stats.xlsx", [ $cast_mean, $cast_stddev ]);

exit;