2014-02-14 96 views
0

我已經從2個載體制成的多項式模型創建在A R數據集的變量:基於兩個矢量

tint <- c(12, 18.8, 25.5, 21.6) 
text <- c(11.4, 17.1, 22.1, 14.0) 
dT <- tint - text 

func.dT2 <- lm(dT[1:3] ~ poly(text[1:3],2,raw=TRUE)) 
現在

,我想使用該模型來從其他variavbels,其預測的dT數量遠大於3:

predict(func.dT2, text=seq(min(text)-1, max(text)+1, 0.01)) 

但我只獲得前三個值。

如何使用func.dT2從seq(min(text)-1,max(text)+1,0.01)中的所有值計算dT?

謝謝!

+0

希望你知道,配件到只有3分是有風險的,而且從相同的推斷是更危險? –

+0

感謝卡爾的評論。我的確很清楚,但這就是我擁有的所有數據,對於這種特殊情況,這將是合理的。不錯的建議,但記住的基礎是非常重要的;) –

回答

2

把你的數據在data.frames:

func.dT2 <- lm(dT ~ poly(text,2,raw=TRUE), 
       data=data.frame(dT, text)[1:3,]) 

predict(func.dT2, 
     newdata=data.frame(text=seq(min(text)-1, max(text)+1, 0.01))) 
+0

感謝羅蘭,工作非常好。 –