L<-function(rho,eta,xi,z,C,n,pn) { L<-0 for(j in 1:nrow(n)) { Lfac<-1 A<-length(z) for(a in 1:A) { if(a==A) shape<-rho*(1-xi)*n[j,a] else shape<-rho*((1-xi)*n[j,a]+xi*n[j,a+1]) if(z[a]