分享一下我的作法(下面)
我的csv数据第二个字段是dependent
Independent是1~length(dependent)的正整数数列
但是算出结果会出现这样讯息
有哪位知道这是什么意思吗?
[1] -1.505057
[1] 4.524962e-08
Warning messages:
1: 读取的项目数必须是行数的倍数
2: 读取的项目数必须是行数的倍数
R-script SOURCE
library(biglm)<br />
path<-choose.dir(default = "", caption = "Select data folder")<br />
ptm <- proc.time()<br />
filelist<-dir(path,full.names = TRUE)<br />
<br />
if(length(filelist) > 0){<br />
ff<-de~ID<br />
if( length(filelist) > 1 ){<br />
df <- scan(filelist[1], what = list("", double(0),double(0),double(0),double(0),double(0),double(0),double(0),double(0),""),sep=",")<br />
de<-df[[2]]<br />
rm(df)<br />
total<-length(de)<br />
ID<-1:total<br />
lmf<-data.frame(lmf=cbind(de,ID))<br />
res<-biglm(ff,lmf)<br />
print(coef(res)[[1]])<br />
print(coef(res)[[2]])<br />
<br />
idcount <- total<br />
for(i in 2:length(filelist))<br />
{<br />
df <- scan(filelist[i], what = list("", double(0),double(0),double(0),double(0),double(0),double(0),double(0),double(0),""),sep=",")<br />
de<-df[[2]]<br />
rm(df)<br />
total<-length(de)<br />
ID<-seq(idcount+1,idcount+total,by=1)<br />
idcount <-idcount+total<br />
lmf<-data.frame(lmf=cbind(de,ID))<br />
res<-update(res,lmf)<br />
print(filelist[i])<br />
print(i)<br />
print(coef(res)[[1]])<br />
print(coef(res)[[2]])<br />
} <br />
}else{<br />
df <- scan(filelist[1], what = list("", double(0),double(0),double(0),double(0),double(0),double(0),double(0),double(0),""),sep=",")<br />
de<-df[[2]]<br />
rm(df)<br />
total<-length(de)<br />
ID<-1:total<br />
lmf<-data.frame(lmf=cbind(de,ID))<br />
res<-biglm(ff,lmf)<br />
}<br />
<br />
print(coef(res)[[1]])<br />
print(coef(res)[[2]])<br />
}<br />
proc.time() - ptm