> x<-read.table("a.txt") > names(x) [1] "V1" "V2" > length(x) [1] 2 > length(V1) 错误: 找不到这个对象"V1" > length(x$V1) [1] 2050 > length(x$V2) [1] 2050 > plot(x$V1,x$V2) > plot(x$V1,x$V2,type="l") > yx<-lm(x$V2~x$V1) > summary(yx) Call: lm(formula = x$V2 ~ x$V1) Residuals: Min 1Q Median 3Q Max -171.85 -50.78 -12.08 47.70 218.90 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) -252.9587 5.6706 -44.61 <2e-16 *** x$V1 0.4014 0.0025 160.53 <2e-16 *** --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Residual standard error: 66.99 on 2048 degrees of freedom Multiple R-squared: 0.9264, Adjusted R-squared: 0.9263 F-statistic: 2.577e+04 on 1 and 2048 DF, p-value: < 2.2e-16 > plot(x$V2~x$V1);abline(yx,col= "blue")