PubMed Health⌕ Search
Home/Archive/PubMed

PubMed · 8970391

On modeling longitudinal pulmonary function data.

D Sherrill·G Viegi
View original source ↗Export metadata ↓

Abstract

In this paper we have discussed how random effects can be included in linear models to accommodate the analyses of longitudinal pulmonary function (REM) and categorical (GEE) types of respiratory data. It was suggested that REM analysis be used for continuous observations that are normally distributed or that can be transformed to have near normal distributions and that GEE be used for categorical or non-normally distributed data. Two methods were reviewed, parallel plots and within-subject regression fitting, which can assist in determining the order of random effects to be included. Using a sample data set of longitudinal FEV1 measures, we outlined the steps that should be taken for selecting the within-subject error structure, the order of random effects, elements of the between-subject covariance matrix, and selecting the most important or predictive fixed effects. Lastly, two different types of residual plots were illustrated, conditional and marginal, which can be used to detect outliers and possible trends from underfitting the observed data. All options discussed herein are not available in all the programs currently available for doing the REM and GEE modeling.

Explore related subjects

Computer SimulationForced Expiratory VolumeHumansLongitudinal StudiesModels, BiologicalRandom AllocationRegression AnalysisRespiratory Function TestsJournal ArticleResearch Support, U.S. Gov't, P.H.S.

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

D Sherrill, G Viegi. 1996. On modeling longitudinal pulmonary function data.. https://doi.org/10.1164/ajrccm%2F154.6_pt_2.s217

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related citations

PubMed2025

Generating correlated data for omics simulation.

Simulation of realistic omics data is a key input for benchmarking studies that help users obtain optimal computational pipelines. Omics data involves large numbers of measured features on each sample and these measures are generally correlated with each other. However, simulation too often ignores these correlations, perhaps due to computational and statistical hurdles of doing so. To alleviate this, we describe three approaches for generating omics-scale data with correlated measures which mimic real datasets. These approaches are all based on a Gaussian copula approach with a covariance matrix that decomposes into a diagonal part and a low-rank part. This decomposition allows for extremely efficient simulation, overcoming a hurdle for adoption of past methods. We use these approaches to demonstrate the importance of including correlation in two benchmarking applications. First, we show that variance of results from the popular DESeq2 method increases when dependence is included. Second, we demonstrate that CYCLOPS, a method for inferring circadian time of collection from transcriptomics, improves in performance when given gene-gene dependencies in some circumstances. We provide an R package, dependentsimr, that has efficient implementations of these methods and can generate dependent data with arbitrary marginal distributions, including discrete (binary, ordered categorical, Poisson, negative binomial), continuous (normal), or with an empirical distribution.

Computer Simulation↗
PubMed2017

Addressing current challenges in cancer immunotherapy with mathematical and computational modelling.

The goal of cancer immunotherapy is to boost a patient's immune response to a tumour. Yet, the design of an effective immunotherapy is complicated by various factors, including a potentially immunosuppressive tumour microenvironment, immune-modulating effects of conventional treatments and therapy-related toxicities. These complexities can be incorporated into mathematical and computational models of cancer immunotherapy that can then be used to aid in rational therapy design. In this review, we survey modelling approaches under the umbrella of the major challenges facing immunotherapy development, which encompass tumour classification, optimal treatment scheduling and combination therapy design. Although overlapping, each challenge has presented unique opportunities for modellers to make contributions using analytical and numerical analysis of model outcomes, as well as optimization algorithms. We discuss several examples of models that have grown in complexity as more biological information has become available, showcasing how model development is a dynamic process interlinked with the rapid advances in tumour-immune biology. We conclude the review with recommendations for modellers both with respect to methodology and biological direction that might help keep modellers at the forefront of cancer immunotherapy development.

Computer Simulation↗
PubMed2001

Protein molecular dynamics with the generalized Born/ACE solvent model.

Implicit solvent models are increasingly important for the study of proteins in aqueous solution. Here, the generalized Born (GB) solvent polarization model as implemented in the analytical ACE potential [Schaefer and Karplus (1996) J Phys Chem 100:1578] is used to perform molecular dynamics simulations of two small, homologous proteins: the immunoglobulin-binding domain of streptococcal protein G and the Ras binding domain of Raf. Several model parameterizations are compared through more than 60 ns of simulation. Results are compared with two simpler solvent models-an accessible surface area model and a distant-dependent dielectric model, with finite-difference Poisson calculations, with existing explicit solvent simulations, and with experimental data. The simpler models yield stable but distorted structures. The best GB/ACE implementation uses a set of atomic Voronoi volumes reported recently, obtained by averaging over a large database of crystallographic protein structures. A 20% reduction is applied to the volumes, compensating in an average sense for an excessive de-screening of individual charges inherent in the ACE self-energy and for an undersolvation of dipolar groups inherent in the GB screening function. This GB/ACE parameterization yields stable trajectories on the 0.5-1-ns time scale that deviate moderately (approximately 1.5-2.5 A) from the X-ray structure, reproduce approximately the surface distribution of charged, polar, and hydrophobic groups, and reproduce accurately backbone flexibility as measured by amide NMR-order parameters. Over longer time scales (1.5-3 ns), some of the protein G runs escape from the native energy basin and deviate strongly (3 A) from the native structure. The conformations sampled during the transition out of the native energy basin are overstabilized by the GB/ACE solvation model, as compared with a numerical treatment of the full dielectric continuum model.

Computer Simulation↗