问我无法定义我的fasta文件，我在R中输出了fasta文件，在终端中使用了"makeblast“命令

EN

library(seqinr)
library(Biostrings)
library(data.table)
library(tidyverse)



#read homo_tabular format
homo_tab = read.csv("proteins_homo.csv", header = TRUE, sep = ",")

homo_tab_1 = homo_tab[,c(7,9:11)]
colnames(homo_tab_1)[2]="ID"


#select longest locus
 son <- homo_tab_1 %>%
  group_by(Locus) %>%
  slice_max(Length, n = 1) %>%
  slice_head(n = 1)


 #read homo protein fasta ile and convert list to df/dt
human_prot <- read.fasta(file= "homo_s.faa", seqtype="AA", as.string =TRUE, set.attributes =TRUE)

human_prot = unlist(human_prot)
human_prot = as.data.frame(human_prot)
human_prot = setDT(human_prot, keep.rownames = "ID")
#rename column
colnames(human_prot)[1] ="ID" 
colnames(human_prot)[2] ="seq"

#merge csv and fasta file
merged = merge(human_prot , son , by="ID", all.x=TRUE)

#remove na rows
library(dplyr)

merged_1 <- na.omit(merged)

#delete column
merged_2 = subset(merged_1, select = -c(3,4,5) )



write.fasta(sequences = merged_2, names = names(merged_3), file.out = "homo_ref.faa")