天使不设防

摘要： do while 循环小练习 1 #include <stdio.h> 2 #include <stdlib.h> 3 4 int main(void) 5 { 6 int a; 7 do{ 8 fprintf(stdout,"%s\n","请输入一个整数:"); 9 scanf("%d",&a); 阅读全文

摘要： 改进了一下，利用zlib可以读取gz格式的压缩文件，也可以直接计算非压缩格式 #include <stdio.h> #include <stdlib.h> #include <string.h> #include <zlib.h> #define buff 1024 typedef unsigned 阅读全文

摘要： C语言小练习：计算非压缩fastq格式的GC含量 1 #include <stdio.h> 2 #include <stdlib.h> 3 #include <string.h> 4 #define buff 1024 5 6 typedef unsigned long long int u_llo 阅读全文

摘要： http://methodspopgen.com/methods-to-infer-populations-history/ 阅读全文

摘要： 使用copy，直接改变原片的值，而不是先创建一个副本。 阅读全文

摘要： getline() 函数无论一行多长，动态分配内存读入行 1 #include <stdio.h> 2 #include <stdlib.h> 3 #include <string.h> 4 5 int main(int argc,const char *argv[]) 6 { 7 FILE *fp 阅读全文

摘要： 需要安装包：ggforce，下面以R自带数据做局部放大演示。 require(ggplot2) require(ggforce) require(reshape2) data(CO2) co2<-melt(CO2,id.vars=c("Plant","Type","Treatment"),varia 阅读全文

摘要： set.seed(1) h1<-hist(rnorm(1000,100,5)) h2<-hist(rnorm(1000,99,5)) plot(h2,col=rgb(255,0,0,50,maxColorValue =255),border = NA) plot(h1,col=rgb(225,225 阅读全文

摘要： 一、芯片数据 此次拿到的illumina芯片数据并不是原始的数据，已经经过GenomeStudio软件处理成了finalreport文件，格式如下： 之前没处理过芯片数据，对于这种编码模式（Forward，top AB）的基因型数据很疑惑，查了很多资料，收效甚微。看过建明大神对芯片这块儿的介绍，发现 阅读全文

摘要： 实现：eval 1 a="indv1" 2 indv1="Sus1" 3 4 eval tmp='$'$a 5 echo $tmp //这里 echo 返回值为Sus1 阅读全文