替换[r]中的lm系数

Question

是否可以替换lm对象中的系数？

我认为以下内容可行

# sample data 
set.seed(2157010)
x1 <- 1998:2011
x2 <- x1 + rnorm(length(x1))
y <- 3*x1 + rnorm(length(x1))
fit <- lm( y ~ x1 + x2)

# view origional coefficeints
coef(fit)

# replace coefficent with new values
fit$coef(fit$coef[2:3]) <- c(5, 1)

# view new coefficents
coef(fit)

非常感谢任何帮助

Answer 1

您的代码无法重现，因为代码中的错误很少。这是更正后的版本，也显示了您的错误：

set.seed(2157010) #forgot set.
x1 <- 1998:2011
x2 <- x1 + rnorm(length(x1))
y <- 3*x2 + rnorm(length(x1)) #you had x, not x1 or x2
fit <- lm( y ~ x1 + x2)

# view original coefficients
coef(fit)
 (Intercept)           x1           x2 
260.55645444  -0.04276353   2.91272272 

# replace coefficients with new values, use whole name which is coefficients:
fit$coefficients[2:3] <- c(5, 1)

# view new coefficents
coef(fit)
(Intercept)          x1          x2 
260.5565      5.0000      1.0000 

问题在于您使用的是fit$coef，尽管lm输出中的组件名称确实是coefficients。缩写版本用于获取值，但不用于设置，因为它创建了名为coef的新组件，coef函数提取了fit$coefficient的值。