data(data)
factor.mat = data[1:342,2:4] # Fama-French 3-factor model
ret.mat = data[1:342,8:ncol(data)] # 25 size-BM portfolio returns
GRS.MLtest(ret.mat,factor.mat) # See column (iv), Table 9C of Fama-French (1993)
Run the code above in your browser using DataLab