Delta Method Standard Error Sas
Contents |
Samples & SAS Notes Focus Areas SUPPORT License Assistance Manage My Software Account Downloads & Hot Fixes TRAINING & BOOKS Books Training Certification SAS Global Academic Program SAS OnDemand For Academics USERS GROUPS Advanced Search support.sas.com Knowledge Base Support Training & Books Store Support Communities Knowledge Base Products & Solutions delta method sas code System Requirements Install Center Third-Party Software Reference Documentation Papers Samples & SAS Notes Browse by Topic standard error sas proc means Search Samples Search Usage Notes Search Installation Notes Search Problem Notes Focus Areas Usage Note 56476: Confidence interval for a ratio of two linear combinations calculate standard error in sas of model parameters The following discusses estimation of a ratio of functions of model parameters. Procedures to estimate other ratios such as risk (prevalence) ratios, odds ratios, hazard ratios, ratios of means are mentioned in the list of Frequently-Asked for
Robust Standard Error Sas
Statistics. Two methods for estimating the confidence limits for a ratio will be discussed using the common example of estimation of the ED50 and the relative potency. In bioassay studies that investigate the effect of various doses of one or more drugs, the ED50 is the dose at which 50% of the subjects respond. The relative potency is the ratio of doses of two drugs that produce the same response. Both the ED50 and the relative potency are estimated as ratios of standard deviation sas model parameters in a binary response model. Quantal bioassay example Stokes et. al. (2000) discuss a study on mice comparing two drugs applied at several doses. The data are presented below. data assay; input drug $ dose ndead total; ldose=log(dose); datalines; N 0.01 0 30 N 0.03 1 30 N 0.10 1 10 N 0.30 1 10 N 1.00 4 10 N 3.00 4 10 N 10.00 5 10 N 30.00 7 10 S 0.30 0 10 S 1.00 0 10 S 3.00 1 10 S 10.00 4 10 S 30.00 5 10 S 100.00 8 10 ; A sequence of logistic models fit to the data show that a model with separate intercepts for the drugs but with a common slope fits well. The following statements fit this final model. Note that with the NOINT option the two DRUG parameters are the intercepts for the two drugs. See this note for more on assessing the need for separate slopes when modeling data from multiple groups. The OUTEST= and COVOUT options are not needed to fit the model but provide a data set containing the estimated model parameters and covariance matrix for use later. The STORE statement is used to save the model information for use later by the NLEstimate macro. This data set is displayed by the PROC PRINT step. proc logistic data=assay outest=parms covout; class drug / param=glm; model ndead/total = drug ldose / noint; store out=assaymod; run; proc print noobs; ru
transformed regression parameters in R using the delta method? The purpose of this page is to introduce estimation of standard errors using the delta method. Examples include manual calculation of standard confidence interval sas errors via the delta method and then confirmation using the function deltamethod so that
Variance Sas
the reader may understand the calculations and know how to use deltamethod. This page uses the following packages Make sure that
T Test Sas
you can load them before trying to run the examples on this page. We will need the msm package to use the deltamethodfunction. If you do not have a package installed, run: install.packages("packagename"), or if you http://support.sas.com/kb/56476.html see the version is out of date, run: update.packages(). library(msm) Version info: Code for this page was tested in R version 3.1.1 (2014-07-10)
On: 2014-08-01
With: pequod 0.0-3; msm 1.4; phia 0.1-5; effects 3.0-0; colorspace 1.2-4; RColorBrewer 1.0-5; xtable 1.7-3; car 2.0-20; foreign 0.8-61; Hmisc 3.14-4; Formula 1.1-2; survival 2.37-7; lattice 0.20-29; mgcv 1.8-1; nlme 3.1-117; png 0.1-7; gridExtra 0.9.1; reshape2 1.4; ggplot2 1.0.0; vcd 1.3-1; rjson 0.2.14; RSQLite 0.11.4; http://www.ats.ucla.edu/stat/r/faq/deltamethod.htm DBI 0.2-7; knitr 1.6 Background to delta method Often in addition to reporting parameters fit by a model, we need to report some transformation of these parameters. The transformation can generate the point estimates of our desired values, but the standard errors of these point estimates are not so easily calculated. They can, however, be well approximated using the delta method. The delta method approximates the standard errors of transformations of random variable using a first-order Taylor approximation. Regression coefficients are themselves random variables, so we can use the delta method to approximate the standard errors of their transformations. Although the delta method is often appropriate to use with large samples, this page is by no means an endorsement of the use of the delta method over other methods to estimate standard errors, such as bootstrapping. Essentially, the delta method involves calculating the variance of the Taylor series approximation of a function. We, thus, first get the Taylor series approximation of the function using the first two terms of the Taylor expansion of the transformation function about the mean of of the random variable. Let \(G\) be the transformation function and \(U\) be the mean vector of random variables \(X=(x1,x2,...)\). The first two terms of the Taylor expansion are then an approximation