R回归分析

1.数据

> data
   YEAR SCALE     PRICE ACREAGE
  2006 82.04   3846.883 1087.25
  2007 82.17   4304.399 1229.74
  2008 82.23   5043.934   811.69
  2009  82.53   5054.121 1265.16
  2010 83.78   5987.330 1206.71
  2011 84.13   8693.126 1031.31
  2012 85.02 11008.840   749.28
  2013 85.27 12114.270 1016.71
  2014 85.43 11793.917   818.50
  2015 85.53 13709.240 1044.89

绘制基本的图层以及散点图

PL=ggplot(data=data,aes(x=YEAR,y=SCALE,group=1))+geom_point()+theme(axis.text.x = element_text(angle=90, hjust=1, vjust=.5))+scale_x_continuous(breaks=seq(2006,2016,1))

在散点图的基础上加上回归线

loess型回归线

PL+stat_smooth(methon='loess')





### R语言回归分析示例教程 回归分析是统计学和数据科学中一种重要的建模方法,用于研究变量之间的关系并进行预测。以下是基于R语言的回归分析示例教程,涵盖了简单线性回归、多项式回归、多元线性回归以及逻辑回归。 #### 1. 简单线性回归 简单线性回归模型旨在揭示一个自变量与一个因变量之间的线性关系。在R语言中,可以使用`lm()`函数实现。 ```r # 示例:简单线性回归 data(mtcars) model <- lm(mpg ~ wt, data = mtcars) # 构建模型 mpg 是因变量,wt 是自变量 summary(model) # 查看模型摘要 ``` 此代码构建了一个以`mpg`为因变量、`wt`为自变量的简单线性回归模型[^2]。 #### 2. 多项式回归 当变量之间的关系不是严格的线性时,可以尝试多项式回归。以下是一个二次多项式回归的例子。 ```r # 示例:多项式回归 data(mtcars) model_poly <- lm(mpg ~ poly(wt, 2), data = mtcars) # 使用 poly() 函数创建多项式项 summary(model_poly) # 查看模型摘要 ``` 通过`poly()`函数,可以轻松地生成多项式项,并将其纳入回归模型中[^2]。 #### 3. 多元线性回归 多元线性回归扩展了简单线性回归的概念,允许同时考虑多个自变量对因变量的影响。 ```r # 示例:多元线性回归 data(mtcars) model_multi <- lm(mpg ~ wt + hp, data = mtcars) # 考虑 wt 和 hp 两个自变量 summary(model_multi) # 查看模型摘要 ``` 此模型评估了`wt`(车辆重量)和`hp`(马力)对`mpg`(每加仑行驶英里数)的联合影响[^2]。 #### 4. 逻辑回归 逻辑回归是一种分类模型,特别适用于二分类问题。以下是使用R语言进行逻辑回归的示例。 ```r # 示例:逻辑回归 library(epiDisplay) data("Ectopic") model_logit <- glm(Pregnant ~ Age + Education, data = Ectopic, family = "binomial") # 构建逻辑回归模型 summary(model_logit) # 查看模型摘要 ``` 在此示例中,`glm()`函数被用来构建逻辑回归模型,`family = "binomial"`指定了模型类型为二项分布[^4]。 #### 5. 模型评估与比较 为了评估不同回归模型的效果,可以使用测试误差或交叉验证等方法。 ```r # 示例:K折交叉验证 library(caret) set.seed(123) train_control <- trainControl(method = "cv", number = 10) # 设置10折交叉验证 model_cv <- train(mpg ~ wt + hp, data = mtcars, method = "lm", trControl = train_control) print(model_cv) # 输出交叉验证结果 ``` 上述代码展示了如何使用`caret`包中的`trainControl()`函数来设置K折交叉验证,并评估模型性能。 --- ###
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值