'."\n".
''."\n".
''.
@@ -385,7 +844,6 @@ sub output_html_stats {
return;
}
-
sub html_preamble {
my $Str='';
$Str .= " ';
+ return $Str;
+}
+
+
+sub sequence_html_output {
+ my ($seq) = @_;
+ my $data = $SeqStat{$seq->symb};
+ my $row = '';
+ foreach my $field (@SeqFields) {
+ next if ($field->{'selected'} ne 'yes');
+ $row .= '{'align'})) {
+ $row .= ' align="'.$field->{'align'}.'"';
+ }
+ $row .= '>';
+ if (exists($field->{'format'})) {
+ $row .= sprintf($field->{'format'},$data->{$field->{'name'}});
+ } else {
+ $row .= $data->{$field->{'name'}};
+ }
+ $row .= ' ';
+ }
+ $row .= ' '."\n";
+ return $row;
+}
+
####################################################
####################################################
##
@@ -480,6 +974,8 @@ sub make_plot {
&sort_data($ENV{'form.sortby'});
if ($plot eq 'degrees') {
°rees_plot($r);
+ } elsif ($plot eq 'tries statistics') {
+ &tries_data_plot($r);
} else {
&make_single_stat_plot($r,$plot);
}
@@ -592,6 +1088,11 @@ sub degrees_plot {
my $diffdata .= ''.join(',',@Labels).' '.$/.
''.join(',',@Diff).' '.$/;
#
+ my $title = 'Degree of Discrimination\nand Degree of Difficulty';
+ if ($xmax > 50) {
+ $title = 'Degree of Discrimination and Degree of Difficulty';
+ }
+ #
$plot=<<"END";
Degree of Discrmination and Degree of Difficulty
+ $title
Problem Number
{'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 {
+ if ($xmax % 5 != 0) {
+ $xmax = 5 * (int($xmax/5)+1);
+ }
+ }
+ $ymax = int($ymax)+1+2;
+ #
+ my $std_data .= ''.join(',',@Labels).' '.$/.
+ ''.join(',',@Mean).' '.$/;
+ #
+ my $std_error_data .= ''.join(',',@Labels).' '.$/.
+ ''.join(',',@Mean).' '.$/.
+ ''.join(',',@STD).' '.$/;
+ #
+ my $title = 'Mean and S.D. of Tries';
+ if ($xmax > 25) {
+ $title = 'Mean and Standard Deviation of Tries';
+ }
+ #
+ $plot=<<"END";
+
+ $title
+ Problem Number
+ Number of Tries
+
+ $std_error_data
+
+
+ $std_data
+
+
+END
+ my $plotresult =
+ ''.&Apache::lonxml::xmlparse($r,'web',$plot).'
'.$/;
+ $r->print($plotresult);
+ return;
+}
+
sub plot_dropdown {
my $current = '';
#
@@ -646,12 +1232,15 @@ sub plot_dropdown {
my @Additional_Plots = (
{ graphable=>'yes',
name => 'degrees',
- title => 'DoDisc and DoDiff' });
+ title => 'Difficulty Indexes' },
+ { graphable=>'yes',
+ name => 'tries statistics',
+ title => 'Tries Statistics' });
#
my $Str= "\n".'';
$Str .= 'none '."\n";
- foreach my $field (@Fields,@Additional_Plots) {
+ foreach my $field (@Additional_Plots,@Fields) {
if (! exists($field->{'graphable'}) ||
$field->{'graphable'} ne 'yes') {
next;
@@ -681,30 +1270,13 @@ sub Excel_output {
&compute_all_statistics($r);
my $c = $r->connection;
return if ($c->aborted());
+ #
+ my ($starttime,$endtime) = &Apache::lonstathelpers::get_time_limits();
##
## Create the excel workbook
- my $filename = '/prtspool/'.
- $ENV{'user.name'}.'_'.$ENV{'user.domain'}.'_'.
- time.'_'.rand(1000000000).'.xls';
- my ($starttime,$endtime) = &Apache::lonstathelpers::get_time_limits();
- #
- # Create sheet
- my $excel_workbook = Spreadsheet::WriteExcel->new('/home/httpd'.$filename);
- #
- # Check for errors
- if (! defined($excel_workbook)) {
- $r->log_error("Error creating excel spreadsheet $filename: $!");
- $r->print(&mt("Problems creating new Excel file. ".
- "This error has been logged. ".
- "Please alert your LON-CAPA administrator."));
- return 0;
- }
- #
- # The excel spreadsheet stores temporary data in files, then put them
- # together. If needed we should be able to disable this (memory only).
- # The temporary directory must be specified before calling 'addworksheet'.
- # File::Temp is used to determine the temporary directory.
- $excel_workbook->set_tempdir($Apache::lonnet::tmpdir);
+ my ($excel_workbook,$filename,$format) =
+ &Apache::loncommon::create_workbook($r);
+ return if (! defined($excel_workbook));
#
# Add a worksheet
my $sheetname = $ENV{'course.'.$ENV{'request.course.id'}.'.description'};
@@ -720,28 +1292,16 @@ sub Excel_output {
#
# Put the course description in the header
$excel_sheet->write($rows_output,$cols_output++,
- $ENV{'course.'.$ENV{'request.course.id'}.'.description'});
+ $ENV{'course.'.$ENV{'request.course.id'}.'.description'},
+ $format->{'h1'});
$cols_output += 3;
#
# Put a description of the sections listed
my $sectionstring = '';
- my @Sections = @Apache::lonstatistics::SelectedSections;
- if (scalar(@Sections) > 1) {
- if (scalar(@Sections) > 2) {
- my $last = pop(@Sections);
- $sectionstring = "Sections ".join(', ',@Sections).', and '.$last;
- } else {
- $sectionstring = "Sections ".join(' and ',@Sections);
- }
- } else {
- if ($Sections[0] eq 'all') {
- $sectionstring = "All sections";
- } else {
- $sectionstring = "Section ".$Sections[0];
- }
- }
- $excel_sheet->write($rows_output,$cols_output++,$sectionstring);
- $cols_output += scalar(@Sections);
+ $excel_sheet->write($rows_output,$cols_output++,
+ &Apache::lonstatistics::section_and_enrollment_description('plaintext'),
+ $format->{'h3'});
+ $cols_output += scalar(&Apache::lonstatistics::get_selected_sections());
#
# Time restrictions
my $time_string;
@@ -757,6 +1317,10 @@ sub Excel_output {
# See note above about lonlocal:locallocaltime
$time_string .= 'Data collected before '.localtime($endtime).'.';
}
+ if (defined($time_string)) {
+ $excel_sheet->write($rows_output,$cols_output++,$time_string);
+ $cols_output+= 5;
+ }
#
# Put the date in there too
$excel_sheet->write($rows_output,$cols_output++,
@@ -764,35 +1328,48 @@ sub Excel_output {
#
$rows_output++;
$cols_output=0;
- #
- # Long Headers
- foreach my $field (@Fields) {
- next if ($field->{'name'} eq 'problem_num');
- if (exists($field->{'long_title'})) {
+ ##
+ ## Sequence Statistics
+ ##
+ &write_headers($excel_sheet,$format,\$rows_output,\$cols_output,
+ \@SeqFields);
+ foreach my $seq (@sequences) {
+ my $data = $SeqStat{$seq->symb};
+ $cols_output=0;
+ foreach my $field (@SeqFields) {
+ next if ($field->{'selected'} ne 'yes');
+ my $fieldformat = undef;
+ if (exists($field->{'excel_format'})) {
+ $fieldformat = $format->{$field->{'excel_format'}};
+ }
$excel_sheet->write($rows_output,$cols_output++,
- $field->{'long_title'});
- } else {
- $excel_sheet->write($rows_output,$cols_output++,'');
+ $data->{$field->{'name'}},$fieldformat);
}
+ $rows_output++;
+ $cols_output=0;
}
+ ##
+ ## Resource Statistics
+ ##
$rows_output++;
$cols_output=0;
- # Brief headers
- foreach my $field (@Fields) {
- 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'});
- }
- $rows_output++;
+ &write_headers($excel_sheet,$format,\$rows_output,\$cols_output,
+ \@Fields);
+ #
foreach my $data (@StatsArray) {
$cols_output=0;
foreach my $field (@Fields) {
+ next if ($field->{'selected'} ne 'yes');
next if ($field->{'name'} eq 'problem_num');
+ my $fieldformat = undef;
+ if (exists($field->{'excel_format'})) {
+ $fieldformat = $format->{$field->{'excel_format'}};
+ }
$excel_sheet->write($rows_output,$cols_output++,
- $data->{$field->{'name'}});
+ $data->{$field->{'name'}},$fieldformat);
}
$rows_output++;
+ $cols_output=0;
}
#
$excel_workbook->close();
@@ -805,6 +1382,44 @@ sub Excel_output {
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;
+}
+
##################################################
##################################################
##
@@ -815,9 +1430,9 @@ sub Excel_output {
sub compute_statistics_on_sequence {
my ($seq) = @_;
my @Data;
- foreach my $res (@{$seq->{'contents'}}) {
- next if ($res->{'type'} ne 'assessment');
- foreach my $part (@{$res->{'parts'}}) {
+ foreach my $res (&Apache::lonstathelpers::get_resources($navmap,$seq)) {
+ foreach my $part (@{$res->parts}) {
+ next if ($res->is_survey($part));
#
# This is where all the work happens
my $data = &get_statistics($seq,$res,$part,scalar(@StatsArray)+1);
@@ -835,9 +1450,9 @@ sub compute_all_statistics {
return;
}
my $c = $r->connection;
- foreach my $seq (&Apache::lonstatistics::Sequences_with_Assess()) {
+ foreach my $seq (@sequences) {
last if ($c->aborted);
- next if ($seq->{'num_assess'} < 1);
+ &compute_sequence_statistics($seq);
&compute_statistics_on_sequence($seq);
}
}
@@ -923,25 +1538,59 @@ sub get_statistics {
my ($sequence,$resource,$part,$problem_num) = @_;
#
my ($starttime,$endtime) = &Apache::lonstathelpers::get_time_limits();
- my $symb = $resource->{'symb'};
+ my $symb = $resource->symb;
my $courseid = $ENV{'request.course.id'};
#
my $data = &Apache::loncoursedata::get_problem_statistics
- (\@Apache::lonstatistics::SelectedSections,
+ ([&Apache::lonstatistics::get_selected_sections()],
$Apache::lonstatistics::enrollment_status,
$symb,$part,$courseid,$starttime,$endtime);
+ $data->{'symb'} = $symb;
$data->{'part'} = $part;
$data->{'problem_num'} = $problem_num;
- $data->{'container'} = $sequence->{'title'};
- $data->{'title'} = $resource->{'title'};
- $data->{'title.link'} = $resource->{'src'}.'?symb='.
- &Apache::lonnet::escape($resource->{'symb'});
- #
- $data->{'deg_of_disc'} = &compute_discrimination_factor($resource,$part,$sequence);
+ $data->{'container'} = $sequence->compTitle;
+ $data->{'title'} = $resource->compTitle;
+ $data->{'title.link'} = $resource->src.'?symb='.
+ &Apache::lonnet::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];
+ $data->{'urlres'}=$urlres;
+ my %storestats =
+ &LONCAPA::lonmetadata::dynamic_metadata_storage($data);
+ my ($dom,$user) = $urlres=~/^(\w+)\/(\w+)/;
+ &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)
+# $data->{'duedate'} =
+# &Apache::lonnet::EXT('resource.'.$part.'.duedate',$symb);
+# $data->{'opendate'} =
+# &Apache::lonnet::EXT('resource.'.$part.'.opendate',$symb);
+# $data->{'maxtries'} =
+# &Apache::lonnet::EXT('resource.'.$part.'.maxtries',$symb);
+# $data->{'hinttries'} =
+# &Apache::lonnet::EXT('resource.'.$part.'.hinttries',$symb);
+ $data->{'weight'} =
+ &Apache::lonnet::EXT('resource.'.$part.'.weight',$symb);
+# $data->{'resptypes'} = join(',',@{$resource->{'partdata'}->{$part}->{'ResponseTypes'}});
return $data;
}
-
###############################################
###############################################
@@ -958,19 +1607,22 @@ Returns: integer between -1 and 1
###############################################
###############################################
sub compute_discrimination_factor {
- my ($resource,$part,$sequence) = @_;
+ my ($resource,$part,$seq) = @_;
+ my $symb = $resource->symb;
my @Resources;
- foreach my $res (@{$sequence->{'contents'}}) {
- next if ($res->{'symb'} eq $resource->{'symb'});
- push (@Resources,$res->{'symb'});
+ 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::SelectedSections,
- $Apache::lonstatistics::enrollment_status,undef);
+ [&Apache::lonstatistics::get_selected_sections()],
+ $Apache::lonstatistics::enrollment_status,undef,
+ $starttime,$endtime);
#
# compute their percent scores on the problems in the sequence,
my $number_to_grab = int(scalar(@{$ranking})/4);
@@ -981,10 +1633,16 @@ sub compute_discrimination_factor {
map {
$_->[&Apache::loncoursedata::RNK_student()];
} @{$ranking}[($num_students-$number_to_grab)..($num_students-1)];
+ 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($resource,$part,\@BottomSet);
+ &Apache::loncoursedata::get_sum_of_scores($symb,$part,\@BottomSet,
+ undef,$starttime,$endtime);
my ($top_sum,$top_max) =
- &Apache::loncoursedata::get_sum_of_scores($resource,$part,\@TopSet);
+ &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';
@@ -998,6 +1656,71 @@ sub compute_discrimination_factor {
###############################################
###############################################
+##
+## 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::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::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;
+}
+
+
=pod
@@ -1058,111 +1781,6 @@ Number of Students had at least one disc
=cut
-
-############################################################
-############################################################
-##
-## How this all works:
-## Statistics are computed by calling &get_statistics with the sequence,
-## resource, and part id to run statistics on. At various places within
-## the loops which compute the statistics, as well as before and after
-## the entire process, subroutines can be called. The subroutines are
-## registered to the following hooks:
-##
-## hook subroutine inputs
-## ----------------------------------------------------------
-## pre $r,$count
-## pre_seq $r,$count,$seq
-## pre_res $r,$count,$seq,$res
-## calc $r,$count,$seq,$res,$data
-## post_res $r,$count,$seq,$res
-## post_seq $r,$count,$seq
-## post $r,$count
-##
-## abort $r
-##
-## subroutines will be called in the order in which they are registered.
-##
-############################################################
-############################################################
-{
-
-my %hooks;
-my $aborted = 0;
-
-sub abort_computation {
- $aborted = 1;
-}
-
-sub clear_hooks {
- $aborted = 0;
- undef(%hooks);
-}
-
-sub register_hook {
- my ($hookname,$subref)=@_;
- if ($hookname !~ /^(pre|pre_seq|pre_res|post|post_seq|post_res|calc)$/){
- return;
- }
- if (ref($subref) ne 'CODE') {
- &Apache::lonnet::logthis('attempt to register hook to non-code: '.
- $hookname,' = '.$subref);
- } else {
- if (exists($hooks{$hookname})) {
- push(@{$hooks{$hookname}},$subref);
- } else {
- $hooks{$hookname} = [$subref];
- }
- }
- return;
-}
-
-sub run_hooks {
- my $context = shift();
- foreach my $hook (@{$hooks{$context}}) {
- if ($aborted && $context ne 'abort') {
- last;
- }
- my $retvalue = $hook->(@_);
- if (defined($retvalue) && $retvalue eq '0') {
- $aborted = 1 if (! $aborted);
- }
- }
-}
-
-sub run_statistics {
- my ($r) = @_;
- my $count = 0;
- &run_hooks('pre',$r,$count);
- foreach my $seq (&Apache::lonstatistics::Sequences_with_Assess()) {
- last if ($aborted);
- next if ($seq->{'num_assess'}<1);
- &run_hooks('pre_seq',$r,$count,$seq);
- foreach my $res (@{$seq->{'contents'}}) {
- last if ($aborted);
- next if ($res->{'type'} ne 'assessment');
- &run_hooks('pre_res',$r,$count,$seq,$res);
- foreach my $part (@{$res->{'parts'}}) {
- last if ($aborted);
- #
- # This is where all the work happens
- my $data = &get_statistics($seq,$res,$part,++$count);
- &run_hooks('calc',$r,$count,$seq,$res,$part,$data);
- }
- &run_hooks('post_res',$r,$count,$seq,$res);
- }
- &run_hooks('post_seq',$r,$count,$seq);
- }
- if ($aborted) {
- &run_hooks('abort',$r);
- } else {
- &run_hooks('post',$r,$count);
- }
- return;
-}
-
-} # End of %hooks scope
-
############################################################
############################################################