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)
