The use of data analysis methods: basic concepts for optometrists

This article explains first, the reasons why a knowledge of statistics is necessary and describes the role that statistics plays in an experimental investigation. Second, the normal distribution is introduced which describes the natural variability shown by many measurements in optometry and vision sciences. Third, the application of the normal distribution to some common statistical problems including how to determine whether an individual observation is a typical member of a population and how to determine the confidence interval for a sample mean is described.

Theory and data for area summation of contrast with and without uncertainty:evidence for a noisy energy model

Contrast sensitivity improves with the area of a sine-wave grating, but why? Here we assess this phenomenon against contemporary models involving spatial summation, probability summation, uncertainty, and stochastic noise. Using a two-interval forced-choice procedure we measured contrast sensitivity for circular patches of sine-wave gratings with various diameters that were blocked or interleaved across trials to produce low and high extrinsic uncertainty, respectively. Summation curves were steep initially, becoming shallower thereafter. For the smaller stimuli, sensitivity was slightly worse for the interleaved design than for the blocked design. Neither area nor blocking affected the slope of the psychometric function. We derived model predictions for noisy mechanisms and extrinsic uncertainty that was either low or high. The contrast transducer was either linear (c1.0) or nonlinear (c2.0), and pooling was either linear or a MAX operation. There was either no intrinsic uncertainty, or it was fixed or proportional to stimulus size. Of these 10 canonical models, only the nonlinear transducer with linear pooling (the noisy energy model) described the main forms of the data for both experimental designs. We also show how a cross-correlator can be modified to fit our results and provide a contemporary presentation of the relation between summation and the slope of the psychometric function.

Short fatigue crack growth behaviour:data analysis effects

Consideration of the influence of test technique and data analysis method is important for data comparison and design purposes. The paper highlights the effects of replication interval, crack growth rate averaging and curve-fitting procedures on crack growth rate results for a Ni-base alloy. It is shown that an upper bound crack growth rate line is not appropriate for use in fatigue design, and that the derivative of a quadratic fit to the a vs N data looks promising. However, this type of averaging, or curve fitting, is not useful in developing an understanding of microstructure/crack tip interactions. For this purpose, simple replica-to-replica growth rate calculations are preferable. © 1988.

Optical data transmission with fluctuating frequency offset and strong phase noise

We examine data transmission during the interval immediately after wavelength switching of a tunable laser and, through simulation, we demonstrate how choice of modulation format can improve the efficacy of an optical burst/packet switched network. © 2013 Optical Society of America.

High-resolution cyclostratigraphy of the Cenomanian-Turonian interval: Paleoceanographic implications, a comparison between deep sea drilling project site 386 and 144

The purpose of this study was to determine the extent to which oceanic anoxic events (OAE's) are recorded in deep-water deposits of the former western Tethyan Sea, by investigating the Cenomanian-Turonian time interval characterized by the worldwide OAE 2 event. The study improved our knowledge of the possible controlling mechanisms that triggered this event at these sites, and furthered our understanding of this global phenomenon. This was examined by high-resolution, multi-proxy analyses of sediments at DSDP Sites 386 and 144, including sedimentology, scanning electron microscopy, stable isotopes, bulk and clay mineralogy, major and trace element geochemistry, biomarkers, and paleontological data. ^ The results provide a better stratigraphic resolution for the Cenomanian-Turonian, which allowed for more precise determination of chronologic boundaries, sedimentation rates at DSDP Site 386, and a more accurate calculation of the frequency of the cycles recorded in the sequence, which fall predominantly within the precession (∼23 kyr) and short eccentricity (∼100 kyr) ranges. The combined proxies allow assessment of the correlation of δ13Corg, and major and trace elements with the predominance of cyanobacteria. These organisms were the main producers of the organic matter during the dysoxic and euxinic conditions of OAE 2 at DSDP Site 386. A huge amount of microcrystalline quartz of eolian origin is also associated with OAE 2. The geochemical proxies further provide evidence that OAE 2 was linked to increased volcanism outside the deep water of the proto-Atlantic Ocean. The clays in the Turonian sediments are terrigenous and derived predominantly from eolian transport. Comparing DSDP Site 386 and 144 with stratotype sections, the δ13C org and TOC data indicate that OAE 2 seems diachronous throughout the proto-Atlantic Ocean. ^ This study concludes that the development of anoxic conditions in the deep water of the Atlantic during the latest Cenomanian-Turonian resulted from a combination of factors related to local oceanic setting and mitigated by global tectonism and climate. The data provide a more comprehensive view of the interacting factors that led to sustained high productivity of the cyanobacteria and photosynthetic protists that produced organic-carbon-rich deposits in the world's oceans. ^

Electrocardiogram breathing: A study of electrocardiogram derived respiration plus multiresolution wavelet analysis of heartbeat interval variability in sleep study patients

This dissertation presents a unique research opportunity by using recordings which provide electrocardiogram (ECG) plus a reference breathing signal (RBS). ECG derived breathing (EDR) is measured and correlated against RBS. Standard deviations of multiresolution wavelet analysis coefficients (SDMW) are obtained from heart rate and classified using RBS. Prior works by others used select patients for sleep apnea scoring with EDR but no RBS. Another prior work classified select heart disease patients with SDMW but no RBS. This study used randomly chosen sleep disorder patient recordings; central and obstructive apneas, with and without heart disease.^ Implementation required creating an application because existing systems were limited in power and scope. A review survey was created to choose a development environment. The survey is presented as a learning tool and teaching resource. Development objectives were rapid development using limited resources (manpower and money). Open Source resources were used exclusively for implementation. ^ Results show: (1) Three groups of patients exist in the study. Grouping RBS correlations shows a response with either ECG interval or amplitude variation. A third group exists where neither ECG intervals nor amplitude variation correlate with breathing. (2) Previous work done by other groups analyzed SDMW. Similar results were found in this study but some subjects had higher SDMW, attributed to a large number of apneas, arousals and/or disconnects. SDMW does not need RBS to show apneic conditions exist within ECG recordings. (3) Results in this study support the assertion that autonomic nervous system variation was measured with SDMW. Measurements using RBS are not corrupted due to breathing even though respiration overlaps the same frequency band.^ Overall, this work becomes an Open Source resource which can be reused, modified and/or expanded. It might fast track additional research. In the future the system could also be used for public domain data. Prerecorded data exist in similar formats in public databases which could provide additional research opportunities. ^

Improved estimation of postmortem interval with Cardiac Troponin T and improved analytical methods

Accurate knowledge of the time since death, or postmortem interval (PMI), has enormous legal, criminological, and psychological impact. In this study, an investigation was made to determine whether the relationship between the degradation of the human cardiac structure protein Cardiac Troponin T and PMI could be used as an indicator of time since death, thus providing a rapid, high resolution, sensitive, and automated methodology for the determination of PMI. ^ The use of Cardiac Troponin T (cTnT), a protein found in heart tissue, as a selective marker for cardiac muscle damage has shown great promise in the determination of PMI. An optimized conventional immunoassay method was developed to quantify intact and fragmented cTnT. A small sample of cardiac tissue, which is less affected than other tissues by external factors, was taken, homogenized, extracted with magnetic microparticles, separated by SDS-PAGE, and visualized with Western blot by probing with monoclonal antibody against cTnT. This step was followed by labeling and available scanners. This conventional immunoassay provides a proper detection and quantitation of cTnT protein in cardiac tissue as a complex matrix; however, this method does not provide the analyst with immediate results. Therefore, a competitive separation method using capillary electrophoresis with laser-induced fluorescence (CE-LIF) was developed to study the interaction between human cTnT protein and monoclonal anti-TroponinT antibody. ^ Analysis of the results revealed a linear relationship between the percent of degraded cTnT and the log of the PMI, indicating that intact cTnT could be detected in human heart tissue up to 10 days postmortem at room temperature and beyond two weeks at 4C. The data presented demonstrates that this technique can provide an extended time range during which PMI can be more accurately estimated as compared to currently used methods. The data demonstrates that this technique represents a major advance in time of death determination through a fast and reliable, semi-quantitative measurement of a biochemical marker from an organ protected from outside factors. ^

Compilation of downward flux observations from sediment trap deployments in the Atlantic Ocean - Contribution to EURO-BASIN's Data integration

We provide a compilation of downward fluxes (total mass, POC, PON, BSiO2, CaCO3, PIC and lithogenic/terrigenous fluxes) from over 6000 sediment trap measurements distributed in the Atlantic Ocean, from 30 degree North to 49 degree South, and covering the period 1982-2011. Data from the Mediterranean Sea are also included. Data were compiled from different sources: data repositories (BCO-DMO, PANGAEA), time series sites (BATS, CARIACO), published scientific papers and/or personal communications from PI's. All sources are specifed in the data set. Data from the World Ocean Atlas 2009 were extracted to provide each flux observation with contextual environmental data, such as temperature, salinity, oxygen (concentration, AOU and percentage saturation), nitrate, phosphate and silicate.

Plant leaf wax (n-alkane) data from sediment core MD96-2048

In this study we investigate Pleistocene vegetation and climate change in southern East Africa by examining plant leaf waxes in a marine sediment core that receives terrestrial runoff from the Limpopo River. The plant leaf wax records are compared to a multi-proxy sea surface temperature (SST) record and pollen assemblage data from the same site. We find that Indian Ocean SST variability, driven by high-latitude obliquity, exerted a strong control on the vegetation of southern East Africa during the past 800,000 yr. Interglacial periods were characterized by relatively wetter and warmer conditions, increased contributions of C3 vegetation, and higher SST, whereas glacial periods were marked by cooler and arid conditions, increased contributions of C4 vegetation, and lower SST. We find that Marine Isotope Stages (MIS) 5e, 11c, 15e and 7a-7c are strongly expressed in the plant leaf wax records but MIS 7e is absent while MIS 9 is rather weak. Our plant leaf wax records also record the climate transition associated with the Mid-Brunhes Event (MBE) suggesting that the pre-MBE interval (430-800 ka) was characterized by higher inputs from grasses in comparison to relatively higher inputs from trees in the post-MBE interval (430 to 0 ka). Differences in vegetation and SST of southern East Africa between the pre- and post-MBE intervals appear to be related to shifts in the location of the Subtropical Front. Comparison with vegetation records from tropical East Africa indicates that the vegetation of southern East Africa, while exhibiting glacial-interglacial variability and notable differences between the pre- and post-MBE portions of the record, likely did not experience such dramatic extremes as occurred to the north at Lake Malawi.

Varve data from a sediment profile of Lake Woserin, Germany

Time-series of varve properties and geochemistry were established from varved sediments of Lake Woserin (north-eastern Germany) covering the recent period AD 2010-1923 and the Mid-Holocene time-window 6400-4950 varve years before present (vyr BP) using microfacies analyses, X-ray fluorescence scanning (µ-XRF), microscopic varve chronology and 14C dating. The microscopic varve chronology was compared to a macroscopic varve chronology for the same sediment interval. Calcite layer thickness during the recent period is significantly correlated to increases in local annual precipitation (r=0.46, p=0.03) and reduced air-pressure (r=-0.72, p<0.0001). Meteorologically consistent with enhanced precipitation at Lake Woserin, a composite 500 hPa anomaly map for years with >1 standard deviation calcite layer thickness depicts a negative wave train air-pressure anomaly centred over southern Europe, with north-eastern Germany at its northern frontal zone. Three centennial-scale intervals of thicker calcite layers around the Mid-Holocene periods 6200-5900, 5750-5400 and 5300-4950 vyr BP might reflect humid conditions favouring calcite precipitation through the transport of Ca2+ ions into Lake Woserin, synchronous to wetter conditions in Europe. Calcite layer thickness oscillations of about 88 and 208 years resemble the solar Gleissberg and Suess cycles suggesting that the recorded hydroclimate changes in north-eastern Germany are modified by solar influences on synoptic-scale atmospheric circulation. However, parts of the periods of thicker calcite layers around 5750-5400 and 5200 vyr BP also coincide to enhanced human catchment activity at Lake Woserin. Therefore, calcite precipitation during these time-windows might have further been favored by anthropogenic deforestation mobilizing Ca2+ ions and/or lake eutrophication.

Dive depth profile and at surface behaviour data of emperor penguins from Pointe Géologie, Adélie Land, Antarctica, from expedition DDU 2005

Adult male and female emperor penguins (Aptenodytes forsteri) were fitted with satellite transmitters at Pointe-Géologie (Adélie Land), Dumont d'Urville Sea coast, in November 2005. Nine of 30 data sets were selected for analyses to investigate the penguins' diving behaviour at high resolution (doi:10.1594/PANGAEA.633708, doi:10.1594/PANGAEA.633709, doi:10.1594/PANGAEA.633710, doi:10.1594/PANGAEA.633711). The profiles are in synchrony with foraging trips of the birds during austral spring (doi:10.1594/PANGAEA.472171, doi:10.1594/PANGAEA.472173, doi:10.1594/PANGAEA.472164, doi:10.1594/PANGAEA.472160, doi:10.1594/PANGAEA.472161). Corresponding high resolution winter data (n = 5; archived elsewhere) were provided by A. Ancel, Centre d'Ecologie et Physiologie Energétiques, CNRS, Strasbourg, France. Air-breathing divers tend to increase their overall dive duration with increasing dive depth. In most penguin species, this occurs due to increasing transit (descent and ascent) durations but also because the duration of the bottom phase of the dive increases with increasing depth. We interpreted the efficiency with which emperor penguins can exploit different diving depths by analysing dive depth profile data of nine birds studied during the early and late chick-rearing period in Adélie Land, Antarctica. Another eight datasets of dive depth and duration frequency recordings (doi:10.1594/PANGAEA.472150, doi:10.1594/PANGAEA.472152, doi:10.1594/PANGAEA.472154, doi:10.1594/PANGAEA.472155, doi:10.1594/PANGAEA.472142, doi:10.1594/PANGAEA.472144, doi:10.1594/PANGAEA.472146, doi:10.1594/PANGAEA.472147), which backup the analysed high resolution depth profile data, and dive depth and duration frequency recordings of another bird (doi:10.1594/PANGAEA.472156, doi:10.1594/PANGAEA.472148) did not match the requirement of high resolution for analyses. Eleven additional data sets provide information on the overall foraging distribution of emperor penguins during the period analysed (doi:10.1594/PANGAEA.472157, doi:10.1594/PANGAEA.472158, doi:10.1594/PANGAEA.472162, doi:10.1594/PANGAEA.472163, doi:10.1594/PANGAEA.472166, doi:10.1594/PANGAEA.472167, doi:10.1594/PANGAEA.472168, doi:10.1594/PANGAEA.472170, doi:10.1594/PANGAEA.472172, doi:10.1594/PANGAEA.472174, doi:10.1594/PANGAEA.472175).

Dinocyst-derived sea-surface characteristics, palynomorph abundances, ice-rafted debris, and planktonic foraminifera abundances and d18O of the 35-41 ka BP interval of core MD99-2285

The last glacial millennial climatic events (i.e. Dansgaard-Oeschger and Heinrich events) constitute outstanding case studies of coupled atmosphere-ocean-cryosphere interactions. Here, we investigate the evolution of sea-surface and subsurface conditions, in terms of temperature, salinity and sea ice cover, at very high-resolution (mean resolution between 55 and 155 years depending on proxies) during the 35-41 ka cal BP interval covering three Dansgaard-Oeschger cycles and including Heinrich event 4, in a new unpublished marine record, i.e. the MD99-2285 core (62.69°N; -3.57s°E). We use a large panel of complementary tools, which notably includes dinocyst-derived sea-ice cover duration quantifications. The high temporal resolution and multiproxy approach of this work allows us to identify the sequence of processes and to assess ocean-cryosphere interactions occurring during these periodic ice-sheet collapse events. Our results evidence a paradoxical hydrological scheme where (i) Greenland interstadials are marked by a homogeneous and cold upper water column, with intensive winter sea ice formation and summer sea ice melting, and (ii) Greenland and Heinrich stadials are characterized by a very warm and low saline surface layer with iceberg calving and reduced sea ice formation, separated by a strong halocline from a less warm and saltier subsurface layer. Our work also suggests that this stadial surface/subsurface warming started before massive iceberg release, in relation with warm Atlantic water advection. These findings thus support the theory that upper ocean warming might have triggered European ice-sheet destabilization. Besides, previous paleoceanographic studies conducted along the Atlantic inflow pathways close to the edge of European ice-sheets suggest that such a feature might have occurred in this whole area. Nonetheless, additional high resolution paleoreconstructions are required to confirm such a regional scheme.

Correction to “Reminder and 2AFC tasks provide similar estimates of the difference limen: A reanalysis of data from Lapid, Ulrich, and Rammsayer (2008) and a discussion of Ulrich and Vorberg (2009)”

We recently published an article (García-Pérez & Alcalá- Quintana, 2010) reanalyzing data presented by Lapid, Ulrich, and Rammsayer (2008) and discussing a theoretical argument developed by Ulrich and Vorberg (2009). The purpose of this note is to correct an error in our study that has some theoretical importance, although it does not affect the conclusion that was raised. The error lies in that asymptote parameters reflecting lapses or finger errors should not enter the constraint relating the psychometric functions that describe performance when the comparison stimulus in a two-alternative forced choice (2AFC) discrimination task is presented in the first or second interval.

Contributions of Dorsal/Ventral Hippocampus and Dorsolateral/Dorsomedial Striatum to Interval Timing

Humans and animals have remarkable capabilities in keeping time and using time as a guide to orient their learning and decision making. Psychophysical models of timing and time perception have been proposed for decades and have received behavioral, anatomical and pharmacological data support. However, despite numerous studies that aimed at delineating the neural underpinnings of interval timing, a complete picture of the neurobiological network of timing in the seconds-to-minutes range remains elusive. Based on classical interval timing protocols and proposing a Timing, Immersive Memory and Emotional Regulation (TIMER) test battery, the author investigates the contributions of the dorsal and ventral hippocampus as well as the dorsolateral and the dorsomedial striatum to interval timing by comparing timing performances in mice after they received cytotoxic lesions in the corresponding brain regions. On the other hand, a timing-based theoretical framework for the emergence of conscious experience that is closely related to the function of the claustrum is proposed so as to serve both biological guidance and the research and evolution of “strong” artificial intelligence. Finally, a new “Double Saturation Model of Interval Timing” that integrates the direct- and indirect- pathways of striatum is proposed to explain the set of empirical findings.

Data-Driven Motion Detection and Characterization in PET Brain Scans Using List Mode

Head motion during a Positron Emission Tomography (PET) brain scan can considerably degrade image quality. External motion-tracking devices have proven successful in minimizing this effect, but the associated time, maintenance, and workflow changes inhibit their widespread clinical use. List-mode PET acquisition allows for the retroactive analysis of coincidence events on any time scale throughout a scan, and therefore potentially offers a data-driven motion detection and characterization technique. An algorithm was developed to parse list-mode data, divide the full acquisition into short scan intervals, and calculate the line-of-response (LOR) midpoint average for each interval. These LOR midpoint averages, known as “radioactivity centroids,” were presumed to represent the center of the radioactivity distribution in the scanner, and it was thought that changes in this metric over time would correspond to intra-scan motion.

Several scans were taken of the 3D Hoffman brain phantom on a GE Discovery IQ PET/CT scanner to test the ability of the radioactivity to indicate intra-scan motion. Each scan incrementally surveyed motion in a different degree of freedom (2 translational and 2 rotational). The radioactivity centroids calculated from these scans correlated linearly to phantom positions/orientations. Centroid measurements over 1-second intervals performed on scans with ~1mCi of activity in the center of the field of view had standard deviations of 0.026 cm in the x- and y-dimensions and 0.020 cm in the z-dimension, which demonstrates high precision and repeatability in this metric. Radioactivity centroids are thus shown to successfully represent discrete motions on the submillimeter scale. It is also shown that while the radioactivity centroid can precisely indicate the amount of motion during an acquisition, it fails to distinguish what type of motion occurred.