Numerical integration with an application in sample size reestimation

Adrian Mander () and Jack Bowden
Additional contact information
Jack Bowden: MRC Biostatistics Unit Hub for Trials Methodology, Cambridge

United Kingdom Stata Users' Group Meetings 2012 from Stata Users Group

Abstract: We introduce a new integrate() function for Mata that evaluates single-dimensional integrals. This function uses three different Gaussian quadrature algorithms: Gauss–Hermite and Gauss–Laguerre for indefinite integrals; and Gauss–Legendre for definite integrals. The algorithms were implemented using the methods of Golub and Welsch (Mathematics of Computation, 1968). The user can specify any integrand by defining a new function in the Mata language. The integrand function is allowed to have two arguments: the first is the variable of integration, and the second is a real scalar. Thus the integrate() function can be used in combination with optimise() to solve for the value of x in an integral f(u,x) du = 0.05. Such calculations are used in the sample size re-estimation methodology introduced by Li, Shih, and Xie (Biostatistics, 2002). We apply these methods to a clinical trial where a single interim analysis is carried out, and the analysis is used to reevaluate the sample size.

Date: 2012-09-22
References: Add references at CitEc
Citations:

Downloads: (external link)
http://repec.org/usug2012/UK12_mander.pdf (application/pdf)

Related works:
This item may be available elsewhere in EconPapers: Search for items with the same title.

Export reference: BibTeX RIS (EndNote, ProCite, RefMan) HTML/Text

Persistent link: https://EconPapers.repec.org/RePEc:boc:usug12:03

Access Statistics for this paper

More papers in United Kingdom Stata Users' Group Meetings 2012 from Stata Users Group Contact information at EDIRC.
Bibliographic data for series maintained by Christopher F Baum ().