[BACK]Return to lonproblemanalysis.pm CVS log [TXT] [DIR] Up to [LON-CAPA] / loncom / interface / statistics

Diff for /loncom/interface/statistics/lonproblemanalysis.pm between versions 1.113 and 1.114

version 1.113, 2005/02/22 05:28:21 version 1.114, 2005/02/22 22:51:16
Line 303  sub numerical_plot_percent { Line 303  sub numerical_plot_percent {
     #      #
     my $total = $stats->{'submission_count'};      my $total = $stats->{'submission_count'};
     return '' if ($total == 0);      return '' if ($total == 0);
     my $min_bin_size = $stats->{'min_percent'};      my $max_bins = 50;
     my $low_bin  = $stats->{'lowest_ans'}-$stats->{'max_bin_size'};      my $lowest_percent = $stats->{'min_percent'};
     my $high_bin = $stats->{'highest_ans'}+$stats->{'max_bin_size'};      my $highest_percent = $stats->{'max_percent'};
     if ($high_bin > 0 && $low_bin > -$high_bin) {      my $percent_spread = $highest_percent - $lowest_percent;
         $low_bin = -$high_bin;      foreach (qw/20 30 40 50 100 200/) {
     } elsif ($low_bin < 0 && $high_bin < -$low_bin) {          if ($percent_spread < $_) {
         $high_bin = -$low_bin;              $highest_percent =$_/2;
               last;
           }
     }      }
     if (($high_bin -$low_bin)/$min_bin_size > 1000) {      my $bin_size = 1;
         $min_bin_size = abs($high_bin - $low_bin) / 1000;      foreach (qw/0.01 0.05 0.1 0.5 1 2 5 10 20 25 50 100/) {
           if ($lowest_percent/2 < $_){
               $bin_size = $_;
               last;
           }
     }      }
     my @bins;      my @bins;
     for (my $num = $low_bin;$num <= $high_bin;$num+=($min_bin_size/2)) {      for (my $bin = -$highest_percent;$bin<0;$bin+=$bin_size) {
         push(@bins,$num);          push (@bins,$bin);
       }
       for (my $bin = 0; $bin<$highest_percent;$bin+=$bin_size) {
           push (@bins,$bin);
     }      }
       push(@bins,$highest_percent);
     #      #
     my @correct;      my @correct;
     my @incorrect;      my @incorrect;
Line 325  sub numerical_plot_percent { Line 335  sub numerical_plot_percent {
     while (my ($ans,$submissions) = each(%$responses)) {      while (my ($ans,$submissions) = each(%$responses)) {
         while (my ($submission,$counts) = each(%$submissions)) {          while (my ($submission,$counts) = each(%$submissions)) {
             my ($correct_count,$incorrect_count) = @$counts;              my ($correct_count,$incorrect_count) = @$counts;
             my $scaled_value = ($submission-$ans)/$ans;              my $scaled_value = 100*($submission-$ans)/abs($ans);
             my $bin=0;              if ($scaled_value < $bins[0]) {
             for ($bin=0;$bin<$#bins;$bin++) { # not <= for a reason                  $bins[0]=$scaled_value -1;
               }
               my $bin;
               for ($bin=0;$bin<$#bins;$bin++) {
                 last if ($bins[$bin]>$scaled_value);                  last if ($bins[$bin]>$scaled_value);
             }              }
             $correct[$bin]+=$correct_count;              $correct[$bin-1]+=$correct_count;
             $incorrect[$bin]+=$incorrect_count;              $incorrect[$bin-1]+=$incorrect_count;
             $count[$bin]+=$correct_count+$incorrect_count;              $count[$bin-1]+=$correct_count+$incorrect_count;
         }          }
     }      }
     #      #
     # Skip empty bins      my @plot_correct   = @correct;
     my (@plot_correct,@plot_incorrect,@new_bins,@new_count);      my @plot_incorrect = @incorrect;
     my $min_skip = 2;      my $max;
     for (my $i=0;$i<=$#bins;$i++) {      for (my $i=0;$i<$#bins;$i++) {
         my $sum=0;  
         for (my $j=-$min_skip;$j<=$min_skip && $i+$j<=$#bins;$j++) {  
             $sum += $correct[$i+$j] + $incorrect[$i+$j];  
         }  
         if ($sum) {  
             push(@new_bins,$bins[$i]);  
             push(@plot_correct,$correct[$i]);  
             push(@plot_incorrect,$incorrect[$i]);  
             push(@new_count,$correct[$i]+$incorrect[$i]);  
         }  
     }  
     @correct   = @plot_correct;  
     @incorrect = @plot_incorrect;  
     @count     = @new_count;  
     @bins      = @new_bins;  
     for (my $i=0;$i<=$#bins;$i++) {  
         $plot_correct[$i]   *= 100/$total;          $plot_correct[$i]   *= 100/$total;
         $plot_incorrect[$i] *= 100/$total;          $plot_incorrect[$i] *= 100/$total;
           if (! defined($max) || 
               $max < $plot_correct[$i]+$plot_incorrect[$i] ) {
               $max = $plot_correct[$i]+$plot_incorrect[$i];
           }
       }
       foreach (qw/1 5 10 15 20 25 30 40 50 75 100/) {
           if ($max <$_) { $max = $_; last; }
     }      }
     #      #
     my $title = &mt('Percent Difference');      my $title = &mt('Percent Difference');
     my @labels = (1..scalar(@bins));      my @labels = (1..scalar(@bins)-1);
     my $graph = &Apache::loncommon::DrawBarGraph      my $graph = &Apache::loncommon::DrawBarGraph
         ($title,'Percent difference from correct','Number of answers',          ($title,'Percent Difference from Correct','Percent of Answers',
          100,['#33FF00','#FF3300'],\@labels,\@plot_correct,\@plot_incorrect,           $max,['#33FF00','#FF3300'],\@labels,\@plot_correct,\@plot_incorrect,
          {xskip=>1});           {xskip=>1});
     #      #
     my $table = $graph.$/.      my $table = $graph.$/.
Line 377  sub numerical_plot_differences { Line 381  sub numerical_plot_differences {
     #      #
     my $total = $stats->{'submission_count'};      my $total = $stats->{'submission_count'};
     return '' if ($total == 0);      return '' if ($total == 0);
     my $max_bins = 50;      my $max_bins = 21;
     my $min_bin_size = $stats->{'min_abs'};      my $min_bin_size = $stats->{'min_abs'};
     my $low_bin  = $stats->{'lowest_ans'}-$stats->{'max_bin_size'};      my $low_bin  = $stats->{'lowest_ans'}-$stats->{'max_bin_size'};
     my $high_bin = $stats->{'highest_ans'}+$stats->{'max_bin_size'};      my $high_bin = $stats->{'highest_ans'}+$stats->{'max_bin_size'};
Line 401  sub numerical_plot_differences { Line 405  sub numerical_plot_differences {
         while (my ($submission,$counts) = each(%$submissions)) {          while (my ($submission,$counts) = each(%$submissions)) {
             my ($correct_count,$incorrect_count) = @$counts;              my ($correct_count,$incorrect_count) = @$counts;
             my $scaled_value = $submission-$ans;              my $scaled_value = $submission-$ans;
               if ($scaled_value < $bins[0]) {
                   $bins[0]=$scaled_value-1;
               }
             my $bin=0;              my $bin=0;
             for ($bin=0;$bin<$#bins;$bin++) { # not <= for a reason              for ($bin=0;$bin<$#bins;$bin++) {
                 last if ($bins[$bin]>$scaled_value);                  last if ($bins[$bin]>$scaled_value);
             }              }
             $correct[$bin]+=$correct_count;              $correct[$bin-1]+=$correct_count;
             $incorrect[$bin]+=$incorrect_count;              $incorrect[$bin-1]+=$incorrect_count;
             $count[$bin]+=$correct_count+$incorrect_count;              $count[$bin-1]+=$correct_count+$incorrect_count;
         }          }
     }      }
     my @plot_correct   = @correct;      my @plot_correct   = @correct;
     my @plot_incorrect = @incorrect;      my @plot_incorrect = @incorrect;
       my $max;
     for (my $i=0;$i<=$#bins;$i++) {      for (my $i=0;$i<=$#bins;$i++) {
         $plot_correct[$i]   *= 100/$total;          $plot_correct[$i]   *= 100/$total;
         $plot_incorrect[$i] *= 100/$total;          $plot_incorrect[$i] *= 100/$total;
           if (! defined($max) || 
               $max < $plot_correct[$i]+$plot_incorrect[$i] ) {
               $max = $plot_correct[$i]+$plot_incorrect[$i];
           }
       }
       foreach (qw/1 5 10 15 20 25 30 40 50 75 100/) {
           if ($max <$_) { $max = $_; last; }
     }      }
     #      #
     my $title = &mt('Difference between submission and correct');      my $title = &mt('Difference between submission and correct');
     my @labels = (1..scalar(@bins));      my @labels = (1..scalar(@bins)-1);
     my $graph = &Apache::loncommon::DrawBarGraph      my $graph = &Apache::loncommon::DrawBarGraph
         ($title,'Difference from Correct','Number of answers',          ($title,'Difference from Correct','Percent of Answers',
          100,['#33FF00','#FF3300'],\@labels,\@plot_correct,\@plot_incorrect,           $max,['#33FF00','#FF3300'],\@labels,\@plot_correct,\@plot_incorrect,
          {xskip=>1});           {xskip=>1});
     #      #
     my $table = $graph.$/.      my $table = $graph.$/.
Line 498  sub numerical_classify_responses { Line 513  sub numerical_classify_responses {
                 &capa::caparesponse_get_real_response($myunit,                  &capa::caparesponse_get_real_response($myunit,
                                                       $mysub,                                                        $mysub,
                                                       \$scaled);                                                        \$scaled);
             &Apache::lonnet::logthis('scaled = '.$scaled.' result ='.$result);  #            &Apache::lonnet::logthis('scaled = '.$scaled.' result ='.$result);
             next if (! defined($scaled));              next if (! defined($scaled));
 #            next if ($result ne '6');  #            next if ($result ne '6');
             my $submission = $scaled;              my $submission = $scaled;
Line 546  sub numerical_bin_table { Line 561  sub numerical_bin_table {
         '</tr>'.$/;          '</tr>'.$/;
     for (my $i=0;$i<scalar(@{$bins}-1);$i++) {      for (my $i=0;$i<scalar(@{$bins}-1);$i++) {
         my $lownum = $bins->[$i];          my $lownum = $bins->[$i];
           if ($i == 0) { $lownum = '-&infin;'; }
         my $highnum = $bins->[$i+1];          my $highnum = $bins->[$i+1];
           if ($i == scalar(@{$bins})-2) { $highnum = '&infin;'; }
         $table .=           $table .= 
             '<tr>'.              '<tr>'.
             '<td>'.$labels->[$i].'</td>'.              '<td>'.$labels->[$i].'</td>'.
Removed from v.1.113  
changed lines
  Added in v.1.114

FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>