Description Usage Arguments Examples
Estimation of all parameters in the the simple linear regression model.
1 | reglin(x, y)
|
x |
|
y |
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))
}
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.