def remove_MGG_unpresent_Augustus( pav_orthofinder, MGG_unpresent_Augustus_unassianed_list_file_name, pav_orthofinder_1574): ''' 用于从pav_orthofinder中删掉1574个出现在unassigned gene中的 input 1: pav_orthofinder input 2: MGG_unpresent_Augustus_unassianed_list output 1: pav_orthofinder_1574 ''' R_code_remove_MGG_unpresent_Augustus = ''' R_remove_MGG_unpresent_Augustus=function( pav_orthofinder_file_name, MGG_unpresent_Augustus_unassianed_list_file_name, pav_orthofinder_1574_file_name ){ require(readxl) require(WriteXLS) require(dplyr) MGG_unpresent_Augustus_unassianed_list=read.table(MGG_unpresent_Augustus_unassianed_list_file_name) pav_orthofinder=read_xlsx(pav_orthofinder_file_name) pav_orthofinder_1574=pav_orthofinder %>% filter(!(protein_id %in% MGG_unpresent_Augustus_unassianed_list$V1)) pav_orthofinder_1574=pav_orthofinder_1574[,-1] pav_orthofinder_1574=pav_orthofinder_1574[,c(1,158,2:157)] WriteXLS::WriteXLS(pav_orthofinder_1574,pav_orthofinder_1574_file_name) } ''' R_remove_MGG_unpresent_Augustus = SignatureTranslatedAnonymousPackage( R_code_remove_MGG_unpresent_Augustus, "R_remove_MGG_unpresent_Augustus") R_remove_MGG_unpresent_Augustus.R_remove_MGG_unpresent_Augustus( str(pav_orthofinder), str(MGG_unpresent_Augustus_unassianed_list_file_name), str(pav_orthofinder_1574))