blob: 028e4369cd9048662a9930397460027aa7b58c57 [file] [log] [blame]
Yingyi Bue3d0ecd2016-02-24 10:56:59 -08001/*
2 * Licensed to the Apache Software Foundation (ASF) under one
3 * or more contributor license agreements. See the NOTICE file
4 * distributed with this work for additional information
5 * regarding copyright ownership. The ASF licenses this file
6 * to you under the Apache License, Version 2.0 (the
7 * "License"); you may not use this file except in compliance
8 * with the License. You may obtain a copy of the License at
9 *
10 * http://www.apache.org/licenses/LICENSE-2.0
11 *
12 * Unless required by applicable law or agreed to in writing,
13 * software distributed under the License is distributed on an
14 * "AS IS" BASIS, WITHOUT WARRANTIES OR CONDITIONS OF ANY
15 * KIND, either express or implied. See the License for the
16 * specific language governing permissions and limitations
17 * under the License.
18 */
19
20drop dataverse test if exists;
21create dataverse test;
22use dataverse test;
23
24create type HRMType as closed {
25 row_id: int32,
26 sid: int32,
27 date: date,
28 day: int32,
29 time: time,
30 bpm: int32,
31 RR: float
32};
33
34create dataset HRM(HRMType) primary key row_id;
35
36declare function median($x)
37{
38 let $c := count($x)
39 return avg(
40 for $i at $p in (for $j in $x order by $j return $j)
41 where $p >= floor(($c-1)/2.0) and $p <= $c/2
42 return $i
43 )
44}
45
46
47declare function MAD($x)
48{
49 let $m := median($x)
50 return median(for $j in $x order by abs($j-$m) return abs($j-$m))
51}
52
53for $i in dataset HRM
54group by $sid := $i.sid, $gdate := $i.date, $gday := $i.day,
55$timebin := interval-bin($i.time, time("00:00:00"), day-time-duration("PT15M")) with $i
56return
57{
58 "sid": $sid,
59 "gdate": $gdate,
60 "gday": $gday,
61 "timebin": $timebin,
62 "stdv": (avg(for $ii in $i return $ii.RR * $ii.RR) -
63 avg(for $ii in $i return $ii.RR) * avg(for $ii in $i return $ii.RR))^(0.5),
64 "MAD": MAD(for $ii in $i return $ii.RR)
65}
66