预测评分的程序

雨恬魅影

1. program the skill
   
2. real obs(160,17),model(160,17),Ps(17)
   
3. integer n1(17),n2(17),n0(17)
   
4. parameter (n=160,f1=5,f2=2)
   
5.       open(10,file='d:\ssx\two\160data.grb',form='binary')
   
6. do it=1,17
   
7. do j=1,160
   
8. read(10)obs(j,it)
   
9. enddo
   
10. do j=1,160
    
11. read(10)model(j,it)
    
12. enddo
    
13. enddo
    
14. do it=1,17
    
15. do j=1,160
    
16. if((abs(obs(j,it))>=50).and.(abs(model(j,it))>=50)) then
    
17. n1(it)=n1(it)+1
    
18. endif
    
19. if((abs(obs(j,it))>=20).and.(abs(model(j,it))>=20)) then
    
20.   if((abs(obs(j,it))<=50).and.(abs(model(j,it))<=50)) then
    
21.    n2(it)=n2(it)+1
    
22.   endif
    
23. endif
    
24. if(obs(j,it)*model(j,it)>0) then
    
25.   n0(it)=n0(it)+1
    
26. else if(obs(j,it)*model(j,it)<0) then
    
27.    if(abs(obs(j,it))<=15.and.abs(model(j,it))<=15) then     
    
28.     n0(it)=n0(it)+1
    
29.    endif
    
30. endif
    
31. enddo;enddo
    
32. open(11,file='d:\ssx\ps.txt')
    
33. do it=1,17
    
34. Ps(it)=100*(n0(it)+n1(it)*f1+n2(it)*f2)/(n+n1(it)*f1+n2(it)*f2)
    
35. write(11,*)'年份  预测评分Ps'
    
36. write(11,"(i4,f8.2)") it+1983,Ps(it)
    
37. enddo
    
38. end program
    

program the skill<br /> real obs(160,17),model(160,17),Ps(17)<br /> integer n1(17),n2(17),n0(17)<br /> parameter (n=160,f1=5,f2=2)<br />       open(10,file='d:\ssx\two\160data.grb',form='binary')<br /> do it=1,17<br /> do j=1,160<br /> read(10)obs(j,it)<br /> enddo<br /> do j=1,160<br /> read(10)model(j,it)<br /> enddo<br /> enddo<br /> do it=1,17<br /> do j=1,160<br /> if((abs(obs(j,it))>=50).and.(abs(model(j,it))>=50)) then <br /> n1(it)=n1(it)+1<br /> endif<br /> if((abs(obs(j,it))>=20).and.(abs(model(j,it))>=20)) then <br />   if((abs(obs(j,it))<=50).and.(abs(model(j,it))<=50)) then<br />    n2(it)=n2(it)+1<br />   endif<br /> endif<br /> if(obs(j,it)*model(j,it)>0) then<br />   n0(it)=n0(it)+1<br /> else if(obs(j,it)*model(j,it)<0) then<br />    if(abs(obs(j,it))<=15.and.abs(model(j,it))<=15) then     <br />     n0(it)=n0(it)+1<br />    endif<br /> endif<br /> enddo;enddo<br /> open(11,file='d:\ssx\ps.txt')<br /> do it=1,17<br /> Ps(it)=100*(n0(it)+n1(it)*f1+n2(it)*f2)/(n+n1(it)*f1+n2(it)*f2)<br /> write(11,*)'年份  预测评分Ps'<br /> write(11,"(i4,f8.2)") it+1983,Ps(it)<br /> enddo<br /> end program<br />

mofangbao

楼主再分享代码的同时，不妨简单介绍一下用的什么方法进行的评分，适用性如何，谢谢了

言深深

刚刚大致浏览了一下，程序大致就是讲read的数据进行了一次分类，比较不同的级别，然后根据级别加权计算。主要应该是此前的grd文件是什么？不知道楼主能够给出来一个test文件，如果涉密的话可以经过处理后上传数据。若能简单介绍一下grd文件的制作过程当真是非常好的
还是很感谢楼主分享的

雨恬魅影

具体说明都在PPT中了。

文字控。洛枫

After the recent sales,Bell & Ross Watches[/ur we have introduced many different colors and designs. Very interested in many of our customers of our blue CC logo bag Chanel Coco Chanel shoulder bag. Blue[url=http://www.replicarolexwatches.org.uk]replica rolex watches CC logo bag Chanel Coco Chanel shoulder bag, good quality in our stores are well known in the world.

rosewei

第一次接触TS评分，在论坛里就找到了楼主的代码，好好学习一下，谢谢楼主分享

雨恬魅影

rosewei 发表于 2012-6-15 11:56

                               
登录/注册后可看大图

第一次接触TS评分，在论坛里就找到了楼主的代码，好好学习一下，谢谢楼主分享

其实原理还比较简单的，加油呀！

兰北

非常感谢楼主，下载了学习一下！！！

shine090

蛮有意思的感觉，看看

鬼王

谢谢楼主，赞一个