Calculates back-transformed means with approximate standard errors and confidence intervals (V.M. Cave).
|Controls printed output (
||The confidence intervals of the back-transformed means to plot (
||Constant c for the
||Parameter k for
||Expression, formed using argument
||Expression, formed using argument
||Probability for the confidence intervals; default 0.95|
||Order in which the back-transformed means are plotted (
||Whether to use the current pen definitions for plotting; (
||Window to use for plot; default 3|
||Supplies the transformed mean(s)|
||Supplies the standard error(s) of the transformed mean(s)|
||Degrees of freedom to construct the confidence intervals; default
||Number of decimal places for printing; default
||Saves the back-transformed means|
||Saves the approximate standard errors for the back-transformed means|
||Saves the approximate confidence intervals for the back-transformed means|
||Saves the back-transformed confidence intervals for the back-transformed means|
||Title for plot; default
||Title for y-axis; default
||Title for x-axis; default
BACKTRANSFORM calculates back-transformed means, with approximate standard errors and confidence intervals. The means and corresponding standard errors, for back-transforming, are supplied using the
SEMEANS parameters, respectively, as either tables, variates or scalars. If
MEANS supplies a table or variate,
SEMEANS can be either of the same type or a scalar, whereas if
MEANS supplies a scalar,
SEMEANS must be a scalar.
The degrees of freedom, used to construct the confidence intervals, can be set using the
DF parameter. If these are not supplied, z-scores are used to form the confidence intervals. The probability for the confidence intervals is specified by the
CIPROBABILITY option; the default 0.95 gives 95% confidence intervals.
The function that was used to transform the data prior to analysis is specified using the
TRANSFORMATION option; the default takes the identity link (i.e. no transformation). The natural logarithm, log10, logit, square root, reciprocal, power, probit, complementary log-log, log-ratio and angular (or arcsine-square root) functions are provided directly by the
You can also define your own transformation by setting
calculated and providing expressions to calculate the inverse and first derivative of the transformation, using the
DERIVATIVE options, respectively. The calculations are specified in terms of the argument
Y. Thus, for example, the
logarithm transformation could be specified by setting options
CLOG option sets the constant c used by the
log10 transformations (default 0), the
EXPONENT option sets the exponent used by the
power transformation (default -2) and the
KLOGRATIO option sets the parameter k used by the
logratio transformation (default 1).
PLOT option allows you to request plots of the results, using settings:
||for back-transformed means and back-transformed confidence intervals (i.e. confidence intervals that maintain exactly the percent coverage on the transformed scale),|
||for back-transformed means with approximate confidence intervals (i.e. derived from the approximate standard errors), and|
||for back-transformed means with both types of confidence interval.|
By default no plots are produced.
XTITLE parameters can supply an overall title, a y-axis title and a x-axis title for the plot, respectively. By default, neither an overall title nor a y-axis title is displayed. The default for the x-axis, when
MEANS supplies the transformed means as a table, is to use the identifiers of the table’s classifying factors to form a title for the x-axis. If
MEANS supplies a variate or scalar, the default is not to display an x-axis title. To omit the x-axis title, you can supply a blank string i.e.
By default, the pen attributes used for plotting are determined automatically within the procedure. However, you can set
yes, to request that the current
THICKNESS pen definitions of pens 1 and 2 are used. Pen 1 controls the colour, symbol, symbol size and line thickness used to plot the back-transformed means with back-transformed confidence intervals, whereas Pen 2 controls these attributes when the back-transformed means are plotted with approximate confidence intervals. The
WINDOW option specifies the window used for plotting; default 3.
By default, the back-transformed means are plotted in order of their ordinal level, however you can use the
DIRECTION option to request that they are plotted in ascending or descending numerical value instead.
SEBACKTRANSFORMEDMEANS parameters can save back-transformed means, approximate standard errors, in data structures of the same type as
CIBACKTRANSFORMED parameters can save pointers of approximate confidence intervals and back-transformed confidence intervals, respectively. The pointers contain two data structures elements, of the same type as
MEANS, storing the lower and the upper confidence limits, respectively.
Printed output is controlled by the
||provides a description of the output,|
||prints the transformed means, with their standard errors and confidence intervals, and|
||prints the back-transformed means, with their approximate standard errors and confidence intervals.|
You can set the number of decimals places to appear in the printed output, using the
BACKTRANSFORM uses a first-order Taylor series expansion, to obtain approximate standard errors for the back-transformed means. The methodology is described in Jørgensen & Pedersen (1998). In brief, let û denote the estimated mean on the transformed scale, seu its standard error, g() the transformation function and g′() the first derivative of the transformation function. On the back-transformed scale the estimated mean (ӯ) and standard error (seӯ) are approximated by g-1(û) and seu / mod(g′(ӯ)), respectively.
A back-transformed confidence interval is given by
(g-1( û – t × seu), g-1( û + t × seu))
where t is the upper (1-
CIPROBABILITY)/2 critical value for the t distribution. Note that, when the degrees of freedom have not been set using the
DF parameter, the z-score (i.e. the Normal distribution) is used to construct the confidence interval instead of the t-value. For the logratio and reciprocal transformations, when the estimated mean and confidence
limit on the transformed scale lie on different sides of the singularity in the inverse function
of the transformation, the back-transformed confidence limit is set to a missing value.
An approximate confidence interval
(ӯ – t × seӯ, ӯ + t × seӯ)
is also provided. However, this should be used only to evaluate the validity of the approximate standard error. If the back-transformed and approximate confidence intervals differ greatly, the approximate standard error is inadequate.
BACKTRANSFORM will work when MEANS and/or SEMEANS supplies a restricted variate; however, if both variates are restricted they must be restricted in the same way. Furthermore, their unrestricted length must be the same.
Jørgensen, E., & Pedersen, A.R. (1998). How to obtain those nasty standard errors from transformed data – and why they shouldn’t be used. Biometry Research Unit, Danish Institute of Agricultural Sciences.
CAPTION 'BACKTRANSFORM examples',\ 'Pup rat data from Weil (1970, Food and Cosmetics Toxicology)';\ STYLE=meta,plain FACTOR [NVALUES=64; LEVELS=32; VALUES=2(1...32)] Litter & [LEVELS=2; LABELS=!t('4','21'); VALUES=(1,2)32] Time & [LEVELS=2; LABELS=!t('control','treated'); VALUES=32(1,2)] Diet VARIATE [NVALUES=64] Pups READ [PRINT=*] Pups 13 13 12 12 9 9 9 9 8 8 8 8 13 12 12 11 10 9 10 9 9 8 13 11 5 4 7 5 10 7 10 7 12 12 11 11 10 10 9 9 11 10 10 9 10 9 9 8 9 8 5 4 9 7 7 4 10 5 6 3 10 3 7 0 : CAPTION 'ANOVA Example: Means and standard errors obtained from AKEEP';\ STYLE=meta BLOCKS Litter TREATMENTS Diet*Time ANOVA [PRINT=aov] log(Pups+1) AKEEP TERMS=Diet.Time; MEANS=meansA; SEMEANS=semsA; DFCEFFECTS=dfA BACKTRANSFORM [TRANSFORMATION=log; CLOG=1; PLOT=both]\ meansA; SEMEANS=semsA; DF=dfA CAPTION 'GLM Example: Means and standard errors obtained from PREDICT';\ STYLE=meta MODEL [DISTRIBUTION=Poisson; LINK=log] Pups FIT [PRINT=summary] Diet*Time RKEEP DF=dfB PREDICT [PRINT=*; PREDICTIONS=meansB; SE=semsB; BACKTRANSFORM=none]\ Diet,Time BACKTRANSFORM [TRANSFORMATION=log; PLOT=both] \ meansB; SEMEANS=semsB; DF=dfB CAPTION 'REML Example: Means and standard errors obtained from VPREDICT';\ STYLE=meta VCOMPONENT [FIXED=Diet*Time] RANDOM=Litter REML log(Pups+1) VKEEP [WMETHOD=add] TERMS=Diet.Time; DDF=dfC VPREDICT [PRINT=*; PREDICTIONS=meansC; SE=semsC] Diet,Time BACKTRANSFORM [TRANSFORMATION=log; CLOG=1; PLOT=both]\ meansC; SEMEANS=semsC; DF=dfC CAPTION 'GLMM Example: Means and standard errors obtained from VPREDICT';\ STYLE=meta GLMM [PRINT=model,components; DISTRIBUTION=Poisson; LINK=log;\ FIXED=Diet*Time; RANDOM=Litter; DISPERSION=*] Pups VKEEP [WMETHOD=add] TERMS=Diet.Time; DDF=dfD VPREDICT [PRINT=*; PREDICTIONS=meansD; SE=semsD] Diet,Time BACKTRANSFORM [TRANSFORMATION=log; PLOT=both]\ meansD; SEMEANS=semsD; DF=dfD