tristanz / OpenIRT / changeset / 4e35255dd116

@@ -1047,6 +1047,7 @@ public:
 	  // Find the log density at the initial point.
     double x0 = parameter_.Value();
     double gx0 = parameter_.LogDensity(x0);
 	 	//cout << parameter_.Label() << " : " << x0 << " = " << gx0 << endl;
     // Determine slice level, in log terms
     double logy = gx0 - myrng.rexp(1);
@@ -1183,6 +1184,7 @@ public:
 	     for(int iter = 0; iter < number_of_iterations; ++iter) {
 	       cout << "Iteration " << iter << endl;
 				 for(int i = 0; i < steps_.size(); ++i) {
 					 //cout << "Step " << i << endl;
 	         steps_[i].DoStep();
+	       }
 	       //++show_progress;

@@ -62,10 +62,12 @@ int main(int argc, char* argv[]) {
 				B2PLParameter b(true, "b"+to_string<int>(items.id(j)), j);
 				sampler.AddStep(new SliceStep<B2PLParameter>(b, 1, -dInf, dInf));
 			} else if (items.type(j) == TYPE_3PL) {
 				// cout << items.id(j) << endl;
 				// Temporarily put some limits on these.  Discrete priors.
 				A3PLParameter a(true, "a"+to_string<int>(items.id(j)), j);
-				sampler.AddStep(new SliceStep<A3PLParameter>(a, 1, 0, dInf));
+				sampler.AddStep(new SliceStep<A3PLParameter>(a, 1, 0, 10));
 				B3PLParameter b(true, "b"+to_string<int>(items.id(j)), j);
-				sampler.AddStep(new SliceStep<B3PLParameter>(b, 1, -dInf, dInf));
+				sampler.AddStep(new SliceStep<B3PLParameter>(b, 1, -10, 10));
 				C3PLParameter c(true, "c"+to_string<int>(items.id(j)), j);
 				sampler.AddStep(new SliceStep<C3PLParameter>(c, 1, 0, 1));
+			}

@@ -101,6 +101,13 @@ void StartUp(int argc, char* argv[]) {
   responses.num_groups = max(responses.group);
   cout << "Responses added: " << responses.num_responses << endl;
   std::cout << "Max int limit: " << std::numeric_limits<int>::max() << std::endl;
 	//   cout << "Items ID:" << endl << items.id << endl;
 	//   cout << "Items Type:" << endl << items.type << endl;
 	//   cout << "Items a:" << endl << items.a << endl;
 	//   cout << "Items b:" << endl << items.b << endl;
 	// cout << "Items c:" << endl << items.c << endl;
+}
 /// 2PL Item Response Functions