Yingyi Bu | e3d0ecd | 2016-02-24 10:56:59 -0800 | [diff] [blame] | 1 | /* |
| 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 | |
| 20 | drop dataverse test if exists; |
| 21 | create dataverse test; |
| 22 | use dataverse test; |
| 23 | |
| 24 | create 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 | |
| 34 | create dataset HRM(HRMType) primary key row_id; |
| 35 | |
| 36 | declare 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 | |
| 47 | declare 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 | |
| 53 | for $i in dataset HRM |
| 54 | group 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 |
| 56 | return |
| 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 | |