We compare phenomenological results from (3 + 1)-dimensional quasiparticle anisotropic hydrodynamics (aHydroQP) with experimental data collected from RHIC 200 GeV Au-Au collisions. We compare identified particle spectra in different centrality classes, charged particle multiplicity versus pseudorapidity, identified particle multiplicity versus centrality across a wide range of particle species, identified particle elliptic flow versus transverse momentum, and charged particle elliptic flow as a function of transverse momentum and rapidity. We use the same aHydroQP and hadronic production codes that were used previously to describe LHC 2.76 TeV data. The aHydroQP hydrodynamic model includes the effects of both shear and bulk viscosities in addition to an infinite number of transport coefficients computed self-consistently in the relaxation-time approximation. To convert to the final-state hadrons, we use anisotropic Cooper-Frye freeze-out performed on a fixed-energy-density hypersurface and compute the production and feed-down by using a customized version of Therminator 2. We find good agreement with many heavy-ion collision observables using only smooth Glauber initial conditions parametrized by an initial central temperature of T-0 = 455 MeV, a constant shear-viscosity-to-entropy-density ratio eta/s = 0.179, and a switching (freeze-out) temperature of T-FO = 130 MeV.