[R] lm fitting with a specified slope

Peter Dalgaard BSA p.dalgaard at biostat.ku.dk
Thu Oct 3 22:53:57 CEST 2002


"Warnes, Gregory R" <gregory_r_warnes at groton.pfizer.com> writes:

> Assume the model:
> 
>   X1 <- rnorm(1000)
>   X2 <- rnorm(1000)
> 
>   B1 = 1.75 # known a-priorori
>   B2 = 0.5  # unknown a-priorori
> 
>   Y  <- 1.75 * X1 + 0.5 * X2  + rnorm(1000,sd=0.25)
> 
> Then wouldn't this do the trick? :
> 
>   reg <- lm( (Y - 1.75 * X1) ~ X2 ) 

I'd prefer

lm(Y ~ X2 + offset(1.75*X1))

This allows you to do things like

>  predict(lm( Y ~ X2 + offset(1.75*X1) ), data.frame(X1=10,X2=10))
[1] 22.41585

whereas

>  predict(lm( (Y - 1.75*X1) ~ X2 ), data.frame(X1=10,X2=10))
[1] 4.915847


-- 
   O__  ---- Peter Dalgaard             Blegdamsvej 3  
  c/ /'_ --- Dept. of Biostatistics     2200 Cph. N   
 (*) \(*) -- University of Copenhagen   Denmark      Ph: (+45) 35327918
~~~~~~~~~~ - (p.dalgaard at biostat.ku.dk)             FAX: (+45) 35327907
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._



More information about the R-help mailing list