def filter_pan_id(pan_id_file_name, length_table_file, filtered_pan_id_file_name): ''' input 1: pan_id_file_name input 2: length_table_file output 1: filtered_pan_id_file_name ''' R_code = ''' filter_pan_id=function(pan_id_file_name,length_table_file,filtered_pan_id_file_name){ require(dplyr) pan_id=read.table(pan_id_file_name) length_table=read.table(length_table_file) length_table_filter=length_table %>% filter(V2>20) pan_id_filter=merge(pan_id,length_table_filter,by.x = 2,by.y = 1,all.y = T) pan_id_filter=pan_id_filter[,-3] pan_id_filter=pan_id_filter[,c(2,1)] write.table(pan_id_filter,filtered_pan_id_file_name,sep = "\t",quote = F,row.names = F,col.names = F) } ''' R_filter_pan_id = SignatureTranslatedAnonymousPackage( R_code, "R_filter_pan_id") R_filter_pan_id.filter_pan_id(pan_id_file_name, length_table_file, filtered_pan_id_file_name)