Re: [R] how to test this

From: <Jin.Li_at_csiro.au>
Date: Thu 04 Aug 2005 - 16:15:51 EST


Thank you all for the reply.
Regards,
Jin

-----Original Message-----
From: Prof Brian Ripley [mailto:ripley@stats.ox.ac.uk] Sent: Wednesday, 3 August 2005 5:20 P
To: Simon Blomberg
Cc: Li, Jin (CSE, Atherton); r-help@stat.math.ethz.ch Subject: Re: [R] how to test this

On Wed, 3 Aug 2005, Simon Blomberg wrote:

> This is two tests: Whether the slope != 1 and whether the intercept !=
0.

Neither model given has an intercept ....

> To do this, include an offset in your model:
>
> fit <- lm(y ~ x + offset(x), data=dat)

but no intercept, so use

summary(lm(y ~ 0 + x + offset(1.05*x), data=dat))

and look if the coefficient of x is significantly different from zero.

E.g.

x <- 1:10
set.seed(1)
y <- 1.05*x + rnorm(10)
summary(lm(y ~ 0 + x + offset(1.05*x)))

Coefficients:

   Estimate Std. Error t value Pr(>|t|)
x 0.03061 0.03910 0.783 0.454

is not.

>
> HTH,
>
> Simon.
>
>
> At 03:44 PM 3/08/2005, Jin.Li@csiro.au wrote:
>>
>> I am wondering how to test whether a simple linear regression model
>> (e.g. y=1.05x) is significantly different from a 1 to 1 line (i.e.
y=x).

-- 
Brian D. Ripley,                  ripley@stats.ox.ac.uk
Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/
University of Oxford,             Tel:  +44 1865 272861 (self)
1 South Parks Road,                     +44 1865 272866 (PA)
Oxford OX1 3TG, UK                Fax:  +44 1865 272595

______________________________________________
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Received on Thu Aug 04 16:21:18 2005

This archive was generated by hypermail 2.1.8 : Sun 23 Oct 2005 - 15:03:45 EST