NAME
Statistics::Descriptive::PDL::Weighted - A close to drop-in replacement for Statistics::Descriptive::Weighted using PDL as the back-end
VERSION
Version 0.11
SYNOPSIS
use Statistics::Descriptive::PDL::Weighted;
my $stats = Statistics::Descriptive::PDL::Weighted->new;
$stats->add_data([1,2,3,4], [1,3,5,6]); # values then weights
my $mean = $stats->mean;
my $var = $stats->variance;
# or you can add data using a hash ref
my %data = (1 => 1, 2 => 3, 3 => 5, 4 => 6);
$stats->add_data(\%data);
# if you want equal weights then you need to supply them yourself
my $data = [1,2,3,4];
$stats->add_data($data, [(1) x scalar @$data]);
DESCRIPTION
This module provides basic functions used in descriptive statistics using weighted values.
METHODS
- new
-
Create a new statistics object. Takes no arguments.
- add_data (\%data)
- add_data ([1,2,3,4], [0.5,1,0.1,2)
-
Add data to the stats object. Appends to any existing data.
If a hash reference is passed then the keys are treated as the numeric data values, with the hash values the weights.
Unlike Statistics::Descriptive::PDL, you cannot pass a single flat array since odd things might happen if we convert it to a hash and the values are multidimensional.
Since the PDL::pdl function is used to process the data and weights you should be able to specify anything pdl accepts as valid.
An exception is raised the weights are <= 0, or are not the same size as the data.
- get_data
-
Returns arrays of the data and the weights.
In scalar context returns an array of arrays, i.e.
[\@data,\@wts]
. - get_data_as_hash
-
Returns the data as a perl hash, with the data values as the hash keys and weights as the hash values. Deduplicates the data if needed, incrementing the weights as appropriate.
Data values are stringified so there is obviously potential for loss of precision.
Returns a hash ref in scalar context.
- values_are_unique
- values_are_unique (1)
-
Flag to indicate if the data have duplicate values. Pass a true value to indicate your data have no duplicate values, making the median and percentile calculations faster (at the risk of you not being correct).
- sum_weights
-
Sum of the weights vector.
- sum_sqr_weights
-
Sum of the squared weights vector. Each weight is squared and the sum of these values then calculated.
- sum_sqr_sample_weights
-
Same as the
sum_sqr_weights
method. - Statistical methods
-
Most of the methods should need no explanation here, except to note that the standard_deviation, skewness and kurtosis use the biased methods. This is because one cannot guarantee the data are sample counts. The same applies to the median and percentiles. The median uses a centre of mass calculation, and the percentiles using analogous approach. This is because the weights are not guaranteed to be integers and so there is no sense interpolating.
Use Statistics::Descriptive::PDL::SampleWeighted when your weights are counts and you need the unbiased methods.
The iqr is the inter-quartile range, calculated as the difference of the 75th and 25th percentiles.
- geometric_mean
- harmonic_mean
- max
- mean
- median
- min
- mode
- sample_range
- standard_deviation
- sum
- variance
- count
- skewness
- kurtosis
- percentile (10)
- percentile (45)
- iqr
Not yet implemented, and possibly won't be.
Any of the trimmed functions, frequency functions and some others.
- least_squares_fit
- trimmed_mean
- quantile
- mindex
- maxdex
AUTHOR
Shawn Laffan, <shawnlaffan at gmail.com>
BUGS
Please report any bugs or feature requests to https://github.com/shawnlaffan/Statistics-Descriptive-PDL/issues
.
ACKNOWLEDGEMENTS
LICENSE AND COPYRIGHT
Copyright 2021 Shawn Laffan.
This program is free software; you can redistribute it and/or modify it under the terms of the the Artistic License (2.0). You may obtain a copy of the full license at:
http://www.perlfoundation.org/artistic_license_2_0
Any use, modification, and distribution of the Standard or Modified Versions is governed by this Artistic License. By using, modifying or distributing the Package, you accept this license. Do not use, modify, or distribute the Package, if you do not accept this license.
If your Modified Version has been derived from a Modified Version made by someone other than you, you are nevertheless required to ensure that your Modified Version complies with the requirements of this license.
This license does not grant you the right to use any trademark, service mark, tradename, or logo of the Copyright Holder.
This license includes the non-exclusive, worldwide, free-of-charge patent license to make, have made, use, offer to sell, sell, import and otherwise transfer the Package with respect to any patent claims licensable by the Copyright Holder that are necessarily infringed by the Package. If you institute patent litigation (including a cross-claim or counterclaim) against any party alleging that the Package constitutes direct or contributory patent infringement, then this Artistic License to you shall terminate on the date that such litigation is filed.
Disclaimer of Warranty: THE PACKAGE IS PROVIDED BY THE COPYRIGHT HOLDER AND CONTRIBUTORS "AS IS' AND WITHOUT ANY EXPRESS OR IMPLIED WARRANTIES. THE IMPLIED WARRANTIES OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE, OR NON-INFRINGEMENT ARE DISCLAIMED TO THE EXTENT PERMITTED BY YOUR LOCAL LAW. UNLESS REQUIRED BY LAW, NO COPYRIGHT HOLDER OR CONTRIBUTOR WILL BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, OR CONSEQUENTIAL DAMAGES ARISING IN ANY WAY OUT OF THE USE OF THE PACKAGE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.