摘要:
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 阅读全文