R语言 时间序列

Python08

R语言 时间序列,第1张

确定时间序列的周期一般用的是谱分析，小波分析方法，这些一般在网上能搜到相关文献！

时间序列是否平稳，ARMA(p,q)中的p,q的确定，这些方法在王文圣，丁晶等著作《随机水文学》中有详细介绍，中国水利水电出版社，第二版，你所提及的内容都有介绍，相信你会搞定！

有现成的包：matchprobes包 里面有个函数basecontent(seq)计算4中碱基每种的含量；

自己做的话：

#List是你的序列

unlist(strsplit(List,""))->sep.letter#把每个字母都单独分开

#遍历所有字母

count_a=0count_g=0count_t=0count_c=0

for(i in 1:length(sep.letter)){

if(sep.letter[i]=="a"){count_a=count_a+1}

if(sep.letter[i]=="g"){count_g=count_g+1}

.........................

}