We analyze the three dimensional anisotropy of the galactic cosmic ray (GCR) intensities observed independently with a muon detector (MD) at Nagoya in Japan and neutron monitors over four solar activity cycles. We clearly see the phase of the free-space diurnal anisotropy shifting toward earlier hours around solar activity minima in A > 0 epochs, due to the reduced anisotropy component parallel to the mean magnetic field. This component is consistent with a rigidity independent spectrum, while the perpendicular anisotropy component increases with GCR rigidity. We suggest that this harder spectrum of the perpendicular component is due to contribution from the drift streaming. We find that the bidirectional latitudinal density gradient is positive in A > 0 epoch, while it is negative in A < 0 epoch, in accord with the drift model prediction. The radial density gradient of GCRs, on the other hand, varies with ∼11-year cycle with maxima (minima) in solar maximum (minimum) periods, but we find no significant difference between the radial gradients in A > 0 and A < 0 epochs.The corresponding parallel mean free path is larger in A < 0 than in A > 0.We also find, however, that parallel mean free path (radial gradient) appears to persistently increase (decreasing) in the last three cycles of weakening solar activity. We suggest that simple differences between these parameters in A > 0 and A < 0 epochs are seriously biased by these long-term trends.