|
求助:y <-scan(file="E:/gt_T2D_ctrl_chr01.txt",list(""),nlines=1)
y<-y[[1]][4:4865]
y<-as.numeric(y)
Name<-character()
pr<-numeric()
for (n in 1:37500)
{x <-scan(file="E:/gt_T2D_ctrl_chr01.txt",list(""),nlines=1,skip = n)
Name[n]<-x[[1]][1]
x<-x[[1]][4:4865]
x<-as.numeric(x)
test<-data.frame(x,y)
mylogit <- glm(y~x,data=test,family = binomial(link = "logit"))
summary(mylogit)
pr[n]<-summary(mylogit)$coefficients[2,4]
Result <- cbind(Name,pr) }
如果这个代码执行到一行SNP的基因型全为0,那么它将会显示出错,终止执行程序。求助,如何在这个代码里嵌入if else 语句,使它能够执行到文件最后。 或者能够删除源文件中全为0的行,, |
|