|
1
|
|
|
<?php |
|
2
|
|
|
/** |
|
3
|
|
|
* webtrees-lib: MyArtJaub library for webtrees |
|
4
|
|
|
* |
|
5
|
|
|
* @package MyArtJaub\Webtrees |
|
6
|
|
|
* @subpackage Sosa |
|
7
|
|
|
* @author Jonathan Jaubart <[email protected]> |
|
8
|
|
|
* @copyright Copyright (c) 2009-2016, Jonathan Jaubart |
|
9
|
|
|
* @license http://www.gnu.org/licenses/gpl.html GNU General Public License, version 3 |
|
10
|
|
|
*/ |
|
11
|
|
|
namespace MyArtJaub\Webtrees\Module\Sosa; |
|
12
|
|
|
|
|
13
|
|
|
use Fisharebest\Webtrees\Auth; |
|
14
|
|
|
use Fisharebest\Webtrees\Controller\PageController; |
|
15
|
|
|
use Fisharebest\Webtrees\Filter; |
|
16
|
|
|
use Fisharebest\Webtrees\I18N; |
|
17
|
|
|
use Fisharebest\Webtrees\Individual; |
|
18
|
|
|
use Fisharebest\Webtrees\Module\AbstractModule; |
|
19
|
|
|
use MyArtJaub\Webtrees\Functions\Functions; |
|
20
|
|
|
use MyArtJaub\Webtrees\Globals; |
|
21
|
|
|
use MyArtJaub\Webtrees\Module\Sosa\Model\SosaProvider; |
|
22
|
|
|
use MyArtJaub\Webtrees\Mvc\Controller\MvcController; |
|
23
|
|
|
use MyArtJaub\Webtrees\Mvc\View\ViewBag; |
|
24
|
|
|
use MyArtJaub\Webtrees\Mvc\View\ViewFactory; |
|
25
|
|
|
|
|
26
|
|
|
/** |
|
27
|
|
|
* Controller for SosaStats |
|
28
|
|
|
*/ |
|
29
|
|
|
class SosaStatsController extends MvcController |
|
30
|
|
|
{ |
|
31
|
|
|
/** |
|
32
|
|
|
* Sosa Provider for the controller |
|
33
|
|
|
* @var SosaProvider $sosa_provider |
|
34
|
|
|
*/ |
|
35
|
|
|
protected $sosa_provider; |
|
36
|
|
|
|
|
37
|
|
|
/** |
|
38
|
|
|
* Constructor for SosaStatsController |
|
39
|
|
|
* @param AbstractModule $module |
|
40
|
|
|
*/ |
|
41
|
|
|
public function __construct(AbstractModule $module) { |
|
42
|
|
|
parent::__construct($module); |
|
43
|
|
|
|
|
44
|
|
|
$this->sosa_provider = new SosaProvider(Globals::getTree(), Auth::user()); |
|
45
|
|
|
} |
|
46
|
|
|
|
|
47
|
|
|
/** |
|
48
|
|
|
* Pages |
|
49
|
|
|
*/ |
|
50
|
|
|
|
|
51
|
|
|
/** |
|
52
|
|
|
* SosaStats@index |
|
53
|
|
|
*/ |
|
54
|
|
|
public function index() { |
|
55
|
|
|
$wt_tree = Globals::getTree(); |
|
56
|
|
|
$controller = new PageController(); |
|
57
|
|
|
$controller |
|
58
|
|
|
->setPageTitle(I18N::translate('Sosa Statistics')) |
|
59
|
|
|
->addInlineJavascript('$(".help_tooltip").tooltip();') |
|
60
|
|
|
; |
|
61
|
|
|
|
|
62
|
|
|
$view_bag = new ViewBag(); |
|
63
|
|
|
$view_bag->set('title', $controller->getPageTitle()); |
|
64
|
|
|
$view_bag->set('is_setup', false); |
|
65
|
|
|
|
|
66
|
|
|
if($this->sosa_provider->isSetup()) { |
|
67
|
|
|
$view_bag->set('is_setup', true); |
|
68
|
|
|
|
|
69
|
|
|
$view_bag->set('root_indi', $this->sosa_provider->getRootIndi()); |
|
70
|
|
|
|
|
71
|
|
|
$sosaCount = $this->sosa_provider->getSosaCount(); |
|
72
|
|
|
$diffSosaCount = $this->sosa_provider->getDifferentSosaCount(); |
|
73
|
|
|
|
|
74
|
|
|
$general_stats = array( |
|
75
|
|
|
'sosa_count' => $sosaCount, |
|
76
|
|
|
'distinct_count' => $diffSosaCount, |
|
77
|
|
|
'sosa_rate' => Functions::safeDivision($diffSosaCount, $this->sosa_provider->getTotalIndividuals()), |
|
78
|
|
|
'pedi_collapse' => 1 - Functions::safeDivision($diffSosaCount, $sosaCount), |
|
79
|
|
|
'mean_gen_time' => $this->sosa_provider->getMeanGenerationTime() |
|
80
|
|
|
); |
|
81
|
|
|
$view_bag->set('general_stats', $general_stats); |
|
82
|
|
|
|
|
83
|
|
|
$stats_gen = $this->sosa_provider->getStatisticsByGeneration(); |
|
84
|
|
|
$view_bag->set('missinganc_url', 'module.php?mod='.$this->module->getName().'&mod_action=SosaList@missing&ged='.$wt_tree->getNameUrl().'&gen='); |
|
85
|
|
|
$view_bag->set('sosaanc_url', 'module.php?mod='.$this->module->getName().'&mod_action=SosaList&ged='.$wt_tree->getNameUrl().'&gen='); |
|
86
|
|
|
|
|
87
|
|
|
$gen_theoretical=1; |
|
88
|
|
|
$total_theoretical=0; |
|
89
|
|
|
$prev_diff=0; |
|
90
|
|
|
$prev_known=0.5; |
|
91
|
|
|
$generation_stats = array(); |
|
92
|
|
|
|
|
93
|
|
|
foreach($stats_gen as $gen => $tab){ |
|
94
|
|
|
$genY1= I18N::translate('-'); |
|
95
|
|
|
$genY2= I18N::translate('-'); |
|
96
|
|
|
if($tab['firstBirth']>0) $genY1=$tab['firstEstimatedBirth']; |
|
97
|
|
|
if($tab['lastBirth']>0) $genY2=$tab['lastEstimatedBirth']; |
|
98
|
|
|
$total_theoretical += $gen_theoretical; |
|
99
|
|
|
$perc_sosa_count_theor = Functions::safeDivision($tab['sosaCount'], $gen_theoretical); |
|
100
|
|
|
$missing=2*$prev_known - $tab['sosaCount']; |
|
101
|
|
|
$gen_diff=$tab['diffSosaTotalCount']-$prev_diff; |
|
102
|
|
|
|
|
103
|
|
|
$generation_stats[$gen] = array( |
|
104
|
|
|
'gen_min_birth' => $genY1, |
|
105
|
|
|
'gen_max_birth' => $genY2, |
|
106
|
|
|
'theoretical' => $gen_theoretical, |
|
107
|
|
|
'known' => $tab['sosaCount'], |
|
108
|
|
|
'perc_known' => $perc_sosa_count_theor, |
|
109
|
|
|
'missing' => $missing, |
|
110
|
|
|
'perc_missing' => 1-Functions::safeDivision($tab['sosaCount'], 2*$prev_known), |
|
111
|
|
|
'total_known' => $tab['sosaTotalCount'], |
|
112
|
|
|
'perc_total_known' => Functions::safeDivision($tab['sosaTotalCount'], $total_theoretical), |
|
113
|
|
|
'different' => $gen_diff, |
|
114
|
|
|
'perc_different' => Functions::safeDivision($gen_diff, $tab['sosaCount']), |
|
115
|
|
|
'total_different' => $tab['diffSosaTotalCount'], |
|
116
|
|
|
'pedi_collapse' => 1 - Functions::safeDivision($tab['diffSosaTotalCount'], $tab['sosaTotalCount']) |
|
117
|
|
|
); |
|
118
|
|
|
|
|
119
|
|
|
$gen_theoretical = $gen_theoretical * 2; |
|
120
|
|
|
$prev_known=$tab['sosaCount']; |
|
121
|
|
|
$prev_diff=$tab['diffSosaTotalCount']; |
|
122
|
|
|
} |
|
123
|
|
|
|
|
124
|
|
|
$view_bag->set('generation_stats', $generation_stats); |
|
125
|
|
|
|
|
126
|
|
|
$gen_depth_stats = $this->sosa_provider->getGenerationDepthStatsAtGen(1); |
|
127
|
|
|
$view_bag->set('mean_gen_depth', count($gen_depth_stats) > 0 ? $gen_depth_stats[1]['mean_gen_depth'] : 0); |
|
128
|
|
|
$view_bag->set('stddev_gen_depth', count($gen_depth_stats) > 0 ? $gen_depth_stats[1]['stddev_gen_depth'] : 0); |
|
129
|
|
|
|
|
130
|
|
|
$view_bag->set('top10multiancestors', $this->getTop10Ancestors()); |
|
131
|
|
|
|
|
132
|
|
|
$view_bag->set('chart_img_g2', $this->htmlAncestorDispersionG2()); |
|
133
|
|
|
$view_bag->set('chart_img_g3', $this->htmlAncestorDispersionG3()); |
|
134
|
|
|
|
|
135
|
|
|
$view_bag->set('chart_img_gendepth3', $this->htmlAncestorGenDepthG3()); |
|
136
|
|
|
//$view_bag->set('mean_gen_depth_gen3_stats', $this->sosa_provider->getMeanGeneratuonDepthAndDeviationAtGen(3)); |
|
137
|
|
|
} |
|
138
|
|
|
|
|
139
|
|
|
ViewFactory::make('SosaStats', $this, $controller, $view_bag)->render(); |
|
140
|
|
|
} |
|
141
|
|
|
|
|
142
|
|
|
/** |
|
143
|
|
|
* Return an array of the top 10 of ancestor appearing multiple times, grouped by the number of occurrences. |
|
144
|
|
|
* The key is the number of appearance, and the value contains the list of Inidvidual |
|
145
|
|
|
* |
|
146
|
|
|
* @return array Top 10 ancestors with count |
|
147
|
|
|
*/ |
|
148
|
|
|
private function getTop10Ancestors() |
|
149
|
|
|
{ |
|
150
|
|
|
$top10multiancestors = $this->sosa_provider->getTopMultiSosaAncestorsNoTies(10); |
|
151
|
|
|
$top10ancestors = array(); |
|
152
|
|
|
if($top10multiancestors !== null && count($top10multiancestors)) { |
|
153
|
|
|
foreach($top10multiancestors as $pid => $count) { |
|
154
|
|
|
$indi = Individual::getInstance($pid, $this->sosa_provider->getTree()); |
|
155
|
|
|
if($indi !== null && $indi->canShowName()) { |
|
156
|
|
|
array_key_exists($count, $top10ancestors) ? |
|
157
|
|
|
$top10ancestors[$count][] = $indi: |
|
158
|
|
|
$top10ancestors[$count] = array($count => $indi); |
|
159
|
|
|
} |
|
160
|
|
|
} |
|
161
|
|
|
} |
|
162
|
|
|
return $top10ancestors; |
|
163
|
|
|
} |
|
164
|
|
|
|
|
165
|
|
|
/** |
|
166
|
|
|
* Returns HTML code for a graph showing the dispersion of ancestors across father & mother |
|
167
|
|
|
* @return string HTML code |
|
168
|
|
|
*/ |
|
169
|
|
|
private function htmlAncestorDispersionG2() |
|
170
|
|
|
{ |
|
171
|
|
|
$ancestorsDispGen2 = $this->sosa_provider->getAncestorDispersionForGen(2); |
|
172
|
|
|
if(count($ancestorsDispGen2) == 0) return; |
|
173
|
|
|
|
|
174
|
|
|
$size = '600x300'; |
|
175
|
|
|
|
|
176
|
|
|
$total = array_sum($ancestorsDispGen2); |
|
177
|
|
|
$father_count = array_key_exists(1, $ancestorsDispGen2) ? $ancestorsDispGen2[1] : 0; |
|
178
|
|
|
$father = array ( |
|
179
|
|
|
'color' => '84beff', |
|
180
|
|
|
'count' => $father_count, |
|
181
|
|
|
'perc' => Functions::safeDivision($father_count, $total), |
|
182
|
|
|
'name' => \Fisharebest\Webtrees\Functions\Functions::getRelationshipNameFromPath('fat') |
|
183
|
|
|
); |
|
184
|
|
|
$mother_count = array_key_exists(2, $ancestorsDispGen2) ? $ancestorsDispGen2[2] : 0; |
|
185
|
|
|
$mother = array ( |
|
186
|
|
|
'color' => 'ffd1dc', |
|
187
|
|
|
'count' => $mother_count, |
|
188
|
|
|
'perc' => Functions::safeDivision($mother_count, $total), |
|
189
|
|
|
'name' => \Fisharebest\Webtrees\Functions\Functions::getRelationshipNameFromPath('mot') |
|
190
|
|
|
); |
|
191
|
|
|
$shared_count = array_key_exists(-1, $ancestorsDispGen2) ? $ancestorsDispGen2[-1] : 0; |
|
192
|
|
|
$shared = array ( |
|
193
|
|
|
'color' => '777777', |
|
194
|
|
|
'count' => $shared_count, |
|
195
|
|
|
'perc' => Functions::safeDivision($shared_count, $total), |
|
196
|
|
|
'name' => I18N::translate('Shared') |
|
197
|
|
|
); |
|
198
|
|
|
|
|
199
|
|
|
$chd = $this->arrayToExtendedEncoding(array(4095 * $father['perc'], 4095 * $shared['perc'], 4095 * $mother['perc'])); |
|
200
|
|
|
$chart_title = I18N::translate('Known Sosa ancestors\' dispersion'); |
|
201
|
|
|
$chl = |
|
202
|
|
|
$father['name'] . ' - ' . I18N::percentage($father['perc'], 1) . '|' . |
|
203
|
|
|
$shared['name'] . ' - ' . I18N::percentage($shared['perc'], 1) . '|' . |
|
204
|
|
|
$mother['name'] . ' - ' . I18N::percentage($mother['perc'], 1); |
|
205
|
|
|
return "<img src=\"https://chart.googleapis.com/chart?cht=p&chp=1.5708&chd=e:{$chd}&chs={$size}&chco={$father['color']},{$shared['color']},{$mother['color']}&chf=bg,s,ffffff00&chl={$chl}\" alt=\"" . $chart_title . "\" title=\"" . $chart_title . "\" />"; |
|
206
|
|
|
} |
|
207
|
|
|
|
|
208
|
|
|
/** |
|
209
|
|
|
* Returns HTML code for a graph showing the dispersion of ancestors across grand-parents |
|
210
|
|
|
* @return string HTML code |
|
211
|
|
|
*/ |
|
212
|
|
|
private function htmlAncestorDispersionG3() |
|
213
|
|
|
{ |
|
214
|
|
|
$ancestorsDispGen2 = $this->sosa_provider->getAncestorDispersionForGen(3); |
|
215
|
|
|
|
|
216
|
|
|
$size = '700x300'; |
|
217
|
|
|
|
|
218
|
|
|
$color_motmot = 'ffd1dc'; |
|
219
|
|
|
$color_motfat = 'b998a0'; |
|
220
|
|
|
$color_fatfat = '577292'; |
|
221
|
|
|
$color_fatmot = '84beff'; |
|
222
|
|
|
$color_shared = '777777'; |
|
223
|
|
|
|
|
224
|
|
|
$total_fatfat = array_key_exists(1, $ancestorsDispGen2) ? $ancestorsDispGen2[1] : 0; |
|
225
|
|
|
$total_fatmot = array_key_exists(2, $ancestorsDispGen2) ? $ancestorsDispGen2[2] : 0; |
|
226
|
|
|
$total_motfat = array_key_exists(4, $ancestorsDispGen2) ? $ancestorsDispGen2[4] : 0; |
|
227
|
|
|
$total_motmot = array_key_exists(8, $ancestorsDispGen2) ? $ancestorsDispGen2[8] : 0; |
|
228
|
|
|
$total_sha = array_key_exists(-1, $ancestorsDispGen2) ? $ancestorsDispGen2[-1] : 0; |
|
229
|
|
|
$total = $total_fatfat + $total_fatmot + $total_motfat+ $total_motmot + $total_sha; |
|
230
|
|
|
|
|
231
|
|
|
$chd = $this->arrayToExtendedEncoding(array( |
|
232
|
|
|
4095 * Functions::safeDivision($total_fatfat, $total), |
|
233
|
|
|
4095 * Functions::safeDivision($total_fatmot, $total), |
|
234
|
|
|
4095 * Functions::safeDivision($total_sha, $total), |
|
235
|
|
|
4095 * Functions::safeDivision($total_motfat, $total), |
|
236
|
|
|
4095 * Functions::safeDivision($total_motmot, $total) |
|
237
|
|
|
)); |
|
238
|
|
|
$chart_title = I18N::translate('Known Sosa ancestors\' dispersion - G3'); |
|
239
|
|
|
$chl = |
|
240
|
|
|
\Fisharebest\Webtrees\Functions\Functions::getRelationshipNameFromPath('fatfat') . ' - ' . I18N::percentage(Functions::safeDivision($total_fatfat, $total), 1) . '|' . |
|
241
|
|
|
\Fisharebest\Webtrees\Functions\Functions::getRelationshipNameFromPath('fatmot') . ' - ' . I18N::percentage(Functions::safeDivision($total_fatmot, $total), 1) . '|' . |
|
242
|
|
|
I18N::translate('Shared') . ' - ' . I18N::percentage(Functions::safeDivision($total_sha, $total), 1) . '|' . |
|
243
|
|
|
\Fisharebest\Webtrees\Functions\Functions::getRelationshipNameFromPath('motfat') . ' - ' . I18N::percentage(Functions::safeDivision($total_motfat, $total), 1) . '|' . |
|
244
|
|
|
\Fisharebest\Webtrees\Functions\Functions::getRelationshipNameFromPath('motmot') . ' - ' . I18N::percentage(Functions::safeDivision($total_motmot, $total), 1); |
|
245
|
|
|
return "<img src=\"https://chart.googleapis.com/chart?cht=p&chp=1.5708&chd=e:{$chd}&chs={$size}&chco={$color_fatfat},{$color_fatmot},{$color_shared},{$color_motfat},{$color_motmot}&chf=bg,s,ffffff00&chl={$chl}\" alt=\"" . $chart_title . "\" title=\"" . $chart_title . "\" />"; |
|
246
|
|
|
} |
|
247
|
|
|
|
|
248
|
|
|
/** |
|
249
|
|
|
* Returns HTML code for a graph showing the grand-parents' mean generation depth and standard deviation |
|
250
|
|
|
* @return string HTML code |
|
251
|
|
|
*/ |
|
252
|
|
|
private function htmlAncestorGenDepthG3() { |
|
253
|
|
|
$ancestorsGenDepth3 = $this->sosa_provider->getGenerationDepthStatsAtGen(3); |
|
254
|
|
|
if(count($ancestorsGenDepth3) == 0) return; |
|
255
|
|
|
|
|
256
|
|
|
$ancestors = array(); |
|
257
|
|
|
$chd_mean = array(); |
|
258
|
|
|
$chd_error_low = array(); |
|
259
|
|
|
$chd_error_high = array(); |
|
260
|
|
|
foreach($ancestorsGenDepth3 as $sosa => $genDepthStat) { |
|
261
|
|
|
$ancestor = Individual::getInstance($genDepthStat['root_ancestor_id'], $this->sosa_provider->getTree()); |
|
262
|
|
|
if($ancestor !== null && $ancestor->canShowName()) { |
|
263
|
|
|
$tmp = $ancestor->getAllNames(); |
|
264
|
|
|
$ancestors[] = Filter::escapeUrl($tmp[$ancestor->getPrimaryName()]['fullNN']); |
|
265
|
|
|
} |
|
266
|
|
|
else { |
|
267
|
|
|
$ancestors[] = I18N::translate('Sosa %s', I18N::number($sosa)); |
|
268
|
|
|
} |
|
269
|
|
|
$chd_mean[] = $genDepthStat['mean_gen_depth']; |
|
270
|
|
|
$chd_error_low[] = $genDepthStat['mean_gen_depth'] - $genDepthStat['stddev_gen_depth']; |
|
271
|
|
|
$chd_error_high[] = $genDepthStat['mean_gen_depth'] + $genDepthStat['stddev_gen_depth']; |
|
272
|
|
|
} |
|
273
|
|
|
|
|
274
|
|
|
$maxChd = ceil(max($chd_error_high)); |
|
275
|
|
|
$chd = implode(',', $chd_mean) . '|' . implode(',', $chd_error_low) . '|' . implode(',', $chd_error_high); |
|
276
|
|
|
$chxl = implode('|', array_reverse($ancestors)); |
|
277
|
|
|
$chbh = 30; |
|
278
|
|
|
$chs = 4 * $chbh + 50; |
|
279
|
|
|
$chxl_title = I18N::translate('Mean generation depth and standard deviation'); |
|
280
|
|
|
|
|
281
|
|
|
return "<img src=\"https://chart.googleapis.com/chart?cht=bhs&chco=84beff&chs=600x{$chs}&chbh={$chbh}&chd=t1:{$chd}&chds=0,{$maxChd}&chxt=x,x,y&chxl=1:|{$chxl_title}|2:|{$chxl}&chxp=1,50&chxr=0,0,{$maxChd},1&chm=E,577292,1:2,,2\" />"; |
|
282
|
|
|
} |
|
283
|
|
|
|
|
284
|
|
|
/** |
|
285
|
|
|
* Convert an array to Google Chart encoding |
|
286
|
|
|
* @param arrat $a Array to encode |
|
287
|
|
|
* @return string |
|
288
|
|
|
*/ |
|
289
|
|
|
private function arrayToExtendedEncoding($a) { |
|
290
|
|
|
$xencoding = WT_GOOGLE_CHART_ENCODING; |
|
291
|
|
|
|
|
292
|
|
|
$encoding = ''; |
|
293
|
|
|
foreach ($a as $value) { |
|
294
|
|
|
if ($value < 0) { |
|
295
|
|
|
$value = 0; |
|
296
|
|
|
} |
|
297
|
|
|
$first = (int) ($value / 64); |
|
298
|
|
|
$second = $value % 64; |
|
299
|
|
|
$encoding .= $xencoding[(int) $first] . $xencoding[(int) $second]; |
|
300
|
|
|
} |
|
301
|
|
|
|
|
302
|
|
|
return $encoding; |
|
303
|
|
|
} |
|
304
|
|
|
|
|
305
|
|
|
|
|
306
|
|
|
} |
jon48 /
webtrees-lib
