Changed ddirich to calculate on log scale rather than logging after. Changed starting tuning values for r's. Fixed lambda_j_prop calculation (didn't work if lambda j not calculated first)