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))
