Diff for /loncom/interface/statistics/lonproblemstatistics.pm between versions 1.16 and 1.122.2.2

version 1.16, 2002/08/06 17:38:26 version 1.122.2.2, 2013/12/30 14:05:21
Line 1 Line 1
 # The LearningOnline Network with CAPA  # The LearningOnline Network with CAPA
 # (Publication Handler  
 #  #
 # $Id$  # $Id$
 #  #
Line 26 Line 25
 # http://www.lon-capa.org/  # http://www.lon-capa.org/
 #  #
 # (Navigate problems for statistical reports  # (Navigate problems for statistical reports
 # YEAR=2001  
 # 5/5,7/9,7/25/1,8/11,9/13,9/26,10/5,10/9,10/22,10/26 Behrouz Minaei  
 # 11/1,11/4,11/16,12/14,12/16,12/18,12/20,12/31 Behrouz Minaei  
 # YEAR=2002  
 # 1/22,2/1,2/6,2/25,3/2,3/6,3/17,3/21,3/22,3/26,4/7,5/6 Behrouz Minaei  
 # 5/12,5/14,5/15,5/19,5/26,7/16,7/25,7/29,8/5  Behrouz Minaei  
 #  #
 ###  ###############################################
   ###############################################
   
 package Apache::lonproblemstatistics;   =pod
   
   =head1 NAME
   
   lonproblemstatistics
   
   =head1 SYNOPSIS
   
   Routines to present problem statistics to instructors via tables,
   Excel files, and plots.
   
   =over 4
   
   =cut
   
   ###############################################
   ###############################################
   
   package Apache::lonproblemstatistics;
   
 use strict;  use strict;
 use Apache::lonnet();  use Apache::lonnet;
   use Apache::loncommon();
 use Apache::lonhtmlcommon;  use Apache::lonhtmlcommon;
 use Apache::loncoursedata;  use Apache::loncoursedata;
 use GDBM_File;  use Apache::lonstatistics;
   use LONCAPA::lonmetadata();
   use Apache::lonlocal;
   use Spreadsheet::WriteExcel;
   use Apache::lonstathelpers();
   use Time::HiRes;
   use LONCAPA;
    
   
 sub BuildProblemStatisticsPage {  my @StatsArray;
     my ($cacheDB, $students, $courseID, $c, $r)=@_;  my %SeqStat;    # keys are symbs, values are hash refs
     my %cache;  
   
     unless(tie(%cache,'GDBM_File',$cacheDB,&GDBM_READER,0640)) {  ##
         return '<html><body>Unable to tie database.</body></html>';  ## Localization notes:
   ##
   ## in @Fields[0]->{'long_title'} is placed in Excel files and is used as the
   ## header for plots created with Graph.pm, both of which more than likely do
   ## not support localization.
   ##
   ## Additional Notes:
   ## Localization can be done and is done before passing the phrases
   ## to the output.
   ## This might conflict with special characters, e.g. German Umlaute or
   ## chinese characters. Do not use such characters in this case.
   ## If this failed, consider that the sequence and folder names
   ## are also passed to the output and would fail the same way.
   ##
   #
   #
   ##
   ## Description of Field attributes
   ##
   ## Attribute     Required   Value       Meaning or Use
   ##
   ## name            yes      any scalar  Used to uniquely identify field
   ## title           yes      any scalar  This is what the user sees to identify
   ##                                      the field.  Passed through &mt().
   ## long_title      yes      any scalar  Used as graph heading and in excel
   ##                                      output.  Passed through &mt().
   ## align           no    (left|right|center)  HTML cell contents alignment
   ## color           yes      html color  HTML cell background color
   ##                                      used to visually group statistics
   ## special         no          (link)   Indicates a link, target is name.link
   ##                                      Currently set in &get_statistics()
   ## graphable       no      (yes|no)     Can a bar graph of the field be 
   ##                                      produced?
   ## sortable        no      (yes|no)     Should a sort link be put in the
   ##                                      column header?
   ## selectable      yes     (yes|no)     Can the column be removed from the
   ##                                      statistics display?
   ## selected        yes     (yes|no)     Is the column selected by default?
   ##
   ## format          no      sprintf format string
   ##
   ## excel_format    no      excel format type 
   ##                               (see &Apache::loncommon::define_excel_formats
   my @Fields = (
              { name => 'problem_num',
                title => 'P#',
                align => 'right',
                color => '#FFFFE6',
                selectable => 'no',
                defaultselected => 'yes',
              },
              { name   => 'container',
                title  => 'Sequence or Folder',
                align  => 'left',
                color  => '#FFFFE6',
                sortable => 'yes',
                selectable => 'no',
                defaultselected => 'yes',
              },
              { name   => 'title',
                title  => 'Title',
                align  => 'left',
                color  => '#FFFFE6',
                special  => 'link',
                sortable => 'yes',
                selectable => 'no',
                defaultselected => 'yes',
              },
              { name   => 'part',
                title  => 'Part',
                align  => 'left',
                color  => '#FFFFE6',
                selectable => 'no',
                defaultselected => 'yes',
              },
              { name   => 'num_students',
                title  => '#Stdnts',
                align  => 'right',
                color  => '#EEFFCC',
                format => '%d',
                sortable  => 'yes',
                graphable => 'yes',
                long_title => 'Number of Students Attempting Problem',
                selectable => 'yes',
                defaultselected => 'yes',
              },
              { name   => 'tries',
                title  => 'Tries',
                align  => 'right',
                color  => '#EEFFCC',
                format => '%d',
                sortable  => 'yes',
                graphable => 'yes',
                long_title => 'Total Number of Tries',
                selectable => 'yes',
                defaultselected => 'yes',
              },
              { name   => 'max_tries',
                title  => 'Max Tries',
                align  => 'right',
                color  => '#DDFFFF',
                format => '%d',
                sortable  => 'yes',
                graphable => 'yes',
                long_title => 'Maximum Number of Tries',
                selectable => 'yes',
                defaultselected => 'yes',
              },
              { name   => 'min_tries',
                title  => 'Min Tries',
                align  => 'right',
                color  => '#DDFFFF',
                format => '%d',
                sortable  => 'yes',
                graphable => 'yes',
                long_title => 'Minumum Number of Tries',
                selectable => 'yes',
                defaultselected => 'yes',
              },
              { name   => 'mean_tries',
                title  => 'Mean Tries',
                align  => 'right',
                color  => '#DDFFFF',
                format => '%5.2f',
                sortable  => 'yes',
                graphable => 'yes',
                long_title => 'Average Number of Tries',
                selectable => 'yes',
                defaultselected => 'yes',
              },
              { name   => 'std_tries',
                title  => 'S.D. tries',
                align  => 'right',
                color  => '#DDFFFF',
                format => '%5.2f',
                sortable  => 'yes',
                graphable => 'yes',
                long_title => 'Standard Deviation of Number of Tries',
                selectable => 'yes',
                defaultselected => 'yes',
              },
              { name   => 'skew_tries',
                title  => 'Skew Tries',
                align  => 'right',
                color  => '#DDFFFF',
                format => '%5.2f',
                sortable  => 'yes',
                graphable => 'yes',
                long_title => 'Skew of Number of Tries',
                selectable => 'yes',
                defaultselected => 'no',
              },
              { name   => 'num_solved',
                title  => '#YES',
                align  => 'right',
                color  => '#FFDDDD',
                format => '%4.1f',#             format => '%d',
                sortable  => 'yes',
                graphable => 'yes',
                long_title => 'Number of Students able to Solve',
                selectable => 'yes',
                defaultselected => 'yes',
              },
              { name   => 'num_override',
                title  => '#yes',
                align  => 'right',
                color  => '#FFDDDD',
                format => '%4.1f',#             format => '%d',
                sortable  => 'yes',
                graphable => 'yes',
                long_title => 'Number of Students given Override',
                selectable => 'yes',
                defaultselected => 'yes',
              },
              { name   => 'tries_per_correct',
                title  => 'tries/correct',
                align  => 'right',
                color  => '#FFDDDD',
                format => '%4.1f',
                sortable  => 'yes',
                graphable => 'yes',
                long_title => 'Tries per Correct Answer',
                selectable => 'yes',
                defaultselected => 'yes',
              },
              { name   => 'num_wrong',
                title  => '#Wrng',
                align  => 'right',
                color  => '#FFDDDD',
                format => '%4.1f',
                sortable  => 'yes',
                graphable => 'yes',
                long_title => 'Number of students whose final answer is wrong',
                selectable => 'yes',
                defaultselected => 'yes',
              },
              { name   => 'per_wrong',
                title  => '%Wrng',
                align  => 'right',
                color  => '#FFDDDD',
                format => '%4.1f',
                sortable  => 'yes',
                graphable => 'yes',
                long_title => 'Percent of students whose final answer is wrong',
                selectable => 'yes',
                defaultselected => 'yes',
              },
              { name   => 'deg_of_diff',
                title  => 'DoDiff',
                align  => 'right',
                color  => '#FFFFE6',
                format => '%5.2f',
                sortable  => 'yes',
                graphable => 'yes',
                long_title => 'Degree of Difficulty'.
                              ' ~[ 1 - ((#YES+#yes) / Tries) ~]',
                selectable => 'yes',
                defaultselected => 'yes',
              },
              { name   => 'deg_of_disc',
                title  => 'DoDisc',
                align  => 'right',
                color  => '#FFFFE6',
                format => '%4.2f',
                sortable  => 'yes',
                graphable => 'yes',
                long_title => 'Degree of Discrimination',
                selectable => 'yes',
                defaultselected => 'yes',
              },
   ##   duedate included for research purposes.  Commented out most of the time.
   #           { name => 'duedate',
   #             title => 'Due Date',
   #             align => 'left',
   #             color => '#FFFFFF',
   #             sortable => 'yes',
   #             graphable => 'no',
   #             long_title => 'Due date of resource for instructor',
   #             selectable => 'no',
   #             defaultselected => 'yes',
   #            },
   ##   opendate included for research purposes.  Commented out most of the time.
   #           { name => 'opendate',
   #             title => 'Open Date',
   #             align => 'left',
   #             color => '#FFFFFF',
   #             sortable => 'yes',
   #             graphable => 'no',
   #             long_title => 'date resource became answerable',
   #             selectable => 'no',
   #             defaultselected => 'yes',
   #            },
   ##   symb included for research purposes.  Commented out most of the time.
   #           { name => 'symb',
   #             title => 'Symb',
   #             align => 'left',
   #             color => '#FFFFFF',
   #             sortable => 'yes',
   #             graphable => 'no',
   #             long_title => 'Unique LON-CAPA identifier for problem',
   #             selectable => 'no',
   #             defaultselected => 'yes',
   #            },
   ##   resptypes included for research purposes.  Commented out most of the time.
   #           { name => 'resptypes',
   #             title => 'Response Types',
   #             align => 'left',
   #             color => '#FFFFFF',
   #             sortable => 'no',
   #             graphable => 'no',
   #             long_title => 'Response Types used in this problem',
   #             selectable => 'no',
   #             defaultselected => 'yes',
   #            },
   ##   maxtries included for research purposes.  Commented out most of the time.
   #           { name => 'maxtries',
   #             title => 'Maxtries',
   #             align => 'left',
   #             color => '#FFFFFF',
   #             sortable => 'no',
   #             graphable => 'no',
   #             long_title => 'Maximum number of tries',
   #             selectable => 'no',
   #             defaultselected => 'yes',
   #            },
   ##   hinttries included for research purposes.  Commented out most of the time.
   #           { name => 'hinttries',
   #             title => 'hinttries',
   #             align => 'left',
   #             color => '#FFFFFF',
   #             sortable => 'no',
   #             graphable => 'no',
   #             long_title => 'Number of tries before a hint appears',
   #             selectable => 'no',
   #             defaultselected => 'yes',
   #            },
   #
   ##   problem weight for instructor
              { name => 'weight',
                title => 'weight',
                align => 'right',
                color => '#FFFFFF',
                sortable => 'no',
                graphable => 'no',
                long_title => 'Problem weight (for instructor)',
                selectable => 'yes',
                defaultselected => 'yes',
               },
   );
   
   my @SeqFields = (
              { name   => 'title',
                title  => 'Sequence',
                align  => 'left',
                color  => '#FFFFE6',
                special  => 'no',
                sortable => 'no',
                selectable => 'yes',
                defaultselected => 'no',
              },
              { name   => 'items',
                title  => '#Items',
                align  => 'right',
                color  => '#FFFFE6',
                format => '%4d',
                sortable  => 'no',
                graphable => 'no',
                long_title => 'Number of Items in Sequence',
                selectable => 'yes',
                defaultselected => 'no',
              },
              { name   => 'scoremean',
                title  => 'Score Mean',
                align  => 'right',
                color  => '#FFFFE6',
                format => '%4.2f',
                sortable  => 'no',
                graphable => 'no',
                long_title => 'Mean Sequence Score',
                selectable => 'yes',
                defaultselected => 'no',
              },
              { name   => 'scorestd',
                title  => 'Score STD',
                align  => 'right',
                color  => '#FFFFE6',
                format => '%4.2f',
                sortable  => 'no',
                graphable => 'no',
                long_title => 'Standard Deviation of Sequence Scores',
                selectable => 'yes',
                defaultselected => 'no',
              },
              { name   => 'scoremax',
                title  => 'Score Max',
                align  => 'right',
                color  => '#FFFFE6',
                format => '%4.2f',
                sortable  => 'no',
                graphable => 'no',
                long_title => 'Maximum Sequence Score',
                selectable => 'yes',
                defaultselected => 'no',
              },
              { name   => 'scoremin',
                title  => 'Score Min',
                align  => 'right',
                color  => '#FFFFE6',
                format => '%4.2f',
                sortable  => 'no',
                graphable => 'no',
                long_title => 'Minumum Sequence Score',
                selectable => 'yes',
                defaultselected => 'no',
              },
              { name   => 'scorecount',
                title  => 'Score N',
                align  => 'right',
                color  => '#FFFFE6',
                format => '%4d',
                sortable  => 'no',
                graphable => 'no',
                long_title => 'Number of Students in score computations',
                selectable => 'yes',
                defaultselected => 'no',
              },
              { name   => 'countmean',
                title  => 'Count Mean',
                align  => 'right',
                color  => '#FFFFFF',
                format => '%4.2f',
                sortable  => 'no',
                graphable => 'no',
                long_title => 'Mean Sequence Score',
                selectable => 'yes',
                defaultselected => 'no',
              },
              { name   => 'countstd',
                title  => 'Count STD',
                align  => 'right',
                color  => '#FFFFFF',
                format => '%4.2f',
                sortable  => 'no',
                graphable => 'no',
                long_title => 'Standard Deviation of Sequence Scores',
                selectable => 'yes',
                defaultselected => 'no',
              },
              { name   => 'countmax',
                title  => 'Count Max',
                align  => 'right',
                color  => '#FFFFFF',
                format => '%4.2f',
                sortable  => 'no',
                graphable => 'no',
                long_title => 'Maximum Number of Correct Problems',
                selectable => 'yes',
                defaultselected => 'no',
              },
              { name   => 'countmin',
                title  => 'Count Min',
                align  => 'right',
                color  => '#FFFFFF',
                format => '%4.2f',
                sortable  => 'no',
                graphable => 'no',
                long_title => 'Minumum Number of Correct Problems',
                selectable => 'yes',
                defaultselected => 'no',
              },
              { name   => 'count',
                title  => 'Count N',
                align  => 'right',
                color  => '#FFFFFF',
                format => '%4d',
                sortable  => 'no',
                graphable => 'no',
                long_title => 'Number of Students in score computations',
                selectable => 'yes',
                defaultselected => 'no',
              },
              { name   => 'KR-21',
                title  => 'KR-21',
                align  => 'right',
                color  => '#FFAAAA',
                format => '%4.2f',
                sortable  => 'no',
                graphable => 'no',
                long_title => 'KR-21 reliability statistic',
                selectable => 'yes',
                defaultselected => 'no',
              },
   );
   
   my %SelectedFields;
   
   sub parse_field_selection {
       #
       # Pull out the defaults
       if (! defined($env{'form.fieldselections'})) {
           $env{'form.fieldselections'} = [];
           foreach my $field (@Fields) {
               next if ($field->{'selectable'} ne 'yes');
               if ($field->{'defaultselected'} eq 'yes') {
                   push(@{$env{'form.fieldselections'}},$field->{'name'});
               }
           }
     }      }
       #
     my $Ptr = '';      # Make sure the data we are plotting is there
     $Ptr .= '<table border="0"><tbody>';      my %NeededFields;
     $Ptr .= '<tr><td align="right"><b>Select Map</b></td>'."\n";      if (exists($env{'form.plot'}) && $env{'form.plot'} ne '' &&
     $Ptr .= '<td align="left">';          $env{'form.plot'} ne 'none') {
     $Ptr .= &Apache::lonhtmlcommon::MapOptions(\%cache, 'ProblemStatistics',          if ($env{'form.plot'} eq 'degrees') {
                                                'Statistics');              $NeededFields{'deg_of_diff'}++;
     $Ptr .= '</td></tr>'."\n";              $NeededFields{'deg_of_disc'}++;
     $Ptr .= '<tr><td align="right"><b>Sorting Type:</b></td>'."\n";          } elsif ($env{'form.plot'} eq 'tries statistics') {
     $Ptr .= '<td align="left">'."\n";              $NeededFields{'mean_tries'}++;
     $Ptr .= &Apache::lonhtmlcommon::AscendOrderOptions(              $NeededFields{'std_tries'}++;
                                             $cache{'ProblemStatisticsAscend'},               $NeededFields{'problem_num'}++;
                                             'ProblemStatistics',          } else {
                                             'Statistics');              $NeededFields{$env{'form.plot'}}++;
     $Ptr .= '</td></tr>'."\n";          }
     $Ptr .= &ProblemStatisticsButtons($cache{'DisplayFormat'});  
     $Ptr .= '</table>';  
     $Ptr .= &ProblemStatisticsLegend();  
     $r->print($Ptr);  
     $r->rflush();  
   
     untie(%cache);  
     &Create_PrgWin($r);  
     my $count=0;  
     foreach (@$students) {  
         &Update_PrgWin(scalar(@$students),$count,$_,$r);  
         my $courseData =   
             &Apache::loncoursedata::DownloadCourseInformation($_, $courseID);  
         last if ($c->aborted());  
         if(tie(%cache,'GDBM_File',$cacheDB,&GDBM_WRCREAT,0640)) {  
             &Apache::loncoursedata::ProcessStudentData(\%cache,   
                                                        $courseData, $_);  
             untie(%cache);  
         }  
         $count++;  
     }  
     &Close_PrgWin($r);  
     if($c->aborted()) { return; }  
   
     unless(tie(%cache,'GDBM_File',$cacheDB,&GDBM_READER,0640)) {  
         return '<html><body>Unable to tie database.</body></html>';  
     }  
   
     my %Header = (0,"Homework Sets Order",1,"#Stdnts",2,"Tries",3,"Mod",  
                   4,"Mean",5,"#YES",6,"#yes",7,"%Wrng",8,"DoDiff",  
                   9,"S.D.",10,"Skew.",11,"D.F.1st",12,"D.F.2nd", 13, "Disc.");  
     my $color=&setbgcolor(0);  
   
     my $state=$ENV{'form.ProblemStatisticsHeading'};   
   
     my $TempCache;  
   
     if ($state) {  
  $TempCache=&CacheStatisticsTable($state,\%cache,\%Header,  
  $r,$color);  
     } else {      
  my %discriminant=();  
  my @list=();  
  my %Discuss=&Apache::loncoursedata::LoadDiscussion($courseID);  
  my $index=0;  
  foreach (@$students) {  
     $index++;  
     &ExtractStudentData(\%cache, $_, \@list,\%Discuss, $r,  
                                 \%discriminant);  
  }  
  my ($upper, $lower) = &Discriminant(\%discriminant,$r);  
  $TempCache= &BuildStatisticsTable(\%cache, $upper, $lower,   
    \@list, \%Header, $students,  
    $r, $color);  
     }  
     untie(%cache);  
   
     foreach (keys %$TempCache) {  
         last if ($c->aborted());  
         if(tie(%cache,'GDBM_File',$cacheDB,&GDBM_WRCREAT,0640)) {  
     $cache{$_}=$TempCache->{$_};  
   
             untie(%cache);  
         }  
     }  
     if($c->aborted()) { return; }  
     untie(%cache);  
 }  
   
   
 #---- Problem Statistics Web Page ---------------------------------------  
   
 sub CreateProblemStatisticsTableHeading {  
     my ($displayFormat,$sequenceSource,$sequenceTitle,$headings,$r)=@_;  
     if($displayFormat eq 'Display CSV Format') {  
         $r->print('<br>"'.$sequenceTitle.'","');  
         $r->print($sequenceSource.'"');  
  return;  
     }      }
     if ($sequenceSource eq 'Sorted by: ') {      #
  $r->print('<br><b>'.$sequenceSource.$sequenceTitle.'</b>');      # This should not happen, but in case it does...
     } else {      if (ref($env{'form.fieldselections'}) ne 'ARRAY') {
         $r->print('<br><a href="'.$sequenceSource.          $env{'form.fieldselections'} = [$env{'form.fieldselections'}];
   '" target="_blank">'.$sequenceTitle.'</a>');      }
       #
       # Set the field data and the selected fields (for easier checking)
       undef(%SelectedFields);
       foreach my $field (@Fields) {
           if ($field->{'selectable'} ne 'yes') {
               $field->{'selected'} = 'yes';
           } else {
               $field->{'selected'} = 'no';
           }
           if (exists($NeededFields{$field->{'name'}})) {
               $field->{'selected'} = 'yes';
               $SelectedFields{$field->{'name'}}++;
           }
           foreach my $selection (@{$env{'form.fieldselections'}}) {
               if ($selection eq $field->{'name'} || $selection eq 'all') {
                   $field->{'selected'} = 'yes';
                   $SelectedFields{$field->{'name'}}++;
               }
           }
     }      }
     my $Result = "\n".'<table border=2><tr><th>P#</th>'."\n";      #
     for(my $nIndex=0; $nIndex < (scalar (keys %$headings)); $nIndex++) {       # Always show all the sequence statistics (for now)
  $Result .= '<th>'.'<input type="submit" name="';      foreach my $field (@SeqFields) {
         $Result .= 'ProblemStatisticsHeading" value="';          $field->{'selected'} = 'yes';
         $Result .= $headings->{$nIndex}.'" />'.'</th>'."\n";  
     }      }
     $Result .= "\n".'</tr>'."\n";          return;
     $r->print($Result);  
     $r->rflush();  
 }  }
   
 sub CloseTable {  sub field_selection_input {
     my ($cache,$r)=@_;      my $Str = '<select name="fieldselections" multiple="multiple" size="5">'."\n";
     if($cache->{'DisplayFormat'} eq 'Display CSV Format') {      $Str .= '<option value="all">all</option>'."\n";
  return;      foreach my $field (@Fields) {
     }              next if ($field->{'selectable'} ne 'yes');
     $r->print("\n".'</table>'."\n");          $Str .= '    <option value="'.$field->{'name'}.'" ';
     $r->rflush();          if ($field->{'selected'} eq 'yes') {
               $Str .= 'selected ';
           }
           $Str .= '>'.$field->{'title'}.'</option>'."\n";
       }
       $Str .= "</select>\n";
 }  }
   
   ###############################################
   ###############################################
   
 # Create progress  =pod 
 sub Create_PrgWin {  
     my ($r)=@_;  
     $r->print(<<ENDPOP);  
     <script>  
     popwin=open('','popwin','width=400,height=100');  
     popwin.document.writeln('<html><body bgcolor="#88DDFF">'+  
       '<title>LON-CAPA Statistics</title>'+  
       '<h4>Computation Progress</h4>'+  
       '<form name=popremain>'+  
       '<input type=text size=35 name=remaining value=Starting></form>'+  
       '</body></html>');  
     popwin.document.close();  
     </script>  
 ENDPOP  
   
     $r->rflush();  =item &CreateInterface()
 }  
   
 # update progress  Create the main intereface for the statistics page.  Allows the user to
 sub Update_PrgWin {  select sections, maps, and output.
     my ($totalStudents,$index,$name,$r)=@_;  
     $r->print('<script>popwin.document.popremain.remaining.value="'.  =cut
               'Computing '.$index.'/'.$totalStudents.': '.  
               $name.'";</script>');  ###############################################
     $r->rflush();  ###############################################
   sub CreateInterface {
       my ($r) = @_;
       #
       &parse_field_selection();
       #
       my $Str = '';
       $Str .= '<p>';
       $Str .= &Apache::loncommon::start_data_table();
       $Str .= &Apache::loncommon::start_data_table_header_row();
       $Str .= '<th>'.&mt('Sections').'</th>';
       $Str .= '<th>'.&mt('Groups').'</th>';
       $Str .= '<th>'.&mt('Access Status').'</th>';
       $Str .= '<th>'.&mt('Sequences and Folders').'</th>';
       $Str .= '<th>'.&mt('Statistics').'</th>';
       $Str .= '<th>'.&mt('Plot Graph').'</th>';
       $Str .= '<th>'.&mt('Time Period').'</th>';
       $Str .= &Apache::loncommon::end_data_table_header_row();
       #
       $Str .= &Apache::loncommon::start_data_table_row();
       $Str .= '<td align="center" valign="top">'."\n";
       $Str .= &Apache::lonstatistics::SectionSelect('Section','multiple',5);
       $Str .= '</td><td align="center" valign="top">';
       $Str .= &Apache::lonstatistics::GroupSelect('Group','multiple',5);
       $Str .= '</td><td align="center" valign="top">';
       $Str .= &Apache::lonhtmlcommon::StatusOptions(undef,undef,5);
       $Str .= '</td><td align="center" valign="top">';
       #
       $Str .= &Apache::lonstatistics::map_select('Maps','multiple,all',5);
       $Str .= '</td><td align="center" valign="top">';
       $Str .= &field_selection_input();
       $Str .= '</td><td align="center" valign="top">';
       $Str .= &plot_dropdown();
       $Str .= '</td>'."\n";
       $Str .= '</td><td align="center" valign="top">';
       $Str .= &Apache::lonstathelpers::limit_by_time_form();
       $Str .= '</td>'."\n";
       $Str .=  &Apache::loncommon::end_data_table_row();
       $Str .= &Apache::loncommon::end_data_table();
       #
       $Str .= '</p>';
       $Str .= '<input type="submit" name="GenerateStatistics" value="'.
           &mt('Generate Statistics').'" />';
       $Str .= ('&nbsp;'x10);
       #
       return $Str;
 }  }
   
 # close Progress Line  ###############################################
 sub Close_PrgWin {  ###############################################
     my ($r)=@_;  
     $r->print('<script>popwin.close()</script>');  
     $r->rflush();   
 }  
   
    =pod 
 # ------ Dump the Student's DB file and handling the data for statistics table   
 sub ExtractStudentData {  
     my ($cache,$name,$list,$Discuss,$r,$discriminant)=@_;  
     my $totalTries = 0;  
     my $totalAwarded = 0;  
     my $spent=0;  
     my $spent_yes=0;  
     my $TotDiscuss=0;  
     my $TotalOpend = 0;  
     my $ProbSolved = 0;  
     my $ProbTot = 0;  
     my $TotFirst = 0;  
     my $TimeTot = 0;  
     my $Discussed=0;  
     my $discrim='';  
     my $tempSequenceOrder=100;  
   
 #$Apache::lonxml::debug=1;  
 #&Apache::lonhomework::showhash(%$cache);  
 #$Apache::lonxml::debug=0;  
   
     foreach my $sequence (split(':', $cache->{'orderedSequences'})) {  
  my $tempProblemOrder=100;  
         if($cache->{'ProblemStatisticsMaps'} ne 'All Maps'  &&  
            $cache->{'ProblemStatisticsMaps'} ne $cache->{$sequence.':title'}) {  
             next;  
         }  
  $tempSequenceOrder++;  
   
         foreach my $problemID (split(':', $cache->{$sequence.':problems'})) {  =item &BuildProblemStatisticsPage()
             my $problem = $cache->{$problemID.':problem'};  
             my $LatestVersion = $cache->{$name.':version:'.$problem};  
     # Output dashes for all the parts of this problem if there  
             # is no version information about the current problem.  
             #if(!$LatestVersion) {  
             #    foreach my $part (split(/\:/,$cache->{$sequence.':'.  
             #                                          $problemID.  
             #                                          ':parts'})) {  
             #        $codes    .= "-,";  
             #        $attempts .= "0,";   
             #    }  
             #    next;  
             #}  
   
             my %partData=undef;  
     $partData{'count'}=0;  
             # Initialize part data, display skips correctly  
             # Skip refers to when a student made no submissions on that  
             # part/problem.  
             foreach my $part (split(/\:/,$cache->{$sequence.':'.  
                                                   $problemID.  
                                                   ':parts'})) {  
  $tempProblemOrder++;  
  $partData{'count'}++;  
                 $partData{$part.':order'}=$tempProblemOrder;  
                 $partData{$part.':tries'}=0;  
                 $partData{$part.':code'}='-';  
             }  
   
             # Looping through all the versions of each part, starting with the  Main interface to problem statistics.
             # oldest version.  Basically, it gets the most recent   
             # set of grade data for each part.  
     for(my $Version=1; $Version<=$LatestVersion; $Version++) {  
  foreach my $part (split(/\:/,$cache->{$sequence.':'.  
                                                       $problemID.  
                                                       ':parts'})) {  
   
                     if(!defined($cache->{$name.":$Version:$problem".  
                                                ":resource.$part.solved"})) {  
                         # No grade for this submission, so skip  
                         next;  
                     }  
   
                     my $tries=0;  =cut
                     my $time=0;  
                     my $awarded=0;  ###############################################
     $Discussed=0;  ###############################################
                     my $code='-';  my $navmap;
   my @sequences;
                     $awarded = $cache->{"$name:$Version:$problem:resource.".  
                                         "$part.awarded"};  sub clean_up {
                     $partData{$part.':awarded'} = ($awarded) ? $awarded : 0;      undef($navmap);
                     $totalAwarded += $awarded;      undef(@sequences);
   }
                     $tries = $cache->{"$name:$Version:$problem".  
                                       ":resource.$part.tries"};  sub BuildProblemStatisticsPage {
                     $partData{$part.':tries'} = ($tries) ? $tries : 0;      my ($r,$c)=@_;
                     $partData{$part.':wrong'} = $partData{$part.':tries'};      undef($navmap);
                     $totalTries += $tries;      undef(@sequences);
       #
                     my $val = $cache->{$name.":$Version:$problem".      my %Saveable_Parameters = ('Status' => 'scalar',
                                        ":resource.$part.solved"};                                 'statsoutputmode' => 'scalar',
                     if    ($val eq 'correct_by_student')   {$code = 'C';}                                  'Section' => 'array',
                     elsif ($val eq 'correct_by_override')  {$code = 'O';}                                 'Groups' => 'array',
                     elsif ($val eq 'incorrect_attempted')  {$code = 'I';}                                  'StudentData' => 'array',
                     elsif ($val eq 'incorrect_by_override'){$code = 'I';}                                 'Maps' => 'array',
                     elsif ($val eq 'excused')              {$code = 'x';}                                 'fieldselections'=> 'array');
                     $partData{$part.':code'}=$code;      &Apache::loncommon::store_course_settings('statistics',
                     if($partData{$part.':wrong'} ne 0 &&                                                 \%Saveable_Parameters);
                        ($code eq 'C' || $code eq 'O')) {      &Apache::loncommon::restore_course_settings('statistics',
                         $partData{$part.':wrong'}--;                                                  \%Saveable_Parameters);
                     }      #
                 }      &Apache::lonstatistics::PrepareClasslist();
       #
       # Clear the package variables
       undef(@StatsArray);
       undef(%SeqStat);
       #
       # Finally let the user know we are here
       $r->print(&Apache::lonhtmlcommon::breadcrumbs('Overall Problem Statistics',
                                                   'Statistics_Overall_Key'));
   
       my $interface = &CreateInterface($r);
       $r->print($interface);
       $r->print('<input type="hidden" name="sortby" value="'.$env{'form.sortby'}.
                 '" />');
       #
       my @CacheButtonHTML =
           &Apache::lonstathelpers::manage_caches($r,'Statistics','stats_status');
       my $Str;
       foreach my $html (@CacheButtonHTML) {
           $Str.=$html.('&nbsp;'x5);
       }
       #
       $r->print($Str);
       if (! exists($env{'form.firstrun'})) {
           $r->print('<p class="LC_info"><b>'.
                     &mt('Press "Generate Statistics" when you are ready.').
                     '</b></p>'.
                     '<p class="LC_info">'.
                     &mt('It may take some time to update the student data '.
                         'for the first analysis. Future analysis this session '.
                         'will not have this delay.').
                     '</p>');
    &clean_up();
           return;
       }
       $r->rflush();
       #
       # This probably does not need to be done each time we are called, but
       # it does not slow things down noticably.
       &Apache::loncoursedata::populate_weight_table();
       #
       ($navmap,@sequences) =
           &Apache::lonstatistics::selected_sequences_with_assessments();
       if (! ref($navmap)) {
           $r->print('<div class="LC_error">'.&mt('A course-wide error occurred.').'</div>'.
                     '<h3>'.$navmap.'</h3>');
    &clean_up();
           return;
       }
       if (exists($env{'form.Excel'})) {
           $r->print('<p>'.
                     &Apache::lonstatistics::section_and_enrollment_description().
                     '</p>');
           &Excel_output($r);
       } else {
           $r->print('<input type="submit" name="Excel" value="'.
                     &mt('Produce Excel Output').'" />'.'&nbsp;'x5);
           $r->rflush();
           $r->print('<p>'.
                     &Apache::lonstatistics::section_and_enrollment_description().
                     '</p>');
           my $count = 0;
           foreach my $seq (@sequences) {
               my @resources =
                   &Apache::lonstathelpers::get_resources($navmap,$seq);
               $count += scalar(@resources);
           }
           if ($count > 10) {
               $r->print('<p>'.
                         &mt('Compiling statistics for [quant,_1,problem]',$count).
                         '</p>');
               if ($count > 30) {
                   $r->print('<p class="LC_info">'.&mt('This will take some time.').'</p>');
               }
               $r->rflush();
           }
           #
           my $sortby = $env{'form.sortby'};
           $sortby = 'container' if (! defined($sortby) || $sortby =~ /^\s*$/);
           my $plot = $env{'form.plot'};
           if ($plot eq '' || $plot eq 'none') {
               undef($plot);
           }
           if ($sortby eq 'container' && ! defined($plot)) {
               &output_sequence_statistics($r);
               &output_html_by_sequence($r);
           } else {
               if (defined($plot)) {
                   &make_plot($r,$plot);
             }              }
               &output_html_stats($r);
               &output_sequence_statistics($r);
           }
       }
       &clean_up();
       return;
   }
   
             # Loop through all the parts for the current problem in the   sub output_sequence_statistics {
             # correct order and prepare the output      my ($r) = @_;
     my $partCounter=0;      my $c=$r->connection();
             foreach (split(/\:/,$cache->{$sequence.':'.$problemID.      $r->print('<h2>'.&mt('Sequence Statistics').
                                          ':parts'})) {        &Apache::loncommon::help_open_topic('Statistics_Sequence').
  $partCounter++;        '</h2>');
                 my $Yes = 0;      $r->print(&Apache::loncommon::start_data_table().&Apache::loncommon::start_data_table_header_row());
                 if($partData{$_.':code'} eq 'C' ||       $r->print(&sequence_html_header());
                    $partData{$_.':code'} eq 'O') {      $r->print(&Apache::loncommon::end_data_table_header_row());
                     $Yes=1;      foreach my $seq (@sequences) {
                 }          last if ($c->aborted);
                 my $pOrder=$partData{$_.':order'};          &compute_sequence_statistics($seq);
                 my $ptr = $tempSequenceOrder.':'.$pOrder.':'.$problemID;          $r->print(&sequence_html_output($seq));
       }
       $r->print(&Apache::loncommon::end_data_table());
       $r->rflush();
       return;
   }
   
                 if($partData{'count'} > 1) {  
                     $ptr .= "*(part $_)";  
                 }  
                 $discrim .= '&';  
   
  my ($pr_no,$dod)=split('&',$ptr);  ##########################################################
 # my $DoDiff=$DoDiff->{$dod};  ##########################################################
 #               $r->print('<br>'.$name.'---'.$ptr.'==='.$DoDiff);  ##
   ## HTML output routines
                 my $Fac = ($partData{$_.':tries'}) ?   ##
                     ($partData{$_.':awarded'}/$partData{$_.':tries'}) : 0;  ##########################################################
                 my $DisF;  ##########################################################
                 if($Fac > 0 &&  $Fac < 1) {   sub output_html_by_sequence {
                     $DisF = sprintf( "%.4f", $Fac );      my ($r) = @_;
                 } else {      my $c = $r->connection();
                     $DisF = $Fac;      $r->print('<br>'.&html_preamble());
                 }      #
       foreach my $seq (@sequences) {
           last if ($c->aborted);
           $r->print("<h3>".$seq->compTitle."</h3>".
                       &Apache::loncommon::start_data_table().
                       &Apache::loncommon::start_data_table_header_row().
                       &statistics_table_header('no container').
                       &Apache::loncommon::end_data_table_header_row()."\n");
           my @Data = &compute_statistics_on_sequence($seq);
           foreach my $data (@Data) {
               $r->print(&Apache::loncommon::start_data_table_row().
                       &statistics_html_table_data($data,'no container').
                       &Apache::loncommon::end_data_table_row()."\n");
           }
           $r->print(&Apache::loncommon::end_data_table()."\n");
           $r->rflush();
       }
       return;
   }
   
                 if ($Discuss->{"$name:$problem"}) {  sub output_html_stats {
     $TotDiscuss++;      my ($r)=@_;
                     $Discussed=1;      &compute_all_statistics($r);
                 }      $r->print(&html_preamble());
                 my $time = $cache->{"$name:$LatestVersion:$problem:timestamp"};      &sort_data($env{'form.sortby'});
                 $discrim .= $tempSequenceOrder.'@'.$pOrder.'='.$DisF.'+'.$Yes;      #
                 $ptr .= '&'.$partData{$_.':tries'}.      my $count=0;
                         '&'.$partData{$_.':wrong'}.      foreach my $data (@StatsArray) {
                         '&'.$partData{$_.':code'};          if ($count++ % 50 == 0) {
                 push (@$list, $ptr."&$Discussed");              $r->print(&Apache::loncommon::end_data_table());
 #                $r->print('<br>'.$_.$name.'---'.$ptr);              $r->print(&Apache::loncommon::start_data_table().
                       &Apache::loncommon::start_data_table_row().
 #### if ($DoDiff>0.85) {                      &statistics_table_header().
                       &Apache::loncommon::end_data_table_row());
           }
           $r->print(&Apache::loncommon::start_data_table_row().
                   &statistics_html_table_data($data).
                   &Apache::loncommon::end_data_table_row());
       }
       $r->print(&Apache::loncommon::end_data_table_row());
       return;
   }
   
                 $TimeTot += $time;  sub html_preamble {
       my $Str='';
       $Str .= "<h2>".
           $env{'course.'.$env{'request.course.id'}.'.description'}.
           "</h2>\n";
       my ($starttime,$endtime) = &Apache::lonstathelpers::get_time_limits();
       if (defined($starttime) || defined($endtime)) {
           # Inform the user what the time limits on the data are.
           $Str .= '<h3>'.&mt('Statistics on submissions from [_1] to [_2]',
                              &Apache::lonlocal::locallocaltime($starttime),
                              &Apache::lonlocal::locallocaltime($endtime)
                              ).'</h3>';
       }
       $Str .= "<p>".&mt('Compiled on [_1]',
                          &Apache::lonlocal::locallocaltime(time))."</p>";
       return $Str;
   }
   
                 if ($Yes==1 && $partData{$_.':tries'}==1) {  
     $TotFirst++;  
                 }  
 # my $Acts= $Activity->{$name.':'.$problem};  
 # if ($Acts) {  
 #    my $Pt=&ProcAct( $Acts, $time );  
     #my ($spe,$beg) = split(/\+/,$Pt);  
 #                    my $spe= $Pt;  
 #    if ($Yes==1) {$spent_yes += $spe;}  
 #    $spent += $spe;  
     #$Beg += $beg;  
 #                   $r->print('<br>'.$name.'---'.$problem.'---'.$spe);  
 # }  
  $TotalOpend++;  
  $ProbTot++;  
   
                 $tempProblemOrder++;  ###############################################
   ###############################################
   ##
   ## Misc HTML output routines
   ##
   ###############################################
   ###############################################
   sub statistics_html_table_data {
       my ($data,$options) = @_;
       my $row = '';
       foreach my $field (@Fields) {
           next if ($options =~ /no $field->{'name'}/);
           next if ($field->{'selected'} ne 'yes');
           $row .= '<td style="background-color:'.$field->{'color'}.'"';
           if (exists($field->{'align'})) {
               $row .= ' align="'.$field->{'align'}.'"';
             }              }
           $row .= '>';
           if (exists($field->{'special'}) && $field->{'special'} eq 'link') {
               $row .= '<a href="'.$data->{$field->{'name'}.'.link'}.'">';
         }          }
           if (exists($field->{'format'}) && $data->{$field->{'name'}} !~ /[A-Z]/i) {
               $row .= sprintf($field->{'format'},$data->{$field->{'name'}});
           } else {
               $row .= $data->{$field->{'name'}};
           }
           if (exists($field->{'special'}) && $field->{'special'} eq 'link') {
               $row.= '</a>';
           }
           $row .= '</td>';
     }      }
     my $pstr;      return $row;
     if($totalTries) {  
  my $DisFac = ($totalAwarded/$totalTries);  
  my $DisFactor = sprintf( "%.4f", $DisFac );  
         my $TS = sprintf( "%.2f", $spent );  
         my $TS_yes = sprintf( "%.2f", $spent_yes );  
  $pstr=$DisFactor.':'.$name.':'.$ProbTot.':'.$TotalOpend.':'.  
               $totalTries.':'.$ProbSolved.':'.$TotFirst.':'.  
               $TS_yes.':'.$TS.':'.$TotDiscuss;  
  (%$discriminant)->{$pstr}=$discrim;  
     }  
 }  }
   
 sub NumericSort {  sub statistics_table_header {
     $a <=> $b;      my ($options) = @_;
 }      my $header_row;
       foreach my $field (@Fields) {
 sub OrderedSort  {          next if ($options =~ /no $field->{'name'}/);
     if ($ENV{'form.order'} eq 'Descending') {          next if ($field->{'selected'} ne 'yes');
  $b <=> $a;          $header_row .= '<th>';
     } else {           if (exists($field->{'sortable'}) && $field->{'sortable'} eq 'yes') {
  $a <=> $b;              $header_row .= '<a href="javascript:'.
                   'document.Statistics.sortby.value='."'".$field->{'name'}."'".
                       ';document.Statistics.submit();">';
           }
           $header_row .= &mt($field->{'title'});
           if ($options =~ /sortable/) {
               $header_row.= '</a>';
           }
           if ($options !~ /no plots/        && 
               exists($field->{'graphable'}) && 
               $field->{'graphable'} eq 'yes') {
               $header_row.=' (';
               $header_row .= '<a href="javascript:'.
                   "document.Statistics.plot.value='$field->{'name'}'".
                       ';document.Statistics.submit();">';
               $header_row .= &mt('plot').'</a>)';
           }
           $header_row .= '</th>';
     }      }
       return $header_row;
 }  }
   
   sub sequence_html_header {
       my $Str .= '<tr>';
       foreach my $field (@SeqFields) {
   #        next if ($field->{'selected'} ne 'yes');
           $Str .= '<th bgcolor="'.$field->{'color'}.'"';
           $Str .= '>'.$field->{'title'}.'</th>';
       }
       $Str .= '</tr>';
       return $Str;
   }
   
   
 sub BuildStatisticsTable {  sub sequence_html_output {
     my ($cache,$upper,$lower,$list,$headings,$students,$r,$color)=@_;      my ($seq) = @_;
     my $NoElements = scalar @$list;      my $data = $SeqStat{$seq->symb};
     my @list=sort(@$list);      my $row = '<tr>';
       foreach my $field (@SeqFields) {
 #6666666          next if ($field->{'selected'} ne 'yes');
 #    my $file="/home/httpd/perl/tmp/183d.txt";          $row .= '<td bgcolor="'.$field->{'color'}.'"';
 #    open(OUT, ">$file");          if (exists($field->{'align'})) {
 #6666666              $row .= ' align="'.$field->{'align'}.'"';
 ##     &Apache::lonstatistics::Create_PrgWin($r);  
 ##777777  
 ##    my (%Activity) = &LoadActivityLog();  
 ##    $r->print('<script>popwin.document.popremain.remaining.value="'.  
 ##              'Loading Discussion...";</script>');  
 ##    my ($doDiffFile) = &LoadDoDiffFile();  
   
 ##777777  
 ##    $Str .= &Classify($discriminantFactor, $students);  
   
     my $p_count = 0;  
     my $dummy;  
     my $p_val;  
     my $ResId;  
     my %TempCache;  
     my $cIdx=0;  
   
     foreach my $sequence (split(':', $cache->{'orderedSequences'})) {  
         if($cache->{'ProblemStatisticsMaps'} ne 'All Maps'  &&  
            $cache->{'ProblemStatisticsMaps'} ne $cache->{$sequence.':title'}) {  
             next;  
         }          }
         &CreateProblemStatisticsTableHeading($cache->{'DisplayFormat'},           $row .= '>';
                                              $cache->{$sequence.':source'},          if (exists($field->{'format'})) {
                                              $cache->{$sequence.':title'},               $row .= sprintf($field->{'format'},$data->{$field->{'name'}});
                                              $headings,$r);          } else {
  my ($tar,$Tries,$Wrongs,$Code,$Disc)=split(/\&/,              $row .= $data->{$field->{'name'}};
                                                    $list[$cIdx]);          }
         my ($SqOrd,$PrOrd,$Prob)=split(/\:/,$tar);          $row .= '</td>';
  $sequence+=100;  
  while ($SqOrd==$sequence && $cIdx<$NoElements) {  
     my %storestats=();  
     my $pOrd=$PrOrd;  
     my $Temp = $Prob;  
     my $MxTries = 0;  
     my $TotalTries = 0;  
     my $YES = 0;  
     my $Incorrect = 0;  
     my $Override = 0;  
     my $StdNo = 0;  
     my $DiscNo=0;  
     my @StdLst;  
     while ($pOrd==$PrOrd && $cIdx<$NoElements)  
     {  
  $cIdx++;  
  $StdNo++;  
  $StdLst[ $StdNo ] = $Tries;  
  $TotalTries += $Tries;  
  if ( $MxTries < $Tries ) { $MxTries = $Tries; }   
  if ( $Code eq 'C' ){ $YES++; }  
  elsif( $Code eq 'I' ) { $Incorrect++; }  
  elsif( $Code eq 'O' ) { $Override++; }  
  elsif( $Code eq '-' ) { $StdNo--; }  
  ($tar,$Tries,$Wrongs,$Code,$Disc)=split(/\&/,  
                                                      $list[$cIdx]);  
         ($SqOrd,$PrOrd,$Prob)=split(/\:/,$tar);  
     }  
   
     $p_count++;  
     my $Dummy;  
     ($ResId,$Dummy)=split(/\*/,$Temp);  
     $Temp = '<a href="'.$cache->{$ResId.':source'}.  
                 '" target="_blank">'.$cache->{$ResId.':title'}.$Dummy.'</a>';  
   
 #    my $urlres = $cache->{$sequence.':source'});  
 #check with Gerd  
     #symb of the problem (already decluttered and cleaned)  
     my $urlres = $cache->{$ResId.':problem'};  
   
 #------------------------ Compute the Average of Tries about one problem  
     my $Average = ($StdNo) ? $TotalTries/$StdNo : 0;  
   
     $storestats{$ENV{'request.course.id'}.'___'.$urlres.'___timestamp'}=time;  
     $storestats{$ENV{'request.course.id'}.'___'.$urlres.'___stdno'}=$StdNo;  
     $storestats{$ENV{'request.course.id'}.'___'.$urlres.'___avetries'}=$Average;  
   
 #-------------------------------- Compute percentage of Wrong tries  
     my $Wrong = ( $StdNo ) ? 100 * ( $Incorrect / $StdNo ) : 0;  
   
 #-------------------------------- Compute Standard Deviation  
     my $StdDev = 0;   
     if ( $StdNo > 1 ) {  
  for ( my $n = 0; $n < $StdNo; $n++ ) {  
     my $Dif = $StdLst[ $n ]-$Average;  
     $StdDev += $Dif*$Dif;  
  }   
  $StdDev /= ( $StdNo - 1 );  
  $StdDev = sqrt( $StdDev );  
     }  
   
 #-------------------------------- Compute Degree of Difficulty  
     my $DoDiff = 0;  
     if( $TotalTries > 0 ) {  
  $DoDiff = 1 - ( ( $YES + $Override ) / $TotalTries );  
 #    $DoDiff =  ($TotalTries)/($YES + $Override+ 0.1);      
     }  
          
     $storestats{$ENV{'request.course.id'}.'___'.$urlres.'___difficulty'}=$DoDiff;  
   
 #-------------------------------- Compute the Skewness  
     my $Skewness = 0;  
     my $Sum = 0;   
     if ( $StdNo > 0 && $StdDev > 0 ) {  
  for ( my $n = 0; $n < $StdNo; $n++ ) {  
     my $Dif = $StdLst[ $n ]-$Average;  
     $Skewness += $Dif*$Dif*$Dif;  
  }   
  $Skewness /= $StdNo;  
  $Skewness /= $StdDev*$StdDev*$StdDev;  
     }  
   
 #--------------------- Compute the Discrimination Factors  
             my ($Up1,$Up2)=split(/\:/,$upper->{$sequence.'@'.$pOrd});  
     my ($Lw1,$Lw2)=split(/\:/,$lower->{$sequence.'@'.$pOrd});  
   
     my $Dis1 = $Up1 - $Lw1;  
     my $Dis2 = $Up2 - $Lw2;  
     my $_D1 = sprintf("%.2f", $Dis1);  
     my $_D2 = sprintf("%.2f", $Dis2);  
   
 #-----------------  Some restition in presenting the float numbers  
     my $Avg = sprintf( "%.2f", $Average );  
     my $Wrng = sprintf( "%.1f", $Wrong );  
     my $SD = sprintf( "%.1f", $StdDev );  
     my $DoD = sprintf( "%.2f", $DoDiff );  
     my $Sk = sprintf( "%.1f", $Skewness );  
     my $join = $sequence.'@'.$pOrd.'&'.$Temp.'&'.$StdNo.'&'.  
                        $TotalTries.'&'.$MxTries.'&'.$Avg.'&'.  
                        $YES.'&'.$Override.'&'.$Wrng.'&'.$DoD.'&'.  
        $SD.'&'.$Sk.'&'.$_D1.'&'.$_D2.'&'.  
                        $DiscNo.'&'.$Prob;  
   
     $TempCache{'CacheTable:'.($p_count-1)}=$join;  
   
 #6666666  
 #    $r->print('<br>'.$out.'&'.$DoD);  
 #            print (OUT $out.'@'.$DoD.'&');  
 #6666666  
   
 #check with Gerd  
     $urlres=~/^(\w+)\/(\w+)/;  
     if ($StdNo) {   
  &Apache::lonnet::put('resevaldata',\%storestats,$1,$2);   
     }  
 #-------------------------------- Row of statistical table  
             &TableRow($cache,$join,$cIdx,($p_count-1),$r,$color,  
                       \%TempCache);  
  }  
  $TempCache{'ProblemCount'}=$p_count;  
  &CloseTable($cache,$r);  
     }  
 ###    &Close_PrgWin();  
 #6666666  
 #    close( OUT );  
 #666666  
     return \%TempCache;  
 }  
   
   
 sub CacheStatisticsTable {  
     my ($state,$cache,$headings,$r,$color)=@_;  
     my @list = ();   
     my %TempCache;  
     my %myHeader = reverse( %$headings );  
     my $pos = $myHeader{$state};  
     if ($pos > 0) {$pos++;}  
     my $p_count = $cache->{'ProblemCount'};  
   
     for ( my $k=0; $k<$p_count;$k++) {  
  my $key=$cache->{'CacheTable:'.$k};  
  my @Temp=split(/\&/,$key);  
  $list[$k]=$Temp[$pos].'+'.$key;  
     }  
       
     if ($pos>0) {  
         @list = sort OrderedSort (@list);  
     } else {  
         @list = sort (@list);  
     }      }
     my $cIdx=0;      $row .= '</tr>'."\n";
       return $row;
     if ( $pos == 0 ) {  }
  foreach my $sequence (split(':', $cache->{'orderedSequences'})) {  
     if($cache->{'ProblemStatisticsMaps'} ne 'All Maps'  &&  
        $cache->{'ProblemStatisticsMaps'} ne $cache->{$sequence.':title'}) {  
  next;  
     }  
     if ($cIdx==$p_count) {  
  return \%TempCache;  
     }  
     &CreateProblemStatisticsTableHeading($cache->{'DisplayFormat'},   
                                              $cache->{$sequence.':source'},  
                                              $cache->{$sequence.':title'},   
                                              $headings,$r);  
   
     my ($tar)=split(/\&/,$list[$cIdx]);  
     $tar=~s/\+//eg;  
     my ($SqOrd)=split(/\@/,$tar);  
     $sequence+=100;  
     while ($SqOrd==$sequence && $cIdx<$p_count) {  
  my($Pre, $Post) = split(/\+/,$list[$cIdx]);   
  &TableRow($cache,$Post,$cIdx,$cIdx,$r,$color,\%TempCache);  
  $cIdx++;  
  my ($tar)=split(/\&/,$list[$cIdx]);  
  $tar=~s/\+//eg;  
  ($SqOrd)=split(/\@/,$tar);  
     }  
     &CloseTable($cache,$r);  
  }  
     }  
     else {  
         &CreateProblemStatisticsTableHeading($cache->{'DisplayFormat'},   
                                              'Sorted by: ',  
      $headings->{$pos-1},  
                                              $headings,$r);  
  for ( my $nIndex = 0; $nIndex < $p_count; $nIndex++ ) {  
     my($Pre, $Post) = split(/\+/,$list[$nIndex]);  
     &TableRow($cache,$Post,$nIndex,$nIndex,$r,$color,\%TempCache);  
  }   
  &CloseTable($cache,$r);  
     }  
   
     return \%TempCache;  
 }  
   
   
 sub TableRow {  
     my ($cache,$Str,$Idx,$RealIdx,$r,$color,$GraphDat)=@_;  
     my($PrOrd,$Temp,$StdNo,$TotalTries,$MxTries,$Avg,$YES,$Override,  
        $Wrng,$DoD,$SD,$Sk,$_D1,$_D2,$DiscNo,$Prob)=split(/\&/,$Str);  
     my $Ptr;  
     if($cache->{'DisplayFormat'} eq 'Display CSV Format') {  
         my ($ResId,$Dummy)=split(/\*/,$Prob);  
         $Ptr="\n".'<br>'.  
              "\n".'"'.($RealIdx+1).'",'.  
              "\n".'"'.$cache->{$ResId.':title'}.$Dummy.'",'.  
              "\n".'"'.$cache->{$ResId.':source'}.'",'.  
              "\n".'"'.$StdNo.'",'.  
              "\n".'"'.$TotalTries.'",'.  
              "\n".'"'.$MxTries.'",'.  
              "\n".'"'.$Avg.'",'.  
              "\n".'"'.$YES.'",'.  
              "\n".'"'.$Override.'",'.  
              "\n".'"'.$Wrng.'",'.  
              "\n".'"'.$DoD.'",'.  
              "\n".'"'.$SD.'",'.  
              "\n".'"'.$Sk.'",'.  
              "\n".'"'.$_D1.'",'.  
      "\n".'"'.$_D2.'"'.  
      "\n".'"'.$DiscNo.'"';  
   
         $r->print("\n".$Ptr);  ####################################################
   ####################################################
   ##
   ##    Plotting Routines
   ##
   ####################################################
   ####################################################
   sub make_plot {
       my ($r,$plot) = @_;
       &compute_all_statistics($r);
       &sort_data($env{'form.sortby'});
       if ($plot eq 'degrees') {
           &degrees_plot($r);
       } elsif ($plot eq 'tries statistics') {
           &tries_data_plot($r);
     } else {      } else {
         $Ptr="\n".'<tr>'.          &make_single_stat_plot($r,$plot);
              "\n".'<td>'.($RealIdx+1).'</td>'.      }
              "\n".'<td>'.$Temp.'</td>'.      return;
              "\n".'<td bgcolor='.$color->{"yellow"}.'> '.$StdNo.'</td>'.  
              "\n".'<td bgcolor='.$color->{"yellow"}.'>'.$TotalTries.'</td>'.  
              "\n".'<td bgcolor='.$color->{"yellow"}.'>'.$MxTries.'</td>'.  
              "\n".'<td bgcolor='.$color->{"gb"}.'>'.$Avg.'</td>'.  
              "\n".'<td bgcolor='.$color->{"gb"}.'> '.$YES.'</td>'.  
              "\n".'<td bgcolor='.$color->{"gb"}.'> '.$Override.'</td>'.  
              "\n".'<td bgcolor='.$color->{"red"}.'> '.$Wrng.'</td>'.  
              "\n".'<td bgcolor='.$color->{"red"}.'> '.$DoD.'</td>'.  
              "\n".'<td bgcolor='.$color->{"green"}.'> '.$SD.'</td>'.  
              "\n".'<td bgcolor='.$color->{"green"}.'> '.$Sk.'</td>'.  
              "\n".'<td bgcolor='.$color->{"purple"}.'> '.$_D1.'</td>'.  
      "\n".'<td bgcolor='.$color->{"purple"}.'> '.$_D2.'</td>'.  
              "\n".'<td bgcolor='.$color->{"yellow"}.'> '.$DiscNo.'</td>';  
         $r->print("\n".$Ptr.'</tr>' );  
     }  
     $GraphDat->{'GraphGif:'.$RealIdx}=$DoD.':'.$Wrng;  
 }  }
   
   sub make_single_stat_plot {
       my ($r,$datafield) = @_;
       #
       my $title; my $yaxis;
       foreach my $field (@Fields) {
           next if ($field->{'name'} ne $datafield);
           $title = &mt($field->{'long_title'});
           $yaxis = &mt($field->{'title'});
           last;
       }
       if ($title eq '' || $yaxis eq '') {
           # datafield is something we do not know enough about to plot
           $r->print('<p class="LC_warning">'.
                     &mt('Unable to plot the requested statistic.').
                     '</p>');
           return;
       }
       #
       # Build up the data sets to plot
       my @Labels;
       my @Data;
       my $max = 1;
       foreach my $data (@StatsArray) {
           push(@Labels,$data->{'problem_num'});
           push(@Data,$data->{$datafield});
           if ($data->{$datafield}>$max) {
               $max = $data->{$datafield};
           }
       }
       foreach (1,2,3,4,5,10,15,20,25,40,50,75,100,150,200,250,300,500,600,750,
                1000,1500,2000,2500,3000,3500,4000,5000,7500,10000,15000,20000) {
           if ($max <= $_) {
               $max = $_;
               last;
           }
       }
       if ($max > 20000) {
           $max = 10000*(int($max/10000)+1);
       }
       #
       $r->print("<p>".&Apache::loncommon::DrawBarGraph($title,
                                                        &mt('Problem Number'),
                                                        $yaxis,
                                                        $max,
                                                        undef, # colors
                                                        \@Labels,
                                                        \@Data)."</p>\n");
       return;
   }
   
 # For loading the colored table for display or un-colored for print  sub degrees_plot {
 sub setbgcolor {      my ($r)=@_;
     my $PrintTable=shift;      my $count = scalar(@StatsArray);
     my %color;      my $width = 50 + 10*$count;
     if ($PrintTable){      $width = 300 if ($width < 300);
  $color{"gb"}="#FFFFFF";      my $height = 300;
  $color{"red"}="#FFFFFF";      my $plot = '';
  $color{"yellow"}="#FFFFFF";      my $ymax = 0;
  $color{"green"}="#FFFFFF";      my $ymin = 0;
  $color{"purple"}="#FFFFFF";      my @Disc; my @Diff; my @Labels;
       foreach my $data (@StatsArray) {
           push(@Labels,$data->{'problem_num'});
           my $disc = $data->{'deg_of_disc'};
           my $diff = $data->{'deg_of_diff'};
           push(@Disc,$disc);
           push(@Diff,$diff);
           #
           $ymin = $disc if ($ymin > $disc);
           $ymin = $diff if ($ymin > $diff);
           $ymax = $disc if ($ymax < $disc);
           $ymax = $diff if ($ymax < $diff);
       }
       #
       # Make sure we show relevant information.
       if ($ymin < 0) {
           if (abs($ymin) < 0.05) {
               $ymin = 0;
           } else {
               $ymin = -1;
           }
       }
       if ($ymax > 0) {
           if (abs($ymax) < 0.05) {
               $ymax = 0;
           } else {
               $ymax = 1;
           }
       }
       #
       my $xmax = $Labels[-1];
       if ($xmax > 50) {
           if ($xmax % 10 != 0) {
               $xmax = 10 * (int($xmax/10)+1);
           }
     } else {      } else {
  $color{"gb"}="#DDFFFF";          if ($xmax % 5 != 0) {
  $color{"red"}="#FFDDDD";              $xmax = 5 * (int($xmax/5)+1);
  $color{"yellow"}="#EEFFCC";          }
  $color{"green"}="#DDFFDD";      }
  $color{"purple"}="#FFDDFF";      #
     }      my $discdata .= '<data>'.join(',',@Labels).'</data>'.$/.
                       '<data>'.join(',',@Disc).'</data>'.$/;
     return \%color;      #
 }      my $diffdata .= '<data>'.join(',',@Labels).'</data>'.$/.
                       '<data>'.join(',',@Diff).'</data>'.$/;
 sub ProblemStatisticsButtons {      #
     my ($displayFormat)=@_;      my $title = &mt('Degree of Discrimination[_1]and Degree of Difficulty','\n');
       if ($xmax > 50) {
     my $Ptr = '<tr><td></td><td align="left">';          $title = &mt('Degree of Discrimination and Degree of Difficulty');
     $Ptr .= '<input type="submit" name="DoDiffGraph" ';      }
     $Ptr .= 'value="DoDiff Graph" />'."\n";      my %lt = &Apache::lonlocal::texthash(
     $Ptr .= '&nbsp;&nbsp;&nbsp;';          'alttag' => 'Degree of Discrimination and Degree of Difficulty Plot',
     $Ptr .= '<input type="submit" name="PercentWrongGraph" ';          'xlabel' => 'Problem Number',
     $Ptr .= 'value="%Wrong Graph" />'."\n";      );
     $Ptr .= '&nbsp;&nbsp;&nbsp;';      #
     $Ptr .= '<input type="submit" name="DisplayCSVFormat" ';      $plot=<<"END";
     if($displayFormat eq 'Display CSV Format') {  <gnuplot 
         $Ptr .= 'value="Display Table Format" />'."\n";      texfont="10"
       fgcolor="x000000"
       plottype="Cartesian"
       font="large"
       grid="on"
       align="center"
       border="on"
       transparent="on"
       alttag="$lt{'alttag'}"
       samples="100"
       bgcolor="xffffff"
       height="$height"
       width="$width">
       <key 
           pos="top right"
           title=""
           box="off" />
       <title>$title</title>
       <axis xmin="0" ymin="$ymin" xmax="$xmax" ymax="$ymax" color="x000000" />
       <xlabel>$lt{'xlabel'}</xlabel>
       <curve 
           linestyle="linespoints" 
           name="DoDisc" 
           pointtype="0" 
           color="x000000">
           $discdata
       </curve>
       <curve 
           linestyle="linespoints" 
           name="DoDiff" 
           pointtype="0" 
           color="xFF0000">
           $diffdata
       </curve>
   </gnuplot>
   END
       my $plotresult =
           '<p>'.&Apache::lonxml::xmlparse($r,'web',$plot).'</p>'.$/;
       $r->print($plotresult);
       return;
   }
   
   sub tries_data_plot {
       my ($r)=@_;
       my $count = scalar(@StatsArray);
       my $width = 50 + 10*$count;
       $width = 300 if ($width < 300);
       my $height = 300;
       my $plot = '';
       my @STD;  my @Mean; my @Max; my @Min;
       my @Labels;
       my $ymax = 5;
       foreach my $data (@StatsArray) {
           my $max = $data->{'mean_tries'} + $data->{'std_tries'};
           $ymax = $max if ($ymax < $max);
           $ymax = $max if ($ymax < $max);
           push(@Labels,$data->{'problem_num'});
           push(@STD,$data->{'std_tries'});
           push(@Mean,$data->{'mean_tries'});
       }
       #
       # Make sure we show relevant information.
       my $xmax = $Labels[-1];
       if ($xmax > 50) {
           if ($xmax % 10 != 0) {
               $xmax = 10 * (int($xmax/10)+1);
           }
     } else {      } else {
         $Ptr .= 'value="Display CSV Format" />'."\n";          if ($xmax % 5 != 0) {
               $xmax = 5 * (int($xmax/5)+1);
           }
     }      }
     $Ptr .= '</td></tr>';      $ymax = int($ymax)+1+2;
       #
       my $std_data .= '<data>'.join(',',@Labels).'</data>'.$/.
                       '<data>'.join(',',@Mean).'</data>'.$/;
       #
       my $std_error_data .= '<data>'.join(',',@Labels).'</data>'.$/.
                             '<data>'.join(',',@Mean).'</data>'.$/.
                             '<data>'.join(',',@STD).'</data>'.$/;
       #
       my $title = &mt('Mean and S.D. of Tries');
       if ($xmax > 30) {
           $title = &mt('Mean and Standard Deviation of Tries');
       }
       #
       my %lt = &Apache::lonlocal::texthash(
           'alttag' => 'Mean and S.D of Tries Plot',
           'xlabel' => 'Problem Number',
           'ylabel' => 'Number of Tries',
       );
       $plot=<<"END";
   <gnuplot 
       texfont="10"
       fgcolor="x000000"
       plottype="Cartesian"
       font="large"
       grid="on"
       align="center"
       border="on"
       transparent="on"
       alttag="$lt{'alttag'}"
       samples="100"
       bgcolor="xffffff"
       height="$height"
       width="$width">
       <title>$title</title>
       <axis xmin="0" ymin="0" xmax="$xmax" ymax="$ymax" color="x000000" />
       <xlabel>$lt{'xlabel'}</xlabel>
       <ylabel>$lt{'ylabel'}</ylabel>
       <curve 
           linestyle="yerrorbars"
           name="S.D. Tries" 
           pointtype="1" 
           color="x666666">
           $std_error_data
       </curve>
       <curve 
           linestyle="points"
           name="Mean Tries" 
           pointtype="1" 
           color="xCC4444">
           $std_data
       </curve>
   </gnuplot>
   END
       my $plotresult =
           '<p>'.&Apache::lonxml::xmlparse($r,'web',$plot).'</p>'.$/;
       $r->print($plotresult);
       return;
   }
   
     return $Ptr;  sub plot_dropdown {
       my $current = '';
       my $title;
       #
       if (defined($env{'form.plot'})) {
           $current = $env{'form.plot'};
       }
       #
       my @Additional_Plots = (
                               { graphable=>'yes',
                                 name => 'degrees',
                                 title => 'Difficulty Indexes' },
                               { graphable=>'yes',
                                 name => 'tries statistics',
                                 title => 'Tries Statistics' });
       #
       my $Str= "\n".'<select name="plot" size="1">';
       $Str .= '<option name="none"></option>'."\n";
       $Str .= '<option name="none2">'.&mt('none').'</option>'."\n";
       foreach my $field (@Additional_Plots,@Fields) {
           if (! exists($field->{'graphable'}) ||
               $field->{'graphable'} ne 'yes') {
               next;
           }
           $Str .= '<option value="'.$field->{'name'}.'"';
           if ($field->{'name'} eq $current) {
               $Str .= ' selected="selected"';
           }
           $title = &mt($field->{'long_title'});
           $title = &mt($field->{'title'}) if (!$title);
           $Str.= '>'.$title.'</option>'."\n";
       }
       $Str .= '</select>'."\n";
       return $Str;
 }  }
   
 sub ProblemStatisticsLegend {  ###############################################
     my $Ptr = '';  ###############################################
     $Ptr = '<table border="0">';  ##
     $Ptr .= '<tr><td>';  ## Excel output routines
     $Ptr .= '<b>#Stdnts</b></td>';  ##
     $Ptr .= '<td>Total Number of Students opened the problem.';  ###############################################
     $Ptr .= '</td></tr><tr><td>';  ###############################################
     $Ptr .= '<b>Tries</b></td>';  sub Excel_output {
     $Ptr .= '<td>Total Number of Tries for solving the problem.';      my ($r) = @_;
     $Ptr .= '</td></tr><tr><td>';      $r->print('<h2>'.&mt('Preparing Excel Spreadsheet').'</h2>');
     $Ptr .= '<b>Mod</b></td>';      ##
     $Ptr .= '<td>Maximunm Number of Tries for solving the problem.';      ## Compute the statistics
     $Ptr .= '</td></tr><tr><td>';      &compute_all_statistics($r);
     $Ptr .= '<b>Mean</b></td>';      my $c = $r->connection;
     $Ptr .= '<td>Average Number of the tries. [ Tries / #Stdnts ]';      return if ($c->aborted());
     $Ptr .= '</td></tr><tr><td>';      #
     $Ptr .= '<b>#YES</b></td>';      my ($starttime,$endtime) = &Apache::lonstathelpers::get_time_limits();
     $Ptr .= '<td>Number of students solved the problem correctly.';      ##
     $Ptr .= '</td></tr><tr><td>';      ## Create the excel workbook
     $Ptr .= '<b>#yes</b></td>';      my ($excel_workbook,$filename,$format) =
     $Ptr .= '<td>Number of students solved the problem by override.';          &Apache::loncommon::create_workbook($r);
     $Ptr .= '</td></tr><tr><td>';      return if (! defined($excel_workbook));
     $Ptr .= '<b>%Wrng</b></td>';      #
     $Ptr .= '<td>Percentage of students tried to solve the problem ';      # Add a worksheet
     $Ptr .= 'but still incorrect. [ 100*((#Stdnts-(#YES+#yes))/#Stdnts) ]';      my $sheetname = $env{'course.'.$env{'request.course.id'}.'.description'};
     $Ptr .= '</td></tr><tr><td>';      if (length($sheetname) > 31) {
     $Ptr .= '<b>DoDiff</b></td>';          $sheetname = substr($sheetname,0,31);
     $Ptr .= '<td>Degree of Difficulty of the problem.  ';      }
     $Ptr .= '[ 1 - ((#YES+#yes) / Tries) ]';      my $excel_sheet = $excel_workbook->addworksheet(
     $Ptr .= '</td></tr><tr><td>';          &Apache::loncommon::clean_excel_name($sheetname));
     $Ptr .= '<b>S.D.</b></td>';      ##
     $Ptr .= '<td>Standard Deviation of the tries.  ';      ## Begin creating excel sheet
     $Ptr .= '[ sqrt(sum((Xi - Mean)^2)) / (#Stdnts-1) ';      ##
     $Ptr .= 'where Xi denotes every student\'s tries ]';      my ($rows_output,$cols_output) = (0,0);
     $Ptr .= '</td></tr><tr><td>';      #
     $Ptr .= '<b>Skew.</b></td>';      # Put the course description in the header
     $Ptr .= '<td>Skewness of the students tries.';      $excel_sheet->write($rows_output,$cols_output++,
     $Ptr .= '[(sqrt( sum((Xi - Mean)^3) / #Stdnts)) / (S.D.^3)]';                     $env{'course.'.$env{'request.course.id'}.'.description'},
     $Ptr .= '</td></tr><tr><td>';                          $format->{'h1'});
     $Ptr .= '<b>Dis.F.</b></td>';      $cols_output += 3;
     $Ptr .= '<td>Discrimination Factor: A Standard for evaluating the ';      #
     $Ptr .= 'problem according to a Criterion<br>';      # Put a description of the sections listed
     $Ptr .= '<b>[Applied Criterion in %27 Upper Students - ';      my $sectionstring = '';
     $Ptr .= 'Applied the same Criterion in %27 Lower Students]</b><br>';      $excel_sheet->write($rows_output,$cols_output++,
     $Ptr .= '<b>1st Criterion</b> for Sorting the Students: ';                          &Apache::lonstatistics::section_and_enrollment_description('plaintext'),
     $Ptr .= '<b>Sum of Partial Credit Awarded / Total Number of Tries</b><br>';                          $format->{'h3'});
     $Ptr .= '<b>2nd Criterion</b> for Sorting the Students: ';      $cols_output += scalar(&Apache::lonstatistics::get_selected_sections());
     $Ptr .= '<b>Total number of Correct Answers / Total Number of Tries</b>';      $cols_output += scalar(&Apache::lonstatistics::get_selected_groups());
     $Ptr .= '</td></tr>';      #
     $Ptr .= '<tr><td><b>Disc.</b></td>';      # Time restrictions
     $Ptr .= '<td>Number of Students had at least one discussion.';      my $time_string;
     $Ptr .= '</td></tr></table>';      if (defined($starttime)) {
           # call localtime but not lonlocal:locallocaltime because excel probably
     return $Ptr;          # cannot handle localized text.  Probably.
 }          $time_string .= 'Data collected from '.localtime($time_string);
           if (defined($endtime)) {
 #------- Processing upperlist and lowerlist according to each problem              $time_string .= ' to '.localtime($endtime);
 sub ProcessDiscriminant {          }
     my ($List,$r) = @_;          $time_string .= '.';
     my @sortedList = sort (@$List);      } elsif (defined($endtime)) {
     my $Count = scalar @sortedList;          # See note above about lonlocal:locallocaltime
     my $Problem;          $time_string .= 'Data collected before '.localtime($endtime).'.';
     my @Dis;      }
     my $Slvd=0;      if (defined($time_string)) {
     my $tmp;          $excel_sheet->write($rows_output,$cols_output++,$time_string);
     my $Sum1=0;          $cols_output+= 5;
     my $Sum2=0;      }
     my $nIndex=0;      #
     my $nStudent=0;      # Put the date in there too
     my %Proc=undef;      $excel_sheet->write($rows_output,$cols_output++,
     while ($nIndex<$Count) {                          'Compiled on '.localtime(time));
 #        $r->print("<br> $nIndex) $sortedList[$nIndex]");      #
  ($Problem,$tmp)=split(/\=/,$sortedList[$nIndex]);      $rows_output++;
  @Dis=split(/\+/,$tmp);      $cols_output=0;
  my $Temp = $Problem;      ##
  do {      ## Sequence Statistics
     $nIndex++;      ##
     $nStudent++;      &write_headers($excel_sheet,$format,\$rows_output,\$cols_output,
     $Sum1 += $Dis[0];                     \@SeqFields);
     $Sum2 += $Dis[1];      foreach my $seq (@sequences) {
     ($Problem,$tmp)=split(/\=/,$sortedList[$nIndex]);          my $data = $SeqStat{$seq->symb};
     @Dis=split(/\+/,$tmp);          $cols_output=0;
  } while ( $Problem eq $Temp && $nIndex < $Count );          foreach my $field (@SeqFields) {
  $Proc{$Temp}=($Sum1/$nStudent).':'.($Sum2/$nStudent);              next if ($field->{'selected'} ne 'yes');
 #        $r->print("<br> $nIndex) $Temp --> ($nStudent) $Proc{$Temp}");              my $fieldformat = undef;
  $Sum1=0;              if (exists($field->{'excel_format'})) {
  $Sum2=0;                  $fieldformat = $format->{$field->{'excel_format'}};
  $nStudent=0;              }
     }              $excel_sheet->write($rows_output,$cols_output++,
                                   $data->{$field->{'name'}},$fieldformat);
     return %Proc;          }
 }          $rows_output++;
           $cols_output=0;
 #------- Creating Discimination factor         }
 sub Discriminant {      ##
     my ($discriminant,$r)=@_;      ## Resource Statistics
     my @discriminantKeys=keys(%$discriminant);      ##
     my $Count = scalar @discriminantKeys;      $rows_output++;
       $cols_output=0;
     my $UpCnt = int(0.27*$Count);      &write_headers($excel_sheet,$format,\$rows_output,\$cols_output,
     my $low=0;                     \@Fields);
     my $up=$Count-$UpCnt;      #
     my @UpList=();      foreach my $data (@StatsArray) {
     my @LowList=();          $cols_output=0;
           foreach my $field (@Fields) {
     $Count=0;              next if ($field->{'selected'} ne 'yes');
     foreach my $key (sort(@discriminantKeys)) {               next if ($field->{'name'} eq 'problem_num');
  $Count++;                  my $fieldformat = undef;
  if($low < $UpCnt || $Count > $up) {              if (exists($field->{'excel_format'})) {
             $low++;                  $fieldformat = $format->{$field->{'excel_format'}};
             my $str=$discriminant->{$key};  
             foreach(split(/\&/,$str)){  
                 if($_) {  
                     if($low<$UpCnt) { push(@LowList,$_); }  
                     else            { push(@UpList,$_);  }  
                 }  
             }              }
               $excel_sheet->write($rows_output,$cols_output++,
                                   $data->{$field->{'name'}},$fieldformat);
           }
           $rows_output++;
           $cols_output=0;
       }
       #
       $excel_workbook->close();
       #
       # Tell the user where to get their excel file
       $r->print('<br />'.
                 '<a href="'.$filename.'">'.
                 &mt('Your Excel Spreadsheet').'</a>'."\n");
       $r->rflush();
       return;
   }
   
   ##
   ## &write_headers
   ##
   sub write_headers {
       my ($excel_sheet,$format,$rows_output,$cols_output,$Fields) = @_;
       ##
       ## First the long titles
       foreach my $field (@{$Fields}) {
           next if ($field->{'name'} eq 'problem_num');
           next if ($field->{'selected'} ne 'yes');
           if (exists($field->{'long_title'})) {
               $excel_sheet->write($$rows_output,${$cols_output},
                                   $field->{'long_title'},
                                   $format->{'bold'});
           } else {
               $excel_sheet->write($$rows_output,${$cols_output},'');
           }
           ${$cols_output}+= 1;
       }
       ${$cols_output} =0;
       ${$rows_output}+=1;
       ##
       ## Then the short titles
       foreach my $field (@{$Fields}) {
           next if ($field->{'selected'} ne 'yes');
           next if ($field->{'name'} eq 'problem_num');
           # Use english for excel as I am not sure how well excel handles
           # other character sets....
           $excel_sheet->write($$rows_output,$$cols_output,
                               $field->{'title'},
                               $format->{'bold'});
           $$cols_output+=1;
       }
       ${$cols_output} =0;
       ${$rows_output}+=1;
       return;
   }
   
   ##################################################
   ##################################################
   ##
   ## Statistics Gathering and Manipulation Routines
   ##
   ##################################################
   ##################################################
   sub compute_statistics_on_sequence {
       my ($seq) = @_;
       my @Data;
       foreach my $res (&Apache::lonstathelpers::get_resources($navmap,$seq)) {
           foreach my $part (@{$res->parts}) {
               next if (($res->is_survey($part)) || ($res->is_anonsurvey($part))) ;
               #
               # This is where all the work happens
               my $data = &get_statistics($seq,$res,$part,scalar(@StatsArray)+1);
               push (@Data,$data);
               push (@StatsArray,$data);
         }          }
     }      }
     my %DisUp =  &ProcessDiscriminant(\@UpList,$r);      return @Data;
     my %DisLow = &ProcessDiscriminant(\@LowList,$r);  }
   
     return (\%DisUp, \%DisLow);  sub compute_all_statistics {
 }         my ($r) = @_;
       if (@StatsArray > 0) {
           # Assume we have already computed the statistics
           return;
       }
       my $c = $r->connection;
       foreach my $seq (@sequences) {
           last if ($c->aborted);
           &compute_sequence_statistics($seq);
           &compute_statistics_on_sequence($seq);
       }
   }
   
 #---- END Problem Statistics Web Page ----------------------------------------  sub sort_data {
       my ($sortkey) = @_;
       return if (! @StatsArray);
       #
       # Sort the data
       my $sortby = undef;
       foreach my $field (@Fields) {
           if ($sortkey eq $field->{'name'}) {
               $sortby = $field->{'name'};
           }
       }
       if (! defined($sortby) || $sortby eq '' || $sortby eq 'problem_num') {
           $sortby = 'container';
       }
       if ($sortby ne 'container') {
           # $sortby is already defined, so we can charge ahead
           if ($sortby =~ /^(title|part)$/i) {
               # Alpha comparison
               @StatsArray = sort {
                   lc($a->{$sortby}) cmp lc($b->{$sortby}) ||
                   lc($a->{'title'}) cmp lc($b->{'title'}) ||
                   lc($a->{'part'}) cmp lc($b->{'part'});
               } @StatsArray;
           } else {
               # Numerical comparison
               @StatsArray = sort {
                   my $retvalue = 0;
                   if ($b->{$sortby} eq 'nan') {
                       if ($a->{$sortby} ne 'nan') {
                           $retvalue = -1;
                       } else {
                           $retvalue = 0;
                       }
                   }
                   if ($a->{$sortby} eq 'nan') {
                       if ($b->{$sortby} ne 'nan') {
                           $retvalue = 1;
                       }
                   }
                   if ($retvalue eq '0') {
                       $retvalue = $b->{$sortby} <=> $a->{$sortby}     ||
                               lc($a->{'title'}) <=> lc($b->{'title'}) ||
                               lc($a->{'part'})  <=> lc($b->{'part'});
                   }
                   $retvalue;
               } @StatsArray;
           }
       }
       #
       # Renumber the data set
       my $count;
       foreach my $data (@StatsArray) {
           $data->{'problem_num'} = ++$count;
       }
       return;
   }
   
 #---- Problem Statistics Graph Web Page --------------------------------------  ########################################################
   ########################################################
   
 # ------------------------------------------- Prepare data for Graphical chart  =pod
   
 sub BuildGraphicChart {  =item &get_statistics()
     my ($ylab,$r,$cacheDB)=@_;  
     my %cache;  
     my $Col;  
     my $data='';  
     my $count = 0;  
     my $Max = 0;  
   
     unless(tie(%cache,'GDBM_File',$cacheDB,&GDBM_READER,0640)) {  Wrapper routine from the call to loncoursedata::get_problem_statistics.
         return '<html><body>Unable to tie database.</body></html>';  Calls lonstathelpers::get_time_limits() to limit the data set by time
     }  and &compute_discrimination_factor
      
     my $p_count = $cache{'ProblemCount'};  Inputs: $sequence, $resource, $part, $problem_num
   
   Returns: Hash reference with statistics data from
   loncoursedata::get_problem_statistics.
   
   =cut
   
   ########################################################
   ########################################################
   sub get_statistics {
       my ($sequence,$resource,$part,$problem_num) = @_;
       #
       my ($starttime,$endtime) = &Apache::lonstathelpers::get_time_limits();
       my $symb = $resource->symb;
       my $courseid = $env{'request.course.id'};
       #
       my $data = &Apache::loncoursedata::get_problem_statistics
                           ([&Apache::lonstatistics::get_selected_sections()],
                            [&Apache::lonstatistics::get_selected_groups()],
                            $Apache::lonstatistics::enrollment_status,
                            $symb,$part,$courseid,$starttime,$endtime);
       $data->{'symb'}        = $symb;
       $data->{'part'}        = $part;
       $data->{'problem_num'} = $problem_num;
       $data->{'container'}   = $sequence->compTitle;
       $data->{'title'}       = $resource->compTitle;
       $data->{'title.link'}  = $resource->src.'?symb='.
           &escape($resource->symb);
       #
       if ($SelectedFields{'deg_of_disc'}) {
           $data->{'deg_of_disc'} =
               &compute_discrimination_factor($resource,$part,$sequence);
       }
       #
       # Store in metadata if computations were done for all students
       if ($data->{'num_students'} > 1) {
           my @Sections = &Apache::lonstatistics::get_selected_sections();
           my $sections = '"'.join(' ',@Sections).'"';
           $sections =~ s/&+/_/g;  # Ensure no special characters
           $data->{'sections'}=$sections;
           $data->{'course'} = $env{'request.course.id'};
           my $urlres=(&Apache::lonnet::decode_symb($resource->symb))[2];
           my %storestats =
               &LONCAPA::lonmetadata::dynamic_metadata_storage($data);
           my ($dom,$user) = ($urlres=~m{^($LONCAPA::domain_re)/($LONCAPA::username_re)});
           &Apache::lonnet::put('nohist_resevaldata',\%storestats,$dom,$user);
       }
       #
       $data->{'tries_per_correct'} = $data->{'tries'} / 
           ($data->{'num_solved'}+0.1);
       #
       # Get the due date for research purposes (commented out most of the time)
   #    my $duedate = &Apache::lonnet::EXT('resource.'.$part.'.duedate',$symb);;
   #    my $opendate = &Apache::lonnet::EXT('resource.'.$part.'.opendate',$symb);
   #    my $maxtries = &Apache::lonnet::EXT('resource.'.$part.'.maxtries',$symb);
   #    my $hinttries = &Apache::lonnet::EXT('resource.'.$part.'.hinttries',$symb);
       my $weight = &Apache::lonnet::EXT('resource.'.$part.'.weight',$symb);
       $data->{'weight'} = $weight;
   #    $data->{'duedate'} = $duedate;
   #    $data->{'opendate'} = $opendate;
   #    $data->{'maxtries'} = $maxtries;
   #    $data->{'hinttries'} = $hinttries;
   #    $data->{'resptypes'} = join(',',@{$resource->{'partdata'}->{$part}->{'ResponseTypes'}});
       return $data;
   }
   
     for ( my $k=0; $k<$p_count;$k++) {  ###############################################
         my @Temp=split(/\:/,$cache{'GraphGif:'.$k});  ###############################################
         my $inf = $Temp[$Col];   
  if ( $Max < $inf ) {$Max = $inf;}  =pod
       $data .= $inf.',';  
         $count++;  
     }  
     untie(%cache);  
 #    $r->print("<br>count=$p_count >>data= $data");  
   
     if ( $Max > 1 ) {   =item &compute_discrimination_factor()
  $Max += (10 - $Max % 10);  
       $Max = int($Max);  
     } else { $Max = 1; }  
   
     my $cid=$ENV{'request.course.id'};  Inputs: $Resource, $Sequence
   
     if ( $ylab eq 'DoDiff Graph' ) {  Returns: integer between -1 and 1
  $ylab = 'Degree-of-Difficulty';  
  $Col = 0;  =cut
   
   ###############################################
   ###############################################
   sub compute_discrimination_factor {
       my ($resource,$part,$seq) = @_;
       my $symb = $resource->symb;
       my @Resources;
       foreach my $res (&Apache::lonstathelpers::get_resources($navmap,$seq)){
           next if ($res->symb eq $symb);
           push (@Resources,$res->symb);
       }
       #
       # rank
       my ($starttime,$endtime) = &Apache::lonstathelpers::get_time_limits();
       my $ranking =
           &Apache::loncoursedata::rank_students_by_scores_on_resources
           (\@Resources,
            [&Apache::lonstatistics::get_selected_sections()],
            [&Apache::lonstatistics::get_selected_groups()],
            $Apache::lonstatistics::enrollment_status,undef,
            $starttime,$endtime, $symb);
       #
       # compute their percent scores on the problems in the sequence,
       my $number_to_grab = int(scalar(@{$ranking})/4);
       my $num_students = scalar(@{$ranking});
       my @BottomSet = map { $_->[&Apache::loncoursedata::RNK_student()];
                         } @{$ranking}[0..$number_to_grab];
       my @TopSet    =
           map {
               $_->[&Apache::loncoursedata::RNK_student()];
             } @{$ranking}[-$number_to_grab..0];
       if (! @BottomSet || (@BottomSet == 1 && $BottomSet[0] eq '') ||
           ! @TopSet    || (@TopSet    == 1 && $TopSet[0]    eq '')) {
           return 'nan';
       }
       my ($bottom_sum,$bottom_max) =
           &Apache::loncoursedata::get_sum_of_scores($symb,$part,\@BottomSet,
                                                     undef,$starttime,$endtime);
       my ($top_sum,$top_max) =
           &Apache::loncoursedata::get_sum_of_scores($symb,$part,\@TopSet,
                                                     undef,$starttime,$endtime);
       my $deg_of_disc;
       if ($top_max == 0 || $bottom_max==0) {
           $deg_of_disc = 'nan';
     } else {      } else {
  $ylab = 'Wrong-Percentage';          $deg_of_disc = ($top_sum/$top_max) - ($bottom_sum/$bottom_max);
  $Col = 1;  
     }      }
     my $Course = $ENV{'course.'.$cid.'.description'};      #&Apache::lonnet::logthis('    '.$top_sum.'/'.$top_max.
     $Course =~ s/\ /"_"/eg;      #                         ' - '.$bottom_sum.'/'.$bottom_max);
     my $GData=$Course.'&'.'Problems#'.'&'.$ylab.'&'.      return $deg_of_disc;
       $Max.'&'.$count.'&'.$data;  }
   
     $r->print('<IMG src="/cgi-bin/graph.gif?'.$GData.'" />');  
   
   ###############################################
   ###############################################
   ##
   ## Compute KR-21
   ##
   ## To compute KR-21, you need the following information:
   ##
   ## K=the number of items in your test
   ## M=the mean score on the test
   ## s=the standard deviation of the scores on your test
   ##
   ## then:
   ##
   ## KR-21 rk= [K/(K-1)] * [1- (M*(K-M))/(K*s^2))]
   ##
   ###############################################
   ###############################################
   sub compute_sequence_statistics {
       my ($seq) = @_;
       my $symb = $seq->symb;
       my @Resources;
       my $part_count;
       foreach my $res (&Apache::lonstathelpers::get_resources($navmap,$seq)) {
           push (@Resources,$res->symb);
           $part_count += scalar(@{$res->parts});
       }
       my ($starttime,$endtime) = &Apache::lonstathelpers::get_time_limits();
       #
       # First compute statistics based on student scores
       my ($smin,$smax,$sMean,$sSTD,$scount,$sMAX) =
           &Apache::loncoursedata::score_stats
                       ([&Apache::lonstatistics::get_selected_sections()],
                        [&Apache::lonstatistics::get_selected_groups()],
                        $Apache::lonstatistics::enrollment_status,
                        \@Resources,$starttime,$endtime,undef);
       $SeqStat{$symb}->{'title'}  = $seq->compTitle;
       $SeqStat{$symb}->{'scoremax'}  = $smax;
       $SeqStat{$symb}->{'scoremin'}  = $smin;
       $SeqStat{$symb}->{'scoremean'} = $sMean;
       $SeqStat{$symb}->{'scorestd'}  = $sSTD;
       $SeqStat{$symb}->{'scorecount'} = $scount;
       $SeqStat{$symb}->{'max_possible'} = $sMAX;
       #
       # Compute statistics based on the number of correct problems
       # 'correct' is taken to mean
       my ($cmin,$cmax,$cMean,$cSTD,$ccount)=
           &Apache::loncoursedata::count_stats
           ([&Apache::lonstatistics::get_selected_sections()],
            [&Apache::lonstatistics::get_selected_groups()],
            $Apache::lonstatistics::enrollment_status,
            \@Resources,$starttime,$endtime,undef);
       my $K = $part_count;
       my $kr_21;
       if ($K > 1 && $cSTD > 0) {
           $kr_21 =  ($K/($K-1)) * (1 - $cMean*($K-$cMean)/($K*$cSTD**2));
       } else {
           $kr_21 = 'nan';
       }
       $SeqStat{$symb}->{'countmax'} = $cmax;
       $SeqStat{$symb}->{'countmin'} = $cmin;
       $SeqStat{$symb}->{'countstd'} = $cSTD;
       $SeqStat{$symb}->{'countmean'} = $cMean;
       $SeqStat{$symb}->{'count'} = $ccount;
       $SeqStat{$symb}->{'items'} = $K;
       $SeqStat{$symb}->{'KR-21'}=$kr_21;
     return;      return;
 }  }
   
   
   
   =pod 
   
   =item ProblemStatisticsLegend
   
   =over 4
   
   =item #Stdnts
   Total number of students attempted the problem.
   
   =item Tries
   Total number of tries for solving the problem.
   
   =item Max Tries
   Largest number of tries for solving the problem by a student.
   
   =item Mean
   Average number of tries. [ Tries / #Stdnts ]
   
   =item #YES
   Number of students solved the problem correctly.
   
   =item #yes
   Number of students solved the problem by override.
   
   =item %Wrong
   Percentage of students who tried to solve the problem 
   but is still incorrect. [ 100*((#Stdnts-(#YES+#yes))/#Stdnts) ]
   
   =item DoDiff
   Degree of Difficulty of the problem.  
   [ 1 - ((#YES+#yes) / Tries) ]
   
   =item S.D.
   Standard Deviation of the tries.  
   [ sqrt(sum((Xi - Mean)^2)) / (#Stdnts-1) 
   where Xi denotes every student\'s tries ]
   
   =item Skew.
   Skewness of the students tries.
   [(sqrt( sum((Xi - Mean)^3) / #Stdnts)) / (S.D.^3)]
   
   =item Dis.F.
   Discrimination Factor: A Standard for evaluating the 
   problem according to a Criterion<br>
   
   =item [Criterion to group students into %27 Upper Students - 
   and %27 Lower Students]
   1st Criterion for Sorting the Students: 
   Sum of Partial Credit Awarded / Total Number of Tries
   2nd Criterion for Sorting the Students: 
   Total number of Correct Answers / Total Number of Tries
   
   =item Disc.
   Number of Students had at least one discussion.
   
   =back
   
   =cut
   
   ############################################################
   ############################################################
   
 1;  1;
 __END__  __END__

Removed from v.1.16  
changed lines
  Added in v.1.122.2.2


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