注册 登录  
 加关注
   显示下一条  |  关闭
温馨提示!由于新浪微博认证机制调整,您的新浪微博帐号绑定已过期,请重新绑定!立即重新绑定新浪微博》  |  关闭

零售创新,创新那些事儿,SPSS,VBA

零售创新

 
 
 

日志

 
 
关于我

新浪微博,零售创新 研究经理,数据分析师 希望和市场研究和零售业的同事共同进步! 本博客发表的都是免费或试用的资料,如果有版权问题请发邮件wangli12a@163.com联系删除。 spss excel vba blog

网易考拉推荐

spss kappa analysis syntax  

2009-02-17 15:41:38|  分类: spss学习 |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |


SPSS SYNTAX

 

Note: For documentation, visit:  ftp://ftp.spss.com/pub/spss/statistics/nichols/macros/

Copy the text that follows into a SPSS syntax window and follow directions below.

 

 

 

**************************************************************

*  MACRO NAME:         MKAPPASC.SPS                          *

*                                                            *

*  README FILE:        MKAPPASC.RM                           *

*                                                            *

*  SPSS REQUIREMENTS:  Release 4.0 or above                  *

*                      Advanced Statistics Module            *

*                                                            *

*  AUTHOR:             David Nichols (nichols@spss.com)      *

*                                                            *

*  LAST UPDATED:       04/08/97                              *

**************************************************************

 

preserve.

set printback=off mprint=off.

save outfile='ka__tmp1.sav'.

define mkappasc (vars=!charend('/')).

set mxloops=1000.

count ms__=!vars (missing).

select if ms__=0.

matrix.

get x /var=!vars.

compute c=mmax(x).

compute y=make(nrow(x),c,0).

loop i=1 to nrow(x).

loop j=1 to ncol(x).

loop k=1 to c.

do if x(i,j)=k.

compute y(i,k)=y(i,k)+1.

end if.

end loop.

end loop.

end loop.

compute pe=msum((csum(y)/msum(y))&**2).

compute k=ncol(x).

compute pa=mssq(y)/(nrow(y)*k*(k-1))-(1/(k-1)).

compute kstat=(pa-pe)/(1-pe).

compute num=2*(pe-(2*k-3)*(pe**2)+2*(k-2)*msum((csum(y)/msum(y))&**3)).

compute den=nrow(y)*k*(k-1)*((1-pe)**2).

compute ase=sqrt(num/den).

compute z=kstat/ase.

compute sig=1-chicdf(z**2,1).

save {kstat,ase,z,sig} /outfile='ka__tmp2.sav'

     /variables=kstat,ase,z,sig.

end matrix.

get file='ka__tmp2.sav'.

formats all (f11.8).

variable labels kstat 'Kappa' /ase 'ASE' /z 'Z-Value' /sig 'P-Value'.

report format=list automatic align(center)

  /variables=kstat ase z sig

  /title "Estimated Kappa, Asymptotic Standard Error,"

         "and Test of Null Hypothesis of 0 Population Value".

get file='ka__tmp1.sav'.

!enddefine.

restore.

 

 

*NOTES: Dataset should be formatted such that

   number of rows = number of items;

   number of columns = number of raters,

   each cell entry represents a single rating.

*Insert column names below in place of rater1, rater2, and rater3.

*Click Run--All.

 

 

 

MKAPPASC VARS=rater1 rater2 rater3.


  评论这张
 
阅读(766)| 评论(0)
推荐 转载

历史上的今天

在LOFTER的更多文章

评论

<#--最新日志,群博日志--> <#--推荐日志--> <#--引用记录--> <#--博主推荐--> <#--随机阅读--> <#--首页推荐--> <#--历史上的今天--> <#--被推荐日志--> <#--上一篇,下一篇--> <#-- 热度 --> <#-- 网易新闻广告 --> <#--右边模块结构--> <#--评论模块结构--> <#--引用模块结构--> <#--博主发起的投票-->
 
 
 
 
 
 
 
 
 
 
 
 
 
 

页脚

网易公司版权所有 ©1997-2017