In this study, ten-minute meteorological data-sets recorded at Burgos, Spain, are used to develop models of Photosynthetic Active Radiation () following two different procedures: multilinear regression and Artificial Neural Networks. Ten Meteorological Indices (MIs) are chosen as inputs to the models: clearness index (), diffuse fraction (), direct fraction (), Perez's clear sky index (ɛ), brightness index (), cloud cover (), air temperature (), pressure (), solar azimuth cosine (), and horizontal global irradiation (). The experimental data are clustered according to the sky conditions, following the CIE standard sky classification. A previous feature selection procedure established the most adequate MIs for modelling in clear, partial and overcast sky conditions. was the common MI used by all models and for all sky conditions. Additional variables were also included: the geometrical parameter, , and three variables related to the sky conditions, , and Both modelling methods, multilinear regression and ANN, yielded very high determination coefficients () with very close results in the models for each of the different sky conditions. Slight improvements can be observed in the ANN models. The results underline the equivalence of multilinear regression models and ANN models of PAR following previous feature selection procedures.Regional Government of Castilla y León, under projects BU021G19 and INVESTUN/19/BU/0004 and the Spanish Ministry of Science & Innovation under the I+D +i state program “Challenges Research Projects” (Ref. RTI2018-098900-B-I00). Diego Granados López expresses his thanks to the Junta de Castilla y León for economic support (PIRTU Program, ORDEN EDU/556/2019)
