fortran程序求教:我有MK检验程序,但不知道怎么把数据导入计算,求帮助

这个程序是执行时直接读入的，有说明语句

C***************************************************

C* N: SAMPLE SIZE *

C* NYEAR: FIRST YEAR OF THE TIME SERIES*

C* Y(N): ORIGINAL TIME SERIES *

C* UF(N): ORIGINAL SERIES OF U(LN) *

C* UB(N): COUNTER SERIES OF U(LN) *

C* A,B: CRITICAL VALUE 1.96 AND -1.96*

C***************************************************

程序中 READ 语句是读入对应数据的，先读入N,NYEAR，然后读入N个Y(N)，但是没有读入UF(N)，看说明语句也应该是原始数据。

Fortran多princeado(站内联系TA)M-K检验是做趋势分析还是突变检验，我编了相关的程序，不是很长的，用R语言编的，不过DPS系统也可以做M-K检验，如果有你可以试一下.guoguo1983(站内联系TA)要是用MK做趋势检验，用那个软件呀princeado(站内联系TA)做趋势分析的话可以用R语言来做，有个包可以做，另外我编了一段的程序，不是很长，DPS系统也可以做，很方便。shzrzheng(站内联系TA)知道原理后，自己编编程序做shhongyuan6331(站内联系TA)DPS:tiger28:leehui789(站内联系TA)现代气候统计诊断与预测技术 气象出版社挺方便的gishuazi(站内联系TA)建议看看陈彦光的书wuyang2468(站内联系TA)用R或者Matlab吧，发文章的话，R的公信度还是比较高的。erqie(站内联系TA)借地问一句，R语言是什么呀？horseechen(站内联系TA)我不是搞地学水文研究的，不过处理离散数据序列时检测突变点也想用这个检验方法，发现网上多数是中国人写的文章，而且错误百出。不这我发现这个方法太受秩序列长度的限制了，对于同样一组数据，整体趋势分析的结果跟分段分析的结果，在突变点检测上几乎没有一致的精度。lijun2011(站内联系TA)R也是一个统计软件，代码公开，不收费的软件，国外用的人多，国内不怎么用geonuist(站内联系TA)EXCEL编个宏，可以看看陈彦光写的《基于EXCEL的地理数据分析》，里面有介绍儿时梦想尚未成(站内联系TA)最近刚刚发现有个EXCEL插件很好用，XLSTATE2013,网上可以下载:victory:geonuist(站内联系TA)EXCEL就可以做，详见本论坛中陈彦光编写的《基于EXCEL的地理数据分析》

---