reglin: Estimation in Simple Linear Regression

Description Usage Arguments Examples

View source: R/reglin.R

Description

Estimation of all parameters in the the simple linear regression model.

Usage

1
reglin(x, y)

Arguments

x
y

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
##---- Should be DIRECTLY executable !! ----
##-- ==>  Define data, use random,
##--	or do  help(data=index)  for the standard data sets.

## The function is currently defined as
function (x, y) 
{
    n <- length(x)
    xbar <- mean(x)
    ybar <- mean(y)
    x2bar <- mean(x^2)
    y2bar <- mean(y^2)
    xybar <- mean(x * y)
    beta1hat <- (xybar - xbar * ybar)/(x2bar - xbar^2)
    beta0hat <- ybar - beta1hat * xbar
    sigmahat <- sqrt(n * (y2bar + beta0hat^2 + beta1hat^2 * x2bar - 
        2 * beta0hat * ybar - 2 * beta1hat * xybar + 2 * beta0hat * 
        beta1hat * xbar)/(n - 2))
    sigmahat.beta1hat <- sigmahat/sqrt(n * (x2bar - xbar^2))
    test.statistic <- beta1hat/sigmahat.beta1hat
    return(list(beta0hat = beta0hat, beta1hat = beta1hat, sigmahat = sigmahat, 
        sigmahat.beta1hat = sigmahat.beta1hat, test.statistic = test.statistic))
  }

D1ng-Ma/MyFirstRPackage documentation built on Oct. 30, 2019, 5:50 a.m.