r语言数据变量分段_【最新】R语言分段回归 数据分析 案例报告(附代码数据)

R

语言分段回归数据数据分析案例报告

#

读取数据

data=

read.csv

(

"artificial-cover.csv"

)

#

查看部分数据

head

(data)

##   tree.cover shurb.grass.cover

## 1       13.2              16.8

## 2       17.2              21.8

## 3       45.4              48.8

## 4       53.6              58.7

## 5       58.5              55.5

## 6       63.3              47.2

#######

先调用

spline

包

library

( splines )

###########

用

lm

拟合，主要注意部分是

bs

(

age,knots=c(...)

)这部分把自变量分

成不同部分

fit =

lm

(tree.cover~

bs

(shurb.grass.cover ,

knots =

c

(

25

,

40

,

60

) ),

data=

da

ta )

############

进行预测，预测数据也要分区

pred=

predict

(fit ,

newdata =

list

(

shurb.grass.cover =

data$shurb.grass.c

over),

se=

T)

#############

然后画图

plot

(fit)