发了个打招呼的邮件过去~

1、对信息安全有非凡兴趣
入侵过学校电脑只为架vpn 0点后不断网... 不过现在没校园网了，电信宽带没法破-。-
2、CET-6
540还是多少来着..
3、参与过github项目
自己写了个小工具算不算... http://v2ex.com/t/68502#reply18（多线程支持已添加
4、爱读书
http://book.douban.com/people/ahxxm/collect 读过不少..烂书也不少..
5、着迷于创造与探索，无比享受它带来的快乐
... 探索能算 创造就差点
6、不拧巴，爱工作，更爱生活
那是必须的

xcode写c++很舒服啊... 虽然我还没入门

为什么要国内的... 免费的blogger就挺好啊

@ling0322 因为不知道怎么用python写这些...

我好像自己搞定了 = = V1一样，V2就都一样的，所以：
merge = merge(x,y, by = "V1")
merge = subset(merge, V9.x = V9.y)
再把行数相除就行了……
其实我都不知道自己在做什么，帮人做作业，抽象成这个过程。这些大概就是实际中需要的技能吧……
求大牛指导

科研应该也算是实际应用吧，我这正好有个问题不知道怎么解决，贴出来给lz看看……
图1 ： http://i.imgur.com/wPa6Fbw.png
pig_sus.blast是用blat处理原始数据后得到的输出文件，C090_1.blast（以及其他blast文件）是用于对比的数据。
pig_sus.blast中，第一列是query，有重复项，取evalue最小的、identity最大的（格式见http://edwards.sdsu.edu/labsite/index.php/ramys/238-blast-output-8），已经处理完了，R代码如下：
/*
x = read.csv("pig_sus.blast", sep = "", header = FALSE)
x = subset(x, subset = !duplicated(x[c("V1")]))
*/
然后合并两组数据发现没对齐，图2：http://i.imgur.com/lY4m7oC.png
合并代码如下：
/*
x = x[1:239,]
c = c("XV2", "XV9", "XV10", "YV2", "YV9", "YV10")
comp = data.frame(x[,2], x[,9], x[,10], y[,2], y[,9], y[,10])
colnames(comp) = c
*/
line 125时，x的V2出现了新项chr11，原本与y[125, ]对应的chr1变成了126。
图3：http://i.imgur.com/2ta3qdG.png
所以想根据V1来进行合并，把类似POR_C090_I10_I10的新项单列成之后，剩余项对齐合并，再进行对比。
对比原则是：两组数据的V2 V9 V10都一致则判定为一致，最后需要得到 一致项/总项数 这个比例。
--------
ps1:看id我们好像在twitter上互fo了哟~
ps2:为了处理数据还重新编译了一遍blat...

@codesaler 给代码让他帮忙运行一下？ 好吧…我暂时没法满足你这个需求… 先记下了>_<

@codesaler 不知道为什么你这两个回复配合头像戳到我笑点……
直接找某个帐号要存档最方便吧，调用api一条条把推文爬下来会不会有点奇怪？
推文存档的js文件里有media_url用来显示官方图床或者什么其他东西，url短链接和expanded_url还原后链接，（也许还得判定一下是不是需要再还原。。）

@codesaler 那么多图床，新手表示有点为难。。。

@AstroProfundis archive不好用吗....