找回密码
 注册
查看: 2225|回复: 0

R语言

[复制链接]
发表于 2012-3-31 13:24:14 | 显示全部楼层 |阅读模式
求助: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的行,,
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

手机版|小黑屋|生物统计家园 网站价格

GMT+8, 2024-11-24 12:43 , Processed in 0.023538 second(s), 15 queries .

Powered by Discuz! X3.5

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表