AN EXACT CONFIDENCE INTERVAL FOR THE RATIO OF MEANS USING NONLINEAR REGRESSION ThponRoy Boehringer Ingelheim PharmaceuticaIs, Inc.
|
|
- Bartholomew Reynolds
- 6 years ago
- Views:
Transcription
1 AN EXACT CONFIDENCE INTERVAL FOR THE RATIO OF MEANS USING NONLINEAR REGRESSION ThponRoy Boehringer Ingelheim PharmaceuticaIs, Inc. ABSTRACT In several statistical applications, the ratio of means is of primary interest. The construction of confidence intervals for such ratios can be done using various techniques, most of which are less than satisfactory. In this paper, a nonlinear regression model will be used to calculate an exact confidence interval for the ratio of means. This result will be compared to a commonly used approximate technique, and to the empirical~istnbution based bootstrap method. A pharmacokinetic example involving the bioequivalence of formulations will illustrate the methodology, and a SAS"" program to perform the calculations will be presented. THE REGRESSION MODEL AND CONFIDENCE INTERVAL The nonlinear regression model can be expressed as: Yi = j(x;, ~ ) + Ei where Yi are responses (i = 1,..., n) taken at design levels Xi (i = 1,...,11), and ~ is ap-element parameter vector. i (i = 1,...,11) are the errors, assumed to be independent, identically distributed normal random variables with mean zero and variance u'. The regression model for the two-sample problem in which the ratio of means is involved can be formulated as: Yi = 9,xi + 9,9,(1 - Xi), where X. is an indicator variable that takes the values o and 1 for means 1 and 2, with n, and n, as the group sample sizes, respectively. This is the Fieller-Creasy model (see Wallace, 1980; and Cook and Witmer, 1985 for a complete discussion). The ratio of means, 9" is of main interest. Following the derivation in Cook and Witmer (1985), the exact marginal confidence region for 9, is: E.(9,) = {9,1 8l (O,-9,t/(1 + Ii, ) ~ c,}, where c, = 202 X 2 (cr; 1)ln. The regression MSE is used to estimate (9,) will be an interval only if o I > c,; otherwise, 0.(9,) will be the complement of an interval or the entire real line. Thus, the exact interval for the estimate of the ratio, 0" is: [ 0, ± {r, (1- r,) +r, iii y'] 1(1- r,), where r,= c,/~ < 1. Aninterestingfeatureofthe Fieller-Creasy model is that though the intrinsic curvature, curvature, f"'/= 0, and the parameter effects r' = ujil,{(o; +n,/n, )U' + IO,IJ In, 10,1, is quite small, indicating low curvature, the linear approximation confidence regions are in poor agreement with the exact regions. This is in contrast to the Bates-Watts criterion guidelines (Bates and Watts, 1980). EXAMPLE An approximate confidence interval for the ratio of means can be represented as: y... ± t(cr/2;n-2)j2'mse / n Y... in the application of testing bioequivalence of pharmaceutical formulations, (std. and test). The MSE is taken from an ANOVA using subject, formulation, and other relevant variables as factors, using an n-subject crossover design. The formulation characteristic of interest could be area under the concentration-time curve (AVe) or a related measure. This type of confidence interval is widely used in bioequivalence testing despite its obvious limitations.' This will be referred to as the conventional method. Effons have been made to develop improved confidence intervals, (see Mandallaz and Mau, 1981; Steinijans and Diletti, 1983; and Locke, 1984). An illustrative example comparing different methods of confidence interval construction is provided in Thble 1 (from Roy, in press). (The data for this example is in the Appendix.) Included is the bootstrap confidence interval which reflects the empirical parameter distribution. Efron (1982) descnbes the bootstrap method and Roy and MacGregor (1991) have applied it to bioequivalence problems. Since in bioequivalence problems, the assumptions of independence and normality may be untenable, the bootstrap method is an attractive alternative to conventional methods. 859
2 The bootstrap technique incorporates distnbutional features, correlation structure, and factor effects when applied appropriately. Since a crossover design was used, the (std., test) pairs of response values were resampled. This ensures that both formulation balance and correlation features are accounted for in the analysis. The distnbution of individual ratios (or their means or medians) can also be characterized using this procedure. For parallel group designs, the resampling should be stratified by formulation. Stratified resampling can be used for any design with multiple factors to guarantee proportionally representative samples. From Thble I it is evident that the Fieller-Creasy and bootstrap intervals are inclose agreement, and the conventional procedure interval is quite different. Since bioequivalence is generally concluded when the 90% confidence interval for the ratio is entirely within the (0.80, 1.20) range «0.80, 1.25) for 10g~tranSformed data), the conventional procedure would indicate bioequivalence while the Fieller~reasy and bootstrap methods do not. 10 examine the validity and accuracy of the bootstrap, bioequivalence study (AUC) data was simulated following the normal, lognormal, and uniform distnbutions. For each distnbution, a two-period crossover (n = 24) design to compare a test formulation with a standard was simulated, incorporating the necessary correlation structure. TABLE I COMPARISON OF CONFIDENCE INTERVALS 90% Confidence Intervals RatiO of Conven- Fieller- Bootstrap, AUC Means tional Creasy B=4, (0.89, (0.85, (0.86, 1.14) 1.21) 1.21) This procedure was repeated 200 times for each distnbution. B = 500 bootstrap replications were performed for each simulated sample (for the lognormal simulations, the 95% confidence intervals are also presented, as this is customary in bioequivalence studies). The conventional procedure for the lognormal simulation data was the usual confidence interval on the difference of the log-transformed means, exponentiating the endpoints to return to the original scale. This construction would have given a correct reference interval if the means were uncorrelated. The results, from Roy and MacGregor (1991), are summarized in Thbles II-Iv. From Thbles II and III, it is evident that the bootstrap intervals are considerably narrower than the conventional intervals. From Thble rv; it is clear that coverage probabilities (the proportion of confidence intervals containing the true values) for the bootstrap analyses are far closer to the nominal levels than those for the conventional intervals, especially for the. lognormal and uniform cases. Theoretically, the assumptions of normality and independence must be tenable for the validity of the exact Fieller-Creasy interval. In Thble V, the stability of the bootstrap intervals is demonstrated for various values of B, the number of bootstrap replications used for normally distnbuted, simulated data. These results indicate that the bootstrap method can be used when the situation defies conventional analysis. In fact, if confidence intervals on the ratio of medians or inferences on individual ratios are desired, it may be the only applicable method. The Fieller-Creasy method is recommended if departures from normality and independence are not excessive. Factor effects can easily be incorporated into the Fieller-Creasy model by the addition of indicator variables, if desired. COMPUTER PROGRAM In the APPendix, a SAS program is given that will calculate the exact Fieller-Creasy confidence intervals. It is written to produce 90% intervals, though this can be changed to any desired level by changing the argument to CINV in the program. The example shown in Thble r is used to demonstrate the input and output for the program. This Version 6.07 program was written for the CMS environment on an IBM 9121 time-sharing system. It should be easily adaptable to any operating system and platform. DISCUSSION In this paper, the computation of an exact confidence interval for the ratio of means has been presented. Comparison of the exact method to a conventional procedure and to the bootstrap technique was performed, and the conditions under which the exact method and bootstrap. technique are valid were examined. Under standard assumptions, the exact method is valid from theoretical considerations. When the usual assumptions do not hold or are difficult to verify, additional factors are present, or other unusual conditions exist, the bootstrap can be used to obtain accurate and reliable confidence intervals. The conventional procedure appears less than satisfactory and should be used with caution. 860
3 TABLE II SIMULATION STUDY RESULTS, 200 TRIALS COMPARISON OF CONFIDENCE INTERVALS FOR RATIOS OF MEANS Conventional Bootstrap, B = % Confidence 95% Confidence 90% Confidence 95% Confidence Distribution Intervals Intervals Intervals Intervals Lower Upper Lower Upper Lower Upper Lower Upper Umlt Umlt Umlt Umlt Umlt Umlt Umlt Umlt Mean Mean Mean Mean Mean Mean Mean Mean (S.D.) (S.D.) (S.D.) (S.D.) (S.D.) (S.D.) (S.D.) (S.D.) Normal ( ) ( ) ( ) ( ) Lognormal ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) Uniform ( ) (0.01S01) ( ) ( ) TABLE III SIMULATION STUDY RESULTS, 200 TRIALS COMPARISON OF CONFIDENCE INTERVAL LENGTHS FOR RATIOS OF MEANS.'.. Conventional Bootstrap, B = % Coi.fldence 95% Confidence 90% Confidence 95% Confidence Distribution Interval Interval Interval Interval Mean Length (S.D.) Mean Length (S.D.) Mean Length (S.D.) Mean Length (S.D.). Normal ( ) ( ) Lognormal , ( ) ( ) ( ) ( ) Uniform (0.0148) ( ) TABLE IV SIMULATION STUDY RESULTS, 200 TRIALS CONFIDENCE INTERVAL COVERAGIO PROBABILITIES FOR RATIOS OF MEANS Conventional Bootstrap, B Distribution Nominal = Normal Lognormal Uniform
4 TABLE V COMPARISON OF RESULTS AT DIFFERENT NUMBERS OF BOOTSTRAP REPLICATIONS: POINT ESTIMATES AND 90% BOOTSTRAP CONADENCE INTERVALS Normal DIStribution B=5OO B=1,000 B=2,000 B=3,000 B=4,000 Ratio of Means (0.8258, ) (0.8254,0.9211) (0.8241, ) (0.8241,0.9199) (0.8241, ) Ratio of Medians (0.7829, ) (0.7829, ) (0.7829, ) (0.7829, ) (0.7829, ) Entries are estimate (5th percentile, 95th percentile) where estimate is the median of the bootstrap distribution. ACknowledgment: The author wishes to thank the BIPI Document Processing Group for their expert word processing. REFERENCES Bates, D.M. and Watts, D.G. (1980). Relative curvature measures of nonlinearity. Journal of the Royal Statistical Society, Ser. B 42, Cook, R.D. and Witmer, J.A. (1985). A note on parameter-effects curvature. Journal of the American Statistical Association 80, Efron, B. (1982). The Jackknife. the Bootstrap, and Other Resampling Plans. Philadelphia: Society for Industrial and Applied Mathematics. Locke, C.S. (1984). An exact confidence interval from untransformed data for the ratio of two formulation means. Journal of Phannacokinetics and Biophannaceutics 12, Mandallaz, D. and Mall, J. (1981). Comparison of different methods for decision-making in bioequivalence assessment. Biometrics 37, Roy, T. and MacGregor, T.R. (1991). Bootstrap confidence intervals for bioequivalence ratios. Paper presented at the 1991 American Association of Pharmaceutical Scientists meeting. Pharmaceutical Research 8, S-63. Roy, T. An exact confidence interval for the ratio of means usfug regression methods: Journal of Statistical Computation and Simulation, in press. Steinijans, v.w. and Diletti, E. (1983). Statistical analysis of bioavaiiability studies:. parametric and nonparametric confidence intervals. European Journal of Clinical Phannacology 24, Wallace, D.L. (1980). The Behrens-Fisher and FieBer-CreasY problems. In RA. FISher: An Appreciation (S.E. Fienberg and D.V. Hinkley, eds.) New York: Springer-Verlag, pp ; SAS is a registered trademark or trademark of SAS Institute, Inc. in the U.S.A. and other countries. indicates U.S.A. registration. Other brand and product names are registered trademarks \>r trademarks of their respective companies. 862
5 APPENDIX: PROGRAM AND OUTPUT.. ", EXACT CONFIDENCE INTERVAL CALCULATION FOR THE RATIO. OF MEANS USING THE FIELLER-CREASY MODEL..:..,, SEE COOK AND WITMER (1985), JABA, 80, AUTHOR: TAPON ROY DATE: 9/10/ ~....~... ': :.:, ; OPTIONS NODATE NONUMBER: READ INPUT DATA IN ; Y VALUES ARE FOR COMPUTATION OF THE TWO ~ANS '00: '00 THE VALUE OF X IS 0 FOR MEAN 1 AND 1 FOR MEAN 2 0": DATA A: SET A.EQVCONF: CALL SYlIPUT( 'N', _N_): PROC PRINT: TITLE 'SAMPLE DATA': PRoe SORT: BY X: PROC MEANs: VAR Y: BY X' TITrE 'GROUP MEANS': '00 IN THE PROCEDURE BELOW, SET THE INITIAL VALUE OF B1 "': TO BE - THE VALUE OF MEAN 2, AND THE INITIAL VALUE OF B2 TO BE - THE RATIO OF MEANl TO MBAH2 ; PROC NLIN BEST=Z5 PLOT METHOD-MARQUARDT OUTEST=POUT: PARMS B1-1Z B2= 1: MODEL Y- (B1 X) + (B1'BZ"(1-X»: DER.BI- X + (BZO(l-X»: DER.B2= B1 (1-X): OUTPUT OUT=D1 P=PD1 R=RD1: TITLE 'NONLINEAR REGRESSION MODEL': TITLE 'PREDICTED VALUES AND RESIDUALS FOR METHOD-MARQUARDT': DATA POUT: SET POUT: IF _fype_='final'; DATA AA1: SET POUT: KEEP _TYPE SSE_ B1 B2: DATA AA1: SET AA1: N=&N: DFE- N-2: C1= (2 (_SSE IDFE) CINV(0.90,l»/N: FOR 90 PCT C.I.: CUTOFF= (Bl T2): R1= Cl/(B1"2): B218= (B2 - SQRT(R1'(1-Rl) (Rl'(82 '2»»/(1-Rl): B2UE= (B2 SQRT(Rl (l-rl) (R1'(82002»»/(1-R1): PROC PRINT: TITLE1 'PARAMETER ESTIMATES AND THE EXACT CONFIDENCE INTERVAL': UTLE2 'FOR THE RATIO OF MEANS':. TITLE3 '(B2UE AND B2LE ARE ONLY VALID IF CUTOFF) C1)':. SAMPLE DATA abs y X
6 VARIABL.E... STANDARD GROUP JiEANS MINIMUM MUIMtJM DEVIATION VALUE VALUE -._. XaO _ STD ERROl OF MEAN.'"' VA.IAHc:E C.V Y S5405UI ; $ X-I y ' & & NOTE: CONVERGENCE CRITERION MET. ITERATION o NDNLINEAR REGRESSION MODEL NON-LINEAR LEAST SQUARES ITERATIVE PHASE DEPENDENT VARIABLE: Y METHOD: MARQUARDT B RESIDUAL SS NONLINEAR REGRESSION MODEL NON-LINEAR LEAST SQUARES SUMMARY STATISTICS DEPENDENT VARIABLE Y SOURCE REGRESSION RESIDUAL OF SUM OF SQUARES MEAN SQUARE UNCORRECTED (CORRECTED TOTP,L TOTAL] PARAMETER Bl B2 ESTIMATE ASYMPTOTIC ASYMPTOTIC 95 % STD. ERROR CONFIDENCE INTERVAL LOWER UPPER USE EXACT C.I. BELOW --- ASYMPTOTIC CORRELATION MATRIX OF THE PARAMETERS CORR B1 B2 B1 B ! l,,, f OBS., PARAMETER ESfIMATES AND THE EXACT CONFIDENCE INTERVAL FOR THE RATIO Dr MEANS (B2UE AND 82LE ARE ONLY VALID IF CUTOFF) ell.2 N DFE C1 CUTOFF R1 12, o EXACT CONFIDENCE INTERVAL B2L.E B2UE
Lecture 1 Linear Regression with One Predictor Variable.p2
Lecture Linear Regression with One Predictor Variablep - Basics - Meaning of regression parameters p - β - the slope of the regression line -it indicates the change in mean of the probability distn of
More informationEXST Regression Techniques Page 1. We can also test the hypothesis H :" œ 0 versus H :"
EXST704 - Regression Techniques Page 1 Using F tests instead of t-tests We can also test the hypothesis H :" œ 0 versus H :" Á 0 with an F test.! " " " F œ MSRegression MSError This test is mathematically
More informationOverview Scatter Plot Example
Overview Topic 22 - Linear Regression and Correlation STAT 5 Professor Bruce Craig Consider one population but two variables For each sampling unit observe X and Y Assume linear relationship between variables
More informationPermutation Tests. Noa Haas Statistics M.Sc. Seminar, Spring 2017 Bootstrap and Resampling Methods
Permutation Tests Noa Haas Statistics M.Sc. Seminar, Spring 2017 Bootstrap and Resampling Methods The Two-Sample Problem We observe two independent random samples: F z = z 1, z 2,, z n independently of
More informationLeast Absolute Value vs. Least Squares Estimation and Inference Procedures in Regression Models with Asymmetric Error Distributions
Journal of Modern Applied Statistical Methods Volume 8 Issue 1 Article 13 5-1-2009 Least Absolute Value vs. Least Squares Estimation and Inference Procedures in Regression Models with Asymmetric Error
More information6 Single Sample Methods for a Location Parameter
6 Single Sample Methods for a Location Parameter If there are serious departures from parametric test assumptions (e.g., normality or symmetry), nonparametric tests on a measure of central tendency (usually
More informationSplineLinear.doc 1 # 9 Last save: Saturday, 9. December 2006
SplineLinear.doc 1 # 9 Problem:... 2 Objective... 2 Reformulate... 2 Wording... 2 Simulating an example... 3 SPSS 13... 4 Substituting the indicator function... 4 SPSS-Syntax... 4 Remark... 4 Result...
More informationGraphical Presentation of a Nonparametric Regression with Bootstrapped Confidence Intervals
Graphical Presentation of a Nonparametric Regression with Bootstrapped Confidence Intervals Mark Nicolich & Gail Jorgensen Exxon Biomedical Science, Inc., East Millstone, NJ INTRODUCTION Parametric regression
More informationThe bootstrap. Patrick Breheny. December 6. The empirical distribution function The bootstrap
Patrick Breheny December 6 Patrick Breheny BST 764: Applied Statistical Modeling 1/21 The empirical distribution function Suppose X F, where F (x) = Pr(X x) is a distribution function, and we wish to estimate
More informationFinite Population Correction Methods
Finite Population Correction Methods Moses Obiri May 5, 2017 Contents 1 Introduction 1 2 Normal-based Confidence Interval 2 3 Bootstrap Confidence Interval 3 4 Finite Population Bootstrap Sampling 5 4.1
More information4 Resampling Methods: The Bootstrap
4 Resampling Methods: The Bootstrap Situation: Let x 1, x 2,..., x n be a SRS of size n taken from a distribution that is unknown. Let θ be a parameter of interest associated with this distribution and
More informationConfidence Interval for the mean response
Week 3: Prediction and Confidence Intervals at specified x. Testing lack of fit with replicates at some x's. Inference for the correlation. Introduction to regression with several explanatory variables.
More informationChapter 6. Estimation of Confidence Intervals for Nodal Maximum Power Consumption per Customer
Chapter 6 Estimation of Confidence Intervals for Nodal Maximum Power Consumption per Customer The aim of this chapter is to calculate confidence intervals for the maximum power consumption per customer
More informationLecture 3: Inference in SLR
Lecture 3: Inference in SLR STAT 51 Spring 011 Background Reading KNNL:.1.6 3-1 Topic Overview This topic will cover: Review of hypothesis testing Inference about 1 Inference about 0 Confidence Intervals
More informationTesting for bioequivalence of highly variable drugs from TR-RT crossover designs with heterogeneous residual variances
Testing for bioequivalence of highly variable drugs from T-T crossover designs with heterogeneous residual variances Christopher I. Vahl, PhD Department of Statistics Kansas State University Qing Kang,
More informationComparison of Different Methods of Sample Size Re-estimation for Therapeutic Equivalence (TE) Studies Protecting the Overall Type 1 Error
Comparison of Different Methods of Sample Size Re-estimation for Therapeutic Equivalence (TE) Studies Protecting the Overall Type 1 Error by Diane Potvin Outline 1. Therapeutic Equivalence Designs 2. Objectives
More informationESTIMATING RHEOLOGICAL PROPERTIES OF YOGURT USING DIFFERENT VERSIONS OF THE FREUNDLICH MODEL AND DESIGN MATRICES
Libraries Conference on Applied Statistics in Agriculture 2004-16th Annual Conference Proceedings ESTIMATING RHEOLOGICAL PROPERTIES OF YOGURT USING DIFFERENT VERSIONS OF THE FREUNDLICH MODEL AND DESIGN
More informationFULL LIKELIHOOD INFERENCES IN THE COX MODEL
October 20, 2007 FULL LIKELIHOOD INFERENCES IN THE COX MODEL BY JIAN-JIAN REN 1 AND MAI ZHOU 2 University of Central Florida and University of Kentucky Abstract We use the empirical likelihood approach
More information1) Answer the following questions as true (T) or false (F) by circling the appropriate letter.
1) Answer the following questions as true (T) or false (F) by circling the appropriate letter. T F T F T F a) Variance estimates should always be positive, but covariance estimates can be either positive
More informationUsing SAS Proc Power to Perform Model-based Power Analysis for Clinical Pharmacology Studies Peng Sun, Merck & Co., Inc.
PharmaSUG00 - Paper SP05 Using SAS Proc Power to Perform Model-based Power Analysis for Clinical Pharmacology Studies Peng Sun, Merck & Co., Inc., North Wales, PA ABSRAC In this paper, we demonstrate that
More informationSupporting Information for Estimating restricted mean. treatment effects with stacked survival models
Supporting Information for Estimating restricted mean treatment effects with stacked survival models Andrew Wey, David Vock, John Connett, and Kyle Rudser Section 1 presents several extensions to the simulation
More informationA Non-parametric bootstrap for multilevel models
A Non-parametric bootstrap for multilevel models By James Carpenter London School of Hygiene and ropical Medicine Harvey Goldstein and Jon asbash Institute of Education 1. Introduction Bootstrapping is
More informationTolerance limits for a ratio of normal random variables
Tolerance limits for a ratio of normal random variables Lanju Zhang 1, Thomas Mathew 2, Harry Yang 1, K. Krishnamoorthy 3 and Iksung Cho 1 1 Department of Biostatistics MedImmune, Inc. One MedImmune Way,
More informationChapter 8 (More on Assumptions for the Simple Linear Regression)
EXST3201 Chapter 8b Geaghan Fall 2005: Page 1 Chapter 8 (More on Assumptions for the Simple Linear Regression) Your textbook considers the following assumptions: Linearity This is not something I usually
More informationChapter 2 Inferences in Simple Linear Regression
STAT 525 SPRING 2018 Chapter 2 Inferences in Simple Linear Regression Professor Min Zhang Testing for Linear Relationship Term β 1 X i defines linear relationship Will then test H 0 : β 1 = 0 Test requires
More informationThe Nonparametric Bootstrap
The Nonparametric Bootstrap The nonparametric bootstrap may involve inferences about a parameter, but we use a nonparametric procedure in approximating the parametric distribution using the ECDF. We use
More information11. Bootstrap Methods
11. Bootstrap Methods c A. Colin Cameron & Pravin K. Trivedi 2006 These transparencies were prepared in 20043. They can be used as an adjunct to Chapter 11 of our subsequent book Microeconometrics: Methods
More informationIntroduction to Linear regression analysis. Part 2. Model comparisons
Introduction to Linear regression analysis Part Model comparisons 1 ANOVA for regression Total variation in Y SS Total = Variation explained by regression with X SS Regression + Residual variation SS Residual
More informationMeasuring relationships among multiple responses
Measuring relationships among multiple responses Linear association (correlation, relatedness, shared information) between pair-wise responses is an important property used in almost all multivariate analyses.
More information2 >1. That is, a parallel study design will require
Cross Over Design Cross over design is commonly used in various type of research for its unique feature of accounting for within subject variability. For studies with short length of treatment time, illness
More informationBusiness Statistics. Lecture 10: Course Review
Business Statistics Lecture 10: Course Review 1 Descriptive Statistics for Continuous Data Numerical Summaries Location: mean, median Spread or variability: variance, standard deviation, range, percentiles,
More information5.3 Three-Stage Nested Design Example
5.3 Three-Stage Nested Design Example A researcher designs an experiment to study the of a metal alloy. A three-stage nested design was conducted that included Two alloy chemistry compositions. Three ovens
More informationDiagnostics. Gad Kimmel
Diagnostics Gad Kimmel Outline Introduction. Bootstrap method. Cross validation. ROC plot. Introduction Motivation Estimating properties of an estimator. Given data samples say the average. x 1, x 2,...,
More informationarxiv: v1 [stat.co] 26 May 2009
MAXIMUM LIKELIHOOD ESTIMATION FOR MARKOV CHAINS arxiv:0905.4131v1 [stat.co] 6 May 009 IULIANA TEODORESCU Abstract. A new approach for optimal estimation of Markov chains with sparse transition matrices
More information- 1 - By H. S Steyn, Statistical Consultation Services, North-West University (Potchefstroom Campus)
- 1 - BIOAVAILABILIY AND BIOEQUIVALENCE By H. S Steyn, Statistical Consultation Services, North-West University (Potchefstroom Campus) 1. Bioavailability (see Westlake, 1988) 1.1 Absorption: he aim is
More informationPractice of SAS Logistic Regression on Binary Pharmacodynamic Data Problems and Solutions. Alan J Xiao, Cognigen Corporation, Buffalo NY
Practice of SAS Logistic Regression on Binary Pharmacodynamic Data Problems and Solutions Alan J Xiao, Cognigen Corporation, Buffalo NY ABSTRACT Logistic regression has been widely applied to population
More informationConfidence Intervals in Ridge Regression using Jackknife and Bootstrap Methods
Chapter 4 Confidence Intervals in Ridge Regression using Jackknife and Bootstrap Methods 4.1 Introduction It is now explicable that ridge regression estimator (here we take ordinary ridge estimator (ORE)
More informationSpecial Topics. Handout #4. Diagnostics. Residual Analysis. Nonlinearity
Special Topics Diagnostics Residual Analysis As with linear regression, an analysis of residuals is necessary to assess the model adequacy. The same techniques may be employed such as plotting residuals
More informationEstimation of AUC from 0 to Infinity in Serial Sacrifice Designs
Estimation of AUC from 0 to Infinity in Serial Sacrifice Designs Martin J. Wolfsegger Department of Biostatistics, Baxter AG, Vienna, Austria Thomas Jaki Department of Statistics, University of South Carolina,
More informationRemedial Measures for Multiple Linear Regression Models
Remedial Measures for Multiple Linear Regression Models Yang Feng http://www.stat.columbia.edu/~yangfeng Yang Feng (Columbia University) Remedial Measures for Multiple Linear Regression Models 1 / 25 Outline
More informationFitting PK Models with SAS NLMIXED Procedure Halimu Haridona, PPD Inc., Beijing
PharmaSUG China 1 st Conference, 2012 Fitting PK Models with SAS NLMIXED Procedure Halimu Haridona, PPD Inc., Beijing ABSTRACT Pharmacokinetic (PK) models are important for new drug development. Statistical
More informationBayesian inference for sample surveys. Roderick Little Module 2: Bayesian models for simple random samples
Bayesian inference for sample surveys Roderick Little Module : Bayesian models for simple random samples Superpopulation Modeling: Estimating parameters Various principles: least squares, method of moments,
More informationUnit 10: Simple Linear Regression and Correlation
Unit 10: Simple Linear Regression and Correlation Statistics 571: Statistical Methods Ramón V. León 6/28/2004 Unit 10 - Stat 571 - Ramón V. León 1 Introductory Remarks Regression analysis is a method for
More informationTHE PRINCIPLES AND PRACTICE OF STATISTICS IN BIOLOGICAL RESEARCH. Robert R. SOKAL and F. James ROHLF. State University of New York at Stony Brook
BIOMETRY THE PRINCIPLES AND PRACTICE OF STATISTICS IN BIOLOGICAL RESEARCH THIRD E D I T I O N Robert R. SOKAL and F. James ROHLF State University of New York at Stony Brook W. H. FREEMAN AND COMPANY New
More informationModel Selection, Estimation, and Bootstrap Smoothing. Bradley Efron Stanford University
Model Selection, Estimation, and Bootstrap Smoothing Bradley Efron Stanford University Estimation After Model Selection Usually: (a) look at data (b) choose model (linear, quad, cubic...?) (c) fit estimates
More informationStatistics - Lecture One. Outline. Charlotte Wickham 1. Basic ideas about estimation
Statistics - Lecture One Charlotte Wickham wickham@stat.berkeley.edu http://www.stat.berkeley.edu/~wickham/ Outline 1. Basic ideas about estimation 2. Method of Moments 3. Maximum Likelihood 4. Confidence
More informationRANDOM and REPEATED statements - How to Use Them to Model the Covariance Structure in Proc Mixed. Charlie Liu, Dachuang Cao, Peiqi Chen, Tony Zagar
Paper S02-2007 RANDOM and REPEATED statements - How to Use Them to Model the Covariance Structure in Proc Mixed Charlie Liu, Dachuang Cao, Peiqi Chen, Tony Zagar Eli Lilly & Company, Indianapolis, IN ABSTRACT
More informationMonte Carlo Study on the Successive Difference Replication Method for Non-Linear Statistics
Monte Carlo Study on the Successive Difference Replication Method for Non-Linear Statistics Amang S. Sukasih, Mathematica Policy Research, Inc. Donsig Jang, Mathematica Policy Research, Inc. Amang S. Sukasih,
More informationPaper Equivalence Tests. Fei Wang and John Amrhein, McDougall Scientific Ltd.
Paper 11683-2016 Equivalence Tests Fei Wang and John Amrhein, McDougall Scientific Ltd. ABSTRACT Motivated by the frequent need for equivalence tests in clinical trials, this paper provides insights into
More informationPaper CD Erika Larsen and Timothy E. O Brien Loyola University Chicago
Abstract: Paper CD-02 2015 SAS Software as an Essential Tool in Statistical Consulting and Research Erika Larsen and Timothy E. O Brien Loyola University Chicago Modelling in bioassay often uses linear,
More informationContents. Acknowledgments. xix
Table of Preface Acknowledgments page xv xix 1 Introduction 1 The Role of the Computer in Data Analysis 1 Statistics: Descriptive and Inferential 2 Variables and Constants 3 The Measurement of Variables
More informationNon-Linear Models. Estimating Parameters from a Non-linear Regression Model
Non-Linear Models Mohammad Ehsanul Karim Institute of Statistical Research and training; University of Dhaka, Dhaka, Bangladesh Estimating Parameters from a Non-linear Regression Model
More informationHigher-Order von Mises Expansions, Bagging and Assumption-Lean Inference
Higher-Order von Mises Expansions, Bagging and Assumption-Lean Inference Andreas Buja joint with: Richard Berk, Lawrence Brown, Linda Zhao, Arun Kuchibhotla, Kai Zhang Werner Stützle, Ed George, Mikhail
More informationBootstrap Confidence Intervals
Bootstrap Confidence Intervals Patrick Breheny September 18 Patrick Breheny STA 621: Nonparametric Statistics 1/22 Introduction Bootstrap confidence intervals So far, we have discussed the idea behind
More informationBootstrap confidence levels for phylogenetic trees B. Efron, E. Halloran, and S. Holmes, 1996
Bootstrap confidence levels for phylogenetic trees B. Efron, E. Halloran, and S. Holmes, 1996 Following Confidence limits on phylogenies: an approach using the bootstrap, J. Felsenstein, 1985 1 I. Short
More informationAppendix D INTRODUCTION TO BOOTSTRAP ESTIMATION D.1 INTRODUCTION
Appendix D INTRODUCTION TO BOOTSTRAP ESTIMATION D.1 INTRODUCTION Bootstrapping is a general, distribution-free method that is used to estimate parameters ofinterest from data collected from studies or
More informationOutline. Confidence intervals More parametric tests More bootstrap and randomization tests. Cohen Empirical Methods CS650
Outline Confidence intervals More parametric tests More bootstrap and randomization tests Parameter Estimation Collect a sample to estimate the value of a population parameter. Example: estimate mean age
More informationFractional Polynomial Regression
Chapter 382 Fractional Polynomial Regression Introduction This program fits fractional polynomial models in situations in which there is one dependent (Y) variable and one independent (X) variable. It
More informationThis is a Randomized Block Design (RBD) with a single factor treatment arrangement (2 levels) which are fixed.
EXST3201 Chapter 13c Geaghan Fall 2005: Page 1 Linear Models Y ij = µ + βi + τ j + βτij + εijk This is a Randomized Block Design (RBD) with a single factor treatment arrangement (2 levels) which are fixed.
More informationInterval Estimation III: Fisher's Information & Bootstrapping
Interval Estimation III: Fisher's Information & Bootstrapping Frequentist Confidence Interval Will consider four approaches to estimating confidence interval Standard Error (+/- 1.96 se) Likelihood Profile
More informationSPECIAL TOPICS IN REGRESSION ANALYSIS
1 SPECIAL TOPICS IN REGRESSION ANALYSIS Representing Nominal Scales in Regression Analysis There are several ways in which a set of G qualitative distinctions on some variable of interest can be represented
More informationCh 2: Simple Linear Regression
Ch 2: Simple Linear Regression 1. Simple Linear Regression Model A simple regression model with a single regressor x is y = β 0 + β 1 x + ɛ, where we assume that the error ɛ is independent random component
More informationCorrelation and the Analysis of Variance Approach to Simple Linear Regression
Correlation and the Analysis of Variance Approach to Simple Linear Regression Biometry 755 Spring 2009 Correlation and the Analysis of Variance Approach to Simple Linear Regression p. 1/35 Correlation
More informationAppendix A Summary of Tasks. Appendix Table of Contents
Appendix A Summary of Tasks Appendix Table of Contents Reporting Tasks...357 ListData...357 Tables...358 Graphical Tasks...358 BarChart...358 PieChart...359 Histogram...359 BoxPlot...360 Probability Plot...360
More informationResampling and the Bootstrap
Resampling and the Bootstrap Axel Benner Biostatistics, German Cancer Research Center INF 280, D-69120 Heidelberg benner@dkfz.de Resampling and the Bootstrap 2 Topics Estimation and Statistical Testing
More informationSTAT Section 2.1: Basic Inference. Basic Definitions
STAT 518 --- Section 2.1: Basic Inference Basic Definitions Population: The collection of all the individuals of interest. This collection may be or even. Sample: A collection of elements of the population.
More informationStatistics for analyzing and modeling precipitation isotope ratios in IsoMAP
Statistics for analyzing and modeling precipitation isotope ratios in IsoMAP The IsoMAP uses the multiple linear regression and geostatistical methods to analyze isotope data Suppose the response variable
More informationTopic 14: Inference in Multiple Regression
Topic 14: Inference in Multiple Regression Outline Review multiple linear regression Inference of regression coefficients Application to book example Inference of mean Application to book example Inference
More informationWhy Do Statisticians Treat Predictors as Fixed? A Conspiracy Theory
Why Do Statisticians Treat Predictors as Fixed? A Conspiracy Theory Andreas Buja joint with the PoSI Group: Richard Berk, Lawrence Brown, Linda Zhao, Kai Zhang Ed George, Mikhail Traskin, Emil Pitkin,
More informationWhat to do today (Nov 22, 2018)?
What to do today (Nov 22, 2018)? Part 1. Introduction and Review (Chp 1-5) Part 2. Basic Statistical Inference (Chp 6-9) Part 3. Important Topics in Statistics (Chp 10-13) Part 4. Further Topics (Selected
More informationBootstrapping, Randomization, 2B-PLS
Bootstrapping, Randomization, 2B-PLS Statistics, Tests, and Bootstrapping Statistic a measure that summarizes some feature of a set of data (e.g., mean, standard deviation, skew, coefficient of variation,
More informationSAS Procedures Inference about the Line ffl model statement in proc reg has many options ffl To construct confidence intervals use alpha=, clm, cli, c
Inference About the Slope ffl As with all estimates, ^fi1 subject to sampling var ffl Because Y jx _ Normal, the estimate ^fi1 _ Normal A linear combination of indep Normals is Normal Simple Linear Regression
More informationOn block bootstrapping areal data Introduction
On block bootstrapping areal data Nicholas Nagle Department of Geography University of Colorado UCB 260 Boulder, CO 80309-0260 Telephone: 303-492-4794 Email: nicholas.nagle@colorado.edu Introduction Inference
More informationReview of Statistics 101
Review of Statistics 101 We review some important themes from the course 1. Introduction Statistics- Set of methods for collecting/analyzing data (the art and science of learning from data). Provides methods
More informationIES 612/STA 4-573/STA Winter 2008 Week 1--IES 612-STA STA doc
IES 612/STA 4-573/STA 4-576 Winter 2008 Week 1--IES 612-STA 4-573-STA 4-576.doc Review Notes: [OL] = Ott & Longnecker Statistical Methods and Data Analysis, 5 th edition. [Handouts based on notes prepared
More information401 Review. 6. Power analysis for one/two-sample hypothesis tests and for correlation analysis.
401 Review Major topics of the course 1. Univariate analysis 2. Bivariate analysis 3. Simple linear regression 4. Linear algebra 5. Multiple regression analysis Major analysis methods 1. Graphical analysis
More informationOutline. Review regression diagnostics Remedial measures Weighted regression Ridge regression Robust regression Bootstrapping
Topic 19: Remedies Outline Review regression diagnostics Remedial measures Weighted regression Ridge regression Robust regression Bootstrapping Regression Diagnostics Summary Check normality of the residuals
More informationThe Steps to Follow in a Multiple Regression Analysis
ABSTRACT The Steps to Follow in a Multiple Regression Analysis Theresa Hoang Diem Ngo, Warner Bros. Home Video, Burbank, CA A multiple regression analysis is the most powerful tool that is widely used,
More informationRatio of Polynomials Fit Many Variables
Chapter 376 Ratio of Polynomials Fit Many Variables Introduction This program fits a model that is the ratio of two polynomials of up to fifth order. Instead of a single independent variable, these polynomials
More informationPackage Rsurrogate. October 20, 2016
Type Package Package Rsurrogate October 20, 2016 Title Robust Estimation of the Proportion of Treatment Effect Explained by Surrogate Marker Information Version 2.0 Date 2016-10-19 Author Layla Parast
More informationGenerating Half-normal Plot for Zero-inflated Binomial Regression
Paper SP05 Generating Half-normal Plot for Zero-inflated Binomial Regression Zhao Yang, Xuezheng Sun Department of Epidemiology & Biostatistics University of South Carolina, Columbia, SC 29208 SUMMARY
More informationMIXED MODELS FOR REPEATED (LONGITUDINAL) DATA PART 2 DAVID C. HOWELL 4/1/2010
MIXED MODELS FOR REPEATED (LONGITUDINAL) DATA PART 2 DAVID C. HOWELL 4/1/2010 Part 1 of this document can be found at http://www.uvm.edu/~dhowell/methods/supplements/mixed Models for Repeated Measures1.pdf
More informationConfidence Intervals of the Simple Difference between the Proportions of a Primary Infection and a Secondary Infection, Given the Primary Infection
Biometrical Journal 42 (2000) 1, 59±69 Confidence Intervals of the Simple Difference between the Proportions of a Primary Infection and a Secondary Infection, Given the Primary Infection Kung-Jong Lui
More informationIn many situations, there is a non-parametric test that corresponds to the standard test, as described below:
There are many standard tests like the t-tests and analyses of variance that are commonly used. They rest on assumptions like normality, which can be hard to assess: for example, if you have small samples,
More informationExact Inference for the Two-Parameter Exponential Distribution Under Type-II Hybrid Censoring
Exact Inference for the Two-Parameter Exponential Distribution Under Type-II Hybrid Censoring A. Ganguly, S. Mitra, D. Samanta, D. Kundu,2 Abstract Epstein [9] introduced the Type-I hybrid censoring scheme
More informationST745: Survival Analysis: Nonparametric methods
ST745: Survival Analysis: Nonparametric methods Eric B. Laber Department of Statistics, North Carolina State University February 5, 2015 The KM estimator is used ubiquitously in medical studies to estimate
More informationMultiple Regression Examples
Multiple Regression Examples Example: Tree data. we have seen that a simple linear regression of usable volume on diameter at chest height is not suitable, but that a quadratic model y = β 0 + β 1 x +
More informationResearch Article A Nonparametric Two-Sample Wald Test of Equality of Variances
Advances in Decision Sciences Volume 211, Article ID 74858, 8 pages doi:1.1155/211/74858 Research Article A Nonparametric Two-Sample Wald Test of Equality of Variances David Allingham 1 andj.c.w.rayner
More informationBootstrap inference for the finite population total under complex sampling designs
Bootstrap inference for the finite population total under complex sampling designs Zhonglei Wang (Joint work with Dr. Jae Kwang Kim) Center for Survey Statistics and Methodology Iowa State University Jan.
More informationInferences for the Ratio: Fieller s Interval, Log Ratio, and Large Sample Based Confidence Intervals
Inferences for the Ratio: Fieller s Interval, Log Ratio, and Large Sample Based Confidence Intervals Michael Sherman Department of Statistics, 3143 TAMU, Texas A&M University, College Station, Texas 77843,
More informationdf=degrees of freedom = n - 1
One sample t-test test of the mean Assumptions: Independent, random samples Approximately normal distribution (from intro class: σ is unknown, need to calculate and use s (sample standard deviation)) Hypotheses:
More informationStat 5101 Lecture Notes
Stat 5101 Lecture Notes Charles J. Geyer Copyright 1998, 1999, 2000, 2001 by Charles J. Geyer May 7, 2001 ii Stat 5101 (Geyer) Course Notes Contents 1 Random Variables and Change of Variables 1 1.1 Random
More informationSTAT 704 Sections IRLS and Bootstrap
STAT 704 Sections 11.4-11.5. IRLS and John Grego Department of Statistics, University of South Carolina Stat 704: Data Analysis I 1 / 14 LOWESS IRLS LOWESS LOWESS (LOcally WEighted Scatterplot Smoothing)
More informationGeneral Linear Model (Chapter 4)
General Linear Model (Chapter 4) Outcome variable is considered continuous Simple linear regression Scatterplots OLS is BLUE under basic assumptions MSE estimates residual variance testing regression coefficients
More informationOPTIMIZATION OF FIRST ORDER MODELS
Chapter 2 OPTIMIZATION OF FIRST ORDER MODELS One should not multiply explanations and causes unless it is strictly necessary William of Bakersville in Umberto Eco s In the Name of the Rose 1 In Response
More informationFiducial Inference and Generalizations
Fiducial Inference and Generalizations Jan Hannig Department of Statistics and Operations Research The University of North Carolina at Chapel Hill Hari Iyer Department of Statistics, Colorado State University
More informationQuantitative Empirical Methods Exam
Quantitative Empirical Methods Exam Yale Department of Political Science, August 2016 You have seven hours to complete the exam. This exam consists of three parts. Back up your assertions with mathematics
More informationHandout 1: Predicting GPA from SAT
Handout 1: Predicting GPA from SAT appsrv01.srv.cquest.utoronto.ca> appsrv01.srv.cquest.utoronto.ca> ls Desktop grades.data grades.sas oldstuff sasuser.800 appsrv01.srv.cquest.utoronto.ca> cat grades.data
More informationThe program for the following sections follows.
Homework 6 nswer sheet Page 31 The program for the following sections follows. dm'log;clear;output;clear'; *************************************************************; *** EXST734 Homework Example 1
More information