In this thesis, we first present a variance estimation technique based on the standardized time series methodology for steady-state simulations. The proposed variance estimator has competitive bias and variance compared to the existing estimators in the literature. We also present the technique of rebatching to further reduce the bias and variance of our variance estimator. Second, we present two fully sequential indifference-zone procedures to select the best system from a number of competing simulated systems when best is defined by the maximum or minimum expected performance. These two procedures have parabola shaped continuation regions rather than the triangular continuation regions employed in several papers. The rocedures we present accommodate unequal and unknown ariances across systems and the use of common random numbers. However, we assume that basic observations are independent and identically normally distributed. Finally, we present procedures for finding a set of feasible or near-feasible systems among a finite number of simulated systems in the presence of multiple stochastic constraints, especially when the number of systems or constraints is large.

In this dissertation, we consider analytic and numeric approaches to the solution of probabilistic steady-state problems with specific applications in simulation and queueing theory. Our first objective on steady-state simulations is to develop new estimators for the variance parameter of a selected output process that have better performance than certain existing variance estimators in the literature. To complete our analysis of these new variance estimators, called linear combinations of overlapping variance estimators, we do the following: establish theoretical asymptotic properties of the new estimators; test the theoretical results on a battery of examples to see how the new estimators perform in practice; and use the estimators for confidence interval estimation for both the mean and the variance parameter. Our theoretical and empirical results indicate the new estimators' potential for improvements in accuracy and computational efficiency. Our second objective on steady-state simulations is to derive the expected values of various competing estimators for the variance parameter. In this research, we do the following: formulate the machinery to calculate the exact expected value of a given estimator for the variance parameter; calculate the exact expected values of various variance estimators in the literature; compute these expected values for certain stochastic processes with complicated covariance functions; and derive expressions for the mean squared error of the estimators studied herein. We find that certain standardized time series estimators outperform their competitors as the sample size becomes large. Our research on queueing theory focuses on pricing of the service provided to individual customers in a queueing system. We find sensitivity results that enable efficient computational procedures for dynamic pricing decisions for maximizing the long-run average reward in a queueing facility with the following properties: there are a fixed number of servers, each with the same constant service rate; the system has a fixed finite capacity; the price charged to a customer entering the system depends on the number of customers in the system; and the customer arrival rate depends on the current price of the service. We show that the sensitivity results considered significantly reduce the computational requirements for finding the optimal pricing policies.