Continuous gravitational wave searches with pulsar timing arrays: Maximization versus marginalization over pulsar phase parameters