r软件求解线性方程组（r软件）

你们好，最近小活发现有诸多的小伙伴们对于r软件求解线性方程组，r软件这个问题都颇为感兴趣的，今天小活为大家梳理了下，一起往下看看吧。

1、打开R软件，文件>新建脚本程序，输入需要进行回归分析的数据。

2、首先绘制出散点图，观察二者的关系是线性的或者非线性的，以建立合适的回归模型。

3、程序代码：

4、#输入数据

5、x=c(1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24)

6、y=c(276,277,283,294,308,330,355,387,422,461,507,558,613,675,741,811,886,966,1053,1142,1239,1340,1445,1556)

7、#绘制散点图

8、plot(x,y,col=1:7)

9、其中plot为绘图命令，col表示散点的颜色。

10、运行后得到下图。

11、从图形可知，y与x成非线性关系，图形比较契合二次曲线，建立回归模型y=a+b*x+c*(x^2)

12、计算回归参数：

13、#回归分析

14、test=lm(y~1+x+I(x^2))

15、test

16、其中lm表示建立模型，其中的公式y=1+x+I(x^2)，1表示常数项，I（x^2）表示x的平方项，注意这里合起来是相当于组成另外一个变量，单独x^2不能得到预想的结果。

17、运行得到参数，模型为：

18、y=280.944-6.909*x +2.501*(x^2)

19、查看显著性检验结果。

20、运行代码summary(test)

21、表示提取模型的计算结果。

22、该结果中，call显示出模型公式，residuals显示出残差的四分位数值，estimate这一列反映了模型的参数估计值，std.Error为标准差，Pr(>|t|)这一列位计算的P值，最下面是解释。可以看出该模型的每个参数的P值的显著性标记均为“***”，表示极为显著。该模型通过回归检验。

23、利用predict()函数求该回归模型的预测值和预测区间。

24、程序代码为：

25、#预测

26、pred=data.frame(x=25)

27、test.pred=predict(test,pred,interval="prediction",level=0.95)

28、test.pred

29、其中pred表示要预测的点，这里必须以数据框的形式输入，predict函数才能计算；interval="prediction"表示给出对应的预测区间，参数level=0.95表示显著性水平。由此得到的预测结果为：Y'=1671.368 ，区间[1669.628,1673.109]。

30、计算全部拟合值，并绘制拟合曲线图。

31、#计算全部拟合值

32、nihe=predict(test)

33、#曲线拟合图

34、plot(x,y,col=1:7)

35、lines(nihe)

36、其中predict计算出全部拟合值。在绘制曲线命令lines()之前，需要先用plot()绘制出散点图。

37、从图中可以看出，拟合效果非常好。

38、在实际数据分析中，还会涉及多元线性回归方程、非线性回归模型等。可以使用命令lm(y~X1+X2+X3+I(X3^2)).....等等建立合适的模型，以得到更好的拟合效果。

以上就是r软件这篇文章的一些介绍，希望对大家有所帮助。

免责声明：本文由用户上传，如有侵权请联系删除！