各位大神好,我是新手一枚,刚用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
下面是我部分数据
| T | 1/T | P | 1/P | cp | 1/cp | u | 1/u | L | 1/L | d | 1/d | h | 1/h | s | 1/s |
| 170.15 | 0.005877167 | 0.40287 | 2.482192 | 1.183824 | 0.84472 | 2.126051 | 0.470356 | 145.0588 | 0.006894 | 1.590316 | 0.628806 | 71.81016 | 0.013926 | 0.414706 | 2.411347 |
| 170.65 | 0.005859947 | 0.425931 | 2.347798 | 1.18385 | 0.844702 | 2.081576 | 0.480405 | 144.7517 | 0.006908 | 1.588997 | 0.629328 | 72.40208 | 0.013812 | 0.41818 | 2.391316 |
| 171.15 | 0.005842828 | 0.450138 | 2.221541 | 1.183893 | 0.844671 | 2.038766 | 0.490493 | 144.4452 | 0.006923 | 1.587678 | 0.629851 | 72.99403 | 0.0137 | 0.421643 | 2.371673 |
| 171.65 | 0.005825808 | 0.475539 | 2.102879 | 1.183953 | 0.844628 | 1.997527 | 0.500619 | 144.1394 | 0.006938 | 1.586358 | 0.630375 | 73.586 | 0.01359 | 0.425097 | 2.352404 |
| 172.15 | 0.005808888 | 0.502182 | 1.991311 | 1.18403 | 0.844574 | 1.957771 | 0.510785 | 143.8341 | 0.006952 | 1.585037 | 0.6309 | 74.17801 | 0.013481 | 0.428541 | 2.3335 |
| 172.65 | 0.005792065 | 0.530119 | 1.886371 | 1.184123 | 0.844507 | 1.919417 | 0.520992 | 143.5294 | 0.006967 | 1.583716 | 0.631426 | 74.77006 | 0.013374 | 0.431975 | 2.31495 |
| 173.15 | 0.005775339 | 0.559401 | 1.787627 | 1.184231 | 0.84443 | 1.882391 | 0.531239 | 143.2254 | 0.006982 | 1.582394 | 0.631954 | 75.36216 | 0.013269 | 0.435399 | 2.296742 |
我想求指导,这是什么回事儿?怎么做才能在结果中显示分式的拟合系数啊?是不是需要加载程序包啊?