有一个数据框，第一列是探针ID，第二列是探针的碱基序列，第三列是芯片平台，模拟数据代码如下；
options(stringsAsFactors = F)
set.seed(123)
ns=sample(1:1000,10)
seqs = do.call(rbind,lapply(ns, function(n){
# n=ns[1]
id=paste0('id',1:n)
seq=rep('aaacccgggtttcccggaaa',n)
gpl=paste0('gpl',n)
df=data.frame(id,seq,gpl)
return(df)
}))
unique(seqs$gpl)
# 这个 seqs 数据框里面，有10个gpl平台，需要拆分成为10个单独的文件
# 把序列fasta格式化
x=as.vector(seqs[1, ])
paste0('>',x[1],'\n',x[2])
# 这个就是fasta序列。
有一个数据框，第一列是探针ID，第二列是探针的碱基序列，第三列是芯片平台，模拟数据代码如下；
下载芯片探针序列并且写成fasta文件 选择在GEO官网的GPL平台下载 : https://www.ncbi....
将一个或多个序列写入FASTA格式的文件中，包括核酸序列和氨基酸序列。在R语言中，需要seqinr package...
有需求将多个单个基因序列.fasta文件合并到一个多序列的.fasta文件里，方便后续的序列比对和进化分析 采用W...
将gbk格式文件转换成fasta格式文件 将多个fasta序列中提取含有'Homo sapiens'的序列 请关注...
有一个蛋白质序列文件： 需要按蛋白将他们分成多个fasta文件。 修改目录和文件即可。
①下载TXT文件和平台文件 若数据平台文件无gene symbol 和探针序列or GBlist，则无法注释。 ②...
fasta文件格式：fasta 文件为一个 ID 对应一个 序列,可以是 转录本序列, 蛋白序列 ''' V350...
如何根据芯片的探针ID找到基因名-基于R语言 使用bioconductor注释包 如果该芯片平台有对应的bioco...
后续代码准备文件夹建立primers.fasta文件，保存接头序列 运行pb的去接头软件lima,（不是基因芯片的...
本文标题：一个R考核题-多个芯片平台的探针序列输出到fasta文件
本文链接：https://www.haomeiwen.com/subject/qpyjzctx.html
网友评论