回归分析
大白羊_Aries
这个作者很懒,什么都没留下…
展开
-
一元线性回归--R实现
R读取数据方法一: 窗口输入 c()为用来创建,矩阵为二维数组,可以通过函数matrix()来创建,而多维数组使用array()函数来创建。 mean()和sd()分别计算变量均值和标准差的函数。方法二: 使用read.table(file,head=,sep=“deliniter”) 其中,file是一个带分割符的文本文件,如.txt和.csv文件;head的取值为TRUE或F...原创 2018-12-23 22:20:59 · 7686 阅读 · 2 评论 -
多项式回归(R语言)
这里,回归采用逐个引入自变量的方式,由此可以清楚次看到各项对回归的贡献,使显著性检验更加明确。依次引入自变量x1,x2,x12,x22,x1x2x_{1},x_{2},x_{1}^{2},x_{2}^{2},x_{1}x_{2}x1,x2,x12,x22,x1x2以查看各变量对回归的贡献。代码实现如下:data9.2<-read.csv("C:/Users/Administr...原创 2019-01-23 15:49:10 · 7943 阅读 · 0 评论 -
可视为线性回归的曲线回归(R语言)
用R软件进行计算,说闲话出原始数据对自变量的散点图,代码实现如下:data9.1<-read.csv("C:/Users/Administrator/Desktop/data9.1.csv",head=TRUE) attach(data9.1) plot(t,y)...原创 2019-01-23 14:41:27 · 3924 阅读 · 2 评论 -
数据中心化、标准化(R语言)
详细参考下面博客https://www.cnblogs.com/ljhdo/archive/2018/08/09/4899086.html转载 2019-01-23 14:29:35 · 3152 阅读 · 0 评论 -
偏最小二乘法(R语言)
偏最小二乘一般针对于下面两个问题:1)自变量的数目较多,但是观测的样本并不是很多;2)回归方程建立后主要用于预测。代码实现如下:data8.2<-read.csv("C:/Users/Administrator/Desktop/data8.2.csv",head=TRUE)datas<-data.frame(scale(data8.2))library(pls)pls1&lt...原创 2019-01-23 14:01:35 · 16754 阅读 · 10 评论 -
主成分分析(R语言)
代码实现如下:data3.3<-read.csv("C:/Users/Administrator/Desktop/data3.3.csv",head=TRUE)datas<-data.frame(scale(data3.3[,1:6])) pr3.3<-princomp(~x1+x2+x3+x4+x5,datas,cor=T)# 对5个变量做主成分分析,其中cor=T表...原创 2019-01-22 19:11:58 · 8777 阅读 · 2 评论 -
岭回归(R语言)
代码实现如下:data3.3<-read.csv("C:/Users/Administrator/Desktop/data3.3.csv",head=TRUE)datas<-data.frame(scale(data3.3[,1:6])) # 对样本数据进行标准化处理并转换为数据框的存储格式library(MASS)ridge3.3<-lm.ridge(y~.-1,da...原创 2019-01-22 15:53:22 · 21018 阅读 · 5 评论 -
多重共线性的诊断(R语言)
多重共线性的诊断1.方差扩大因子法经验表明,VIFj≥10VIF_{j}\geq10VIFj≥10时,就说明自变量xjx_{j}xj与其余变量之间有严重的多重共线性,且这种多重共线性可能会过度地影响最小二乘估计值。代码实现如下:data5.5&amp;lt;-read.csv(&quot;C:/Users/Administrator/Desktop/data5.5.csv&quot;,head=TRUE)lm5...原创 2019-01-22 14:48:22 · 32898 阅读 · 0 评论 -
逐步回归(R语言)
R软件提供了非常方便地进行逐步回归分析的计算函数step(),它是以AIC1.前进法代码实现如下:data3.1&amp;lt;-read.csv(&quot;C:/Users/Administrator/Desktop/data3.1.csv&quot;,head=TRUE)lmo3.1&amp;lt;-lm(y~1,data=data3.1)lm3.1.for&amp;lt;-step(lmo3.1,sc原创 2019-01-22 14:02:04 · 54783 阅读 · 6 评论 -
自变量选择(R语言)
1.自有度调整复决定系数Ra2R_{a}^{2}Ra2设Ra2R_{a}^{2}Ra2为调整的复决定系数,nnn为样本量,ppp为自变量的个数,则 Ra2=1−n−1n−p−1(1−R2)R_{a}^{2}=1-\frac{n-1}{n-p-1}(1-R^{2})Ra2=1−n−p−1n−1(1−R2)在一个实际问题的回归建模中,自由度调整复决定系数Ra2R_{a}^{2...原创 2019-01-20 22:15:58 · 7510 阅读 · 0 评论 -
异常点与强影响点(R语言)
1.关于因变量y的异常点我们可以直接用R软件的rstudent()函数计算出删除学生化残差SRE(i)SRE_{(i)}SRE(i)的数值,∣SRE(i)∣&gt;3\left | SRE_{(i)} \right |&gt;3∣∣SRE(i)∣∣>3的观测值即判定为异常值。代码实现如下:data2.2<-read.csv("C:/Users/Admini...原创 2019-01-20 16:34:40 · 9393 阅读 · 0 评论 -
BOX-COX变换(R语言)
1.消除异方差代码实现如下:data3.2&lt;-read.csv("C:/Users/Administrator/Desktop/data3.2.csv",head=TRUE)library(MASS)bc3.2&lt;-boxcox(y~x1+x2,data=data3.2,lambda=seq(-2,2,0.01))# λ的取值为区间[-2,2]上步长为0.01的值,bc3.2中...原创 2019-01-20 15:48:36 · 16031 阅读 · 6 评论 -
违背基本假设的几种情况——自相关性(R语言)
自相关性的诊断DW检验DW检验的代码实现如下:data2.2&lt;-read.csv("C:/Users/Administrator/Desktop/data2.2.csv",head=TRUE)lm2.2&lt;-lm(y~x,data=data2.2)library(lmtest)dwtest(lm2.2,alternative="two.sided")输出结果如下:原创 2019-01-20 15:01:42 · 5756 阅读 · 0 评论 -
违背基本假设的几种情况——异方差性(R语言)
在建立实际问题的回归模型时,经常存在于此假设想违背的情况,一种是计量经济模型中常说的异方差性,即 var(εi)≠var(εj),当i≠j时var(\varepsilon _{i})\neq var(\varepsilon _{j}),当i\neq j时var(εi)̸=var(εj),当i̸=j时诊断方法1.残差图分析法data4.3&amp;amp;lt;-read.csv(&amp;quot;C:/User...原创 2019-01-20 14:20:56 · 25628 阅读 · 11 评论 -
多元线性回归--R实现
原创 2019-01-19 20:15:47 · 9168 阅读 · 1 评论 -
因变量顺序类别的回归(R语言)
本例使用MASS包中的housing数据集,该数据集是关于哥本哈根住房情况的调查数据。其中包括5个变量,分别为:房主对他们目前住房的满意度(高、中、低),记为Sat是有序变量;房主认为物业管理的影响程度(高、中、低),记为Infl;租赁住房的类型(塔式、中庭、公寓、露宿),记为Type;与其他住户的沟通程度(低、高),记为Cont;每组对应的居民人数,记为Freq,其中共3∗3∗4∗2=723*3...原创 2020-01-02 11:14:40 · 7546 阅读 · 0 评论 -
Logistic回归模型(R语言)
1.分组数据的Logistic回归模型下面我们以一道例题来说明,R软件中实现分组数据的logistics回归模型:代码实现如下:data10.4<-read.csv("C:/Users/Administrator/Desktop/data10.4.csv",head=TRUE)# data10.4中保留的p1变量为逻辑变换后的变量lm10.4<-lm(p1~x,weigh...原创 2019-01-24 16:04:19 · 54511 阅读 · 5 评论 -
非线性模型(R语言)
实例一首先,绘制出原数据的散点图。代码试下如下:data9.3<-read.csv("C:/Users/Administrator/Desktop/data9.3.csv",head=TRUE)attach(data9.3)plot(x,y)输出结果为: 可以看出,这时yyy与xxx之间呈现出非线性,因此需要对数据进行非线性回归分析。代码实现如下:nls9.3<-...原创 2019-01-23 16:50:36 · 10628 阅读 · 1 评论