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

version 1.111, 2005/02/16 17:43:20 version 1.113, 2005/02/22 05:28:21
Line 260  sub numerical_response_analysis { Line 260  sub numerical_response_analysis {
                                          $respid,$students);                                           $respid,$students);
         if ($c->aborted()) { return; };          if ($c->aborted()) { return; };
         #          #
         my $responses = &numerical_classify_responses($response_data,$correct,          my ($responses,$stats) = 
                                                       $restriction_function);              &numerical_classify_responses($response_data,$correct,
         if ($responses->{'_count'} == 0) {                                            $restriction_function);
           if ($stats->{'submission_count'} == 0) {
             $analysis_html.=               $analysis_html.= 
                 '<tr><td colspan="2"><font size="+1"><b>'.                  '<tr><td colspan="2"><font size="+1"><b>'.
                 &mt($no_data_message,$plot_num,@extra_data).                  &mt($no_data_message,$plot_num,@extra_data).
Line 274  sub numerical_response_analysis { Line 275  sub numerical_response_analysis {
                 '</b></font></td></tr>'.                  '</b></font></td></tr>'.
                 '<tr><td colspan="2" align="center">'.                  '<tr><td colspan="2" align="center">'.
                 &mt($stats_message,                  &mt($stats_message,
                     $responses->{'_count'},                      $stats->{'submission_count'},
                     $responses->{'_correct'},                      $stats->{'correct_count'},
                     $responses->{'_count'}-$responses->{'_correct'},                      $stats->{'incorrect_count'},
                     $responses->{'_students'},                      $stats->{'students'},
                     @extra_data).                      @extra_data).
                     '</td></tr>'.                      '</td></tr>'.
                     '<tr>'.'<td valign="top" align="center">'.                      '<tr>'.'<td valign="top" align="center">'.
                     &numerical_plot_percent($r,$responses).'</td>'.                      &numerical_plot_percent($r,$responses,$stats).'</td>'.
                     '<td align="center" valign="top">'.                      '<td align="center" valign="top">'.
                     &numerical_plot_differences($r,$responses).'</td>'.                      &numerical_plot_differences($r,$responses,$stats).'</td>'.
                     '</tr>';                      '</tr>';
         }          }
         if ($post_message ne '') {          if ($post_message ne '') {
Line 298  sub numerical_response_analysis { Line 299  sub numerical_response_analysis {
 }  }
   
 sub numerical_plot_percent {  sub numerical_plot_percent {
     my ($r,$responses) = @_;      my ($r,$responses,$stats) = @_;
     #      #
     my $total = $responses->{'_count'};      my $total = $stats->{'submission_count'};
     return '' if ($total == 0);      return '' if ($total == 0);
     my $minbin = 5;      my $min_bin_size = $stats->{'min_percent'};
     while (my ($interval,$submissions) = each(%$responses)) {      my $low_bin  = $stats->{'lowest_ans'}-$stats->{'max_bin_size'};
         next if ($interval =~ /^_/);      my $high_bin = $stats->{'highest_ans'}+$stats->{'max_bin_size'};
         my ($ans,$ans_low,$ans_high) = split(" ",$interval);      if ($high_bin > 0 && $low_bin > -$high_bin) {
         my $low_percent  = abs(100*($ans-$ans_low)/$ans);          $low_bin = -$high_bin;
         my $high_percent = abs(100*($ans_high-$ans)/$ans);      } elsif ($low_bin < 0 && $high_bin < -$low_bin) {
         if ($minbin > $high_percent) { $minbin = $high_percent; }          $high_bin = -$low_bin;
         if ($minbin > $low_percent) { $minbin = $low_percent; }      }
       if (($high_bin -$low_bin)/$min_bin_size > 1000) {
           $min_bin_size = abs($high_bin - $low_bin) / 1000;
     }      }
     #  
     my @bins;      my @bins;
     if ($minbin < 1) {      for (my $num = $low_bin;$num <= $high_bin;$num+=($min_bin_size/2)) {
         @bins = ('0.1','0.5','1.0','1.5','2.0','2.5','3.0','4.0','5.0',10,20,50,100);                  push(@bins,$num);
     } elsif ($minbin < 2) {  
         @bins = ('0.5','1.0','1.5','2.0','2.5','3.0','4.0','5.0',10,20,50,100);  
     } elsif ($minbin < 5) {  
         @bins = (1,2,3,4,5,10,25,50,75,100,200);  
     } elsif ($minbin < 10) {  
         @bins = (2,4,6,8,10,12,15,20,25,30,50,75,100,200);  
     } else {  
         @bins = (5,10,15,20,25,30,50,75,100,200);  
     }      }
     my @labels = (1..scalar(@bins));  
     #      #
     my @correct;      my @correct;
     my @incorrect;      my @incorrect;
     my @count;      my @count;
     while (my ($interval,$submissions) = each(%$responses)) {      while (my ($ans,$submissions) = each(%$responses)) {
         next if ($interval =~ /^_/);  
         my ($ans,$ans_low,$ans_high) = split(" ",$interval);  
         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 = abs(($submission-$ans)/$ans);              my $scaled_value = ($submission-$ans)/$ans;
             my $bin=0;              my $bin=0;
             for ($bin=0;$bin<$#bins;$bin++) { # not <= for a reason              for ($bin=0;$bin<$#bins;$bin++) { # not <= for a reason
                 last if ($bins[$bin]>$scaled_value);                  last if ($bins[$bin]>$scaled_value);
Line 345  sub numerical_plot_percent { Line 336  sub numerical_plot_percent {
         }          }
     }      }
     #      #
     my @plot_correct;      # Skip empty bins
     my @plot_incorrect;      my (@plot_correct,@plot_incorrect,@new_bins,@new_count);
       my $min_skip = 2;
       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++) {      for (my $i=0;$i<=$#bins;$i++) {
         $plot_correct[$i] = $correct[$i]*100/$total;          $plot_correct[$i]   *= 100/$total;
         $plot_incorrect[$i] = $incorrect[$i]*100/$total;          $plot_incorrect[$i] *= 100/$total;
     }      }
     my $title = &mt('Distribution by Percent');      #
       my $title = &mt('Percent Difference');
       my @labels = (1..scalar(@bins));
     my $graph = &Apache::loncommon::DrawBarGraph      my $graph = &Apache::loncommon::DrawBarGraph
         ($title,'Percent difference from correct','Number of answers',          ($title,'Percent difference from correct','Number of answers',
          100,['#33FF00','#FF3300'],\@labels,\@plot_correct,\@plot_incorrect,           100,['#33FF00','#FF3300'],\@labels,\@plot_correct,\@plot_incorrect,
Line 363  sub numerical_plot_percent { Line 373  sub numerical_plot_percent {
 }  }
   
 sub numerical_plot_differences {  sub numerical_plot_differences {
     my ($r,$responses) = @_;      my ($r,$responses,$stats) = @_;
     #      #
     my $total = $responses->{'_count'};      my $total = $stats->{'submission_count'};
     return '' if ($total == 0);      return '' if ($total == 0);
     my $minbin = undef;      my $max_bins = 50;
     my $maxbin = undef;      my $min_bin_size = $stats->{'min_abs'};
     while (my ($interval,$submissions) = each(%$responses)) {      my $low_bin  = $stats->{'lowest_ans'}-$stats->{'max_bin_size'};
         next if ($interval =~ /^_/);      my $high_bin = $stats->{'highest_ans'}+$stats->{'max_bin_size'};
         my ($ans,$ans_low,$ans_high) = split(" ",$interval);      if ($high_bin > 0 && $low_bin > -$high_bin) {
         my $low_diff  = abs($ans-$ans_low);          $low_bin = -$high_bin;
         my $high_diff = abs($ans_high-$ans);      } elsif ($low_bin < 0 && $high_bin < -$low_bin) {
         if (! defined($maxbin)) { $maxbin = $low_diff;}          $high_bin = -$low_bin;
         if (! defined($minbin)) { $minbin = $low_diff;}      }
         #      if (($high_bin -$low_bin)/$min_bin_size * 2 > $max_bins) {
         if ($minbin > $high_diff) { $minbin = $high_diff; }          $min_bin_size = abs($high_bin - $low_bin) / $max_bins * 2;
         if ($minbin > $low_diff ) { $minbin = $low_diff; }      }
         #  
         if ($maxbin < $high_diff) { $maxbin = $high_diff; }  
         if ($maxbin < $low_diff ) { $maxbin = $low_diff; }  
     }      
     #  
     my @bins;      my @bins;
     my @labels;      for (my $num = $low_bin;$num <= $high_bin;$num+=($min_bin_size/2)) {
     # Hmmmm, should switch to absolute difference          push(@bins,$num);
     for (my $i=1;$i<=20;$i++) {  
         push(@bins,$i*$minbin/2);  
         push(@labels,$i);  
     }      }
     #      #
     my @correct;      my @correct;
     my @incorrect;      my @incorrect;
     my @count;      my @count;
     while (my ($interval,$submissions) = each(%$responses)) {      while (my ($ans,$submissions) = each(%$responses)) {
         next if ($interval =~ /^_/);  
         my ($ans,$ans_low,$ans_high) = split(" ",$interval);  
         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 $value = abs($submission-$ans);              my $scaled_value = $submission-$ans;
             my $bin=0;              my $bin=0;
             for ($bin=0;$bin<$#bins;$bin++) { # not <= for a reason              for ($bin=0;$bin<$#bins;$bin++) { # not <= for a reason
                 last if ($bins[$bin]>$value);                  last if ($bins[$bin]>$scaled_value);
             }              }
             $correct[$bin]+=$correct_count;              $correct[$bin]+=$correct_count;
             $incorrect[$bin]+=$incorrect_count;              $incorrect[$bin]+=$incorrect_count;
             $count[$bin]+=$correct_count+$incorrect_count;              $count[$bin]+=$correct_count+$incorrect_count;
         }          }
     }      }
     #      my @plot_correct   = @correct;
     my @plot_correct;      my @plot_incorrect = @incorrect;
     my @plot_incorrect;  
     for (my $i=0;$i<=$#bins;$i++) {      for (my $i=0;$i<=$#bins;$i++) {
         $plot_correct[$i]   =   $correct[$i]*100/$total;          $plot_correct[$i]   *= 100/$total;
         $plot_incorrect[$i] = $incorrect[$i]*100/$total;          $plot_incorrect[$i] *= 100/$total;
     }      }
     my $title = &mt('Distribution by Magnitude');      #
       my $title = &mt('Difference between submission and correct');
       my @labels = (1..scalar(@bins));
     my $graph = &Apache::loncommon::DrawBarGraph      my $graph = &Apache::loncommon::DrawBarGraph
         ($title,'magnitude difference from correct','Number of answers',          ($title,'Difference from Correct','Number of answers',
          100,['#33FF00','#FF3300'],\@labels,\@plot_correct,\@plot_incorrect,           100,['#33FF00','#FF3300'],\@labels,\@plot_correct,\@plot_incorrect,
          {xskip=>1});           {xskip=>1});
     #      #
Line 429  sub numerical_plot_differences { Line 430  sub numerical_plot_differences {
 }  }
   
 sub numerical_classify_responses {  sub numerical_classify_responses {
     &Apache::lonnet::logthis('--------------');  
     my ($full_row_data,$correct,$function) = @_;      my ($full_row_data,$correct,$function) = @_;
     my %submission_data;      my %submission_data;
     my %students;      my %students;
       my %stats;
     my $max=0;      my $max=0;
     foreach my $row (@$full_row_data) {      foreach my $row (@$full_row_data) {
 #        &Apache::lonnet::logthis(' row = '.join(',',@$row));  
         my %subm = &hashify_attempt($row);          my %subm = &hashify_attempt($row);
         if (ref($correct) eq 'HASH') {          if (ref($correct) eq 'HASH') {
             $subm{'correct'} = $correct->{$subm{'student'}}->{'answer'};              my $s_correct = $correct->{$subm{'student'}};
             $subm{'unit'} = $correct->{$subm{'student'}}->{'unit'};              $subm{'correct'} = $s_correct->{'answer'};
               foreach my $item ('unit','ans_low','ans_high') {
                   $subm{$item} = $s_correct->{$item};
               }
         } else { # This probably never happens....          } else { # This probably never happens....
             $subm{'correct'} = $correct->{'answer'};              $subm{'correct'} = $correct->{'answer'};
             $subm{'unit'} = $correct->{'unit'};              $subm{'unit'} = $correct->{'unit'};
         }          }
           # 
           my $abs_low =abs($subm{'correct'}-$subm{'ans_low'});
           my $abs_high=abs($subm{'correct'}-$subm{'ans_high'});
           if (! defined($stats{'min_abs'}) ||
               $stats{'min_abs'} > $abs_low) {
               $stats{'min_abs'} = $abs_low;
           }
           if ($stats{'min_abs'} > $abs_high) {
               $stats{'min_abs'} = $abs_high;
           }
           if (! defined($stats{'max_abs'}) ||
               $stats{'max_abs'} < $abs_low) {
               $stats{'max_abs'} = $abs_low;
           }
           if ($stats{'max_abs'} < $abs_high) {
               $stats{'max_abs'} = $abs_high;
           }
           my $low_percent  = 100 * abs($abs_low  / $subm{'correct'});
           my $high_percent = 100 * abs($abs_high / $subm{'correct'});
           if (! defined($stats{'min_percent'}) ||
               $stats{'min_percent'} > $low_percent) {
               $stats{'min_percent'} = $low_percent;
           }
           if ($stats{'min_percent'} > $high_percent) {
               $stats{'min_percent'} = $high_percent;
           }
           if (! defined($stats{'max_percent'}) ||
               $stats{'max_percent'} < $low_percent) {
               $stats{'max_percent'} = $low_percent;
           }
           if ($stats{'max_percent'} < $high_percent) {
               $stats{'max_percent'} = $high_percent;
           }
           if (! defined($stats{'lowest_ans'}) ||
               $stats{'lowest_ans'} > $subm{'correct'}) {
               $stats{'lowest_ans'} = $subm{'correct'};
           }
           if (! defined($stats{'highest_ans'}) ||
               $stats{'highest_ans'} < $subm{'correct'}) {
               $stats{'highest_ans'} = $subm{'correct'};
           }
           # 
         $subm{'submission'} =~ s/=\d+\s*$//;          $subm{'submission'} =~ s/=\d+\s*$//;
         if (&$function(\%subm)) {          if (&$function(\%subm)) {
             my $scaled = '1';              my $scaled = '1';
             my ($sname,$sdom) = split(':',$subm{'student'});              my ($sname,$sdom) = split(':',$subm{'student'});
             # Note that $subm{'unit'} is modified by the following call  
             # We do not use it again but you should be aware just in case.  
             my ($myunit,$mysub) = ($subm{'unit'},$subm{'submission'});              my ($myunit,$mysub) = ($subm{'unit'},$subm{'submission'});
             my $result =               my $result = 
                 &capa::caparesponse_get_real_response($myunit,                  &capa::caparesponse_get_real_response($myunit,
                                                       $mysub,                                                        $mysub,
                                                       \$scaled);                                                        \$scaled);
               &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;
             $students{$subm{'student'}}++;              $students{$subm{'student'}}++;
               $stats{'submission_count'}++;
             if (&numerical_submission_is_correct($subm{'award'})) {               if (&numerical_submission_is_correct($subm{'award'})) { 
                 &Apache::lonnet::logthis('correct:'.$submission.':'.$subm{'correct'});                  $stats{'correct_count'}++;
                 $submission_data{'_correct'}++;  
                 $submission_data{'_count'}++;  
                 $submission_data{$subm{'correct'}}->{$submission}->[0]++;                  $submission_data{$subm{'correct'}}->{$submission}->[0]++;
             } elsif (&numerical_submission_is_incorrect($subm{'award'})) {               } elsif (&numerical_submission_is_incorrect($subm{'award'})) { 
                 &Apache::lonnet::logthis('incorrect:'.$submission.':'.$subm{'correct'});                  $stats{'incorrect_count'}++;
                 $submission_data{'_count'}++;  
                 $submission_data{$subm{'correct'}}->{$submission}->[1]++;                  $submission_data{$subm{'correct'}}->{$submission}->[1]++;
             }              }
             my $value =   
                 $submission_data{$subm{'correct'}}->{$submission}->[0]+   
                 $submission_data{$subm{'correct'}}->{$submission}->[1];  
             if ($max < $value) { $max = $value; }  
         }          }
     }      }
     $submission_data{'_max'} = $max;      $stats{'students'}=scalar(keys(%students));
     $submission_data{'_students'}=scalar(keys(%students));      return (\%submission_data,\%stats);
     return \%submission_data;  
 }  }
   
 sub numerical_submission_is_correct {  sub numerical_submission_is_correct {
     my ($award) = @_;      my ($award) = @_;
     &Apache::lonnet::logthis('award = "'.$award.'"');  
     if ($award =~ /^(APPROX_ANS|EXACT_ANS)$/) {      if ($award =~ /^(APPROX_ANS|EXACT_ANS)$/) {
         return 1;          return 1;
     } else {      } else {
Line 508  sub numerical_bin_table { Line 544  sub numerical_bin_table {
         '<th>'.&mt('Correct').'</th>'.          '<th>'.&mt('Correct').'</th>'.
         '<th>'.&mt('Count').'</th>'.          '<th>'.&mt('Count').'</th>'.
         '</tr>'.$/;          '</tr>'.$/;
     for (my $i=0;$i<scalar(@{$bins});$i++) {      for (my $i=0;$i<scalar(@{$bins}-1);$i++) {
         my $lownum;          my $lownum = $bins->[$i];
         if ($i == 0) {          my $highnum = $bins->[$i+1];
             $lownum = 0;  
         } else {  
             $lownum = $bins->[$i-1];  
         }  
         my $highnum = $bins->[$i];  
         $table .=           $table .= 
             '<tr>'.              '<tr>'.
             '<td>'.$labels->[$i].'</td>'.              '<td>'.$labels->[$i].'</td>'.
Line 535  sub numerical_determine_answers { Line 566  sub numerical_determine_answers {
     my ($r,$resource,$partid,$respid,$students)=@_;      my ($r,$resource,$partid,$respid,$students)=@_;
     my $c = $r->connection();      my $c = $r->connection();
     #      #
     # FIX ME: May need progress dialog updates      my %prog_state=&Apache::lonhtmlcommon::Create_PrgWin
           ($r,'Student Answer Compilation Status',
            'Student Answer Compilation Progress', scalar(@$students),
            'inline',undef,'Statistics','stats_status');
     #      #
     # Read in the cache (if it exists) before we start timing things.      # Read in the cache (if it exists) before we start timing things.
     &Apache::lonstathelpers::ensure_proper_cache($resource->{'symb'});      &Apache::lonstathelpers::ensure_proper_cache($resource->{'symb'});
Line 553  sub numerical_determine_answers { Line 587  sub numerical_determine_answers {
                                                                 $sdom);                                                                  $sdom);
         # make the key          # make the key
         my $key = $partid.'.'.$respid;          my $key = $partid.'.'.$respid;
         $correct->{$sname.':'.$sdom}->{'answer'} =           foreach my $item ('answer','unit','ans_high','ans_low') {
             $analysis->{$key.'.answer'}->[0];              $correct->{$sname.':'.$sdom}->{$item} = 
         $correct->{$sname.':'.$sdom}->{'unit'} =                   $analysis->{$key.'.'.$item}->[0];
             $analysis->{$key.'.unit'}->[0];          }
         $answers{$analysis->{$key.'.answer'}->[0]}++;          $answers{$analysis->{$key.'.answer'}->[0]}++;
           &Apache::lonhtmlcommon::Increment_PrgWin($r,\%prog_state,
                                                    &mt('last student'));
     }      }
     &Apache::lonstathelpers::write_analysis_cache();      &Apache::lonstathelpers::write_analysis_cache();
       &Apache::lonhtmlcommon::Close_PrgWin($r,\%prog_state);
     return ($correct,\%answers);      return ($correct,\%answers);
 }  }
   
Line 888  sub ensure_start_end_times { Line 925  sub ensure_start_end_times {
   
 sub RR_concept_plot {  sub RR_concept_plot {
     my ($concepts,$foil_data,$title) = @_;      my ($concepts,$foil_data,$title) = @_;
     &Apache::lonnet::logthis('got to RR_concept_plot');  
     #      #
     my %correct_by_concept;      my %correct_by_concept;
     my %incorrect_by_concept;      my %incorrect_by_concept;
Line 936  sub RR_concept_plot { Line 972  sub RR_concept_plot {
   
 sub RR_create_percent_selected_plot {  sub RR_create_percent_selected_plot {
     my ($concepts,$foils,$foil_data,$title) = @_;      my ($concepts,$foils,$foil_data,$title) = @_;
     &Apache::lonnet::logthis('got to RR_create_percent_selected_plot');  
     #      #
     if ($foil_data->{'_count'} == 0) { return ''; };      if ($foil_data->{'_count'} == 0) { return ''; };
     my %correct_selections;      my %correct_selections;
Line 1008  sub RR_create_percent_selected_plot { Line 1043  sub RR_create_percent_selected_plot {
 sub RR_create_stacked_selection_plot {  sub RR_create_stacked_selection_plot {
     my ($foils,$foil_data,$title,$true_foils)=@_;      my ($foils,$foil_data,$title,$true_foils)=@_;
     #      #
     &Apache::lonnet::logthis('got to RR_create_stacked_selection_plot');  
     my @dataset; # array of array refs - multicolor rows $datasets[row]->[col]      my @dataset; # array of array refs - multicolor rows $datasets[row]->[col]
     my @labels;      my @labels;
     my $count;      my $count;

Removed from v.1.111  
changed lines
  Added in v.1.113


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