Don't just run simple time-series regression to see if you get statistically significant betas. This procedure will not tell you if the factors are actually priced. You run a high risk of finding spurious correlations.
There is a fairly well established standard program to test factor models, called the Fama-MacBeth method. It is based on two sets of regressions. A factor model essentially gives you two (linear) equations that you can estimate using regression analysis:
$$ r_i = a_i + \beta_{i1} f_1 + … + \beta_{iK} f_K + \epsilon_i \ (1)$$
$$\mu = \lambda_0 + \beta_{i1} \lambda_1 + … + \beta_{iK} \lambda_K \ (2)$$
where $r_i$ is the (realized) return of security $i$, $a_i$ is the intercept of security I (often taken to be the risk-free rate), $\beta_{ij}$ are the factor sensitivities of security $i$ to factor $j$, and $f_j$ are the factor realizations of the relevant factors.
Given this structure, the APT implies that given absence of arbitrage, the expected return of asset i, $\mu_i$, is given by the second equation where the $\lambda_j$ are the factor risk premia and $\lambda_0$ is the model’s zero beta parameter.
The standard way to test the factors is to first run time-series regression of equation (1) using rolling windows to obtain the beta parameters. With monthly data, for example, one would usually estimate the regression using 5 years of data and take the estimated betas of that regression as the beta observation for the last date of those 5 years. Then you repeat the same setup by moving the estimation window one month further, and so on. This will give you a time series of beta estimates for each stock.
In the next step, you run cross-sectional regressions for each time $t$ using the returns on the securities at time $t$ as the dependent variable and the beta estimates as the independent variables. This will give you estimates of the factor risk premia for each time. Finally, you take the average of the factor risk premia over time to see whether they are statistically different from zero.
You will have to adjust your standard errors appropriately because of the fact that the betas in the cross sectional regressions are estimated. You can read up on the technical details of the method in chapter 12 of Cochrane’s book on asset pricing. More details on some other methods are given in Campbell, Lo, MacKinley’s book on econometrics. You should also check out the paper by Peterson which contains a detailed discussion of the statistical issues with this method. He also published STATA and R codes for this paper on his website.
Campbell, J., Lo, A., MacKinley, A. (1997) The Econometrics of Financial Markets. Princeton University Press
Cochrane, J. (2001) Asset Pricing (Revised Edition). Princeton University Press
Fama, Eugene F.; MacBeth, James D. (1973). "Risk, Return, and Equilibrium: Empirical Tests". Journal of Political Economy 81 (3): 607–636.
Petersen, Mitchell (2009). "Estimating Standard Errors in Finance Panel Data Sets: Comparing Approaches". Review of Financial Studies 22 (1): 435–480.