Diff for /loncom/interface/statistics/lonproblemanalysis.pm between versions 1.103 and 1.112

version 1.103, 2004/11/10 16:56:39 version 1.112, 2005/02/22 04:32:46
Line 37  use Apache::lonstathelpers(); Line 37  use Apache::lonstathelpers();
 use Apache::lonstudentsubmissions();  use Apache::lonstudentsubmissions();
 use HTML::Entities();  use HTML::Entities();
 use Time::Local();  use Time::Local();
 use Spreadsheet::WriteExcel();  use capa;
   
 my $plotcolors = ['#33ff00',   my $plotcolors = ['#33ff00', 
                   '#0033cc', '#990000', '#aaaa66', '#663399', '#ff9933',                    '#0033cc', '#990000', '#aaaa66', '#663399', '#ff9933',
Line 52  my @SubmitButtons = ({ name => 'PrevProb Line 52  my @SubmitButtons = ({ name => 'PrevProb
                        text => 'Next Problem' },                         text => 'Next Problem' },
                      { name => 'break'},                       { name => 'break'},
                      { name => 'SelectAnother',                       { name => 'SelectAnother',
                        text => 'Choose a different Problem' },                         text => 'Choose a different Problem' });
                      { name => 'ExcelOutput',  
                        text => 'Produce Excel Output' });  
   
 sub BuildProblemAnalysisPage {  sub BuildProblemAnalysisPage {
     my ($r,$c)=@_;      my ($r,$c)=@_;
Line 83  sub BuildProblemAnalysisPage { Line 81  sub BuildProblemAnalysisPage {
         &Apache::lonstathelpers::manage_caches($r,'Statistics','stats_status');          &Apache::lonstathelpers::manage_caches($r,'Statistics','stats_status');
     $r->rflush();      $r->rflush();
     #      #
     # Support for numerical and radio response isn't complete enough to      my $problem_types = '(option|radiobutton|numerical)';
     # include in 1.2 release.  
     # my $problem_types = '(option|radiobutton|numerical)';  
     my $problem_types = '.';#(option)';  
     if (exists($ENV{'form.problemchoice'}) &&       if (exists($ENV{'form.problemchoice'}) && 
         ! exists($ENV{'form.SelectAnother'})) {          ! exists($ENV{'form.SelectAnother'})) {
         foreach my $button (@SubmitButtons) {          foreach my $button (@SubmitButtons) {
Line 193  sub numerical_response_analysis { Line 188  sub numerical_response_analysis {
          $Apache::lonstatistics::enrollment_status,           $Apache::lonstatistics::enrollment_status,
          $resource->{'symb'},$respid);           $resource->{'symb'},$respid);
     #      #
       $problem_analysis->{'answercomputed'} = 1;
       if ($problem_analysis->{'answercomputed'}) {
           my $answers = 
               &Apache::lonstathelpers::GetStudentAnswers($r,$problem,$students,
                                                          'Statistics',
                                                          'stats_status');
           $r->print(&numerical_one_dimensional_plot($r,600,150,$answers));
       }
       #
       if (ref($response_data) ne 'ARRAY') {
           $r->print('<h2>'.
                     &mt('There is no submission data for this resource').
                     '</h2>');
           return;
       }
       my $analysis_html = '<table>';
     for (my $plot_num = 1;$plot_num<=$ENV{'form.NumPlots'};$plot_num++) {      for (my $plot_num = 1;$plot_num<=$ENV{'form.NumPlots'};$plot_num++) {
         my $restriction_function;          my $restriction_function;
         $restriction_function = sub {($_[0]->{'tries'} == $plot_num?1:0)};          my $header_message;
           my $stats_message;
           my $post_message; # passed through &mt sooner rather than later
           my $no_data_message;
           my @extra_data;
           if ($ENV{'form.AnalyzeOver'} eq 'tries') {
               $restriction_function = sub {($_[0]->{'tries'} == $plot_num?1:0)};
               $header_message = 'Attempt [_1]';
               $stats_message = 
                   '[_1] submissions, [_2] correct, [_3] incorrect';
               $post_message = '';
               $no_data_message = 'No data exists for attempt [_1]';
           } else {
               my $starttime = &Apache::lonhtmlcommon::get_date_from_form
                   ('startdate_'.$plot_num);
               my $endtime = &Apache::lonhtmlcommon::get_date_from_form
                   ('enddate_'.$plot_num);
               ($starttime,$endtime) = &ensure_start_end_times
                   ($starttime,$endtime,
                    &get_time_from_row($response_data->[0]),
                    &get_time_from_row($response_data->[-1]),
                    $plot_num);
               $header_message = 'Data from [_2] to [_3]';
               $extra_data[0] = &Apache::lonlocal::locallocaltime($starttime);
               $extra_data[1] = &Apache::lonlocal::locallocaltime($endtime);
               #
               $stats_message = 
                   '[_1] submissions from [_4] students, [_2] correct, [_3] incorrect';
               #
               $post_message = 
                   &mt('Start time: [_1]',
                       &Apache::lonhtmlcommon::date_setter
                       ('Statistics','startdate_'.$plot_num,$starttime)).
                   '<br />'.
                   &mt('End time: [_1]',
                       &Apache::lonhtmlcommon::date_setter
                       ('Statistics','enddate_'.$plot_num,$endtime));
               $restriction_function = 
                   sub { 
                       my $t = $_[0]->{'timestamp'};
                       if ($t >= $starttime && $t < $endtime) {
                           return 1;
                       } else { 
                           return 0;
                       }
                   };
               $no_data_message = 'No data for [_2] to [_3]';
           }
         #          #
         my ($correct,$intervals,$answers) =           my ($correct,$answers) = 
             &numerical_response_determine_intervals($r,$resource,$partid,              &numerical_determine_answers($r,$resource,$partid,
                                                     $respid,$students);                                           $respid,$students);
         if ($c->aborted()) { return; };          if ($c->aborted()) { return; };
         #          #
         my $responses = &classify_response_data($response_data,          my ($responses,$stats) = 
                                                 $correct,              &numerical_classify_responses($response_data,$correct,
                                                 $restriction_function);                                            $restriction_function);
         my $student_count = $responses->{'_students'};          if ($stats->{'submission_count'} == 0) {
         my $correct_count = $responses->{'_correct'};              $analysis_html.= 
         my $total_count   = $responses->{'_count'};                  '<tr><td colspan="2"><font size="+1"><b>'.
         my $max           = $responses->{'_max'};                  &mt($no_data_message,$plot_num,@extra_data).
         #                  '</b></font></td></tr>';
         my $analysis_html = '<h2>'.&mt('Attempt [_1]',$plot_num).'</h2>';          } else {
         $analysis_html.= &numerical_one_dimensional_plot($r,600,150,              $analysis_html.= 
                                                          scalar(@$students),                  '<tr><td colspan="2" align="center"><font size="+1"><b>'.
                                                          $answers);                  &mt($header_message,$plot_num,@extra_data).
         $analysis_html.= '<table><tr>'.                  '</b></font></td></tr>'.
             '<td valign="top" align="center">'.                  '<tr><td colspan="2" align="center">'.
             &numerical_plot_percent($r,$responses).'</td>'.                  &mt($stats_message,
             '<td align="center" valign="top">'.                      $stats->{'submission_count'},
             &numerical_plot_differences($r,$responses).'</td>'.                      $stats->{'correct_count'},
             '</tr>'.$/.'</table>';                      $stats->{'incorrect_count'},
         $r->print($analysis_html);                      $stats->{'students'},
                       @extra_data).
                       '</td></tr>'.
                       '<tr>'.'<td valign="top" align="center">'.
                       &numerical_plot_percent($r,$responses,$stats).'</td>'.
                       '<td align="center" valign="top">'.
                       &numerical_plot_differences($r,$responses,$stats).'</td>'.
                       '</tr>';
           }
           if ($post_message ne '') {
               $analysis_html .=
                   '<tr><td colspan="2">'.$post_message.'</td></tr>';
           }
     }      }
       $analysis_html.='</table>';
       $r->print($analysis_html);
     #      #
     return;      return;
 }  }
   
 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 = 0.5;      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);  
         my $low_percent  = abs(($ans-$ans_low)/$ans);  
         my $high_percent = abs(($ans_high-$ans)/$ans);  
         if ($minbin > $high_percent) { $minbin = $high_percent; }  
         if ($minbin > $low_percent) { $minbin = $low_percent; }  
     }      
     #  
     my @bins;      my @bins;
     if ($minbin < 1) {      for (my $num = $low_bin;$num <= $high_bin;$num+=($min_bin_size/2)) {
         @bins = (0.1, ".5",1,1.5,2,2.5,3,4,5,10,20,50,100);                  push(@bins,$num);
     } elsif ($minbin < 2) {  
         @bins = (0.5,1,1.5,2,2.5,3,4,5,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;
     while (my ($interval,$submissions) = each(%$responses)) {      my @count;
         next if ($interval =~ /^_/);      while (my ($ans,$submissions) = each(%$responses)) {
         my ($ans,$ans_low,$ans_high) = split(" ",$interval);          while (my ($submission,$counts) = each(%$submissions)) {
         while (my ($submission,$count) = each(%$submissions)) {              my ($correct_count,$incorrect_count) = @$counts;
             my $scaled_value = abs(($submission-$ans)/$ans);              my $scaled_value = abs(($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);
             }              }
             if (($submission <= $ans_high) && ($submission >= $ans_low)) {              $correct[$bin]+=$correct_count;
                 $correct[$bin]+=$count;              $incorrect[$bin]+=$incorrect_count;
             } else {              $count[$bin]+=$correct_count+$incorrect_count;
                 $incorrect[$bin]+=$count;  
             }  
         }          }
     }      }
     #      #
       # Skip empty bins
       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++) {
         if (! defined(  $correct[$i])) {   $correct[$i]=0;}          $plot_correct[$i]   *= 100/$total;
         if (! defined($incorrect[$i])) { $incorrect[$i]=0;}          $plot_incorrect[$i] *= 100/$total;
         $correct[$i]*=100/$total;  
         $incorrect[$i]*=100/$total;  
     }      }
       #
     my $title = &mt('Distribution by Percent');      my $title = &mt('Distribution by Percent');
       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,\@correct,\@incorrect,           100,['#33FF00','#FF3300'],\@labels,\@plot_correct,\@plot_incorrect,
          {xskip=>1});           {xskip=>1});
     #      #
     my $table = $graph.'<table><tr><th>'.&mt('Bar').'</th>'.      my $table = $graph.$/.
         '<th colspan="3">'.&mt('Range').'</th>';          &numerical_bin_table(\@bins,\@labels,\@incorrect,\@correct,\@count).$/;
     for (my $i=0;$i<=$#bins;$i++) {  
         my $lownum;  
         if ($i == 0) {  
             $lownum = 0;  
         } else {  
             $lownum = $bins[$i-1];  
         }  
         my $highnum = $bins[$i];  
         $table .=   
             '<tr>'.  
             '<td>'.$labels[$i].'</td>'.  
             '<td align="right">'.$lownum.'</td>'.  
             '<td>&nbsp;-&nbsp;</td>'.  
             '<td align="right">'.$highnum.'</td>'.'</tr>'.$/;  
     }  
     $table.= '</table>'.$/;  
     return $table;      return $table;
 }  }
   
 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 $maxbin = undef;  
     while (my ($interval,$submissions) = each(%$responses)) {  
         next if ($interval =~ /^_/);  
         my ($ans,$ans_low,$ans_high) = split(" ",$interval);  
         my $low_diff  = abs($ans-$ans_low);  
         my $high_diff = abs($ans_high-$ans);  
         if (! defined($maxbin)) { $maxbin = $low_diff;}  
         if (! defined($minbin)) { $minbin = $low_diff;}  
         #  
         if ($minbin > $high_diff) { $minbin = $high_diff; }  
         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;      my @labels;
     # Hmmmm, should switch to absolute difference      # Hmmmm, should switch to absolute difference
     for (my $i=1;$i<=20;$i++) {      for (my $i=1;$i<=20;$i++) {
         push(@bins,$i*$minbin/2);          push(@bins,$i/2);
         push(@labels,$i);          push(@labels,$i);
     }      }
     #      #
     my @correct;      my @correct;
     my @incorrect;      my @incorrect;
     while (my ($interval,$submissions) = each(%$responses)) {      my @count;
         next if ($interval =~ /^_/);      while (my ($ans,$submissions) = each(%$responses)) {
         my ($ans,$ans_low,$ans_high) = split(" ",$interval);          next if ($ans =~ /^_/);
         while (my ($submission,$count) = each(%$submissions)) {          while (my ($submission,$counts) = each(%$submissions)) {
               my ($correct_count,$incorrect_count) = @$counts;
             my $value = abs($submission-$ans);              my $value = abs($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]>$value);
             }              }
             if (($submission <= $ans_high) && ($submission >= $ans_low)) {              $correct[$bin]+=$correct_count;
                 $correct[$bin]+=$count;              $incorrect[$bin]+=$incorrect_count;
             } else {              $count[$bin]+=$correct_count+$incorrect_count;
                 $incorrect[$bin]+=$count;  
             }  
         }          }
     }      }
     #      #
       my @plot_correct;
       my @plot_incorrect;
     for (my $i=0;$i<=$#bins;$i++) {      for (my $i=0;$i<=$#bins;$i++) {
         if (! defined(  $correct[$i])) {   $correct[$i]=0;}          $plot_correct[$i]   =   $correct[$i]*100/$total;
         if (! defined($incorrect[$i])) { $incorrect[$i]=0;}          $plot_incorrect[$i] = $incorrect[$i]*100/$total;
         $correct[$i]*=100/$total;  
         $incorrect[$i]*=100/$total;  
     }      }
     my $title = &mt('Distribution by Magnitude');      my $title = &mt('Distribution by Magnitude');
     my $graph = &Apache::loncommon::DrawBarGraph      my $graph = &Apache::loncommon::DrawBarGraph
         ($title,'magnitude difference from correct','Number of answers',          ($title,'magnitude difference from correct','Number of answers',
          100,['#33FF00','#FF3300'],\@labels,\@correct,\@incorrect,{xskip=>1});           100,['#33FF00','#FF3300'],\@labels,\@plot_correct,\@plot_incorrect,
            {xskip=>1});
     #      #
     my $table = $graph.'<table><tr><th>'.&mt('Bar').'</th>'.      my $table = $graph.$/.
         '<th colspan="3">'.&mt('Range').'</th>';          &numerical_bin_table(\@bins,\@labels,\@incorrect,\@correct,\@count).$/;
     for (my $i=0;$i<=$#bins;$i++) {      return $table;
   }
   
   sub numerical_classify_responses {
       my ($full_row_data,$correct,$function) = @_;
       my %submission_data;
       my %students;
       my %stats;
       my $max=0;
       foreach my $row (@$full_row_data) {
           my %subm = &hashify_attempt($row);
           if (ref($correct) eq 'HASH') {
               my $s_correct = $correct->{$subm{'student'}};
               $subm{'correct'} = $s_correct->{'answer'};
               foreach my $item ('unit','ans_low','ans_high') {
                   $subm{$item} = $s_correct->{$item};
               }
           } else { # This probably never happens....
               $subm{'correct'} = $correct->{'answer'};
               $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*$//;
           if (&$function(\%subm)) {
               my $scaled = '1';
               my ($sname,$sdom) = split(':',$subm{'student'});
               my ($myunit,$mysub) = ($subm{'unit'},$subm{'submission'});
               my $result = 
                   &capa::caparesponse_get_real_response($myunit,
                                                         $mysub,
                                                         \$scaled);
               &Apache::lonnet::logthis('scaled = '.$scaled.' result ='.$result);
               next if (! defined($scaled));
   #            next if ($result ne '6');
               my $submission = $scaled;
               $students{$subm{'student'}}++;
               $stats{'submission_count'}++;
               if (&numerical_submission_is_correct($subm{'award'})) { 
                   $stats{'correct_count'}++;
                   $submission_data{$subm{'correct'}}->{$submission}->[0]++;
               } elsif (&numerical_submission_is_incorrect($subm{'award'})) { 
                   $stats{'incorrect_count'}++;
                   $submission_data{$subm{'correct'}}->{$submission}->[1]++;
               }
           }
       }
       $stats{'students'}=scalar(keys(%students));
       return (\%submission_data,\%stats);
   }
   
   sub numerical_submission_is_correct {
       my ($award) = @_;
       if ($award =~ /^(APPROX_ANS|EXACT_ANS)$/) {
           return 1;
       } else {
           return 0;
       }
   }
   
   sub numerical_submission_is_incorrect {
       my ($award) = @_;
       if ($award =~ /^(INCORRECT)$/) {
           return 1;
       } else {
           return 0;
       }
   }
   
   sub numerical_bin_table {
       my ($bins,$labels,$incorrect,$correct,$count)=@_;
       my $table = 
           '<table><tr><th>'.&mt('Bar').'</th>'.
           '<th colspan="3">'.&mt('Range').'</th>'.
           '<th>'.&mt('Incorrect').'</th>'.
           '<th>'.&mt('Correct').'</th>'.
           '<th>'.&mt('Count').'</th>'.
           '</tr>'.$/;
       for (my $i=0;$i<scalar(@{$bins});$i++) {
         my $lownum;          my $lownum;
         if ($i == 0) {          if ($i == 0) {
             $lownum = 0;              $lownum = 0;
         } else {          } else {
             $lownum = $bins[$i-1];              $lownum = $bins->[$i-1];
         }          }
         my $highnum = $bins[$i];          my $highnum = $bins->[$i];
         $table .=           $table .= 
             '<tr>'.              '<tr>'.
             '<td>'.$labels[$i].'</td>'.              '<td>'.$labels->[$i].'</td>'.
             '<td align="right">'.$lownum.'</td>'.              '<td align="right">'.$lownum.'</td>'.
             '<td>&nbsp;-&nbsp;</td>'.              '<td>&nbsp;-&nbsp;</td>'.
             '<td align="right">'.$highnum.'</td>'.'</tr>'.$/;              '<td align="right">'.$highnum.'</td>'.
               '<td align="right">'.$incorrect->[$i].'</td>'.
               '<td align="right">'.$correct->[$i].'</td>'.
               '<td align="right">'.$count->[$i].'</td>'.
               '</tr>'.$/;
     }      }
     $table.= '</table>'.$/;      $table.= '</table>';
     return $table;      return $table;
 }  }
   
 sub numerical_response_determine_intervals {  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: 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'});
     #      #
     my $correct;      my $correct;
     # %intervals differs from %answers because it may be possible for two  
     # students to have the same correct answer but different intervals.  
     my %intervals;   
     my %answers;      my %answers;
     foreach my $student (@$students) {      foreach my $student (@$students) {
         last if ($c->aborted());          last if ($c->aborted());
Line 413  sub numerical_response_determine_interva Line 576  sub numerical_response_determine_interva
                                                                 $sdom);                                                                  $sdom);
         # make the key          # make the key
         my $key = $partid.'.'.$respid;          my $key = $partid.'.'.$respid;
         my $interval = join(' ',($analysis->{$key.'.answer'}->[0],          foreach my $item ('answer','unit','ans_high','ans_low') {
                                  $analysis->{$key.'.ans_low'}->[0],              $correct->{$sname.':'.$sdom}->{$item} = 
                                  $analysis->{$key.'.ans_high'}->[0]));                  $analysis->{$key.'.'.$item}->[0];
         $correct->{$sname.':'.$sdom} = $interval;          }
         $intervals{$interval}++;  
         $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();
     return ($correct,\%intervals,\%answers);      &Apache::lonhtmlcommon::Close_PrgWin($r,\%prog_state);
       return ($correct,\%answers);
 }  }
   
   #
   # Inputs: $r, $width, $height, $data
   #         $n = number of students
   #         $data = hashref of $answer => $frequency pairs
 sub numerical_one_dimensional_plot {  sub numerical_one_dimensional_plot {
     my ($r,$width,$height,$n,$data)=@_;      my ($r,$width,$height,$data)=@_;
     #      #
     # Compute data -> image scaling factors      # Compute data -> image scaling factors
     my $max_y = 0;      my $max_y = 0;
     my $min_x = undef;      my $min_x = undef;
     my $max_x = undef;      my $max_x = undef;
       my $n = 0;
     while (my ($answer,$count) = each(%$data)) {      while (my ($answer,$count) = each(%$data)) {
           $n+=$count;
         $max_y = $count if ($max_y < $count);          $max_y = $count if ($max_y < $count);
         if (! defined($min_x) || $answer < $min_x) {          if (! defined($min_x) || $answer < $min_x) {
             $min_x = $answer;              $min_x = $answer;
Line 441  sub numerical_one_dimensional_plot { Line 612  sub numerical_one_dimensional_plot {
         }          }
     }      }
     #      #
     my $h_scale = ($width-10)/($max_x-$min_x);      my $min_max_difference = $max_x - $min_x;
       if (! defined($min_max_difference) || $min_max_difference == 0) {
           $min_max_difference = 1;
       }
       my $h_scale = ($width-10)/$min_max_difference;
     #      #
     my $ticscale = 5;      my $ticscale = 5;
     if ($max_y * $ticscale > $height/2) {      if ($max_y * $ticscale > $height/2) {
Line 471  sub numerical_one_dimensional_plot { Line 646  sub numerical_one_dimensional_plot {
     }      }
     $plot .= '</drawimage>';      $plot .= '</drawimage>';
     my $plotresult =  &Apache::lonxml::xmlparse($r,'web',$plot);      my $plotresult =  &Apache::lonxml::xmlparse($r,'web',$plot);
       
     my $result = '<table>'.      my $result = '<table>'.
         '<tr><td colspan="3" align="center">'.          '<tr><td colspan="3" align="center">'.
         '<font size="+2">'.&mt('Distribution of correct answers').'</font>'.          '<font size="+2">'.&mt('Distribution of correct answers').'</font>'.
Line 522  sub circle { Line 696  sub circle {
     return &arc($x,$y,$radius,$radius,0,360,$color,$thickness,$filled);      return &arc($x,$y,$radius,$radius,0,360,$color,$thickness,$filled);
 }  }
   
 sub build_student_data_worksheet {  
     my ($workbook,$format) = @_;  
     my $rows_output = 3;  
     my $cols_output = 0;  
     my $worksheet  = $workbook->addworksheet('Student Data');  
     $worksheet->write($rows_output++,0,'Student Data',$format->{'h3'});  
     my @Headers = ('full name','username','domain','section',  
                    "student\nnumber",'identifier');  
     $worksheet->write_row($rows_output++,0,\@Headers,$format->{'header'});  
     my @Students = @Apache::lonstatistics::Students;  
     my $studentrows = &Apache::loncoursedata::get_student_data(\@Students);  
     my %ids;  
     foreach my $row (@$studentrows) {  
         my ($mysqlid,$student) = @$row;  
         $ids{$student}=$mysqlid;  
     }  
     foreach my $student (@Students) {  
         my $name_domain = $student->{'username'}.':'.$student->{'domain'};  
         $worksheet->write_row($rows_output++,0,  
                           [$student->{'fullname'},  
                            $student->{'username'},$student->{'domain'},  
                            $student->{'section'},$student->{'id'},  
                            $ids{$name_domain}]);  
     }  
     return $worksheet;  
 }  
   
 #########################################################  #########################################################
 #########################################################  #########################################################
 ##  ##
Line 570  sub radio_response_analysis { Line 717  sub radio_response_analysis {
     my $analysis_html;      my $analysis_html;
     my $foildata = $problem_analysis->{'_Foils'};      my $foildata = $problem_analysis->{'_Foils'};
     my ($table,$foils,$concepts) = &build_foil_index($problem_analysis);      my ($table,$foils,$concepts) = &build_foil_index($problem_analysis);
       if (! defined($concepts)) {
           $concepts = [];
       }
     #      #
     my %true_foils;      my %true_foils;
     my $num_true = 0;      my $num_true = 0;
Line 778  sub RR_concept_plot { Line 928  sub RR_concept_plot {
                 } else {                  } else {
                     $incorrect_by_concept{$concept->{'name'}} +=                      $incorrect_by_concept{$concept->{'name'}} +=
                         $foil_data->{$foil}->{$choice};                          $foil_data->{$foil}->{$choice};
                       
                 }                  }
             }              }
         }          }
Line 965  sub classify_response_data { Line 1114  sub classify_response_data {
         if (&$function(\%subm)) {          if (&$function(\%subm)) {
             $students{$subm{'student'}}++;              $students{$subm{'student'}}++;
             $submission_data{'_count'}++;              $submission_data{'_count'}++;
               
             if (&submission_is_correct($subm{'award'})) {               if (&submission_is_correct($subm{'award'})) { 
                 $submission_data{'_correct'}++;                  $submission_data{'_correct'}++;
             }              }
Line 1004  sub OptionResponseAnalysis { Line 1152  sub OptionResponseAnalysis {
                   '</h2>');                    '</h2>');
     }  else {      }  else {
         $r->rflush();          $r->rflush();
         if (exists($ENV{'form.ExcelOutput'})) {          if ($ENV{'form.AnalyzeOver'} eq 'tries') {
             my $result = &OR_excel_sheet($r,$resource,              my $analysis_html = &OR_tries_analysis($r,
                                          $PerformanceData,                                                     $PerformanceData,
                                          $problem_data);                                                     $problem_data);
             $r->print($result);              $r->print($analysis_html);
               $r->rflush();
           } elsif ($ENV{'form.AnalyzeOver'} eq 'time') {
               my $analysis_html = &OR_time_analysis($PerformanceData,
                                                     $problem_data);
               $r->print($analysis_html);
             $r->rflush();              $r->rflush();
         } else {          } else {
             if ($ENV{'form.AnalyzeOver'} eq 'tries') {              $r->print('<h2>'.
                 my $analysis_html = &OR_tries_analysis($r,                        &mt('The analysis you have selected is '.
                                                     $PerformanceData,                            'not supported at this time').
                                                     $problem_data);                        '</h2>');
                 $r->print($analysis_html);  
                 $r->rflush();  
             } elsif ($ENV{'form.AnalyzeOver'} eq 'time') {  
                 my $analysis_html = &OR_time_analysis($PerformanceData,  
                                                    $problem_data);  
                 $r->print($analysis_html);  
                 $r->rflush();  
             } else {  
                 $r->print('<h2>'.  
                           &mt('The analysis you have selected is '.  
                               'not supported at this time').  
                           '</h2>');  
             }  
         }          }
     }      }
 }  }
Line 1042  sub OR_tries_analysis { Line 1182  sub OR_tries_analysis {
     my $mintries = 1;      my $mintries = 1;
     my $maxtries = $ENV{'form.NumPlots'};      my $maxtries = $ENV{'form.NumPlots'};
     my ($table,$Foils,$Concepts) = &build_foil_index($ORdata);      my ($table,$Foils,$Concepts) = &build_foil_index($ORdata);
       if (! defined($Concepts)) {
           $Concepts = [];
       }
     my %response_data = &OR_analyze_by_tries($r,$PerformanceData,      my %response_data = &OR_analyze_by_tries($r,$PerformanceData,
                                                      $mintries,$maxtries);                                                       $mintries,$maxtries);
     my $analysis = '';      my $analysis = '';
Line 1129  sub OR_tries_analysis { Line 1272  sub OR_tries_analysis {
             #              #
             $concept_graph = &Apache::loncommon::DrawBarGraph              $concept_graph = &Apache::loncommon::DrawBarGraph
                 ('Correct Concepts','Concept Number','Percent Correct',                  ('Correct Concepts','Concept Number','Percent Correct',
                  100,$plotcolors,undef,\@concept_plot_data);                   100,$plotcolors,undef,\@concept_plot_data,{xskip=>1});
         }          }
         #          #
         # Create Foil Plots          # Create Foil Plots
Line 1167  sub OR_tries_analysis { Line 1310  sub OR_tries_analysis {
         #          #
         my $correct_graph = &Apache::loncommon::DrawBarGraph          my $correct_graph = &Apache::loncommon::DrawBarGraph
             ('Correct Statements','Statement','% Answered Correct',              ('Correct Statements','Statement','% Answered Correct',
              100,$plotcolors,\@Labels,$Datasets[0]);               100,$plotcolors,\@Labels,$Datasets[0],{xskip=>1});
                   
         #          #
         #          #
Line 1179  sub OR_tries_analysis { Line 1322  sub OR_tries_analysis {
                                            $response_data{'_correct'}->[$try];                                             $response_data{'_correct'}->[$try];
         my $incorrect_graph = &Apache::loncommon::DrawBarGraph          my $incorrect_graph = &Apache::loncommon::DrawBarGraph
             ('Incorrect Statements','Statement','% Chosen Incorrectly',              ('Incorrect Statements','Statement','% Chosen Incorrectly',
              100,$plotcolors,\@Labels,@Datasets);               100,$plotcolors,\@Labels,@Datasets,{xskip=>1});
         $analysis_html.=           $analysis_html.= 
             '<tr><td colspan="4" align="center">'.              '<tr><td colspan="4" align="center">'.
             '<font size="+1">'.              '<font size="+1">'.
Line 1373  sub OR_Foil_Time_Analysis { Line 1516  sub OR_Foil_Time_Analysis {
                                                         100,                                                          100,
                                                         $plotcolors,                                                          $plotcolors,
                                                         undef,                                                          undef,
                                                         $plotdata[0]);                                                          $plotdata[0],
                                                           {xskip=>1});
     for (my $j=0; $j< scalar(@{$plotdata[0]});$j++) {      for (my $j=0; $j< scalar(@{$plotdata[0]});$j++) {
         $plotdata[0]->[$j]=0;          $plotdata[0]->[$j]=0;
     }      }
Line 1384  sub OR_Foil_Time_Analysis { Line 1528  sub OR_Foil_Time_Analysis {
                                          100,                                           100,
                                          $plotcolors,                                           $plotcolors,
                                          undef,                                           undef,
                                          @plotdata);                                           @plotdata,{xskip=>1});
     return ($correct_plot,$incorrect_plot);      return ($correct_plot,$incorrect_plot);
 }  }
   
Line 1415  sub OR_Concept_Time_Analysis { Line 1559  sub OR_Concept_Time_Analysis {
                                             100,                                              100,
                                             $plotcolors,                                              $plotcolors,
                                             undef,                                              undef,
                                             \@plotdata);                                              \@plotdata,{xskip=>1});
 }  }
   
 sub OR_time_process_data {  sub OR_time_process_data {
Line 1456  sub OR_time_process_data { Line 1600  sub OR_time_process_data {
             scalar(keys(%distinct_students)));              scalar(keys(%distinct_students)));
 }  }
   
 #########################################################  
 #########################################################  
 ##  
 ##             Excel output   
 ##  
 #########################################################  
 #########################################################  
 sub OR_excel_sheet {  
     my ($r,$resource,$performance_data,$ORdata) = @_;  
     my $response = '';  
     my (undef,$Foils,$Concepts) = &build_foil_index($ORdata);  
     #  
     # Create excel worksheet  
     my $filename = '/prtspool/'.  
         $ENV{'user.name'}.'_'.$ENV{'user.domain'}.'_'.  
         time.'_'.rand(1000000000).'.xls';  
     my $workbook  = Spreadsheet::WriteExcel->new('/home/httpd'.$filename);  
     if (! defined($workbook)) {  
         $r->log_error("Error creating excel spreadsheet $filename: $!");  
         $r->print('<p>'.&mt("Unable to create new Excel file.  ".  
                             "This error has been logged.  ".  
                             "Please alert your LON-CAPA administrator").  
                   '</p>');  
         return undef;  
     }  
     #  
     $workbook->set_tempdir('/home/httpd/perl/tmp');  
     my $format = &Apache::loncommon::define_excel_formats($workbook);  
     #  
     # Create and populate main worksheets  
     my $problem_data_sheet  = $workbook->addworksheet('Problem Data');  
     my $student_data_sheet = &build_student_data_worksheet($workbook,$format);  
     my $response_data_sheet = $workbook->addworksheet('Response Data');  
     foreach my $sheet ($problem_data_sheet,$student_data_sheet,  
                        $response_data_sheet) {  
         $sheet->write(0,0,$resource->{'title'},$format->{'h2'});  
         $sheet->write(1,0,$resource->{'src'},$format->{'h3'});  
     }  
     #  
     my $result;  
     $result = &OR_build_problem_data_worksheet($problem_data_sheet,$format,  
                                             $Concepts,$ORdata);  
     if ($result ne 'okay') {  
         # Do something useful  
     }  
     $result = &OR_build_response_data_worksheet($response_data_sheet,$format,  
                                              $performance_data,$Foils,  
                                              $ORdata);  
     if ($result ne 'okay') {  
         # Do something useful  
     }  
     $response_data_sheet->activate();  
     #  
     # Close the excel file  
     $workbook->close();  
     #  
     # Write a link to allow them to download it  
     $result .= '<h2>'.&mt('Excel Raw Data Output').'</h2>'.  
               '<p><a href="'.$filename.'">'.  
               &mt('Your Excel spreadsheet.').  
               '</a></p>'."\n";  
     return $result;  
 }  
   
 sub OR_build_problem_data_worksheet {  
     my ($worksheet,$format,$Concepts,$ORdata) = @_;  
     my $rows_output = 3;  
     my $cols_output = 0;  
     $worksheet->write($rows_output++,0,'Problem Structure',$format->{'h3'});  
     ##  
     ##  
     my @Headers;  
     if (@$Concepts > 1) {  
         @Headers = ("Concept\nNumber",'Concept',"Foil\nNumber",  
                     'Foil Name','Foil Text','Correct value');  
     } else {  
         @Headers = ('Foil Number','FoilName','Foil Text','Correct value');  
     }  
     $worksheet->write_row($rows_output++,0,\@Headers,$format->{'header'});  
     my %Foildata = %{$ORdata->{'_Foils'}};  
     my $conceptindex = 1;  
     my $foilindex = 1;  
     foreach my $concept (@$Concepts) {  
         my @FoilsInConcept = @{$concept->{'foils'}};  
         my $firstfoil = shift(@FoilsInConcept);  
         if (@$Concepts > 1) {  
             $worksheet->write_row($rows_output++,0,  
                                   [$conceptindex,  
                                    $concept->{'name'},  
                                    $foilindex++,  
                                    $Foildata{$firstfoil}->{'name'},  
                                    $Foildata{$firstfoil}->{'text'},  
                                    $Foildata{$firstfoil}->{'value'},]);  
         } else {  
             $worksheet->write_row($rows_output++,0,  
                                   [ $foilindex++,  
                                     $Foildata{$firstfoil}->{'name'},  
                                     $Foildata{$firstfoil}->{'text'},  
                                     $Foildata{$firstfoil}->{'value'},]);  
         }  
         foreach my $foilid (@FoilsInConcept) {  
             if (@$Concepts > 1) {  
                 $worksheet->write_row($rows_output++,0,  
                                       ['',  
                                        '',  
                                        $foilindex,  
                                        $Foildata{$foilid}->{'name'},  
                                        $Foildata{$foilid}->{'text'},  
                                        $Foildata{$foilid}->{'value'},]);  
             } else {  
                 $worksheet->write_row($rows_output++,0,                  
                                       [$foilindex,  
                                        $Foildata{$foilid}->{'name'},  
                                        $Foildata{$foilid}->{'text'},  
                                        $Foildata{$foilid}->{'value'},]);  
             }                  
         } continue {  
             $foilindex++;  
         }  
     } continue {  
         $conceptindex++;  
     }  
     $rows_output++;  
     $rows_output++;  
     ##  
     ## Option data output  
     $worksheet->write($rows_output++,0,'Options',$format->{'header'});  
     foreach my $string (@{$ORdata->{'_Options'}}) {  
         $worksheet->write($rows_output++,0,$string);  
     }  
     return 'okay';  
 }  
   
 sub OR_build_response_data_worksheet {  
     my ($worksheet,$format,$performance_data,$Foils,$ORdata)=@_;  
     my $rows_output = 3;  
     my $cols_output = 0;  
     $worksheet->write($rows_output++,0,'Response Data',$format->{'h3'});  
     $worksheet->set_column(1,1,20);  
     $worksheet->set_column(2,2,13);  
     my @Headers = ('identifier','time','award detail','attempt');  
     foreach my $foil (@$Foils) {  
         push (@Headers,$foil.' submission');  
         push (@Headers,$foil.' grading');  
     }  
     $worksheet->write_row($rows_output++,0,\@Headers,$format->{'header'});  
     #  
     foreach my $row (@$performance_data) {  
         next if (! defined($row));  
         my ($student,$award,$grading,$submission,$time,$tries) = @$row;  
         my @Foilgrades = split('&',$grading);  
         my @Foilsubs   = split('&',$submission);  
         my %response_data;  
         for (my $j=0;$j<=$#Foilgrades;$j++) {  
             my ($foilid,$correct)  = split('=',$Foilgrades[$j]);  
             my (undef,$submission) = split('=',$Foilsubs[$j]);  
             $submission = &Apache::lonnet::unescape($submission);  
             $response_data{$foilid.' submission'}=$submission;  
             $response_data{$foilid.' award'}=$correct;  
         }  
         $worksheet->write($rows_output,$cols_output++,$student);  
         $worksheet->write($rows_output,$cols_output++,  
              &Apache::lonstathelpers::calc_serial($time),$format->{'date'});  
         $worksheet->write($rows_output,$cols_output++,$award);  
         $worksheet->write($rows_output,$cols_output++,$tries);  
         foreach my $foilid (@$Foils) {  
             $worksheet->write($rows_output,$cols_output++,  
                               $response_data{$foilid.' submission'});  
             $worksheet->write($rows_output,$cols_output++,  
                               $response_data{$foilid.' award'});  
         }  
         $rows_output++;  
         $cols_output = 0;  
     }  
     return;  
 }  
   
 sub build_foil_index {  sub build_foil_index {
     my ($ORdata) = @_;      my ($ORdata) = @_;
     return if (! exists($ORdata->{'_Foils'}));      return if (! exists($ORdata->{'_Foils'}));

Removed from v.1.103  
changed lines
  Added in v.1.112


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