Many Polyakov loop models can be written in a dual formulation which is free of sign problem even when a non-vanishing baryon chemical potential is introduced in the action. Here, results of numerical simulations of a dual representation of one such effective Polyakov loop model at finite baryon density are presented. We compute various local observables such as energy density, baryon density, quark condensate and describe in details the phase diagram of the model. The regions of the first order phase transition and the crossover, as well as the line of the second order phase transition, are established. We also compute several correlation functions of the Polyakov loops.