2013-08-27 86 views
0

我試圖將三次曲線擬合到我的散點圖中。我能夠在minitab中做到這一點,沒有問題,但我發現很難將立方非線性迴歸擬合到我的數據中。我的數據設置如下:R中的立方迴歸

AGE Value 
3 10 
4 10 
5 11 
5 13 
6 10 
7 9 
8 8 

回答

2

?poly信息(默認)在x值的範圍規定開度的設計矩陣就緒正交多項式建立你的手吧:

df <- read.table(text="AGE Value 
3 10 
4 10 
5 11 
5 13 
6 10 
7 9 
8 8", header=TRUE) 

lm(Value ~ poly(AGE, degree=3), data=df) 

# Call: 
# lm(formula = Value ~ poly(AGE, 3), data = df) 
# 
# Coefficients: 
# (Intercept) poly(AGE, 3)1 poly(AGE, 3)2 poly(AGE, 3)3 
#  10.1429  -2.0026  -2.3908   0.6019 
+0

或使用'polyfit' tho'授予,不返回有關適合質量的信息。 –

+0

或者,如果您使用的是(B-)樣條擬合,則可以在公式中使用'bs(AGE,degree = 3)'(來自** splines **包),將總數'df' [作爲demo'd在這裏](http://stackoverflow.com/questions/12925533/quadratic-spline/12938310#12938310)。 –

+0

@CarlWitthoft - 我在哪裏可以找到'polyfit'? –