R語言行篩選的方法之filter函數(shù)詳解
更新時間：2022年07月01日 15:39:04 作者：育種數(shù)據(jù)分析之放飛自我
最近工作中要處理一大批數(shù)據(jù),所以這篇文章主要給大家介紹了關(guān)于R語言行篩選的方法之filter函數(shù)的相關(guān)資料,文中通過實例代碼介紹的非常詳細(xì),需要的朋友可以參考下
下面介紹一下R語言中行篩選的方法，主要介紹filter函數(shù)
1. 數(shù)據(jù)
這里，使用asreml分析中的BLUP值為例，相關(guān)的模型為：
m1 = asreml(Phen ~ G , random = ~ vm(Progeny,ainv) + vm(Dam,ainv) + vm(Progeny,dinv),
workspace = "10Gb", residual = ~ idv(units),data = dat)
summary(m1)$varcomp
計算育種值：
blup = coef(m1)$random head(blup) tail(blup)
數(shù)據(jù)特點：
- 沒有ID列，rownames的前綴為類型，比如
vm(Progeny, ainv)為加性效應(yīng)的BLUP值，vm(Progeny,dinv)為顯性效應(yīng)的BLUP值。
提取目的：
- 提取加性效應(yīng)的BLUP值，顯性效應(yīng)的BLUP值和母體效應(yīng)的BLUP
- 值提取BLUP值大于0.1的個體
2. 生成ID列和類型
首先，把rowname提取，作為新的一列
blup1 = blup %>% as.data.frame() %>% mutate(ID = rownames(.)) head(blup1)
根據(jù)下劃線，進行分列：
blup2 = blup1 %>% separate(ID,into = c("Type","IID"),sep = "_",remove = F)
head(blup2)
3. 提取effect大于0.1的行
re1 = blup2 %>% filter(effect>0.1) head(re1)
4. 提取加性效應(yīng)，且effect小于0的行
這里，條件之間，默認(rèn)是并集，如果想用交集，用|間隔。
re2 = blup2 %>% filter(Type == "vm(Progeny, ainv)",effect <0) head(re2)
5. 根據(jù)部分行名刪選
select函數(shù)，可以根據(jù)開頭，中間，結(jié)尾，進行列的刪選。
filter結(jié)合其它函數(shù)，也可以進行行的篩選。
如果想對ID中，包含ainv的行，進行篩選，可以這樣操作：
re3 = blup2 %>% filter(str_detect(ID,"ainv")) %>% arrange(-effect) head(re3)
注意，這里str_detect的pattern是正則表達式。如果直接用原始的字符：
re3 = blup2 %>% filter(str_detect(ID,"vm(Progeny, ainv)")) %>% arrange(-effect) head(re3)
可以看到，報錯，如果想要支持，需要對括號用兩個反斜線進行轉(zhuǎn)義。
轉(zhuǎn)義后的代碼：
re3 = blup2 %>% filter(str_detect(ID,"vm\\(Progeny, ainv\\)")) %>% arrange(-effect) head(re3)
6. 固定字符特征進行行篩選
str_detect沒有fixed = T的選項，如果想固定字符匹配，可以用fixed()函數(shù)：
re3 = blup2 %>% filter(str_detect(ID,fixed("vm(Progeny, ainv)"))) %>% arrange(-effect)
head(re3)
總結(jié)
到此這篇關(guān)于R語言行篩選的方法之filter函數(shù)詳解的文章就介紹到這了,更多相關(guān)R語言行篩選filter函數(shù)內(nèi)容請搜索腳本之家以前的文章或繼續(xù)瀏覽下面的相關(guān)文章希望大家以后多多支持腳本之家！
相關(guān)文章
R語言讀取柵格數(shù)據(jù)的方法(raster包讀取)
這篇文章主要介紹了R語言讀取柵格數(shù)據(jù)的方法,本文就只是對R語言raster包讀取、處理柵格數(shù)據(jù)加以基本的方法介紹，需要的朋友可以參考下2023-05-05
基于R語言時間序列的平穩(wěn)時間序列模型預(yù)測圖文詳解
時間序列是將統(tǒng)一統(tǒng)計值按照時間發(fā)生的先后順序來進行排列,時間序列分析的主要目的是根據(jù)已有數(shù)據(jù)對未來進行預(yù)測,下面這篇文章主要給大家介紹了基于R語言時間序列的平穩(wěn)時間序列模型預(yù)測的相關(guān)資料,需要的朋友可以參考下2022-12-12
這篇文章主要介紹了Java調(diào)用R語言實例講解，文中代碼講解的很清晰，有需要的同學(xué)可以研究下2021-03-03
這篇文章主要介紹了R語言時間序列中時間年、月、季、日的處理操作，具有很好的參考價值，希望對大家有所幫助。一起跟隨小編過來看看吧2021-03-03
這篇文章主要介紹了R語言讀取csv文件出錯的解決方案，具有很好的參考價值，希望對大家有所幫助。一起跟隨小編過來看看吧2021-04-04
這篇文章主要介紹了R語言導(dǎo)入CSV數(shù)據(jù)的簡單方法，文中通過示例代碼介紹的非常詳細(xì)，對大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價值，需要的朋友們下面隨著小編來一起學(xué)習(xí)學(xué)習(xí)吧2021-03-03
這篇文章主要介紹了R語言中Fisher判別的使用方法，文中2021-03-03
這篇文章主要介紹了R語言數(shù)據(jù)建模流程分析,本篇中包含了數(shù)據(jù)導(dǎo)入，清洗，可視化，特征工程，建模的代碼，大家可以選擇需要的去參考2022-01-01
這篇文章主要介紹了R語言ARMA模型的參數(shù)選擇說明，具有很好的參考價值，希望對大家有所幫助。一起跟隨小編過來看看吧2021-04-04
這篇文章主要介紹了R語言作圖:坐標(biāo)軸的設(shè)置方式，具有很好的參考價值，希望對大家有所幫助。一起跟隨小編過來看看吧2021-03-03