my $infile = $ARGV[0]; open (IFILE, "<", $infile) || die "Cannot open : $!\n"; ## read the input file my @data; my $COLUMNTOPARSE; my $local; while () { chomp; my @local = split ('\t',$_); $COLUMNTOPARSE = $local[$ARGV[1]]; push @data; $COLUMNTOPARSE; } close IFILE; my $count = @data; #### my $variance; my $var; foreach my $data[$index] (@data) { $var = 1/($count-1); $variance = $var*(($data[$index]-$Average)**2); } print "$variance\n";