| Literature DB >> 34294827 |
Amir Baghdadi1, Sanju Lama1, Rahul Singh1,2, Hamidreza Hoshyarmanesh1, Mohammadsaleh Razmi1, Garnette R Sutherland3.
Abstract
Surgical error and resulting complication have significant patient and economic consequences. Inappropriate exertion of tool-tissue force is a common variable for such error, that can be objectively monitored by sensorized tools. The rich digital output establishes a powerful skill assessment and sharing platform for surgical performance and training. Here we present SmartForceps data app incorporating an Expert Room environment for tracking and analysing the objective performance and surgical finesse through multiple interfaces specific for surgeons and data scientists. The app is enriched by incoming geospatial information, data distribution for engineered features, performance dashboard compared to expert surgeon, and interactive skill prediction and task recognition tools to develop artificial intelligence models. The study launches the concept of democratizing surgical data through a connectivity interface between surgeons with a broad and deep capability of geographic reach through mobile devices with highly interactive infographics and tools for performance monitoring, comparison, and improvement.Entities:
Year: 2021 PMID: 34294827 PMCID: PMC8298519 DOI: 10.1038/s41598-021-94487-9
Source DB: PubMed Journal: Sci Rep ISSN: 2045-2322 Impact factor: 4.379
Figure 1SmartForceps timeseries data of the Right prong across the 5 surgical tasks of Retracting, Manipulation, Dissecting, Pulling, and Coagulation overlaid for 50 cases. Differences in the range and duration of force are shown in the overlaid data profiles. Please refer to our Supplementary Materials for both left and right prong data. These charts have fully interactive capability including zoom, pan, download, etc. Figure created by R Plotly library version 2.0: https://plotly.com/r/.
Figure 2Design interface of the SmartForceps surgical data monitoring application on a mobile device. The user can click to visit the general dashboard without login or login with exclusive credentials to visit their own reports in surgeon or data scientist views. Visualization created through framing the PWA created in JavaScript inside a phone view using MockuPhone mock-up generator: https://mockuphone.com.
Figure 3SmartForceps Data Analytics Dashboard in "General" view shown in a desktop mode. The current view includes two tabs of "Geospatial Information" (a) and “Surgical Force Data” (b). These charts have fully interactive capability including zoom, pan, download, etc. The figure for dashboard was created by Python Dash library version 0.43.0: https://plotly.com/dash/.
Figure 4SmartForceps Data Analytics Dashboard in "Surgeon" view. The current view includes three tabs of "Geospatial Information", “Surgical Force Data”, and “Performance Comparison Dashboard”. These charts have fully interactive capability including zoom, pan, download, etc. This figure shows the overtime performance report (with the slide bar at the top to select range of cases) for a Novice surgeon with PGY > 4. The name is deidentified for privacy reasons. The gauge charts show the performance (purple bar) compared to the Expert surgeon (mean and standard deviation indicated as red mark and green area, respectively). In this graph, the representative surgeon gauge starts from zero as the baseline with the goal of reaching to the expert level values denoted by a red bar and green area. The figure for dashboard was created by Python Dash library version 0.43.0: https://plotly.com/dash/.
Figure 5SmartForceps Data Analytics Dashboard in "Data Scientist" view. The current view includes four tabs of "Geospatial Information", “Surgical Force Data”, “Skill Prediction Tool” (a), and “Task Recognition Tool” (b). The figure for dashboard was created by Python Dash library version 0.43.0: https://plotly.com/dash/.
Various features were extracted from the force signal data of each prong in SmartForceps to be used in the monitoring models for surgeon skill assessment.
| Derived force signal features | Description |
|---|---|
| Force duration | Duration of force application in one task segment |
| Force average | Average of force values for one task segment |
| Force max | Maximum of force values in one task segment |
| Force min | Minimum of force values in one task segment |
| Force range | Range of force values across one task segment |
| Force median | Median of force values for one task segment |
| Force SD | Standard deviation of force values in one task segment |
| Force CV | Coefficient of variation of force values in one task segment |
| Force mean CI (0.95) | Confidence interval on the mean with 95% probability |
| Force data skewness | The extent to which the force data distribution deviates from a normal distribution |
| Force data skewness 2SE | The significance of skewness in force data based on dividing by 2 standards errors (significant when > 1) |
| Force data kurtosis | The extent to which the force data distribution is tailed in a normal distribution |
| Force data kurtosis 2SE | The significance of kurtosis in force data based on dividing by 2 standards errors (significant when > 1) |
| Force data normality | Shapiro–Wilk test of normality in force data distribution |
| Force data significance of normality | Significance of Shapiro–Wilk test of normality |
| Force peak value | Peak force value in one task segment |
| Force peak counts | Number of force peaks in one task segment |
| 1st derivative SD | Standard deviation for the first derivative of the force signal in one task segment |
| Force signal flat spots | Maximum run length for each section of force time-series when divided into ten equal-sized intervals |
| Force signal frequency | Dominant time-series harmonics extracted from Fast Fourier Transform (FFT) of force value in one task segment |
| Force cycle length | Average time length of force cycles in one task segment |
| Force signal trend | Force time-series trend in one task segment |
| Force signal fluctuations | Force time-series fluctuation index in one task segment |
| Force signal spikiness | Force time series spikiness index (variance of the leave-one-out variances of the remainder component) in one task segment |
| Force signal linearity | Force time-series linearity index (from Teräsvirta’s nonlinearity test) in one task segment |
| Force signal stability | Force time-series stability index (variance of the means) in one task segment |
| Force signal lumpiness | Force time-series lumpiness index (variance of the variances) in one task segment |
| Force signal curvature | Force time-series curvature index in one task segment (calculated based on the coefficients of an orthogonal quadratic regression) |
| Force signal mean shift | Force time-series largest mean shift between two consecutive windows in one task segment |
| Force signal variance shift | Force time-series largest variance shift between two consecutive windows in one task segment |
| Force signal divergence | Force time-series divergence index in one task segment (largest shift in Kulback-Leibler divergence between two consecutive windows) |
| Force signal stationary index | Force time-series stationary index around a deterministic trend in one task segment (based on Kwiatkowski-Phillips-Schmidt-Shin (KPSS) unit root test with linear trend and lag one) |
| Force signal entropy | Force time-series forcastabilty in one task segment (low values indicate a high signal-to-noise ratio) |
| First autocorrelation minimum | Time of first minimum of the autocorrelation function in force time-series signal from one task segment |
| First autocorrelation zero | Time of first zero crossing of the autocorrelation function in force time-series signal from one task segment |
| Autocorrelation function E1 | First autocorrelation coefficient from force time-series signal in one task segment |
| Autocorrelation function E10 | Sum of the first ten squared autocorrelation coefficients from force time-series signal in one task segment |
This table provides a full list of these time-series and statistical features along with their detailed definitions.
Patient demographics, pathology, location, tumour size and surgeon experience level.
| Mean (SD) age | 54.7 (14.1) |
| Gender (male/female) | 30 M/20 F |
| Disease type (count) | Hemangioblastoma (3) |
| Glioma (10) | |
| Vestibular schwannoma (15) | |
| Meningioma (10) | |
| Cavernous angioma (2) | |
| Trigeminal neuralgia/hemifacial spasm (4) | |
| Chiari malformation (2) | |
| Cervical spondylosis (2) | |
| Arteriovenous malformation (1) | |
| Paraganglioma (1) | |
| Location | Frontal (10), temporal (4), parietal (4), occipital (1), brainstem (2), posterior fossa (27), cervical spine (2) |
| Mean (SD) tumour size (cm, max. diameter R1 × R2 × R3) | 3.3 (1.8) × 2.8 (1.5) × 2.8 (1.5) |
| Surgeon experience year (count) | 30 + (1), fellow (1), PGY 5–6 (4), PGY 3–4 (3), PGY 1–2 (4) |
Figure 6Aggregative data distribution of both Expert (green violin plots) and Novice (purple violin plots) surgeons across the surgical tasks for each time-series extracted feature (Force Range in this sample graph). Detailed statistical information including min, max, median, mean, q1, and q3 are available to view on mouse hover in the original app and the Supplementary Material. Figure created by R Plotly library version 2.0: https://plotly.com/r/.
Workflow diagram of SmartForceps data management and analysis pipeline. Forces of tool-tissue interaction along with de-identified case information is uploaded to a HIPAA-compliant data storage and analytics architecture, i.e., Microsoft Azure. Force data were manually segmented and labelled by listening to the operating room voice recordings where each surgeon name, surgical tasks, and important incidents were properly narrated. Data Pre-processing was performed for noise filtering (Butterworth low-pass filter) and outlier removal (1st and 99th percentiles of either maximum force, minimum force, or task completion time). To generate a holistic information from tool-tissue interaction force profiles, 37 hand-crafted time-series features were extracted in Feature Engineering phase. In Data Analytics phase, two-way ANOVA tests were examined to monitor the representation power of each feature set for different surgeon skill and task categories and a subset of 25 features was selected. The force profiles and selected features were used in Data Analytics Dashboard for performance monitoring and machine learning modelling tools to perform skill prediction and task recognition. The visualization was created in Microsoft PowerPoint version 16.49 with the icons obtained from a Google search: e.g., https://www.iconfinder.com.