# for the first 15 colomns ========== dat=read.table("data1.txt",head=TRUE) na=colnames(dat) dat1=dat[,na[-1]] nc=dim(dat1)[2] # 列数 nr=dim(dat1)[1] # 行数 logr=-dat1