function [samp]=impsamp(pars,L,sk)

samp.hmp=(1./rand(L,1) -1);
if sk==1
    samp.ga=gamrnd(pars.a_ga,1/pars.b_ga,L,1);
end
samp.df=gamrnd(pars.a_df,1/pars.b_df,L,1);
samp.la=gamrnd(pars.a_la,1/pars.b_la,L,1);
samp.ht=gamrnd(pars.a_ht,1/pars.b_ht,L,1);
samp.aal=gamrnd(pars.a_aal,1/pars.b_aal,L,1);
samp.bal=gamrnd(pars.a_bal,1/pars.b_bal,L,1);
samp.al=2*(betarnd(ones(L,1)*pars.a_al, ones(L,1)*pars.b_al))-1; 
samp.hb=pars.mhb + sqrt(pars.vhb)*randn(L,1);% samp.hb=reshape(samp.hb,L,[]);
samp.be=Mnormrnd(pars.mbe,pars.vbe,L);
samp.mu=Mnormrnd(pars.mmu,pars.vmu,L);
samp.hmu=Mnormrnd(pars.mhmu,pars.vhmu,L);