*log using clogg2.log, replace

input f1-f4
44  11  38  62
59  41 147 293
23  11  13  27
27   8  16  27
258 57 105 110   
end

input r1-r4
1 1 1 1  
2 2 2 2  
3 3 3 3  
4 4 4 4
5 5 5 5

input c1-c4
1 2 3 4  
1 2 3 4  
1 2 3 4  
1 2 3 4  
1 2 3 4  

stack f1 r1 c1 f2 r2 c2 f3 r3 c3 f4 r4 c4, into(freq row column) clear

tab row column, sum(freq) nost nofreq

gen u_i = row
gen v_j = column

recode u_i (1=-0.0753) (2=-0.7763) (3=0.0983) (4=0.1550) (5=0.5983)
recode v_j (1=-0.7432) (2=-0.1265) (3=0.2713) (4=0.5984) 

gen inter=u_i*v_j  

xi:glm freq i.row i.column , family(poisson)


xi:glm freq i.row i.column inter, family(poisson)