前言
(nature重磅综述|关于rna-seq你想知道的全在这)()()()ggseqlogo,只要你根据此包要求的数据格式上传一堆dna序列或者氨基酸序列,再根据现成的命令流程就能画出logo图。
安装到作图的代码如下:
安装方式有两种
从github中安装devtools::install.github("omarwagih/ggseqlogo")seqs_aahead(seqs_aa) $akt1 "vvgarrsswrvvssi" "gprsrsrsrdrrrke" "llclrrsslkaygng" "terprpntfiirclq" "lsrervfsedrarfy" "pstsrrfsppssslq"也可以用自己的数据集,支持两种格式,矩阵。
长度为5的motif矩阵示例,每一列代表一个位置,及碱基在该位置的出现次数。共4行,每一行代表一种碱基matrix - "base 1 2 3 4 5a 10 2 0 8 1c 1 12 1 2 3g 4 0 9 1 1t 0 0 0 1 9"matrix_input - read.table(matrix, header=t, row.names=1)matrix_input - as.matrix(matrix_input)ggplot()+geom_logo(seqs_dna$ma0001.1)+theme_logo()
ggseqlogo支持以下几种类型数据输入:
下面是使用数据中的位置频率矩阵生成的seqlogo
ggseqlogo(pfms_dna$ma0018.2)

ggseqlogo提供函数make_col_scheme
来自定义离散或者连续配色方案csl - make_col_scheme(chars = c("a","t", "c", "g"), groups = c("gr1","gr1", "gr2","gr2"), cols = c("purple","purple","blue","blue"))ggseqlogo(seqs_dna$ma0001.1,col_scheme=csl)
cs2 - make_col_scheme(chars = c("a", "t", "c", "g"), values = 1:4)ggseqlogo(seqs_dna$ma0001.1, col_scheme=cs2)
ggseqlogo(seqs_dna, ncol = 4)
上述命令实际上等同于
ggplot()+geom_logo(seqs_dna)+theme_logo()+ facet_wrap(~seq_group,ncol = 4,scales = "free_x")通过创建矩阵可以生成每个标志的高度,还可以有负值高度
set.seed(1234)custom_mat - matrix(rnorm(20), nrow = 4, dimnames = list(c("a","t","c", "g")))ggseqlogo(custom_mat,method="custom",seq_type="dna")+ ylab("my custom height")
可以通过font
参数来设置字体,具体可参考?list_fontsfonts - list_fonts(f)p_list - lapply(fonts, function(f){ ggseqlogo(seqs_dna$ma0001.1,font=f)+ggtitle(f)})do.call(gridextra::grid.arrange,c(p_list, ncol=4))
注释的话跟ggseqlogo生成的图形与阅读原文跳转其博客。
更多数据可视化可参考在线工具:
1、路易·威登lv (1854年法国,世界品牌500强,箱包和皮具领域的全球第一品牌,)lv 经典款2、香奈儿chanel (于1913年法国,世界品牌500强,世界知名奢侈品品牌,世界十大皮具品牌,法国香奈儿(chanel)公司)channel经典款3、古琦欧·古琦gucci (于1923年意大利,意大利最大的时装集团之一,世界十大皮具品牌,世界著名品牌,意大利gucci公司)
快递已经慢慢开始恢复了,爱美的小仙女们最想买的春日单品是什么呢?对于在家憋了这么久的我们来说,每一次出门都要盛装出席,用最精致的样子去迎接第一缕春光。俗话说得好:“包”治百病。衣服可以少买一件,但是小包包可不能少,平时能装装手机、化妆品,还能成为搭配小心机,所以包包就成为了复工后必买的第一单品。编辑在时髦商店“女王节”大促来临之际,给大家挑选了一些非常好搭配的小众包包,能充分凸显个人风格不说,而且
下班后摆地摊卖什么赚钱?这是一位网友问我的问题,相信这也是很多朋友都想了解的情况,那么今天我们就来聊聊这个话题,希望能帮到大家,也希望大家帮忙转发一下,谢谢大家!最近几年由于口罩的影响,全国各地的经济都大受影响,各地也开始鼓励人们,允许居民下班后摆个地摊创业了,于是很多网友都开始利用下班时间去摆起地摊赚钱了。
最近很多朋友都在问我,今年秋冬流行什么包?街拍博主和秀场模特们仿佛在按头安利尺寸超大的包包,女明星们也纷纷提上了复古手提袋。那么,今年秋冬哪些新包值得买?我看上了哪些包?今天这篇大家可以参考~(一起攒钱)1.saint laurent「icare」今年秋冬包包关键词:黑金。好多品牌都在出!近期最火的,应该就是ysl这只超大购物袋了吧?
入圈指南:搜索推荐:头条搜索缝纫之家,发现更多精彩缝纫视频!搜更多精彩内容来缝纫之家学缝纫技巧牛仔裤改造成随身小包包,这个旧物改造的小妙招你一定不能错过!01材料准备牛仔裤1条(拆下兜,标牌及裤腿部分),30cm花边1条。02缝纫制作裁正1片宽17cm,高16cm的布片,将布边向内折2cm,正面固定好花边,车明线;底边向上4cm处居中固定车缝标牌。
适合新手的包包教程又来了,这个系列的包包教程钩法相对都比较简单,也是入门级的钩包教程,用到的基本都是最基础的针法,没什么繁琐的花样,但钩出来的包包外形还是很靓丽的,所以说,钩包最重要的还是包包的型,并不是说钩得多复杂就好看。而且有很多人就喜欢朴素大方的包包,搭配衣服也不显突兀,太花哨的有时反而不好搭配。
你以为我的手柄只是边油坏了,其实我还在偷偷地笑。你以为我没法修好了,但我转瞬间就满血复活了。这个包包的手柄不仅边油老化、脱落,表面的皮质也开裂破损了,修复的意义不大。所以我决定直接为它定制一条新的。首先,我要进行一比一的制版和剪裁,将皮料和衬料刷上胶水,然后进行粘合和压实。接着,我要对手柄两端进行均匀的打孔,方便后续的缝制。
诊治项目:五金养护年龄:1年至3年不等科别:养护门诊养护样本号:001症状说明:金属氧化变色、磨损、划痕、腐蚀01金属钩环镀色养护 五金掉色、无光泽、氧化,是奢侈品包包最常见的问题之一,其原因主要是:五金长期与空气接触易发生氧化,日常使用的磨损,以及平时与化妆品、香水等化学物质的接触,再加上使用过程中手部汗渍的侵蚀,个人汗渍的酸碱值不同,对金属表面的伤害程度也不一样。
这两天,走在时尚前沿欢喜轧闹猛的上海人又开始排队了!哪能回事体?原来,继喜茶fendi之后,lv也开始和咖啡店合作了在上海上线了3家lv限时书店。分别为:薄荷绿manner咖啡馆奉贤路300号落日黄plusone咖啡馆武康路286号经典红metalhands咖啡馆永嘉路37号在这3家快闪店,你可以花20元买到
丝巾可以直接系在包侧或者缠绕在包带上,增添包包的活泼感,也给炎热的夏日带来了色彩。看似不经意的点缀就会让你的包装焕然一新。百变的丝巾系法与包包的的碰撞,绑对了丝巾会让你的包包看起来很贵。今天,龙翔纺织分享给大家四种包包丝巾的系法。1.麻花辫系法丝巾在包柄中间对折,把丝带交叉绕过包柄,注意一定要绷紧丝带,系出来的辫子才会整齐不容易散,绕到底部的时候,把丝带穿过金属扣,然后打结。