R操作建立回归模型

《R操作建立回归模型》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、R操作建立回归模型操作一：载入数据创建一个记事本，输入相关数据，并以"toothpaste.txt"命名，保存在桌面上。修改路径，在R软件中打开file---change.dir,修改在上述已经保存的路径中。tp=read.table("toothpaste.txt",header=T);tp价格均价广告价差销量13.853.805.50-0.057.3823.754.006.750.258.5133.704.307.250.609.5243.703.705.500.007.5053.603.857.000.259.3363.6

2、03.806.500.208.2873.603.756.750.158.7583.803.855.250.057.8793.803.655.25-0.157.10103.854.006.000.158.00113.904.106.500.207.89123.904.006.250.108.15133.704.107.000.409.10143.754.206.900.458.86153.754.106.800.358.90163.804.106.800.308.87173.704.207.100.509.26183.804.30

3、7.000.509.00193.704.106.800.408.75203.803.756.50-0.057.95213.803.756.25-0.057.65223.753.656.00-0.107.27233.703.906.500.208.00243.553.657.000.108.50253.604.106.800.508.75263.654.256.800.609.21273.703.656.50-0.058.27283.753.755.750.007.67293.803.855.800.057.93303.704.2

4、56.800.559.26操作二：分析数据>x1=tp[,4];>x2=tp[,3];>y=tp[,5];>plot(x1,y)abline(lm(y~x1))#y~x1的一次回归曲线操作三：建立模型3.1plot(x2,y)#画出散点图，建立模型sol=lm(y~x1+x2+I(x2^2));summary(sol)Call:lm(formula=y~x1+x2+I(x2^2))Residuals:Min1QMedian3QMax-0.40330-0.14509-0.030350.154880.46602Coefficient

5、s:EstimateStd.ErrortvaluePr(>

6、t

7、)(Intercept)17.32445.64153.0710.00495**x11.30700.30364.3050.00021***x2-3.69561.8503-1.9970.05635.I(x2^2)0.34860.15122.3060.02934*---Signif.codes:0‘***’0.001‘**’0.01‘*’0.05‘.’0.1‘’1Residualstandarderror:0.2213on26degreesoffreedomMultipl

8、eR-squared:0.9054,（你建立的模型对销售量的贡献）AdjustedR-squared:0.8945F-statistic:82.94on3and26DF,p-value:1.944e-133.2sol1=lm(y~x1+x2+I(x1*x2)+I(x2^2));summary(sol1)Call:lm(formula=y~x1+x2+I(x1*x2)+I(x2^2))Residuals:Min1QMedian3QMax-0.43725-0.117540.004890.122630.38410Coefficient

9、s:EstimateStd.ErrortvaluePr(>

10、t

11、)(Intercept)29.11337.48323.8900.000656***x111.13424.44592.5040.019153*x2-7.60802.4691-3.0810.004963**I(x1*x2)-1.47770.6672-2.2150.036105*I(x2^2)0.67120.20273.3120.002824**---Signif.codes:0‘***’0.001‘**’0.01‘*’0.05‘.’0.1‘’1Residualstand

12、arderror:0.2063on25degreesoffreedomMultipleR-squared:0.9209,AdjustedR-squared:0.9083F-statistic:72.78（稳定量）on4and25DF,p-value:2.107e