Performs inverse prediction given a fitted simple linear regression model.
SUMWTF — Sum of
products of weights with frequencies from the fitted regression.
(Input, if INTCEP
= 1)
In the ordinary case when weights and frequencies are all one, SUMWTF equals the
number of observations used in the fit of the model.
DFS2 — Degrees of
freedom for estimate of σ2.
(Input)
If IY0 = 1, DFS2 is the degrees of
freedom for error from the fitted regression. If IY0 = 0,
DFS2 is the
pooled degrees of freedom from the estimate of sigma-squared based on the fitted
regression and the additional responses used to compute the mean Y0.
B — Vector of length INTCEP + 1 containing a least-squares solution for the intercept and slope. (Input)
INTCEP |
Intercept |
Slope |
0 |
|
B(1) |
1 |
B(1) |
B(2) |
XYMEAN — Vector
of length 2 with the mean of the independent and dependent variables,
respectively. (Input, if INTCEP
= 1)
If INTCEP
= 0, XYMEAN is
not referenced and can be a vector of length 1.
SSX — Sum of
squares for x. (Input)
If INTCEP
= 1, SSX
is the sum of squares of deviations of x from its mean. If INTCEP = 0,
SSX
must not be corrected for the mean.
S2 —
s2, the estimate of the variance of the
error in the model. (Input)
If IY0= 1, S2 is the estimate of
σ2 from the
fitted regression. If IY0= 0, S2 is the pooled
estimate of σ2 based on the
fitted regression, and the additional responses used to compute the mean Y0.
Y0 — Value of the response variable for which an interval estimate of the corresponding independent variable value is desired. (Input)
X0HAT — Point estimate of the independent variable. (Output)
XLOWER — Lower limit of the interval estimate for the independent variable. (Output)
XUPPER — Upper limit of the interval estimate for the independent variable. (Output)
INTCEP —
Intercept option. (Input)
Default: INTCEP = 1.
INTCEP Action
0 An intercept is not in the model.
1 An intercept is in the model.
CONPER —
Confidence level for the interval estimation. (Input)
CONPER must be
expressed as a percentage between 0.0 and 100.0. CONPER often will be
90.0, 95.0, 99.0. For one-sided confidence intervals with confidence level ONECL, set CONPER = 100.0 − 2.0 * (100.0 − ONECL).
Default:
CONPER =
95.0.
IY0 — Option for
Y0.
(Input)
Default: IY0 = 1.
IY0 Meaning
0 Y0 is a sample mean of one or more responses.
1 Y0 is the true mean response.
SWTFY0 — Sum of
products of weights with frequencies for Y0.
(Input, if IY0 = 0)
In the ordinary case, when weights and frequencies are all one, SWTFY0 is the number
of observations used to obtain the mean Y0. If IY0 = 1, SWTFY0 is not
referenced.
Generic: CALL RINPF (SUMWTF,DFS2 , B , XYMEAN , SSX , S2 , Y0 , X0HAT, XLOWER , XUPPER [,…])
Specific: The specific interface names are S_RINPF and D_RINPF.
Single: CALL RINPF (SUMWTF,DFS2 , INTCEP, B , XYMEAN , SSX , S2 , CONPER , IY0, SWTFY0, Y0, X0HAT, XLOWER , XUPPER)
Double: The double precision name is DRINPF.
Routine RINPF computes a confidence interval on the independent variable setting x0 for a given response y0 from the results of a straight line fit. Here, y0 may represent the mean of k responses or the true mean response. The results of routine RLINE or RONE can be used for input into RINPF. The simple linear regression model is assumed,
where the ɛ i’s are
independently distributed normal errors with mean zero and variance σ2/wi. Here, n
is the total number of observations used in the fit of the line, i.e., n = DFE + INTCEP
+ 1 where DFE is
the degrees of freedom from the fitted regression. Also, k is the number
of additional responses used to determine y0. The
wi’s
are the weights that must be used in the fit of the model. The methodology is
discussed by Graybill (1976, pages 280−283). For the
case when IY0 =
1, the estimate of σ2,
s2
(stored in S2),
is the usual estimate of σ2 from the
fitted regression based on the first n observations. If IY0 = 0, the estimate
of σ2 is a pooled
estimator based on the fitted regression and the k responses that produce
.
This pooled estimator (stored in S2) is given by
where (n − 2) + (k − 1) (stored in DFS2) is the pooled degrees of freedom for s2.
First, a point estimate (stored in X0HAT)
is computed by
Then, a test of the hypothesis H0 : β1 = 0 vs. Ha : β1 ≠ 0 is performed. If H0 is accepted, the model becomes yi = β0 + ɛ i, and therefore no confidence interval exists for x0 because it is no longer in the model. In this case, a type 3 warning error is issued. If H0 is rejected, a confidence interval exists and is computed for the case IY0 = 1 by
where
and t is the 50 + CONPER/2 percentile of the t distribution with DFS2 degrees of freedom. The interval actually has a confidence coefficient less than that specified by CONPER.
In the weighted case, which was discussed earlier, the means (stored in XYMEAN) and the sum of squares for x (stored in SSX) are all weighted. When the variances of the ɛ i’s are all equal, ordinary least squares must be used, this corresponds to all wi= 1.
Modifications are necessary to the preceding discussion for
other cases. For the case when an intercept is not in the model, let the pooled degrees of freedom of s2 equal to
(n − 1) + (k − 1), and replace
the first term under the square root symbol with zero, with zero, and
with zero.
For the case of the true response mean, i.e, when IY0 = 1, replace the second term under the square root symbol by zero.
Informational errors
Type Code
3 2 The slope is not significant at the (100 − CONPER)% level. Confidence limits XLOWER and XUPPER cannot be obtained.
This example fits a line to a set of data discussed by Draper and Smith (1981, Table 1.1, page 9). The response y is the amount of steam used per month (in pounds), and the independent variable x is the average atmospheric temperature (in degrees Fahrenheit). A 95% confidence interval for the temperature x0 is computed given a single response of y0 = 10.
USE RINPF_INT
USE RLINE_INT
USE UMACH_INT
IMPLICIT NONE
INTEGER NOBS
PARAMETER (NOBS=25)
!
INTEGER INTCEP, IY0, NOUT
REAL B(2), B0, B1, CONPER, DFS2, S2, SSX, STAT(12), &
SUMWTF, SWTFY0, X0HAT, XDATA(NOBS), XLOWER, XUPPER, &
XYMEAN(2), Y0, YDATA(NOBS)
!
DATA XDATA/35.3, 29.7, 30.8, 58.8, 61.4, 71.3, 74.4, 76.7, 70.7, &
57.5, 46.4, 28.9, 28.1, 39.1, 46.8, 48.5, 59.3, 70.0, 70.0, &
74.5, 72.1, 58.1, 44.6, 33.4, 28.6/
DATA YDATA/10.98, 11.13, 12.51, 8.4, 9.27, 8.73, 6.36, 8.5, &
7.82, 9.14, 8.24, 12.19, 11.88, 9.57, 10.94, 9.58, 10.09, &
8.11, 6.83, 8.88, 7.68, 8.47, 8.86, 10.36, 11.08/
!
CALL RLINE (XDATA, YDATA, B0, B1, STAT=STAT)
SUMWTF = NOBS
DFS2 = STAT(10)
INTCEP = 1
B(1) = B0
B(2) = B1
XYMEAN(1) = STAT(1)
XYMEAN(2) = STAT(2)
SSX = STAT(3)*(NOBS-1)
S2 = STAT(11)/STAT(10)
CONPER = 95.0
IY0 = 0
SWTFY0 = 1.0
Y0 = 10.0
CALL RINPF (SUMWTF, DFS2, B, XYMEAN, SSX, S2, Y0, X0HAT, XLOWER, &
XUPPER, IY0=IY0, SWTFY0=SWTFY0)
CALL UMACH (2, NOUT)
WRITE (NOUT,*) 'X0HAT = ', X0HAT
WRITE (NOUT,*) '(XLOWER,XUPPER) = (', XLOWER, ',', XUPPER, ')'
END
X0HAT = 45.3846
(XLOWER,XUPPER) = (20.2627,69.347)
PHONE: 713.784.3131 FAX:713.781.9260 |