编辑代码

lambda2=function(dat){
    my=max(apply(dat,1,sum))
    mx=max(apply(dat,2,sum))
    sumy=sum(apply(dat,2,max))
    sumx=sum(apply(dat,1,max))
    n=sum(dat)
    lam=(sumx+sumy-mx-my)/(2*n-mx-my)
    return(lam)
}

c1=c(10,40,10)
c2=c(30,10,0)
dat65=cbind(c1,c2)
lambda2(dat65)