溫馨提示×
您好，登錄后才能下訂單哦！
點擊 登錄注冊 即表示同意《億速云用戶服務條款》
您好，登錄后才能下訂單哦！
有三種方法 !is.na , na.omit, complete.cases
> d <- read.table("GWAS_s2.qassoc", header=T, stringsAsFactors=F)
// 文件行數
> nrow(d)
[1] 431493
> d1 <- subset(d, select=c("CHR", "SNP", "BP", "P"))
// 計算非NA 的行數
> num.bool <- complete.cases(d1)
> head(num.bool)
[1] FALSE TRUE TRUE FALSE TRUE TRUE
> sum(num.bool)
[1] 363836
> dn1 <- d1[which(!is.na(d1$P)),]
> nrow(dn1)
[1] 363836
> dn2 <- na.omit(d1)
> nrow(dn2)
[1] 363836
> dn3 <-d1[complete.cases(d1[,4]),]
> nrow(dn3)
[1] 363836
> dn4 <-d1[complete.cases(d1),]
> nrow(dn4)
[1] 363836方法三和方法四， 一個是根據第四列是否為NA判斷的， 一個是根據所有列。
免責聲明：本站發布的內容（圖片、視頻和文字）以原創、轉載和分享為主，文章觀點不代表本網站立場，如果涉及侵權請聯系站長郵箱：is@yisu.com進行舉報，并提供相關證據，一經查實，將立刻刪除涉嫌侵權內容。
