Re: runs test for evaluation of model fit

*To*: mathgroup at smc.vnet.net*Subject*: [mg57227] Re: [mg57173] runs test for evaluation of model fit*From*: Clifford Martin <camartin at snet.net>*Date*: Sat, 21 May 2005 02:39:56 -0400 (EDT)*Sender*: owner-wri-mathgroup at wolfram.com

Why don't you just simply do differences between the second and third column such as: assume your data function name is data1. cols2= data1[[All,2]; cols3= data2[[All,3]]; diffs= cols2-cols3; maxdiff=Max[diffs]; mindiff=Min[diffs]; you can then do standard deviations and standard statistics stuff on you differences. I know it's kind of simple but it's a start. Cliff --- Csukas Attila <attila at biking.taiiku.tsukuba.ac.jp> wrote: > Dear Everybody, > > I am facing to a new problem and if it is possible > would like to ask > some help from experts as you are. > I have observed values (second in the brackets) for > three ids (first in > the brackets) and also have predicted values (third > in the brackets) > estimated by a model. I would like to use runs test > to prove that the > model fitted well, that is there is no significant > difference between > observed and predicted values for each id. > > Does anybody know how can it be done? Any help is > appreciated! Thanks > in advance! > > Out[10]= > {{id,obs,pred},{2,116.9,116.486},{2,122.1,122.073},{2,126.1,127.074},{2, > > 131.1,131.598},{2,137.1,135.899},{2,141.1,140.88},{2,148.3,149.053},{2, > 161.2,160.338},{2,165.9,166.697},{2, > > 167.8,168.316},{2,168.,168.617},{2,170.1,168.67},{4,120.8,121.477},{4, > > 128.2,127.612},{4,134.5,133.438},{4,138.9,139.369},{4,145.2,146.045},{4, > 153.7,154.016},{4,163.7,162.616},{4, > > 170.1,169.62},{4,172.1,173.697},{4,174.4,175.536},{4,177.3,176.255},{4, > > 177.3,176.518},{7,111.8,110.578},{7,115.5,116.887},{7,122.1,122.456},{ > > 7,126.8,127.377},{7,132.4,131.743},{7,136.4,135.698},{7,139.1,139.616},{ > > 7,145.3,144.669},{7,152.1,152.685},{7,161.,160.428},{7,163.2,163.488},{ > 7,164.1,164.176}} > > One of the mathgroup contributors had comments on > Wald-Wolfowitz test > but I have difficulties with the application for the > above data. Hints > I have found on W-W are the next: > > In[1]:= > x = Table[Random[Integer],{n = 25}] > > Out[1]= > {1,1,0,0,1,1,1,0,0,1,1,1,1,0,0,0,1,1,0,1,1,0,1,0,0} > > In[2]:= > {{n1 = Tr@x, n0 = n - n1}, r = Length@Split@x} > > Out[2]= > {{14,11},12} > > In[3]:= > N@{m = 1 + 2*n0*n1/n, sd = Sqrt[(m-1)(m-2)/(n-1)]} > > Out[3]= > {13.32,2.41059} > > In[6]:= > f[r_,n0_,n1_] := If[EvenQ@r, > With[{k = r/2 - 1},2*Binomial[n0-1,k > ]*Binomial[n1-1,k ]], > With[{k = (r-1)/2}, Binomial[n0-1,k > ]*Binomial[n1-1,k-1] + > (**) > Binomial[n0-1,k-1]*Binomial[n1-1,k ]]]; > rf=Table[{r,f[r,n0,n1]},{r,2,2Min[n0,n1]+Boole[n0??n1]}] > > Out[7]= > {{2,2},{3,23},{4,260},{5,1365},{6, > > 7020},{7,22230},{8,68640},{9,145860},{10,300300},{11,450450},{12,648648} > > ,{?_ > 13,702702},{ > > 14,720720},{15,566280},{16,411840},{17,231660},{18,115830},{19,45045},{2 > > 0,?_ > 14300},{21,3575},{22,572},{23,78}} > > In[8]:= > {#, # === Tr@rf[[All,2]]}&@Binomial[n,n1] > > Out[8]= > {4457400,True} > >