lambda=function(dat){
    my=max(apply(dat,1,sum))
    sumy=sum(apply(dat,2,max))
    n=sum(dat)
    lam=(sumy-my)/(n-my)
    return(lam)
}

c1=c(18,20);c2=c(29,13)
c3=c(18,16);c4=c(12,20)
dat64=cbind(c1,c2,c3,c4)
lambda(dat64)