全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
3326 1
2013-08-07
各位大神好,我是新手一枚,刚用R语言没多久。我在做线性回归分析时,拟合关系式中含有分式,但是在输出的结果中却没有显示分式的拟合系数,下面是我的代码:
cp=read.table("clipboard",header=T)
fm=lm(cp~T+1/T+P+1/P+u+1/u+L+1/L+d+1/d+h+1/h+s+1/s,data=cp)
summary(fm)
anova(fm)
library(leaps)
varsel=regsubsets(cp~T+1/T+P+1/P+u+1/u+L+1/L+d+1/d+h+1/h+s+1/s,data=cp)
result=summary(varsel)
data.frame(result$outmat,RSS=result$rss,R2=result$rsq,adjR2=result$adjr2,Cp=result$cp,BIC=result$bic)
fm.step=step(fm,direction="forward")
fm.step=step(fm,direction="backward")
fm.step=step(fm,direction="both")

下面是我输出的结果:
Call:
lm(formula = cp ~ T + 1/T + P + 1/P + u + 1/u + L + 1/L + d +
    1/d + h + 1/h + s + 1/s, data = cp)
Residuals:
       Min         1Q     Median         3Q        Max
-11.525279  -1.062020   0.004545   1.076789   4.112983
Coefficients:
              Estimate Std. Error t value Pr(>|t|)   
(Intercept)  8.773e+03  5.170e+02   16.97   <2e-16 ***
T            2.023e+02  9.798e+00   20.65   <2e-16 ***
P            3.921e-01  1.680e-02   23.34   <2e-16 ***
u           -1.124e+02  6.513e+00  -17.26   <2e-16 ***
L            7.185e+00  9.142e-02   78.60   <2e-16 ***
d           -1.802e+04  8.957e+02  -20.12   <2e-16 ***
h           -1.659e+02  7.996e+00  -20.75   <2e-16 ***
s           -8.248e+03  4.460e+02  -18.50   <2e-16 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Residual standard error: 1.389 on 401 degrees of freedom
Multiple R-squared: 0.9954,     Adjusted R-squared: 0.9953
F-statistic: 1.233e+04 on 7 and 401 DF,  p-value: < 2.2e-16
下面是我部分数据
T1/TP1/Pcp1/cpu1/uL1/Ld1/dh1/hs1/s
170.150.0058771670.402872.4821921.1838240.844722.1260510.470356145.05880.0068941.5903160.62880671.810160.0139260.4147062.411347
170.650.0058599470.4259312.3477981.183850.8447022.0815760.480405144.75170.0069081.5889970.62932872.402080.0138120.418182.391316
171.150.0058428280.4501382.2215411.1838930.8446712.0387660.490493144.44520.0069231.5876780.62985172.994030.01370.4216432.371673
171.650.0058258080.4755392.1028791.1839530.8446281.9975270.500619144.13940.0069381.5863580.63037573.5860.013590.4250972.352404
172.150.0058088880.5021821.9913111.184030.8445741.9577710.510785143.83410.0069521.5850370.630974.178010.0134810.4285412.3335
172.650.0057920650.5301191.8863711.1841230.8445071.9194170.520992143.52940.0069671.5837160.63142674.770060.0133740.4319752.31495
173.150.0057753390.5594011.7876271.1842310.844431.8823910.531239143.22540.0069821.5823940.63195475.362160.0132690.4353992.296742
我想求指导,这是什么回事儿?怎么做才能在结果中显示分式的拟合系数啊?是不是需要加载程序包啊?






二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

全部回复
2013-8-7 17:10:48
我觉得在回归拟合中,R并没有把1/T、1/P等作为一个预测变量,你把1/T 改成T1,1/P改P1,其它几个做同样的改变就应该可以了,你试试
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群