我有一個數據集,包含幾個時間序列,我想在這些時間序列內計算估計值以替換NA值。我熟悉ddply()和LM(),但我有一點與這個麻煩.... 這裏是前5個站從一個數據幀 x <- read.table(text="
LOC YEAR JAN FEB MAR APR MAY JUN JUL AUG SEP OCT NOV DEC
A 1980 0 0 0 8 104 NA 534 399 20
是否有一個單一的函數,類似於「runif」,「rnorm」等會產生線性模型的模擬預測?我可以自己編寫代碼,但代碼很難看,而且我認爲這是以前做過的事情。 slope = 1.5
intercept = 0
x = as.numeric(1:10)
e = rnorm(10, mean=0, sd = 1)
y = slope * x + intercept + e
fit = lm(y ~